Changeset 496 for branches/GRISLIv3/SOURCES/ablation_mod.f90

Timestamp:

05/03/24 15:49:57 (2 months ago)

Author:

dumas

Message:

ablation_mod : PDD compatible with negative accumulation due to sublimation

File:

• branches/GRISLIv3/SOURCES/ablation_mod.f90 (modified) (6 diffs)

branches/GRISLIv3/SOURCES/ablation_mod.f90

@@ -191 +191 @@
   IF (methode_abl.EQ.0) THEN
      IF (pdd_type.EQ.1) THEN
-        PDS(:,:)=ACC(:,:)/Csnow_2D(:,:)
-        SIMAX(:,:)=ACC(:,:)*CSI_2D(:,:)
+        PDS(:,:)=max(ACC(:,:),0.)/Csnow_2D(:,:)
+        SIMAX(:,:)=max(ACC(:,:),0.)*CSI_2D(:,:)
         PDSI(:,:)=SIMAX(:,:)/Cice_2D(:,:)
 ! avec regel de 60% puis fonte (2 premiers where) :
@@ -216 +216 @@
         endwhere
      ELSEIF (PDD_type.EQ.2) THEN  ! PDD Tarasov
-        PDS(:,:)=ACC(:,:)/Csnow_2D(:,:)
-        pr_ice_eq(:,:) = amax1(0.,((PRECIP(:,:)/DICE)-ACC(:,:)))  ! precipe liquide (ice equivalent)
+        PDS(:,:)=max(ACC(:,:),0.)/Csnow_2D(:,:)
+        pr_ice_eq(:,:) = amax1(0.,((PRECIP(:,:)/DICE)-max(ACC(:,:),0.)))  ! precipe liquide (ice equivalent)
 
         WHERE (PDD(:,:).LE.PDS(:,:))
@@ -225 +225 @@
         ENDWHERE
 
-        snowmelt(:,:) = amin1(totmelt(:,:),ACC(:,:))
+        snowmelt(:,:) = amin1(totmelt(:,:),max(ACC(:,:),0.))
 
 ! Deux formules possibles pour la capacité calorifique (en J/kg.K):
@@ -234 +234 @@
 
 !        where(snowmelt(:,:).lt.ACC(:,:))
-        refr2(:,:) = 2.2*(ACC(:,:)-snowmelt(:,:))-(cpsurf(:,:)/CL)*amin1(TANN(:,:),0.)
+        refr2(:,:) = 2.2*(max(ACC(:,:),0.)-snowmelt(:,:))-(cpsurf(:,:)/CL)*amin1(TANN(:,:),0.)
         refreezed_ice(:,:) = amin1(pr_ice_eq(:,:)+snowmelt(:,:),refr2(:,:))
 !        elsewhere
@@ -246 +246 @@
      ELSE ! pdd standard reeh
 !       (* Positive degrees required to melt the snow layer *)
-        PDS(:,:)=ACC(:,:)/Csnow_2D(:,:)
+        PDS(:,:)=max(ACC(:,:),0.)/Csnow_2D(:,:)
 !       (* Maximum amount of super. ice that can be formed *)
-        SIMAX(:,:)=ACC(:,:)*CSI_2D(:,:)
+        SIMAX(:,:)=max(ACC(:,:),0.)*CSI_2D(:,:)
 !       (* Pos. degrees required to melt the superimposed ice *)
         PDSI(:,:)=SIMAX(:,:)/Cice_2D(:,:)
@@ -274 +274 @@
      ENDIF
   ELSEIF ( methode_abl.EQ.1 ) THEN ! Insolation Temperature Melt equation, van den Berg, 2008
-     SIMAX(:,:)=ACC(:,:)*CSI
+     SIMAX(:,:)=max(ACC(:,:),0.)*CSI
      SIF(:,:)=SIMAX(:,:)
   ELSE