source: LMDZ6/trunk/libf/phylmd/wake_ini_mod.F90 @ 4161

Last change on this file since 4161 was 4085, checked in by fhourdin, 2 years ago

Reecriture de wake sur de nouveaux standards
Permet de rejouer les calculs (fonctionalité replay) et
pave le chemin vers les GPU.
Frédéric, Jean-Yves et Etienne
+ corrections makelmdz/create_make_gcm.
File size: 6.9 KB
Line 
1MODULE wake_ini_mod
2IMPLICIT NONE
3
4
5  ! ============================================================================
6
7
8  ! But : Decrire le comportement des poches froides apparaissant dans les
9  ! grands systemes convectifs, et fournir l'energie disponible pour
10  ! le declenchement de nouvelles colonnes convectives.
11
12  ! State variables :
13  ! deltatw    : temperature difference between wake and off-wake regions
14  ! deltaqw    : specific humidity difference between wake and off-wake regions
15  ! sigmaw     : fractional area covered by wakes.
16  ! wdens      : number of wakes per unit area
17
18  ! -------------------------------------------------------------------------
19  ! Déclaration de variables
20  ! -------------------------------------------------------------------------
21
22  ! Variables à fixer
23!jyg<
24!!  REAL, SAVE                                            :: stark, wdens_ref, coefgw, alpk
25  INTEGER, SAVE, PROTECTED                                         :: prt_level
26  REAL, SAVE, PROTECTED, DIMENSION(2)                              :: wdens_ref
27  REAL, SAVE, PROTECTED                                            :: stark, coefgw, alpk
28!>jyg
29  REAL, SAVE, PROTECTED                                            :: crep_upper, crep_sol 
30  !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, crep_upper, crep_sol)
31
32  REAL, SAVE, PROTECTED                                            :: tau_cv
33  !$OMP THREADPRIVATE(tau_cv)
34  REAL, SAVE, PROTECTED                                            :: rzero, aa0 ! minimal wake radius and area
35  !$OMP THREADPRIVATE(rzero, aa0)
36
37  LOGICAL, SAVE, PROTECTED                                         :: flag_wk_check_trgl
38  !$OMP THREADPRIVATE(flag_wk_check_trgl)
39  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_act
40  !$OMP THREADPRIVATE(iflag_wk_act)
41
42  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_check_trgl
43  !$OMP THREADPRIVATE(iflag_wk_check_trgl)
44  INTEGER, SAVE, PROTECTED                                         :: iflag_wk_pop_dyn
45  !$OMP THREADPRIVATE(iflag_wk_pop_dyn)
46
47  REAL, SAVE, PROTECTED                                            :: wdensmin
48  !$OMP THREADPRIVATE(wdensmin)
49  REAL, SAVE, PROTECTED                                            :: sigmad, hwmin, wapecut, cstart
50  !$OMP THREADPRIVATE(sigmad, hwmin, wapecut, cstart)
51  REAL, SAVE, PROTECTED                                            :: sigmaw_max
52  !$OMP THREADPRIVATE(sigmaw_max) 
53  REAL, SAVE, PROTECTED                                            :: dens_rate
54  !$OMP THREADPRIVATE(dens_rate)
55  REAL, SAVE, PROTECTED                                            :: epsilon_loc
56  !$OMP THREADPRIVATE(epsilon_loc)
57  REAL, SAVE, PROTECTED                                            :: epsim1,RG,RD
58  !$OMP THREADPRIVATE(epsim1,RG,RD)
59
60
61
62CONTAINS
63
64  ! =========================================================================
65  SUBROUTINE wake_ini(rg_in,rd_in,rv_in,prt_lev)
66  ! =========================================================================
67
68  ! **************************************************************
69  ! *
70  ! WAKE                                                        *
71  ! retour a un Pupper fixe                                *
72  ! *
73  ! written by   :  GRANDPEIX Jean-Yves   09/03/2000            *
74  ! modified by :   ROEHRIG Romain        01/29/2007            *
75  ! **************************************************************
76
77  ! -------------------------------------------------------------------------
78  ! Initialisations
79  ! -------------------------------------------------------------------------
80
81  USE ioipsl_getin_p_mod, ONLY : getin_p
82  real eps
83  integer, intent(in) :: prt_lev
84  real, intent(in) :: rg_in,rd_in,rv_in
85
86  prt_level=prt_lev
87  epsilon_loc=1.E-15
88  wapecut=1. ! previously 5.
89  ! Essais d'initialisation avec sigmaw = 0.02 et hw = 10.
90  sigmad=0.02
91  hwmin=10.
92!!  DATA wdensmin/1.e-12/
93  wdensmin=1.e-14
94  ! cc nrlmd
95  sigmaw_max=0.4
96  dens_rate=0.1
97
98  eps = rd_in/rv_in
99  epsim1 = 1.0/eps - 1.0
100  RG=rg_in
101  RD=rd_in
102
103
104  ! cc
105  ! Longueur de maille (en m)
106  ! -------------------------------------------------------------------------
107
108  ! ALON = 3.e5
109  ! alon = 1.E6
110
111  ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei)
112
113  ! coefgw : Coefficient pour les ondes de gravité
114  ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE)
115  ! wdens : Densité surfacique de poche froide
116  ! -------------------------------------------------------------------------
117
118  ! cc nrlmd      coefgw=10
119  ! coefgw=1
120  ! wdens0 = 1.0/(alon**2)
121  ! cc nrlmd      wdens = 1.0/(alon**2)
122  ! cc nrlmd      stark = 0.50
123  ! CRtest
124  ! cc nrlmd      alpk=0.1
125  ! alpk = 1.0
126  ! alpk = 0.5
127  ! alpk = 0.05
128
129
130
131  crep_upper = 0.9
132  crep_sol = 1.0
133
134  ! Get wapecut from parameter file
135  wapecut = 1.
136
137print*,'wapecut',wapecut
138  CALL getin_p('wapecut', wapecut)
139print*,'wapecut',wapecut
140
141  ! cc nrlmd Lecture du fichier wake_param.data
142
143
144  ! cc nrlmd Lecture du fichier wake_param.data
145  stark=0.33
146  CALL getin_p('stark',stark)
147  cstart = stark*sqrt(2.*wapecut)
148
149  alpk=0.25
150  CALL getin_p('alpk',alpk)
151!jyg<
152!!  wdens_ref=8.E-12
153!!  CALL getin_p('wdens_ref',wdens_ref)
154  wdens_ref(1)=8.E-12
155  wdens_ref(2)=8.E-12
156  CALL getin_p('wdens_ref_o',wdens_ref(1))    !wake number per unit area ; ocean
157  CALL getin_p('wdens_ref_l',wdens_ref(2))    !wake number per unit area ; land
158!>jyg
159!
160!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
161!!!!!!!!!  Population dynamics parameters    !!!!!!!!!!!!!!!!!!!!!!!!!!!!
162!------------------------------------------------------------------------
163
164  iflag_wk_pop_dyn = 0
165  CALL getin_p('iflag_wk_pop_dyn',iflag_wk_pop_dyn) ! switch between wdens prescribed
166                                                    ! and wdens prognostic
167  iflag_wk_act = 0
168  CALL getin_p('iflag_wk_act',iflag_wk_act) ! 0: act(:)=0.
169                                            ! 1: act(:)=1.
170                                            ! 2: act(:)=f(Wape)
171
172  rzero = 5000.
173  CALL getin_p('rzero_wk', rzero)
174  aa0 = 3.14*rzero*rzero
175!
176  tau_cv = 4000.
177  CALL getin_p('tau_cv', tau_cv)
178
179!------------------------------------------------------------------------
180
181  coefgw=4.
182  CALL getin_p('coefgw',coefgw)
183
184  WRITE(*,*) 'stark=', stark
185  WRITE(*,*) 'alpk=', alpk
186!jyg<
187!!  WRITE(*,*) 'wdens_ref=', wdens_ref
188  WRITE(*,*) 'wdens_ref_o=', wdens_ref(1)
189  WRITE(*,*) 'wdens_ref_l=', wdens_ref(2)
190!>jyg
191  WRITE(*,*) 'iflag_wk_pop_dyn=',iflag_wk_pop_dyn
192  WRITE(*,*) 'iflag_wk_act',iflag_wk_act
193  WRITE(*,*) 'coefgw=', coefgw
194
195  flag_wk_check_trgl=.false.
196  CALL getin_p('flag_wk_check_trgl ', flag_wk_check_trgl)
197  WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
198  WRITE(*,*) 'flag_wk_check_trgl OBSOLETE. Utilisr iflag_wk_check_trgl plutot'
199  iflag_wk_check_trgl=0 ; IF (flag_wk_check_trgl) iflag_wk_check_trgl=1
200  CALL getin_p('iflag_wk_check_trgl ', iflag_wk_check_trgl)
201  WRITE(*,*) 'iflag_wk_check_trgl=', iflag_wk_check_trgl
202
203 RETURN
204
205END SUBROUTINE wake_ini
206
207END MODULE wake_ini_mod
Note: See TracBrowser for help on using the repository browser.