Changeset 4099

Timestamp:

Mar 15, 2022, 2:06:03 PM (3 years ago)

Author:

oboucher

Message:

update to ice_sursat module

Location:

LMDZ6/trunk/libf/phylmd

Files:

• ice_sursat_mod.F90 (modified) (4 diffs)
• lscp_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/ice_sursat_mod.F90

@@ -22 +22 @@
 REAL, PARAMETER :: erf_lim = 5., eps = 1.e-10
 !
-!--Parametres de tuning
-REAL, PARAMETER :: chi=1.1, l_turb=50.0, tun_N=1.3, tun_ratqs=3.0
+!--Tuning parameters (and their default values)
+!
+!--chi gère la répartition statistique de la longueur des frontières
+!  entre les zones nuages et ISSR/ciel clair sous-saturé. Gamme de valeur :
+!  chi > 1, je n'ai pas regardé de limite max (pour chi = 1, la longueur de
+!  la frontière entre ne nuage et l'ISSR est proportionnelle à la
+!  répartition ISSR/ciel clair sous-sat dans la maille, i.e. il n'y a pas
+!  de favorisation de la localisation de l'ISSR près de nuage. Pour chi = inf, 
+!  le nuage n'est en contact qu'avec de l'ISSR, quelle que soit la taille 
+!  de l'ISSR dans la maille.)
+!
+!--l_turb est la longueur de mélange pour la turbulence. 
+!  dans les tests, ça n'a jamais été modifié pour l'instant. 
+!
+!--tun_N est le paramètre qui contrôle l'importance relative de N_2 par rapport à N_1. 
+!  La valeur est comprise entre 1 et 2 (tun_N = 1 => N_1 = N_2)
+!
+!--tun_ratqs : paramètre qui modifie ratqs en fonction de la valeur de
+!  alpha_cld selon la formule ratqs_new = ratqs_old / ( 1 + tun_ratqs *
+!  alpha_cld ). Dans le rapport il est appelé beta. Il varie entre 0 et 5
+!  (tun_ratqs = 0 => pas de modification de ratqs).
+
+REAL, SAVE :: chi=1.1, l_turb=50.0, tun_N=1.3, tun_ratqs=3.0
+!$OMP THREADPRIVATE(chi,l_turb,tun_N,tun_ratqs)
 !
 !--contrail cross section, typical value found in Freudenthaler et al, GRL, 22, 3501-3504, 1995
@@ -32 +54 @@
 
 !*******************************************************************
-
+!
+SUBROUTINE ice_sursat_init()
+
+  USE print_control_mod, ONLY: lunout
+  USE ioipsl_getin_p_mod, ONLY : getin_p
+
+  IMPLICIT NONE
+
+  CALL getin_p('flag_chi',chi)
+  CALL getin_p('flag_l_turb',l_turb)
+  CALL getin_p('flag_tun_N',tun_N)
+  CALL getin_p('flag_tun_ratqs',tun_ratqs)
+
+  WRITE(lunout,*) 'flag_chi = ', chi
+  WRITE(lunout,*) 'flag_l_turb = ', l_turb
+  WRITE(lunout,*) 'flag_tun_N = ', tun_N
+  WRITE(lunout,*) 'flag_tun_ratqs = ', tun_ratqs
+
+END SUBROUTINE ice_sursat_init
+
+!*******************************************************************
+!
 SUBROUTINE airplane(debut,pphis,pplay,paprs,t_seri)
 
@@ -97 +140 @@
       !
       IF ( klona /= klon_glo ) THEN
-        WRITE(lunout,*), 'klona & klon_glo =', klona, klon_glo
+        WRITE(lunout,*) 'klona & klon_glo =', klona, klon_glo
         CALL abort_physic(modname,'problem klon in aircraft_phy.nc file',1)
       ENDIF
       !
       IF ( ntimea /= 12 ) THEN
-        WRITE(lunout,*), 'ntimea=', ntimea
+        WRITE(lunout,*) 'ntimea=', ntimea
         CALL abort_physic(modname,'problem ntime<>12 in aircraft_phy.nc file',1)
       ENDIF
@@ -207 +250 @@
 
 !********************************************************************
-! temporary routine to initialise flight_m and test a flight corridor
+! simple routine to initialise flight_m and test a flight corridor
+!--Olivier Boucher - 2021 
+!
 SUBROUTINE flight_init()
   USE dimphy

LMDZ6/trunk/libf/phylmd/lscp_mod.F90

@@ -279 +279 @@
   !$OMP THREADPRIVATE(itap)
 
-
   LOGICAL lognormale(klon)
   LOGICAL convergence(klon)
@@ -286 +285 @@
   !$OMP THREADPRIVATE(appel1er)
 
-
   CHARACTER (len = 20) :: modname = 'lscp'
   CHARACTER (len = 80) :: abort_message
-
   
   DATA appel1er /.TRUE./
-
-
-
 
 !===============================================================================
@@ -370 +364 @@
     ENDDO
 
+    IF (ok_ice_sursat) CALL ice_sursat_init()
+
     appel1er = .FALSE.
 
@@ -378 +374 @@
 zfrac_lessi = 0.
 
-  
 ! Initialisation of output variables (JYG):
 prfl(:,:) = 0.0
@@ -1198 +1193 @@
     ! --------------------------------
 
-
     ! Outputs:
     ! Precipitation fluxes at layer interfaces
@@ -1210 +1204 @@
         d_t(i,k) = zt(i) - t(i,k)
     ENDDO
-