Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/mass_redistribution_mod.F90

Timestamp:

Mar 19, 2026, 2:35:46 PM (11 days ago)

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/mass_redistribution_mod.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/mass_redistribution_mod.F90

@@ -13 +13 @@
        USE planete_mod, only: bp
        use comcstfi_mod, only: g
-       USE callkeys_mod, ONLY: water
+       USE callkeys_mod, ONLY: water, generic_condensation
+       use generic_tracer_index_mod, only: generic_tracer_index
+       use generic_cloud_common_h
        
        IMPLICIT NONE
@@ -96 +98 @@
       REAL zqm(nlayer+1,nq),zqm1(nlayer+1)
       REAL sigma(nlayer+1)
-
-!   local saved variables
-      LOGICAL, SAVE :: firstcall=.true.
-!$OMP THREADPRIVATE(firstcall)
-
-!----------------------------------------------------------------------
-
-!   Initialisation
-!   --------------
-!
-      IF (firstcall) THEN
-         firstcall=.false.       
-      ENDIF
-!
+      
+!    Generic condensation      
+      integer igcm_generic_vap, igcm_generic_ice
+      logical call_ice_vap_generic
+
 !======================================================================
 !    Calcul of h2o condensation  
@@ -159 +152 @@
          enddo
       endif
+      
+      if (generic_condensation) then
+      
+         do iq=1,nq
+
+           call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
+
+         end do ! do iq=1,nq loop on tracers
+      
+         do ig=1,ngrid
+            call Tsat_generic(pplev(ig,1)+zdmass_sum(ig,1)*g*ptimestep,Tsat(ig))
+         enddo
+#ifndef MESOSCALE
+         call writediagfi(ngrid,'Tsat','saturation temperature at surface','',2,Tsat)
+#endif
+         
+         do ig=1,ngrid
+            if (ztsrf(ig).gt.Tsat(ig)) then
+               zmassboil(ig)=(ptsrf(ig)-Tsat(ig))*pcapcal(ig)/RLVTT_generic/ptimestep
+               if ((zmassboil(ig)*ptimestep.gt.pqs(ig,igcm_generic_vap)).and.(nint(rnat(ig)).eq.1)) then
+                  zmassboil(ig)=pqs(ig,igcm_generic_vap)/ptimestep
+               endif
+               zmassboil(ig)=zmassboil(ig)*0.0 !momentary, should be 1. JL12
+               pdqsmr(ig,igcm_generic_vap)=-zmassboil(ig)
+               pdtsrfmr(ig)=-zmassboil(ig)*RLVTT_generic/pcapcal(ig)
+               ztsrf(ig)=ptsrf(ig)+pdtsrfmr(ig)*ptimestep
+            else
+               zmassboil(ig)=0.
+               pdtsrfmr(ig)=0.
+            endif
+         enddo
+      endif
 
 !     *************************
@@ -225 +250 @@
          zqm(1,1:nq)=0. ! most tracer do not condense !
          if (water) zqm(1,igcm_h2o_vap)=1. ! flux is 100% h2o at surface
+         if (generic_condensation) zqm(1,igcm_generic_vap)=1. ! flux is 100% condensable at surface
          
 !        Van Leer scheme: