Changeset 4392

Timestamp:

Jan 23, 2023, 2:36:19 PM (16 months ago)

Author:

evignon

Message:

super session de commentaires pour preparer linclusion des isotopes
dans lscp. Camille, Cecile, Niels, Sebastien et Etienne

File:

• LMDZ6/trunk/libf/phylmd/lscp_mod.F90 (modified) (20 diffs)

LMDZ6/trunk/libf/phylmd/lscp_mod.F90

@@ -329 +329 @@
 fcontrP(:,:) = 0.0
 
+!c_iso: variable initialisation for iso 
+
+
 !===============================================================================
 !              BEGINNING OF VERTICAL LOOP FROM TOP TO BOTTOM
@@ -341 +344 @@
         zt(i)=t(i,k)
         zq(i)=q(i,k)
+        !c_iso init of iso
     ENDDO
     
@@ -390 +394 @@
     ! --------------------------------------------------------------------
 
-    IF (iflag_evap_prec.GE.1) THEN
+    IF (iflag_evap_prec.GE.1) THEN ! useless -> to be removed
 
         ! Calculation of saturation specific humidity
@@ -406 +410 @@
 
             ! LTP: we only account for precipitation evaporation in the clear-sky (iflag_evap_prec=4).
+            ! c_iso: likely important to distinguish cs from neb isotope precipitation
+
                 IF (iflag_evap_prec.EQ.4) THEN
                     zrfl(i) = zrflclr(i)
@@ -478 +484 @@
 
 
-                ! New solid and liquid precipitation fluxes
+                ! New solid and liquid precipitation fluxes after evap and sublimation
                 zrfln(i) = Max(0.,zrfl(i) - zqev*(paprs(i,k)-paprs(i,k+1))  &
                          /RG/dtime)
                 zifln(i) = Max(0.,zifl(i) - zqevi*(paprs(i,k)-paprs(i,k+1)) &
                          /RG/dtime)
+
 
                 ! vapor, temperature, precip fluxes update
@@ -501 +508 @@
                 zifl(i) = zifln(i)
 
+                ! c_iso here call_reevap that updates isotopic zrfl, zifl (in inout)
+                !                         due to evap + sublim                                     
+                                           
+
                 IF (iflag_evap_prec.EQ.4) THEN
                     zrflclr(i) = zrfl(i)
@@ -511 +522 @@
                 ENDIF
 
-
-                ! CR Be careful: ice melting has been moved
+                ! c_iso duplicate for isotopes or loop on isotopes
+
+                ! Melting:
                 zmelt = ((zt(i)-273.15)/(ztfondue-273.15)) ! JYG
                 ! precip fraction that is melted
                 zmelt = MIN(MAX(zmelt,0.),1.)
+
                 ! Ice melting
                 IF (iflag_evap_prec.EQ.4) THEN
@@ -524 +537 @@
                     zrfl(i)=zrfl(i)+zmelt*zifl(i)
                 ENDIF
+                ! c_iso: melting of isotopic precipi with zmelt (no fractionation)
 
                 ! Latent heat of melting with precipitation thermalization
@@ -536 +550 @@
                     zifl(i)=zifl(i)*(1.-zmelt)
                 ENDIF
+                !c_iso: same for isotopic precip.
 
             ELSE
@@ -563 +578 @@
                 zq(i)=1.e-15
             ENDIF
+            ! c_iso: do something similar for isotopes
 
      ENDDO
@@ -826 +842 @@
             ENDIF
 
+            ! c_iso : routine that computes in-cloud supersaturation
+            ! c_iso condensation of isotopes (zcond, zsursat, zfice, zq in input) 
+
             ! temperature update due to phase change
             zt(i) = zt(i) + (1.-zfice(i))*zcond(i) & 
@@ -885 +904 @@
                 ziflclr(i) = ziflclr(i) + d_zifl_cld_clr(i) - d_zifl_clr_cld(i)
 
+                ! c_iso  do the same thing for isotopes precip
         ENDDO
     ENDIF
@@ -894 +914 @@
             zoliqi(i)= zoliq(i)*zfice(i)
             zoliql(i)= zoliq(i)*(1.-zfice(i))
+            ! c_iso : initialisation of zoliq* also for isotopes
             zrho(i,k)  = pplay(i,k) / zt(i) / RD
             zdz(i)   = (paprs(i,k)-paprs(i,k+1)) / (zrho(i,k)*RG)
@@ -925 +946 @@
                 zsnow=0.
                 zprecip=0.
+                ! c_iso same init for isotopes. Externalisation?
 
                 IF (zneb(i).EQ.seuil_neb) THEN
@@ -949 +971 @@
                     ! reduces the in-cloud water).
 
+
+                    ! warning: here should be zoliql?
                     IF ((iflag_cloudth_vert.GE.3).AND.(iflag_rain_incloud_vol.EQ.1)) THEN
                         zchau = zct   *dtime/REAL(ninter) * zoliq(i) &
@@ -966 +990 @@
                 zoliqi(i) = MAX(zoliq(i)*zfice(i)-1.*zsnow  , 0.0)
                 zoliq(i)  = MAX(zoliq(i)-zprecip , 0.0)
+
+                ! c_iso: call isotope_conversion (for readibility) or duplicate 
 
                 radliq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1)
@@ -1021 +1047 @@
                     zqpreci(i) = max (zqpreci(i) - Deltaqprecl - Deltaq, 0.)
 
+                    ! c_iso : mv here condensation of isotopes + redispatchage en precip 
+
                 IF (iflag_evap_prec.EQ.4) THEN
                     zrflcld(i) = zrflcld(i)+zqprecl(i) &
@@ -1035 +1063 @@
                    *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
                 ENDIF
+                ! c_iso : same for isotopes
  
            ENDIF ! rneb>0
@@ -1076 +1105 @@
         d_qi(i,k) = zfice(i)*zoliq(i)
         d_q(i,k) = zq(i) - q(i,k)
+        ! c_iso: same for isotopes
         d_t(i,k) = zt(i) - t(i,k)
     ENDDO
@@ -1191 +1221 @@
       snow(i) = zifl(i)
       rain(i) = zrfl(i)
+      ! c_iso final output
   ENDDO