Changeset 4062 for LMDZ6/trunk/libf/phylmd/ice_sursat_mod.F90

Timestamp:

Jan 21, 2022, 6:50:11 PM (2 years ago)

Author:

oboucher

Message:

more keys for plane effects to avoid pb on outputs

File:

• LMDZ6/trunk/libf/phylmd/ice_sursat_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/ice_sursat_mod.F90

@@ -24 +24 @@
 !--Parametres de tuning
 REAL, PARAMETER :: chi=1.1, l_turb=50.0, tun_N=1.3, tun_ratqs=3.0
-!
-! Contrail cirrus
-LOGICAL, PARAMETER :: contrail_cirrus=.FALSE., h2o_emissions=.FALSE.
 !
 !--contrail cross section, typical value found in Freudenthaler et al, GRL, 22, 3501-3504, 1995
@@ -257 +254 @@
   INCLUDE "FCTTRE.h"
   INCLUDE "fisrtilp.h"
+  INCLUDE "clesphys.h"
+
   !
   ! Input
@@ -335 +334 @@
      ! flight_h2O is in kg H2O / s / cell
      !
-     IF (h2o_emissions) THEN 
+     IF (ok_plane_h2o) THEN 
        q = ( M_cell*q + flight_h2o(i,k)*dtime*(1.-q) ) / (M_cell + flight_h2o(i,k)*dtime*(1.-q) ) 
      ENDIF
@@ -634 +633 @@
 
      !--add a source of cirrus from aviation contrails
-     IF (contrail_cirrus) THEN
+     IF (ok_plane_contrail) THEN
        drneb_avi(i,k) = rnebss*flight_m(i,k)*contrail_cross_section/V_cell    !--tendency rneb due to aviation [s-1]
        drneb_avi(i,k) = MIN(drneb_avi(i,k), rnebss/dtime)                     !--majoration