Changeset 2954

Timestamp:

May 2, 2023, 4:58:20 PM (20 months ago)

Author:

emillour

Message:

Generic PCM:
Improve diagnostics given by callcorrk (get min/max over whole globe, not
just over the master domain).
While at it also moved "Nmix_co2" and "radfixed" from callkeys_mod to radii_mod
since they are only used in that module.
EM

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/callcorrk.F90 (modified) (4 diffs)
• libf/phystd/callkeys_mod.F90 (modified) (2 diffs)
• libf/phystd/inifis_mod.F90 (modified) (1 diff)
• libf/phystd/radii_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1786 +1786 @@
 More code tidying: turn aeropacity, aeroptproperties, gfluxi, gfluxv,
 sfluxi and sfluxv into modules.
+
+== 02/05/2023 == EM
+Improve diagnostics given by callcorrk (get min/max over whole globe, not
+just over the master domain).
+While at it also moved "Nmix_co2" and "radfixed" from callkeys_mod to radii_mod
+since they are only used in that module.

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

@@ -48 +48 @@
       use generic_cloud_common_h, only: Psat_generic, epsi_generic
       use generic_tracer_index_mod, only: generic_tracer_index
+      use planetwide_mod, only: planetwide_maxval, planetwide_minval
       implicit none
 
@@ -220 +221 @@
 !$OMP THREADPRIVATE(qvap_deep)
 
+      REAL :: maxvalue,minvalue
+
 !===============================================================
 !           I.a Initialization on first call
@@ -449 +452 @@
          if ((iaer.eq.iaero_co2).and.tracer.and.(igcm_co2_ice.gt.0)) then ! Treat condensed co2 particles.
             call co2_reffrad(ngrid,nlayer,nq,pq,reffrad(1,1,iaero_co2))
+
+            call planetwide_maxval(reffrad(:,:,iaero_co2),maxvalue)
+            call planetwide_minval(reffrad(:,:,iaero_co2),minvalue)
             if (is_master) then
-               print*,'Max. CO2 ice particle size = ',maxval(reffrad(1:ngrid,1:nlayer,iaer))/1.e-6,' um'
-               print*,'Min. CO2 ice particle size = ',minval(reffrad(1:ngrid,1:nlayer,iaer))/1.e-6,' um'
+               print*,'Max. CO2 ice particle size = ',maxvalue/1.e-6,' um'
+               print*,'Min. CO2 ice particle size = ',minvalue/1.e-6,' um'
             end if
          end if
@@ -458 +464 @@
             call h2o_reffrad(ngrid,nlayer,pq(1,1,igcm_h2o_ice),pt, &
                              reffrad(1,1,iaero_h2o),nueffrad(1,1,iaero_h2o))
+            
+            call planetwide_maxval(reffrad(:,:,iaero_h2o),maxvalue)
+            call planetwide_minval(reffrad(:,:,iaero_h2o),minvalue)
             if (is_master) then
-               print*,'Max. H2O cloud particle size = ',maxval(reffrad(1:ngrid,1:nlayer,iaer))/1.e-6,' um'
-               print*,'Min. H2O cloud particle size = ',minval(reffrad(1:ngrid,1:nlayer,iaer))/1.e-6,' um'
+               print*,'Max. H2O cloud particle size = ',maxvalue/1.e-6,' um'
+               print*,'Min. H2O cloud particle size = ',minvalue/1.e-6,' um'
+            end if
+
+            call planetwide_maxval(nueffrad(:,:,iaero_h2o),maxvalue)
+            call planetwide_minval(nueffrad(:,:,iaero_h2o),minvalue)
+            if (is_master) then
+               print*,'Max. H2O cloud particle variance = ',maxvalue
+               print*,'Min. H2O cloud particle variance = ',minvalue
             end if
          endif

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

@@ -37 +37 @@
       logical,save :: varactive
       logical,save :: varfixed
-      logical,save :: radfixed
       logical,save :: sedimentation
       logical,save :: generic_condensation
       logical,save :: generic_rain
-!$OMP THREADPRIVATE(varactive,varfixed,radfixed,sedimentation,generic_condensation,generic_rain)
+!$OMP THREADPRIVATE(varactive,varfixed,sedimentation,generic_condensation,generic_rain)
       logical,save :: water ,watercond, waterrain, moistadjustment
 !$OMP THREADPRIVATE(water, watercond, waterrain, moistadjustment)
@@ -104 +103 @@
 !$OMP THREADPRIVATE(tplanckmin,tplanckmax,dtplanck)
       real,save :: topdustref
-      real,save :: Nmix_co2
       real,save :: dusttau
       real,save :: Fat1AU
       real,save :: stelTbb
-!$OMP THREADPRIVATE(topdustref,Nmix_co2,dusttau,Fat1AU,stelTbb)
+!$OMP THREADPRIVATE(topdustref,dusttau,Fat1AU,stelTbb)
       real,save :: Tstrat
       real,save :: tplanet

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

@@ -11 +11 @@
   use init_print_control_mod, only: init_print_control
   use radinc_h, only: ini_radinc_h, naerkind
+  use radii_mod, only: radfixed, Nmix_co2
   use datafile_mod, only: datadir
   use comdiurn_h, only: sinlat, coslat, sinlon, coslon

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

@@ -3 +3 @@
 !==================================================================
 !  module to centralize the radii calculations for aerosols
-! OK for water but should be extended to other aerosols (CO2,...)
 !==================================================================
       
-!     water cloud optical properties
-
-      use callkeys_mod, only: radfixed,Nmix_co2
-      
+!     CO2 cloud properties (initialized in inifis)
+      real,save :: Nmix_co2 ! Number mixing ratio of CO2 ice particles
+
+      ! flag to specify if we assume a constant fixed radius for particles
+      logical,save :: radfixed ! initialized in inifis
+
+!     water cloud optical properties (initialized in su_aer_radii below)
       real, save ::  rad_h2o
       real, save ::  rad_h2o_ice