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

Timestamp:

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

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/condensation_generic_mod.F90 (modified) (9 diffs)

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

@@ -6 +6 @@
     subroutine condensation_generic(ngrid,nlayer,nq,ptimestep, pplev, pplay,    &
                 pt, pq, pdt, pdq, pdtlsc, pdqvaplsc, pdqliqlsc, rneb)
-        use ioipsl_getin_p_mod, only: getin_p !-> to get the metallicity 
+        use callkeys_mod, only: thermo_phy,metallicity,qvap_deep,qvap_top,align_strato_cold_trap
         use generic_cloud_common_h
         USE tracer_h
         USE mod_phys_lmdz_para, only: is_master
         use generic_tracer_index_mod, only: generic_tracer_index
+        use thermo_mod
+        use comcstfi_mod
         IMPLICIT none
 
@@ -45 +47 @@
         REAL, intent(out) :: rneb(ngrid,nlayer,nq)    ! fraction nuageuse
 
-!       Options : 
-        real, save :: metallicity !metallicity of planet
-        REAL, SAVE :: qvap_deep   ! deep mixing ratio of vapor when simulating bottom less planets
-        REAL, SAVE :: qvap_top   ! top mixing ratio of vapor when simulating bottom less planets
-        logical, save :: align_strato_cold_trap
-!$OMP THREADPRIVATE(metallicity, qvap_deep, qvap_top, align_strato_cold_trap)
-
 !       Local variables
 
@@ -66 +61 @@
         DOUBLE PRECISION zq(ngrid)
         DOUBLE PRECISION zcond(ngrid),zcond_iter
-        DOUBLE PRECISION zqs(ngrid)
-        real zt(ngrid),local_p,psat_tmp,dlnpsat_tmp,Lcp,zqs_temp,zdqs
+        DOUBLE PRECISION zqs(ngrid), zqs_tmp(ngrid), zt_temp(ngrid)
+        real Lcp(ngrid,nlayer)
+        real zt(ngrid),local_p,psat_tmp,dlnpsat_tmp,zqs_temp,zdqs
         real zqs_temp_1, zqs_temp_2
         integer igcm_generic_vap, igcm_generic_ice ! index of the vap and ice of generic_tracer 
@@ -77 +73 @@
         DOUBLE PRECISION zx_q(ngrid)
         DOUBLE PRECISION zy_q(ngrid)
-        LOGICAL,SAVE :: firstcall=.true.
-!$OMP THREADPRIVATE(firstcall)
-        IF (firstcall) THEN
-                write(*,*) "value for metallicity? "
-                metallicity=0.0 ! default value
-                call getin_p("metallicity",metallicity)
-                write(*,*) " metallicity = ",metallicity
-
-                write(*,*) "Deep generic tracer vapor mixing ratio ? (no effect if negative) "
-                qvap_deep=-1. ! default value
-                call getin_p("qvap_deep",qvap_deep)
-                write(*,*) " qvap_deep = ",qvap_deep   
-
-                write(*,*) "top generic tracer vapor mixing ratio ? (no effect if negative) "
-                qvap_top=-1. ! default value
-                call getin_p("qvap_top",qvap_top)
-                write(*,*) " qvap_top = ",qvap_top 
-
-                write(*,*) " align_strato_cold_trap ? "
-                align_strato_cold_trap=.false. ! default value
-                call getin_p("align_strato_cold_trap",align_strato_cold_trap)
-                write(*,*) " align_strato_cold_trap = ",align_strato_cold_trap
-
-                firstcall = .false.
-        ENDIF
+        DOUBLE PRECISION zzx_q(ngrid)
+        DOUBLE PRECISION zzy_q(ngrid)
+        DOUBLE PRECISION zzt(ngrid)
+
 !       Initialisation of outputs and local variables
         pdtlsc(1:ngrid,1:nlayer)  = 0.0
@@ -125 +100 @@
                         metallicity_coeff=constants_metallicity_coeff(iq)
 
-                        Lcp=RLVTT_generic/cpp ! need to be init here
+                        Lcp(:,:)=RLVTT_generic/cpp(:,:) ! need to be init here
 
                         !  Vertical loop (from top to bottom)
                         DO k = nlayer, 1, -1
                                 zt(1:ngrid)=pt(1:ngrid,k)+pdt(1:ngrid,k)*ptimestep
+                                zzt(1:ngrid)=pt(1:ngrid,k)+pdt(1:ngrid,k)*ptimestep
 
                                 ! Computes Psat and the partial condensation
@@ -139 +115 @@
                                         endif
                                         zx_q(i) = pq(i,k,igcm_generic_vap)+pdq(i,k,igcm_generic_vap)*ptimestep
+                                        zzx_q(i) = zx_q(i)
 
                                         call Psat_generic(zt(i),local_p,metallicity,psat_tmp,zqs_temp)
                                         zy_q(i) = pq(i,k,igcm_generic_ice)+pdq(i,k,igcm_generic_ice)*ptimestep
+                                        zzy_q(i) = zy_q(i)
                                         
                                         if ((zx_q(i) .le. zqs_temp) .and. (zy_q(i) .eq. 0.)) then
@@ -153 +131 @@
 
                                                 ! iterative process to stabilize the scheme when large water amounts JL12
+                                                
+                                                SELECT CASE (TRIM(thermo_phy))
+                                                
+                                                CASE ('thermo_uni_ideal')
+                                                
                                                 zcond(i) = 0.0d0
                                                 Do nn=1,nitermax  
                                                         call Psat_generic(zt(i),local_p,metallicity,psat_tmp,zqs_temp)
                                                         zqs(i)=zqs_temp
-                                                        call Lcpdqsat_generic(zt(i),local_p,psat_tmp,zqs_temp,zdqs,dlnpsat_tmp)
+                                                        call Lcpdqsat_generic(zt(i),local_p,cpp(i,k),psat_tmp,zqs_temp,zdqs,dlnpsat_tmp)
                                                         zcond_iter = alpha*(zx_q(i)-zqs(i))/(1.d0+zdqs)
                                                         !zcond can be negative here
                                                         zx_q(i) = zx_q(i) - zcond_iter
                                                         zcond(i) = zcond(i) + zcond_iter
-                                                        zt(i) = zt(i) + zcond_iter*Lcp
+                                                        zt(i) = zt(i) + zcond_iter*Lcp(i,k)
                                                         if (ABS(zcond_iter/alpha/zqs(i)).lt.qthreshold) exit
                                                         if (nn.eq.nitermax) print*,'itermax in largescale'
@@ -173 +156 @@
                                                 !                       so the line below is to check how much we can evaporate.
                                                 !                       If there is no ice available, zcond(i) will be 0. (first condidition of "if")
-
+                                                
                                                 zcond(i)=MAX(zcond(i),-(pq(i,k,igcm_generic_ice)+pdq(i,k,igcm_generic_ice)*ptimestep))
+                                                zt(i) = zzt(i)+Lcp(i,k)*zcond(i)
+                                                END SELECT
                                         
                                         endif
@@ -190 +175 @@
                                 pdqvaplsc(1:ngrid,k,igcm_generic_vap)  = - zcond(1:ngrid) 
                                 pdqliqlsc(1:ngrid,k,igcm_generic_ice) = - pdqvaplsc(1:ngrid,k,igcm_generic_vap)
-                                pdtlsc(1:ngrid,k)  = pdtlsc(1:ngrid,k) + pdqliqlsc(1:ngrid,k,igcm_generic_ice)*Lcp
+                                pdtlsc(1:ngrid,k)  = pdtlsc(1:ngrid,k) + (zt(1:ngrid)-zzt(1:ngrid))/ptimestep
 
                         Enddo ! k= nlayer, 1, -1