source: LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.F90 @ 5069

Last change on this file since 5069 was 4908, checked in by jyg, 8 months ago

New wake height formulation in the wake parametrization

File size: 9.4 KB
Line 
1MODULE lmdz_wake_ini
2IMPLICIT NONE
3
4  ! ============================================================================
5
6
7  ! But : Decrire le comportement des poches froides apparaissant dans les
8  ! grands systemes convectifs, et fournir l'energie disponible pour
9  ! le declenchement de nouvelles colonnes convectives.
10
11  ! State variables :
12  ! deltatw    : temperature difference between wake and off-wake regions
13  ! deltaqw    : specific humidity difference between wake and off-wake regions
14  ! sigmaw     : fractional area covered by wakes.
15  ! wdens      : number of wakes per unit area
16
17  ! -------------------------------------------------------------------------
18  ! Declaration de variables
19  ! -------------------------------------------------------------------------
20
21  ! Variables a fixer
22!jyg<
23!!  REAL, SAVE                                            :: stark, wdens_ref, coefgw, alpk
24  INTEGER, SAVE, PROTECTED                                    :: prt_level
25  REAL, SAVE, PROTECTED, DIMENSION(2)                         :: wdens_ref
26  REAL, SAVE, PROTECTED                                       :: stark, coefgw, alpk, wk_pupper
27!>jyg
28  REAL, SAVE, PROTECTED                                       :: crep_upper, crep_sol 
29  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, wk_pupper, crep_upper, crep_sol)
30
31  REAL, SAVE, PROTECTED                                       :: tau_cv
32  !$OMP THREADPRIVATE(tau_cv)
33
34   REAL, SAVE, PROTECTED                                       :: wk_delta_t_min
35  !$OMP THREADPRIVATE(wk_delta_t_min)
36
37   REAL, SAVE, PROTECTED                                       :: wk_frac_int_delta_t
38  !$OMP THREADPRIVATE(wk_frac_int_delta_t)
39
40  REAL, SAVE, PROTECTED                                       :: rzero, aa0 ! minimal wake radius and area
41  !$OMP THREADPRIVATE(rzero, aa0)
42
43  LOGICAL, SAVE, PROTECTED                                    :: ok_bug_gfl
44  !$OMP THREADPRIVATE(ok_bug_gfl)
45  LOGICAL, SAVE, PROTECTED                                    :: flag_wk_check_trgl
46  !$OMP THREADPRIVATE(flag_wk_check_trgl)
47  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_act
48  !$OMP THREADPRIVATE(iflag_wk_act)
49
50  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_check_trgl
51  !$OMP THREADPRIVATE(iflag_wk_check_trgl)
52  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_pop_dyn
53  !$OMP THREADPRIVATE(iflag_wk_pop_dyn)
54
55  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_profile
56  !$OMP THREADPRIVATE(iflag_wk_profile)
57
58  INTEGER, SAVE, PROTECTED                                    :: wk_nsub
59  !$OMP THREADPRIVATE(wk_nsub)
60
61  INTEGER, SAVE, PROTECTED                                    :: iflag_wk_new_ptop
62  !$OMP THREADPRIVATE(iflag_wk_new_ptop)
63
64  REAL, SAVE, PROTECTED                                       :: wdensinit ! Minimum wake density used to restart wakes from a wake-free state
65  !$OMP THREADPRIVATE(wdensinit)
66  REAL, SAVE, PROTECTED                                       :: wdensthreshold ! Threshold wake density below which wakes are killed
67  !$OMP THREADPRIVATE(wdensthreshold)
68  REAL, SAVE, PROTECTED                                       :: sigmad, hwmin, wapecut, cstart
69  !$OMP THREADPRIVATE(sigmad, hwmin, wapecut, cstart)
70  REAL, SAVE, PROTECTED                                       :: sigmaw_max
71  !$OMP THREADPRIVATE(sigmaw_max) 
72  REAL, SAVE, PROTECTED                                       :: dens_rate
73  !$OMP THREADPRIVATE(dens_rate)
74  REAL, SAVE, PROTECTED                                       :: epsilon_loc
75  !$OMP THREADPRIVATE(epsilon_loc)
76  REAL, SAVE, PROTECTED                                       :: epsim1,RG,RD
77  !$OMP THREADPRIVATE(epsim1,RG,RD)
78  REAL, SAVE, PROTECTED                                        ::smallestreal
79  !$OMP THREADPRIVATE(smallestreal)
80  REAL, SAVE, PROTECTED                                        :: wk_int_delta_t_min
81  !$OMP THREADPRIVATE(wk_int_delta_t_min)
82
83
84
85CONTAINS
86
87  ! =========================================================================
88  SUBROUTINE wake_ini(rg_in,rd_in,rv_in,prt_lev)
89  ! =========================================================================
90
91  ! **************************************************************
92  ! *
93  ! WAKE                                                        *
94  ! retour a un Pupper fixe                                *
95  ! *
96  ! written by   :  GRANDPEIX Jean-Yves   09/03/2000            *
97  ! modified by :   ROEHRIG Romain        01/29/2007            *
98  ! **************************************************************
99
100  ! -------------------------------------------------------------------------
101  ! Initialisations
102  ! -------------------------------------------------------------------------
103
104  USE ioipsl_getin_p_mod, ONLY : getin_p
105  real eps
106  integer, intent(in) :: prt_lev
107  real, intent(in) :: rg_in,rd_in,rv_in
108
109  smallestreal=tiny(smallestreal)
110!
111  prt_level=prt_lev
112  epsilon_loc=1.E-15
113  wapecut=1. ! previously 5.
114!
115  rzero = 5000.
116  CALL getin_p('rzero_wk', rzero)
117  aa0 = 3.14*rzero*rzero
118!
119  ! Essais d'initialisation avec sigmaw = 0.02 et hw = 10.
120!!  sigmad=0.005
121  sigmad=0.02
122  CALL getin_p('sigmad', sigmad)
123  hwmin=10.
124!
125!!wdensthreshold=1.e-12
126  wdensthreshold=1.e-14
127  wdensthreshold=2.e-14
128  CALL getin_p('wdensthreshold', wdensthreshold)
129!
130  IF (sigmad < 0.) THEN
131    sigmad = abs(sigmad)
132!!    wdensmin=sigmad/(3.14*rzero**2)
133    wdensinit=sigmad/(3.14*rzero**2)
134  ELSE
135    wdensinit = wdensthreshold/2.
136  ENDIF
137!
138!
139  ! cc nrlmd
140  sigmaw_max=0.4
141  dens_rate=0.1
142
143  eps = rd_in/rv_in
144  epsim1 = 1.0/eps - 1.0
145  RG=rg_in
146  RD=rd_in
147
148
149  ! cc
150  ! Longueur de maille (en m)
151  ! -------------------------------------------------------------------------
152
153  ! ALON = 3.e5
154  ! alon = 1.E6
155
156  ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei)
157
158  ! coefgw : Coefficient pour les ondes de gravite
159  ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE)
160  ! wdens : Densite surfacique de poche froide
161  ! -------------------------------------------------------------------------
162
163  ! cc nrlmd      coefgw=10
164  ! coefgw=1
165  ! wdens0 = 1.0/(alon**2)
166  ! cc nrlmd      wdens = 1.0/(alon**2)
167  ! cc nrlmd      stark = 0.50
168  ! CRtest
169  ! cc nrlmd      alpk=0.1
170  ! alpk = 1.0
171  ! alpk = 0.5
172  ! alpk = 0.05
173
174
175
176  crep_upper = 0.9
177  crep_sol = 1.0
178
179  ! Flag concerning the bug in gfl computation
180  ok_bug_gfl = .True.
181  call getin_p('ok_bug_gfl', ok_bug_gfl)
182
183  ! Get wapecut from parameter file
184  wapecut = 1.
185
186print*,'wapecut',wapecut
187  CALL getin_p('wapecut', wapecut)
188print*,'wapecut',wapecut
189
190  ! cc nrlmd Lecture du fichier wake_param.data
191
192
193  ! cc nrlmd Lecture du fichier wake_param.data
194  stark=0.33
195  CALL getin_p('stark',stark)
196  cstart = stark*sqrt(2.*wapecut)
197
198  alpk=0.25
199  CALL getin_p('alpk',alpk)
200 
201  wk_pupper=0.6
202  CALL getin_p('wk_pupper',wk_pupper)
203
204
205!jyg<
206!!  wdens_ref=8.E-12
207!!  CALL getin_p('wdens_ref',wdens_ref)
208  wdens_ref(1)=8.E-12
209  wdens_ref(2)=8.E-12
210  CALL getin_p('wdens_ref_o',wdens_ref(1))    !wake number per unit area ; ocean
211  CALL getin_p('wdens_ref_l',wdens_ref(2))    !wake number per unit area ; land
212!>jyg
213!
214!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
215!!!!!!!!!  Population dynamics parameters    !!!!!!!!!!!!!!!!!!!!!!!!!!!!
216!------------------------------------------------------------------------
217
218  iflag_wk_pop_dyn = 0
219  CALL getin_p('iflag_wk_pop_dyn',iflag_wk_pop_dyn) ! switch between wdens prescribed
220                                                    ! and wdens prognostic
221  iflag_wk_act = 0
222  CALL getin_p('iflag_wk_act',iflag_wk_act) ! 0: act(:)=0.
223                                            ! 1: act(:)=1.
224                                            ! 2: act(:)=f(Wape)
225
226  iflag_wk_profile = 0
227  CALL getin_p('iflag_wk_profile',iflag_wk_profile) ! switch between wdens prescribed
228                                                    ! and wdens prognostic
229 ! iflag_wk_profile = 0
230  iflag_wk_new_ptop = 0
231  CALL getin_p('iflag_wk_new_ptop',iflag_wk_new_ptop)
232
233  wk_nsub = 10
234  CALL getin_p('wk_nsub',wk_nsub)
235
236  tau_cv = 4000.
237  CALL getin_p('tau_cv', tau_cv)
238 
239  wk_delta_t_min = 0.
240  CALL getin_p('wk_delta_t_min', wk_delta_t_min)
241 
242  wk_int_delta_t_min = 10.
243  CALL getin_p('wk_int_delta_t_min', wk_int_delta_t_min)
244 
245  wk_frac_int_delta_t = 0.9
246  CALL getin_p('wk_frac_int_delta_t', wk_frac_int_delta_t)
247
248
249!------------------------------------------------------------------------
250
251  coefgw=4.
252  CALL getin_p('coefgw',coefgw)
253
254  WRITE(*,*) 'stark=', stark
255  WRITE(*,*) 'alpk=', alpk
256  WRITE(*,*) 'wk_pupper=', wk_pupper
257!jyg<
258!!  WRITE(*,*) 'wdens_ref=', wdens_ref
259  WRITE(*,*) 'wdens_ref_o=', wdens_ref(1)
260  WRITE(*,*) 'wdens_ref_l=', wdens_ref(2)
261!>jyg
262  WRITE(*,*) 'iflag_wk_pop_dyn=',iflag_wk_pop_dyn
263  WRITE(*,*) 'iflag_wk_act',iflag_wk_act
264  WRITE(*,*) 'coefgw=', coefgw
265
266  flag_wk_check_trgl=.false.
267  CALL getin_p('flag_wk_check_trgl ', flag_wk_check_trgl)
268  WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
269  WRITE(*,*) 'flag_wk_check_trgl OBSOLETE. Utilisr iflag_wk_check_trgl plutot'
270  iflag_wk_check_trgl=0 ; IF (flag_wk_check_trgl) iflag_wk_check_trgl=1
271  CALL getin_p('iflag_wk_check_trgl ', iflag_wk_check_trgl)
272  WRITE(*,*) 'iflag_wk_check_trgl=', iflag_wk_check_trgl
273 
274  WRITE(*,*) 'wk_delta_t_min=', wk_delta_t_min
275  WRITE(*,*) 'wk_int_delta_t_min=', wk_int_delta_t_min
276  WRITE(*,*) 'wk_frac_int_delta_t=', wk_frac_int_delta_t
277  WRITE(*,*) 'iflag_wk_new_ptop=', iflag_wk_new_ptop
278  WRITE(*,*) 'wk_nsub=', wk_nsub
279
280 RETURN
281
282END SUBROUTINE wake_ini
283
284
285
286END MODULE lmdz_wake_ini
Note: See TracBrowser for help on using the repository browser.