Changeset 4908 for LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.F90

Timestamp:

Apr 15, 2024, 7:30:55 PM (5 weeks ago)

Author:

jyg

Message:

New wake height formulation in the wake parametrization

File:

• LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.F90

@@ -55 +55 @@
   INTEGER, SAVE, PROTECTED                                    :: iflag_wk_profile
   !$OMP THREADPRIVATE(iflag_wk_profile)
+
+  INTEGER, SAVE, PROTECTED                                    :: wk_nsub
+  !$OMP THREADPRIVATE(wk_nsub)
 
   INTEGER, SAVE, PROTECTED                                    :: iflag_wk_new_ptop
@@ -75 +78 @@
   REAL, SAVE, PROTECTED                                        ::smallestreal
   !$OMP THREADPRIVATE(smallestreal)
+  REAL, SAVE, PROTECTED                                        :: wk_int_delta_t_min
+  !$OMP THREADPRIVATE(wk_int_delta_t_min)
 
 
@@ -226 +231 @@
   CALL getin_p('iflag_wk_new_ptop',iflag_wk_new_ptop)
 
+  wk_nsub = 10
+  CALL getin_p('wk_nsub',wk_nsub)
+
   tau_cv = 4000.
   CALL getin_p('tau_cv', tau_cv)
@@ -231 +239 @@
   wk_delta_t_min = 0.
   CALL getin_p('wk_delta_t_min', wk_delta_t_min)
+  
+  wk_int_delta_t_min = 10.
+  CALL getin_p('wk_int_delta_t_min', wk_int_delta_t_min)
   
   wk_frac_int_delta_t = 0.9
@@ -262 +273 @@
   
   WRITE(*,*) 'wk_delta_t_min=', wk_delta_t_min
+  WRITE(*,*) 'wk_int_delta_t_min=', wk_int_delta_t_min
   WRITE(*,*) 'wk_frac_int_delta_t=', wk_frac_int_delta_t
   WRITE(*,*) 'iflag_wk_new_ptop=', iflag_wk_new_ptop
+  WRITE(*,*) 'wk_nsub=', wk_nsub
 
  RETURN