Changeset 3280 for trunk/LMDZ.GENERIC/libf

Timestamp:

Mar 20, 2024, 4:29:25 PM (7 months ago)

Author:

jleconte

Message:

new flag align_strato_cold_trap (Noe C)

File:

• trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90 (modified) (4 diffs)

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

@@ -49 +49 @@
         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 :: perfect_vap_profile
-!$OMP THREADPRIVATE(metallicity, qvap_deep, qvap_top, perfect_vap_profile)
+        logical, save :: align_strato_cold_trap
+!$OMP THREADPRIVATE(metallicity, qvap_deep, qvap_top, align_strato_cold_trap)
 
 !       Local variables
@@ -68 +68 @@
         DOUBLE PRECISION zqs(ngrid)
         real zt(ngrid),local_p,psat_tmp,dlnpsat_tmp,Lcp,zqs_temp,zdqs
-        real zqs_temp_1, zqs_temp_2, zqs_temp_3
-        integer igcm_generic_vap, igcm_generic_ice! index of the vap and ice of generic_tracer
-        ! CHARACTER(len=*) :: tname_ice
+        real zqs_temp_1, zqs_temp_2
+        integer igcm_generic_vap, igcm_generic_ice ! index of the vap and ice of generic_tracer 
         ! evaporation calculations
         REAL dqevap(ngrid,nlayer),dtevap(ngrid,nlayer)     
@@ -95 +94 @@
                 call getin_p("qvap_top",qvap_top)
                 write(*,*) " qvap_top = ",qvap_top 
-                
-                write(*,*) " perfect_vap_profile ? "
-                perfect_vap_profile=.false. ! default value
-                call getin_p("perfect_vap_profile",perfect_vap_profile)
-                write(*,*) " perfect_vap_profile = ",perfect_vap_profile
+
+                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.
@@ -195 +194 @@
                         Enddo ! k= nlayer, 1, -1
 
-                        if ((perfect_vap_profile) .and. (ngrid.eq.1)) then
-                                ! perfect_vap_profile is a mode that should a priori only be used in 1D:
-                                ! as it is written below, it aims to force the vap profile:
-                                ! - below condensation, it is forced to qvap_deep
-                                ! - at condensation levels, it is forced to 99% of sat
-                                ! - above the cold trap, it is forced to the value allowed by the cold trap
-
-                                ! perfect_vap_profile can be customed as you want
-                                
-                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
-
-                                k_cold_trap = 2
+                        if (align_strato_cold_trap .and. (ngrid.eq.1)) then
+                                ! aligns GCS vapor in the stratosphere to the quantity fixed by the cold trap
+                                ! works only in 1D
+                                tau = 10. * ptimestep
+                                ! tau is in seconds and must not be lower than the physical step time.
+
+                                ! Look for the cold trap
+                                k_cold_trap = 1 ! initialization of cold trap layer
+                                zqs_temp_2  = 1 ! initialization of minimum of saturation
                                 do k=2,nlayer-1
-
-                                        zt(1)=pt(1,k-1)+pdt(1,k-1)*ptimestep
-                                        call Psat_generic(zt(1),pplay(1,k-1),metallicity,psat_tmp,zqs_temp_1)
-                                        zt(1)=pt(1,k)+pdt(1,k)*ptimestep
-                                        call Psat_generic(zt(1),pplay(1,k),metallicity,psat_tmp,zqs_temp_2)
-                                        zt(1)=pt(1,k+1)+pdt(1,k+1)*ptimestep
-                                        call Psat_generic(zt(1),pplay(1,k+1),metallicity,psat_tmp,zqs_temp_3)
-
-                                        if ((zqs_temp_1 .gt. zqs_temp_2) .and. (zqs_temp_3 .gt. zqs_temp_2)) then
+                                        zt(1) = pt(1,k)+pdt(1,k)*ptimestep
+                                        call Psat_generic(zt(1),pplay(1,k),metallicity,psat_tmp,zqs_temp_1)
+                                        if (zqs_temp_1 .lt. zqs_temp_2) then
                                                 k_cold_trap = k
-                                                exit
+                                                zqs_temp_2 = zqs_temp_1
                                         endif
                                 enddo
-                                if (k_cold_trap .lt. nlayer) then
-                                        do k=k_cold_trap+1,nlayer
-                                                pdqvaplsc(1,k,igcm_generic_vap) = (pq(1,k_cold_trap,igcm_generic_vap) - pq(1,k,igcm_generic_vap))/tau - pdq(1,k,igcm_generic_vap)
-                                        enddo
-                                endif
-
-                                do k=1,k_cold_trap
-                                        zt(1)=pt(1,k)+pdt(1,k)*ptimestep
-                                        call Psat_generic(zt(1),pplay(1,k),metallicity,psat_tmp,zqs_temp)
-                                        if (zqs_temp .gt. qvap_deep) then
-                                                pdqvaplsc(1,k,igcm_generic_vap)  = (qvap_deep - pq(1,k,igcm_generic_vap))/tau - pdq(1,k,igcm_generic_vap)
-                                        endif
-                                        if (zqs_temp .lt. qvap_deep) then
-                                                pdqvaplsc(1,k,igcm_generic_vap)  = (0.99*zqs_temp - pq(1,k,igcm_generic_vap))/tau - pdq(1,k,igcm_generic_vap)
-                                        endif
+
+                                ! aligning 
+                                do k=k_cold_trap,nlayer
+                                        pdqvaplsc(1,k,igcm_generic_vap) = (zqs_temp_2 - pq(1,k,igcm_generic_vap))/tau - pdq(1,k,igcm_generic_vap)
+                                        pdtlsc(1,k) = 0.
                                 enddo
-
-                                pdqliqlsc(1:ngrid,:,igcm_generic_ice) = 0.
-                                pdtlsc(1:ngrid,:)  = 0.
                         endif
-
+                        
                         if (qvap_deep >= 0.) then
-                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
+                                tau = 10. * ptimestep
                                 ! brings lower generic vapor ratio to a fixed value.
                                 ! tau is in seconds and must not be lower than the physical step time.
                                 pdqvaplsc(1:ngrid,1,igcm_generic_vap) = (qvap_deep - pq(1:ngrid,1,igcm_generic_vap))/tau - pdq(1:ngrid,1,igcm_generic_vap)
                         endif
+
                         if (qvap_top >= 0.) then
-                                tau = 10. * ptimestep ! tau must not be lower than the physical step time. 
+                                tau = 10. * ptimestep
                                 ! brings lower generic vapor ratio to a fixed value.
                                 ! tau is in seconds and must not be lower than the physical step time.
                                 pdqvaplsc(1:ngrid,nlayer,igcm_generic_vap) = (qvap_top - pq(1:ngrid,nlayer,igcm_generic_vap))/tau - pdq(1:ngrid,nlayer,igcm_generic_vap)
                         endif
+                        
                 endif !if(call_ice_vap_generic)
         enddo ! iq=1,nq