Changeset 4828 for LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90

Timestamp:

Feb 17, 2024, 9:25:57 PM (3 months ago)

Author:

evignon

Message:

petites modification neige soufflee

File:

• LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90

@@ -110 +110 @@
     DO i = 1, ngrid
 
+          rhoair  = pplay(i,k) / zt(i) / RD
+          dz      = (paprs(i,k)-paprs(i,k+1)) / (rhoair*RG)
+          ! BS fall velocity assumed to be constant for now
+          zvelo(i) = fallv_bs
+
+
           IF (zt(i) .GT. RTT) THEN
              
@@ -117 +123 @@
            
              taumeltbs=taumeltbs0*exp(-max(0.,(zt(i)-RTT)/(tbsmelt-RTT)))
-             dqmelt=min(zqbs(i),zqbs(i)/taumeltbs*dtime)
+             dqmelt=min(zqbs(i),-1.*zqbs(i)*(exp(-dtime/taumeltbs)-1.))
              maxdqmelt= max(RCPD*(1.0+RVTMP2*(zq(i)+zqbs(i)))*(zt(i)-RTT)/(RLMLT+RLVTT),0.)
-             dqmelt=min(max(dqmelt,0.),maxdqmelt)
+             dqmelt=min(dqmelt,maxdqmelt)
              ! q update, melting + evaporation
              zq(i) = zq(i) + dqmelt
@@ -130 +136 @@
               ! negative celcius temperature     
               ! Sublimation scheme   
-
+             
               rhoair=pplay(i,k)/zt(i)/RD 
               dqsub = 1./rhoair*coef_eva_bs*(1.0-zq(i)/qsi(i))*((rhoair*zvelo(i)*zqbs(i))**expo_eva_bs)*dtime
@@ -148 +154 @@
 
         ! Sedimentation scheme
-           
-        rhoair  = pplay(i,k) / zt(i) / RD
-        dz      = (paprs(i,k)-paprs(i,k+1)) / (rhoair*RG)
-        ! BS fall velocity assumed to be constant for now
-        zvelo(i) = fallv_bs
+        !----------------------
 
         sedimn(i) = rhoair*zqbs(i)*zvelo(i)
@@ -173 +175 @@
                  zt(i) = zt(i) - zqbs(i) * RLSTT/RCPD/(1.0+RVTMP2*(zq(i)+zqbs(i)))
               ELSE
-               zt(i) = zt(i) - zqbs(i) * (RLVTT+RLMLT)/RCPD/(1.0+RVTMP2*(zq(i)+zqbs(i)))
+                 zt(i) = zt(i) - zqbs(i) * (RLVTT+RLMLT)/RCPD/(1.0+RVTMP2*(zq(i)+zqbs(i)))
               ENDIF
               zqbs(i)=0.