Changeset 5796 for LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

Timestamp:

Aug 4, 2025, 3:03:07 PM (11 days ago)

Author:

aborella

Message:

Additional diags for contrails + simplified coupling between deep conv and cirrus clouds + small modifsin RRTM for RF of contrails alone

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90 (modified) (4 diffs)

LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

@@ -12 +12 @@
   LOGICAL, PROTECTED :: ok_icefra_lscp, ok_new_lscp
   LOGICAL, PROTECTED :: ok_plane_contrail
-  LOGICAL, PROTECTED :: ok_higher_cirrus_cover
   REAL, PROTECTED :: bl95_b0, bl95_b1 ! Parameter in B&L 95-Formula
   REAL, ALLOCATABLE :: latitude_deg(:)
@@ -26 +25 @@
   REAL, PROTECTED :: zepsec
   REAL, PROTECTED :: eff2vol_radius_contrails=0.7
+  REAL, PROTECTED :: rei_min_contrails=10.
   REAL, PROTECTED :: rho_ice=920. ! Ice crystal density (assuming spherical geometry) [kg/m3]
   REAL, PARAMETER :: thres_tau=0.3, thres_neb=0.001
@@ -45 +45 @@
 !$OMP THREADPRIVATE(rei_coef, rei_min_temp)
 !$OMP THREADPRIVATE(zepsec) 
-!$OMP THREADPRIVATE(eff2vol_radius_contrails, rho_ice)
-!$OMP THREADPRIVATE(ok_plane_contrail, ok_higher_cirrus_cover)
+!$OMP THREADPRIVATE(eff2vol_radius_contrails, rei_min_contrails, rho_ice)
+!$OMP THREADPRIVATE(ok_plane_contrail)
 
   
@@ -117 +117 @@
     CALL getin_p('eff2vol_radius_contrails', eff2vol_radius_contrails)
     write(lunout,*)'eff2vol_radius_contrails=',eff2vol_radius_contrails
-    CALL getin_p('ok_higher_cirrus_cover', ok_higher_cirrus_cover)
-    write(lunout,*)'ok_higher_cirrus_cover=',ok_higher_cirrus_cover
-
-    IF ( ok_higher_cirrus_cover .AND. iflag_rei .GT. 0 ) THEN
-      abort_message = 'in cloud_optics, ok_higher_cirrus_cover is not implemented for iflag_rei > 0'
-      CALL abort_physic (modname,abort_message,1)
-    ENDIF
+    CALL getin_p('rei_min_contrails', rei_min_contrails)
+    write(lunout,*)'rei_min_contrails=',rei_min_contrails
    
   END SUBROUTINE cloud_optics_prop_ini