Changeset 4013 for LMDZ6/branches/Ocean_skin/libf/phylmd/inlandsis/sisvat_qsn.F

Timestamp:

Nov 19, 2021, 4:58:59 PM (3 years ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/inlandsis/sisvat_qsn.F (modified) (6 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/inlandsis/sisvat_qsn.F

@@ -61 +61 @@
       use VARxSV
       use VARySV
+      use surface_data, only: is_ok_slush,opt_runoff_ac
+
 
       IMPLICIT NONE
@@ -235 +237 @@
  
       DO ikl=1,knonv
-       DO isn=min(nsno,isnoSV(ikl)+1),1,-1
+
+      DO isn=min(nsno,isnoSV(ikl)+1),1,-1
 ! EV          DO isn=nsno,1,-1
 C +--Energy, store Previous Content
@@ -243 +246 @@
      .                + ro__SV(ikl,isn) * Cn_dSV * dTSnow
      .                                           * dzsnSV(ikl,isn)
-
-          Tsave       = TsisSV(ikl,isn)
-
           TsisSV(ikl,isn) =                        TfSnow
  
@@ -312 +312 @@
           rdzNEW      = WaFrez + rdzsno
           ro__SV(ikl,isn) =      rdzNEW /max(epsi, dzsnSV(ikl,isn))
-
-! EV: condition on Enfrez
-!          if (EnFrez .eq. 0.) then
-          
-          TsisSV(ikl,isn) = Tsave
-!          else
           TsisSV(ikl,isn) =      TfSnow
      .                + EnFrez /(Cn_dSV *max(epsi, rdzNEW)        )
-!          end if
           EExcsv(ikl) =          EExcsv(ikl)     - EnFrez
           wer_SV(ikl) = WaFrez
@@ -499 +492 @@
           rusnew      = rusnSV(ikl) * SWf_SV(ikl)
  
-          if(isnoSV(ikl)<=1) rusnew = 0.
+          if(isnoSV(ikl)<=1 .OR. opt_runoff_ac) rusnew = 0.
           !if(ivgtSV(ikl)>=1) rusnew = 0.
  
 c #EU                        rusnew = 0.
-c #AC                        rusnew = 0.
+c #AC               rusnew = 0.
+
           RnofSV(ikl) = RnofSV(ikl)
      .                +(rusnSV(ikl) - rusnew     ) / dt__SV
@@ -545 +539 @@
         ENDDO
  
-C +--Slush Formation (CAUTION: ADD RunOff Possibility before Activation)
+C +--Slush Formation (Activated. CAUTION: ADD RunOff Possibility before Activation)
 C +  ---------------  ^^^^^^^  ^^^
  
- 
-c #SU DO  ikl=1,knonv
-c #SU  DO isn=1,isnoSV(ikl)
-c #SU     kSlush = min(1,max(0,isn+1-ispiSV(ikl)))        ! Slush Switch
+      IF (is_ok_slush) THEN 
+
+      DO  ikl=1,knonv
+       DO isn=1,isnoSV(ikl)
+          kSlush = min(1,max(0,isn+1-ispiSV(ikl)))        ! Slush Switch
  
 C +--Available Additional Pore   Volume [-]
 C +  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
-c #SU     PorVol = 1. - ro__SV(ikl,isn)                    ! [--]
-c #SU.           *(1. - eta_SV(ikl,isn))/ ro_Ice           !
-c #SU.           -      eta_SV(ikl,isn)                    !
-c #SU.                 *ro__SV(ikl,isn) / ro_Wat           !
-c #SU     PorVol =  max(PorVol          , zero  )          !
-c #SU     zWater =      dzsnSV(ikl,isn) * PorVol * 1000.   ! [mm] OR [kg/m2]
-c #SU.           * (1. -SWS_SV(ikl)                        ! 0 <=> freezing
-c #SU.                *(1 -min(1,iabs(isn-isnoSV(ikl)))))  ! 1 <=> isn=isnoSV
-c #SU     zSlush =  min(rusnSV(ikl)     , zWater)          ! [mm] OR [kg/m2]
-c #SU     ro_new      =(dzsnSV(ikl,isn) * ro__SV(ikl,isn)  !
-c #SU.                 +zSlush                           ) !
-c #SU.            / max(dzsnSV(ikl,isn) , epsi           ) !
-c #SU     if(ro_new<ro_Ice+20) then ! MAX 940kg/m3         !
-c #SU      rusnSV(ikl)  = rusnSV(ikl)          - zSlush    ! [mm] OR [kg/m2]
-c #SU      RuofSV(ikl,4)= max(0.,RuofSV(ikl,4) - zSlush/dt__SV)
-c #SU      eta_SV(ikl,isn) =(ro_new - ro__SV(ikl,isn)      !
-c #SU.                     *(1.     - eta_SV(ikl,isn)))    !
-c #SU.                / max (ro_new , epsi            )    !
-c #SU      ro__SV(ikl,isn) =      ro_new                   !
-c #SU     endif
-c #SU   END DO
-c #SU END DO
- 
+          PorVol = 1. - ro__SV(ikl,isn)                    ! [--]
+     .           *(1. - eta_SV(ikl,isn))/ ro_Ice           !
+     .           -      eta_SV(ikl,isn)                    !
+     .                 *ro__SV(ikl,isn) / ro_Wat           !
+          PorVol =  max(PorVol          , zero  )          !
+          zWater =      dzsnSV(ikl,isn) * PorVol * 1000.   ! [mm] OR [kg/m2]
+     .           * (1. -SWS_SV(ikl)                        ! 0 <=> freezing
+     .                *(1 -min(1,iabs(isn-isnoSV(ikl)))))  ! 1 <=> isn=isnoSV
+          zSlush =  min(rusnSV(ikl)     , zWater)          ! [mm] OR [kg/m2]
+          ro_new      =(dzsnSV(ikl,isn) * ro__SV(ikl,isn)  !
+     .                 +zSlush                           ) !
+     .            / max(dzsnSV(ikl,isn) , epsi           ) !
+          if(ro_new<ro_Ice+20) then ! MAX 940kg/m3         !
+           rusnSV(ikl)  = rusnSV(ikl)          - zSlush    ! [mm] OR [kg/m2]
+           RuofSV(ikl,4)= max(0.,RuofSV(ikl,4) - zSlush/dt__SV)
+           eta_SV(ikl,isn) =(ro_new - ro__SV(ikl,isn)      !
+     .                     *(1.     - eta_SV(ikl,isn)))    !
+     .                / max (ro_new , epsi            )    !
+           ro__SV(ikl,isn) =      ro_new                   !
+          endif
+        END DO
+      END DO
+      END IF
  
 C +--Impact of the Sublimation/Deposition on the Surface Mass Balance