Changeset 2726 for trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

Timestamp:

Jun 21, 2022, 11:05:47 AM (3 years ago)

Author:

aslmd

Message:

Some comments and metallicity now generalised

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (5 diffs)

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

@@ -421 +421 @@
       real rneb_generic(ngrid,nlayer,nq) ! GCS cloud fraction (generic condensation).
       real psat_tmp_generic
-      real metallicity
-      
+      real, save :: metallicity ! metallicity of planet --- is not used here, but necessary to call function Psat_generic
+!$OMP THREADPRIVATE(metallicity)
+
       logical clearsky ! For double radiative transfer call. By BC
       
@@ -716 +717 @@
 
          endif
+
+!        Set metallicity for GCS
+!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+         metallicity=0.0 ! default value --- is not used here but necessary to call function Psat_generic
+         call getin_p("metallicity",metallicity) ! --- is not used here but necessary to call function Psat_generic
          
 !        Set some parameters for the thermal plume model
@@ -1581 +1587 @@
             end if
 
-            if (.not. water) then ! Compute GCS cloud fraction. For now we can not have both water cloud fraction and GCS cloud fraction 
+            if (.not. water) then 
+               ! Compute GCS (Generic Condensable Specie) cloud fraction. For now we can not have both water cloud fraction and GCS cloud fraction
+               ! Water is the priority
+               ! If you have set water and generic_condensation, then cloudfrac will be water cloudfrac
+               !
+               ! If you have set generic_condensation (and not water) and you have set several GCS
+               ! then cloudfrac will be only the cloudfrac of the last generic tracer
+               ! (Because it is rewritten every tracer in the loop)
+               !
+               ! Maybe one should create a cloudfrac_generic(ngrid,nlayer,nq) with 3 dimensions, the last one for tracers
 
                ! Let's loop on tracers
@@ -1592 +1607 @@
                         do ig=1,ngrid
                            cloudfrac(ig,l)=MAX(rneb_generic(ig,l,iq),rneb_man(ig,l)) 
-                           ! To be improved : for now, cloudfrac is only the cloudfrac of the last generic tracer
-                           ! Because it is rewritten every tracer
-                           ! One should create a cloudfrac_generic(ngrid,nlayer,nq)
                         enddo
                      enddo
@@ -2121 +2133 @@
                do l = 1, nlayer
                   do ig=1,ngrid
-                     metallicity = 0.
                      call Psat_generic(zt(ig,l),pplay(ig,l),metallicity,psat_tmp_generic,qsat_generic(ig,l,iq))
                      RH_generic(ig,l,iq) = zq(ig,l,igcm_generic_vap) / qsat_generic(ig,l,iq)