source: LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90 @ 3150

Last change on this file since 3150 was 3150, checked in by jyg, 6 years ago

Introducing upper bounds on the convective
tendencies above which convection scheme is called
every time step.
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 28.0 KB
Line 
1!
2! $Id: phys_state_var_mod.F90 3150 2018-01-15 17:25:50Z jyg $
3!
4      MODULE phys_state_var_mod
5! Variables sauvegardees pour le startphy.nc
6!======================================================================
7!
8!
9!======================================================================
10! Declaration des variables
11      USE dimphy
12      USE netcdf, only: nf90_fill_real
13      INTEGER, PARAMETER :: nlevSTD=17
14      INTEGER, PARAMETER :: nlevSTD8=8
15      INTEGER, PARAMETER :: nlevSTD3=3
16      INTEGER, PARAMETER :: nout=3
17      INTEGER, PARAMETER :: napisccp=1
18      INTEGER, SAVE :: radpas  ! radiation is called every "radpas" step
19      INTEGER, SAVE :: cvpas   ! convection is called every "cvpas" step
20      INTEGER, SAVE :: cvpas_0 ! reference value for cvpas
21      INTEGER, SAVE :: wkpas   ! wake scheme is called every "wkpas" step
22      REAL, PARAMETER :: missing_val_nf90=nf90_fill_real
23!$OMP THREADPRIVATE(radpas)
24!$OMP THREADPRIVATE(cvpas)
25!$OMP THREADPRIVATE(cvpas_0)
26!$OMP THREADPRIVATE(wkpas)
27      REAL, SAVE :: dtime, solaire_etat0
28!$OMP THREADPRIVATE(dtime, solaire_etat0)
29
30      REAL, ALLOCATABLE, SAVE :: pctsrf(:,:)
31!$OMP THREADPRIVATE(pctsrf)
32      REAL, ALLOCATABLE, SAVE :: ftsol(:,:)
33!$OMP THREADPRIVATE(ftsol)
34      REAL,ALLOCATABLE,SAVE :: qsol(:),fevap(:,:),z0m(:,:),z0h(:,:),agesno(:,:)
35!$OMP THREADPRIVATE(qsol,fevap,z0m,z0h,agesno)
36!FC drag des arbres
37      REAL, ALLOCATABLE, SAVE :: treedrg(:,:,:)
38!$OMP THREADPRIVATE(treedrg)
39
40!      character(len=6), SAVE :: ocean
41!!!!!!$OMP THREADPRIVATE(ocean)
42!      logical, SAVE :: ok_veget
43!!!!!!$OMP THREADPRIVATE(ok_veget)
44      REAL, ALLOCATABLE, SAVE :: falb1(:,:), falb2(:,:)
45!$OMP THREADPRIVATE(falb1, falb2)
46
47!albedo SB >>>
48      REAL, ALLOCATABLE, SAVE :: falb_dif(:,:,:), falb_dir(:,:,:)
49      real, allocatable, save :: chl_con(:)
50!$OMP THREADPRIVATE(falb_dir,falb_dif,chl_con)
51!albedo SB <<<
52
53
54      REAL, ALLOCATABLE, SAVE :: rain_fall(:), snow_fall(:)
55!$OMP THREADPRIVATE( rain_fall, snow_fall)
56      REAL, ALLOCATABLE, SAVE :: solsw(:), sollw(:)
57!$OMP THREADPRIVATE(solsw, sollw)
58      REAL, ALLOCATABLE, SAVE :: radsol(:)
59!$OMP THREADPRIVATE(radsol)
60      REAL, ALLOCATABLE, SAVE :: swradcorr(:)
61!$OMP THREADPRIVATE(swradcorr)
62
63!clesphy0 param physiq
64!
65! Parametres de l'Orographie a l'Echelle Sous-Maille (OESM):
66!
67      REAL, ALLOCATABLE, SAVE :: zmea(:), zstd(:), zsig(:), zgam(:)
68!$OMP THREADPRIVATE(zmea, zstd, zsig, zgam)
69      REAL, ALLOCATABLE, SAVE :: zthe(:), zpic(:), zval(:)
70!$OMP THREADPRIVATE(zthe, zpic, zval)
71!     REAL tabcntr0(100)
72      REAL, ALLOCATABLE, SAVE :: rugoro(:)
73!$OMP THREADPRIVATE(rugoro)
74      REAL, ALLOCATABLE, SAVE :: t_ancien(:,:), q_ancien(:,:)
75!$OMP THREADPRIVATE(t_ancien, q_ancien)
76      REAL, ALLOCATABLE, SAVE :: ql_ancien(:,:), qs_ancien(:,:)
77!$OMP THREADPRIVATE(ql_ancien, qs_ancien)
78      REAL, ALLOCATABLE, SAVE :: prw_ancien(:), prlw_ancien(:), prsw_ancien(:)
79!$OMP THREADPRIVATE(prw_ancien, prlw_ancien, prsw_ancien)
80      REAL, ALLOCATABLE, SAVE :: u_ancien(:,:), v_ancien(:,:)
81!$OMP THREADPRIVATE(u_ancien, v_ancien)
82!!! RomP >>>
83      REAL, ALLOCATABLE, SAVE :: tr_ancien(:,:,:)
84!$OMP THREADPRIVATE(tr_ancien)
85!!! RomP <<<
86      LOGICAL, SAVE :: ancien_ok
87!$OMP THREADPRIVATE(ancien_ok)
88      REAL, ALLOCATABLE, SAVE :: clwcon(:,:),rnebcon(:,:)
89!$OMP THREADPRIVATE(clwcon,rnebcon)
90      REAL, ALLOCATABLE, SAVE :: qtc_cv(:,:),sigt_cv(:,:)
91!$OMP THREADPRIVATE(qtc_cv,sigt_cv)
92      REAL, ALLOCATABLE, SAVE :: ratqs(:,:)
93!$OMP THREADPRIVATE(ratqs)
94      REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy
95      REAL, ALLOCATABLE, SAVE :: coefh(:,:,:) ! Kz enthalpie
96      REAL, ALLOCATABLE, SAVE :: coefm(:,:,:) ! Kz momentum
97!$OMP THREADPRIVATE(pbl_tke, coefh,coefm)
98!nrlmd<
99      REAL, ALLOCATABLE, SAVE :: delta_tsurf(:,:) ! Surface temperature difference inside-outside cold pool
100!$OMP THREADPRIVATE(delta_tsurf)
101!>nrlmd
102      REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) !
103!$OMP THREADPRIVATE(zmax0,f0)
104      REAL, ALLOCATABLE, SAVE :: sig1(:,:), w01(:,:)
105!$OMP THREADPRIVATE(sig1,w01)
106      REAL, ALLOCATABLE, SAVE :: entr_therm(:,:), fm_therm(:,:)
107!$OMP THREADPRIVATE(entr_therm,fm_therm)
108      REAL, ALLOCATABLE, SAVE :: detr_therm(:,:)
109!$OMP THREADPRIVATE(detr_therm)
110!IM 150408
111!     pour phsystoke avec thermiques
112      REAL,ALLOCATABLE,SAVE :: clwcon0th(:,:),rnebcon0th(:,:)
113!$OMP THREADPRIVATE(clwcon0th,rnebcon0th)
114! radiation outputs
115      REAL,ALLOCATABLE,SAVE :: swdnc0(:,:), swdn0(:,:), swdn(:,:)
116!$OMP THREADPRIVATE(swdnc0,swdn0,swdn)
117      REAL,ALLOCATABLE,SAVE :: swupc0(:,:), swup0(:,:), swup(:,:)
118!$OMP THREADPRIVATE(swupc0, swup0,swup)
119      REAL,ALLOCATABLE,SAVE :: SWdn200clr(:), SWdn200(:)
120!$OMP THREADPRIVATE(SWdn200clr,SWdn200)
121      REAL,ALLOCATABLE,SAVE :: SWup200clr(:), SWup200(:)
122!$OMP THREADPRIVATE(SWup200clr,SWup200)
123      REAL,ALLOCATABLE,SAVE :: lwdnc0(:,:), lwdn0(:,:), lwdn(:,:)
124!$OMP THREADPRIVATE(lwdnc0,lwdn0,lwdn)
125      REAL,ALLOCATABLE,SAVE :: lwupc0(:,:), lwup0(:,:), lwup(:,:)
126!$OMP THREADPRIVATE(lwupc0,lwup0,lwup)
127      REAL,ALLOCATABLE,SAVE :: LWdn200clr(:), LWdn200(:)
128!$OMP THREADPRIVATE(LWdn200clr,LWdn200)
129      REAL,ALLOCATABLE,SAVE :: LWup200clr(:), LWup200(:)
130!$OMP THREADPRIVATE(LWup200clr,LWup200)
131      REAL,ALLOCATABLE,SAVE :: LWdnTOA(:), LWdnTOAclr(:)
132!$OMP THREADPRIVATE(LWdnTOA,LWdnTOAclr)
133! pressure level
134      REAL,ALLOCATABLE,SAVE :: tsumSTD(:,:,:)
135!$OMP THREADPRIVATE(tsumSTD)
136      REAL,ALLOCATABLE,SAVE :: usumSTD(:,:,:), vsumSTD(:,:,:)
137!$OMP THREADPRIVATE(usumSTD,vsumSTD)
138      REAL,ALLOCATABLE,SAVE :: wsumSTD(:,:,:), phisumSTD(:,:,:)
139!$OMP THREADPRIVATE(wsumSTD,phisumSTD)
140      REAL,ALLOCATABLE,SAVE :: qsumSTD(:,:,:), rhsumSTD(:,:,:)
141!$OMP THREADPRIVATE(qsumSTD,rhsumSTD)
142      REAL,ALLOCATABLE,SAVE :: tnondef(:,:,:)
143!$OMP THREADPRIVATE(tnondef)
144      REAL,ALLOCATABLE,SAVE :: uvsumSTD(:,:,:)
145!$OMP THREADPRIVATE(uvsumSTD)
146      REAL,ALLOCATABLE,SAVE :: vqsumSTD(:,:,:)
147!$OMP THREADPRIVATE(vqsumSTD)
148      REAL,ALLOCATABLE,SAVE :: vTsumSTD(:,:,:)
149!$OMP THREADPRIVATE(vTsumSTD)
150      REAL,ALLOCATABLE,SAVE :: wqsumSTD(:,:,:)
151!$OMP THREADPRIVATE(wqsumSTD)
152      REAL,ALLOCATABLE,SAVE :: vphisumSTD(:,:,:)
153!$OMP THREADPRIVATE(vphisumSTD)
154      REAL,ALLOCATABLE,SAVE :: wTsumSTD(:,:,:)
155!$OMP THREADPRIVATE(wTsumSTD)
156      REAL,ALLOCATABLE,SAVE :: u2sumSTD(:,:,:)
157!$OMP THREADPRIVATE(u2sumSTD)
158      REAL,ALLOCATABLE,SAVE :: v2sumSTD(:,:,:)
159!$OMP THREADPRIVATE(v2sumSTD)
160      REAL,ALLOCATABLE,SAVE :: T2sumSTD(:,:,:)
161!$OMP THREADPRIVATE(T2sumSTD)
162      REAL,ALLOCATABLE,SAVE :: O3sumSTD(:,:,:), O3daysumSTD(:,:,:)
163!$OMP THREADPRIVATE(O3sumSTD,O3daysumSTD)
164!IM begin
165      REAL,ALLOCATABLE,SAVE :: wlevSTD(:,:), ulevSTD(:,:), vlevSTD(:,:)
166!$OMP THREADPRIVATE(wlevSTD,ulevSTD,vlevSTD)
167      REAL,ALLOCATABLE,SAVE :: tlevSTD(:,:), qlevSTD(:,:), rhlevSTD(:,:)
168!$OMP THREADPRIVATE(tlevSTD,qlevSTD,rhlevSTD)
169      REAL,ALLOCATABLE,SAVE :: philevSTD(:,:)
170!$OMP THREADPRIVATE(philevSTD)
171      REAL,ALLOCATABLE,SAVE :: uvSTD(:,:)
172!$OMP THREADPRIVATE(uvSTD)
173      REAL,ALLOCATABLE,SAVE :: vqSTD(:,:)
174!$OMP THREADPRIVATE(vqSTD)
175      REAL,ALLOCATABLE,SAVE :: vTSTD(:,:)
176!$OMP THREADPRIVATE(vTSTD)
177      REAL,ALLOCATABLE,SAVE :: wqSTD(:,:)
178!$OMP THREADPRIVATE(wqSTD)
179      REAL,ALLOCATABLE,SAVE :: vphiSTD(:,:)
180!$OMP THREADPRIVATE(vphiSTD)
181      REAL,ALLOCATABLE,SAVE :: wTSTD(:,:)
182!$OMP THREADPRIVATE(wTSTD)
183      REAL,ALLOCATABLE,SAVE :: u2STD(:,:)
184!$OMP THREADPRIVATE(u2STD)
185      REAL,ALLOCATABLE,SAVE :: v2STD(:,:)
186!$OMP THREADPRIVATE(v2STD)
187      REAL,ALLOCATABLE,SAVE :: T2STD(:,:)
188!$OMP THREADPRIVATE(T2STD)
189      REAL,ALLOCATABLE,SAVE :: O3STD(:,:), O3daySTD(:,:)
190!$OMP THREADPRIVATE(O3STD,O3daySTD)
191!IM end
192      INTEGER,ALLOCATABLE,SAVE :: seed_old(:,:)
193!$OMP THREADPRIVATE(seed_old)
194      REAL,ALLOCATABLE,SAVE :: zuthe(:),zvthe(:)
195!$OMP THREADPRIVATE(zuthe,zvthe)
196      REAL,ALLOCATABLE,SAVE :: alb_neig(:)
197!$OMP THREADPRIVATE(alb_neig)
198!cloud base mass flux
199      REAL,ALLOCATABLE,SAVE :: ema_cbmf(:)
200!$OMP THREADPRIVATE(ema_cbmf)
201!cloud base pressure & cloud top pressure
202      REAL,ALLOCATABLE,SAVE :: ema_pcb(:), ema_pct(:)
203!$OMP THREADPRIVATE(ema_pcb,ema_pct)
204      REAL,ALLOCATABLE,SAVE :: Ma(:,:)        ! undilute upward mass flux
205!$OMP THREADPRIVATE(Ma)
206      REAL,ALLOCATABLE,SAVE :: qcondc(:,:)    ! in-cld water content from convect
207!$OMP THREADPRIVATE(qcondc)
208      REAL,ALLOCATABLE,SAVE :: wd(:) ! sb
209!$OMP THREADPRIVATE(wd)
210      REAL,ALLOCATABLE,SAVE :: sigd(:)
211!$OMP THREADPRIVATE(sigd)
212!
213      REAL,ALLOCATABLE,SAVE :: cin(:)
214!$OMP THREADPRIVATE(cin)
215! ftd : convective heating due to unsaturated downdraughts
216      REAL,ALLOCATABLE,SAVE :: ftd(:,:)
217!$OMP THREADPRIVATE(ftd)
218! fqd : convective moistening due to unsaturated downdraughts
219      REAL,ALLOCATABLE,SAVE :: fqd(:,:)     
220!$OMP THREADPRIVATE(fqd)
221!34EK
222! -- Variables de controle de ALE et ALP
223!ALE : Energie disponible pour soulevement : utilisee par la
224!      convection d'Emanuel pour le declenchement et la regulation
225      REAL,ALLOCATABLE,SAVE :: ALE(:)
226!$OMP THREADPRIVATE(ALE)
227!ALP : Puissance  disponible pour soulevement
228      REAL,ALLOCATABLE,SAVE :: ALP(:)
229!$OMP THREADPRIVATE(ALP)
230!
231! nouvelles variables pour le couplage convection-couche limite
232      REAL,ALLOCATABLE,SAVE :: Ale_bl(:)
233!$OMP THREADPRIVATE(Ale_bl)
234      REAL,ALLOCATABLE,SAVE :: Alp_bl(:)
235!$OMP THREADPRIVATE(Alp_bl)
236      INTEGER,ALLOCATABLE,SAVE :: lalim_conv(:)
237!$OMP THREADPRIVATE(lalim_conv)
238      REAL,ALLOCATABLE,SAVE :: wght_th(:,:)
239!$OMP THREADPRIVATE(wght_th)
240      REAL,ALLOCATABLE,SAVE    :: ale_wake(:)
241!$OMP THREADPRIVATE(ale_wake)
242      REAL,ALLOCATABLE,SAVE    :: ale_bl_stat(:)
243!$OMP THREADPRIVATE(ale_bl_stat)
244!
245! variables de la wake
246! wake_deltat : ecart de temperature avec la zone non perturbee
247! wake_deltaq : ecart d'humidite avec la zone non perturbee
248! wake_s      : fraction surfacique occupee par la poche froide
249! wake_dens   : number of wakes per unit area
250! wake_occ    : occurence of wakes (= 1 if wakes occur, =0 otherwise)
251! wake_Cstar  : vitesse d'etalement de la poche
252! wake_pe     : wake potential energy - WAPE
253! wake_fip    : Gust Front Impinging power - ALP
254      REAL,ALLOCATABLE,SAVE :: wake_deltat(:,:)
255!$OMP THREADPRIVATE(wake_deltat)
256      REAL,ALLOCATABLE,SAVE :: wake_deltaq(:,:)
257!$OMP THREADPRIVATE(wake_deltaq)
258      REAL,ALLOCATABLE,SAVE :: wake_s(:)
259!$OMP THREADPRIVATE(wake_s)
260      REAL,ALLOCATABLE,SAVE :: wake_dens(:)
261!$OMP THREADPRIVATE(wake_dens)
262      REAL,ALLOCATABLE,SAVE :: wake_Cstar(:)
263!$OMP THREADPRIVATE(wake_Cstar)
264      REAL,ALLOCATABLE,SAVE :: wake_pe(:)
265!$OMP THREADPRIVATE(wake_pe)
266      REAL,ALLOCATABLE,SAVE :: wake_fip(:)
267!$OMP THREADPRIVATE(wake_fip)
268!
269!jyg<
270! variables related to the spitting of the PBL between wake and
271! off-wake regions.
272! wake_delta_pbl_TKE : difference TKE_w - TKE_x
273      REAL,ALLOCATABLE,SAVE :: wake_delta_pbl_TKE(:,:,:)
274!$OMP THREADPRIVATE(wake_delta_pbl_TKE)
275!>jyg
276!
277! pfrac_impa : Produits des coefs lessivage impaction
278! pfrac_nucl : Produits des coefs lessivage nucleation
279! pfrac_1nucl: Produits des coefs lessi nucl (alpha = 1)
280      REAL,ALLOCATABLE,SAVE :: pfrac_impa(:,:), pfrac_nucl(:,:)
281!$OMP THREADPRIVATE(pfrac_impa,pfrac_nucl)
282      REAL,ALLOCATABLE,SAVE :: pfrac_1nucl(:,:)
283!$OMP THREADPRIVATE(pfrac_1nucl)
284!
285      REAL,ALLOCATABLE,SAVE :: total_rain(:), nday_rain(:) 
286!$OMP THREADPRIVATE(total_rain,nday_rain)
287! albsol1: albedo du sol total pour SW visible
288! albsol2: albedo du sol total pour SW proche IR
289      REAL,ALLOCATABLE,SAVE :: albsol1(:), albsol2(:)
290!$OMP THREADPRIVATE(albsol1,albsol2)
291
292!albedo SB >>>
293      REAL,ALLOCATABLE,SAVE :: albsol_dif(:,:),albsol_dir(:,:)
294!$OMP THREADPRIVATE(albsol_dif,albsol_dir)
295!albedo SB <<<
296
297
298      REAL, ALLOCATABLE, SAVE:: wo(:, :, :)
299      ! column-density of ozone in a layer, in kilo-Dobsons
300      ! Third dimension has size 1 or 2.
301      ! "wo(:, :, 1)" is for the average day-night field,
302      ! "wo(:, :, 2)" is for daylight time.
303      !$OMP THREADPRIVATE(wo)
304
305! heat : chauffage solaire
306! heat0: chauffage solaire ciel clair
307! cool : refroidissement infrarouge
308! cool0 : refroidissement infrarouge ciel clair
309! sollwdown : downward LW flux at surface
310! sollwdownclr : downward CS LW flux at surface
311! toplwdown : downward CS LW flux at TOA
312! toplwdownclr : downward CS LW flux at TOA
313      REAL,ALLOCATABLE,SAVE :: clwcon0(:,:),rnebcon0(:,:)
314!$OMP THREADPRIVATE(clwcon0,rnebcon0)
315      REAL,ALLOCATABLE,SAVE :: heat(:,:)   
316!$OMP THREADPRIVATE(heat)
317      REAL,ALLOCATABLE,SAVE :: heat0(:,:)
318!$OMP THREADPRIVATE(heat0)
319      REAL,ALLOCATABLE,SAVE :: cool(:,:)
320!$OMP THREADPRIVATE(cool)
321      REAL,ALLOCATABLE,SAVE :: cool0(:,:)
322!$OMP THREADPRIVATE(cool0)
323      REAL,ALLOCATABLE,SAVE :: topsw(:), toplw(:)
324!$OMP THREADPRIVATE(topsw,toplw)
325      REAL,ALLOCATABLE,SAVE :: sollwdown(:)
326!$OMP THREADPRIVATE(sollwdown)
327      REAL,ALLOCATABLE,SAVE :: gustiness(:)
328!$OMP THREADPRIVATE(gustiness)
329      REAL,ALLOCATABLE,SAVE :: sollwdownclr(:)
330!$OMP THREADPRIVATE(sollwdownclr)
331      REAL,ALLOCATABLE,SAVE :: toplwdown(:)
332!$OMP THREADPRIVATE(toplwdown)
333      REAL,ALLOCATABLE,SAVE :: toplwdownclr(:)
334!$OMP THREADPRIVATE(toplwdownclr)
335      REAL,ALLOCATABLE,SAVE :: topsw0(:),toplw0(:),solsw0(:),sollw0(:)
336!$OMP THREADPRIVATE(topsw0,toplw0,solsw0,sollw0)
337      REAL,ALLOCATABLE,SAVE :: albpla(:)
338!$OMP THREADPRIVATE(albpla)
339
340!IM ajout variables CFMIP2/CMIP5
341      REAL,ALLOCATABLE,SAVE :: heatp(:,:), coolp(:,:)
342!$OMP THREADPRIVATE(heatp, coolp)
343      REAL,ALLOCATABLE,SAVE :: heat0p(:,:), cool0p(:,:)
344!$OMP THREADPRIVATE(heat0p, cool0p)
345      REAL,ALLOCATABLE,SAVE :: radsolp(:), topswp(:), toplwp(:)
346!$OMP THREADPRIVATE(radsolp, topswp, toplwp)
347      REAL,ALLOCATABLE,SAVE :: albplap(:)
348!$OMP THREADPRIVATE(albplap)
349      REAL,ALLOCATABLE,SAVE :: solswp(:), sollwp(:)
350!$OMP THREADPRIVATE(solswp, sollwp)
351      REAL,ALLOCATABLE,SAVE :: sollwdownp(:)
352!$OMP THREADPRIVATE(sollwdownp)
353      REAL,ALLOCATABLE,SAVE :: topsw0p(:),toplw0p(:)
354      REAL,ALLOCATABLE,SAVE :: solsw0p(:),sollw0p(:)
355!$OMP THREADPRIVATE(topsw0p,toplw0p,solsw0p,sollw0p)
356      REAL,ALLOCATABLE,SAVE :: lwdnc0p(:,:), lwdn0p(:,:), lwdnp(:,:)
357      REAL,ALLOCATABLE,SAVE :: lwupc0p(:,:), lwup0p(:,:), lwupp(:,:)
358!$OMP THREADPRIVATE(lwdnc0p, lwdn0p, lwdnp, lwupc0p, lwup0p, lwupp)
359      REAL,ALLOCATABLE,SAVE :: swdnc0p(:,:), swdn0p(:,:), swdnp(:,:)
360      REAL,ALLOCATABLE,SAVE :: swupc0p(:,:), swup0p(:,:), swupp(:,:)
361!$OMP THREADPRIVATE(swdnc0p, swdn0p, swdnp, swupc0p, swup0p, swupp)
362
363! pbase : cloud base pressure
364! bbase : cloud base buoyancy
365      REAL,ALLOCATABLE,SAVE :: cape(:)
366!$OMP THREADPRIVATE(cape)
367      REAL,ALLOCATABLE,SAVE :: pbase(:)
368!$OMP THREADPRIVATE(pbase)
369      REAL,ALLOCATABLE,SAVE :: bbase(:)
370!$OMP THREADPRIVATE(bbase)
371!
372      REAL,SAVE,ALLOCATABLE :: zqasc(:,:)
373!$OMP THREADPRIVATE( zqasc)
374      INTEGER,ALLOCATABLE,SAVE :: ibas_con(:), itop_con(:)
375!$OMP THREADPRIVATE(ibas_con,itop_con)
376      REAL,SAVE,ALLOCATABLE :: rain_con(:)
377!$OMP THREADPRIVATE(rain_con)
378      REAL,SAVE,ALLOCATABLE :: snow_con(:)
379!$OMP THREADPRIVATE(snow_con)
380!
381      REAL,SAVE,ALLOCATABLE :: rlonPOS(:)
382!$OMP THREADPRIVATE(rlonPOS)
383      REAL,SAVE,ALLOCATABLE :: newsst(:)
384!$OMP THREADPRIVATE(newsst)
385      REAL,SAVE,ALLOCATABLE :: ustar(:,:),u10m(:,:), v10m(:,:),wstar(:,:)
386!$OMP THREADPRIVATE(ustar,u10m,v10m,wstar)
387!
388! ok_ade=T -ADE=topswad-topsw
389! ok_aie=T ->
390!       ok_ade=T -AIE=topswai-topswad
391!       ok_ade=F -AIE=topswai-topsw
392!
393!topswad, solswad : Aerosol direct effect
394      REAL,SAVE,ALLOCATABLE :: topswad(:), solswad(:)
395!$OMP THREADPRIVATE(topswad,solswad)
396!topswai, solswai : Aerosol indirect effect
397      REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:)
398!$OMP THREADPRIVATE(topswai,solswai)
399
400      REAL,SAVE,ALLOCATABLE :: tau_aero(:,:,:,:), piz_aero(:,:,:,:), cg_aero(:,:,:,:)
401!$OMP THREADPRIVATE(tau_aero, piz_aero, cg_aero)
402      REAL,SAVE,ALLOCATABLE :: tau_aero_sw_rrtm(:,:,:,:), piz_aero_sw_rrtm(:,:,:,:), cg_aero_sw_rrtm(:,:,:,:)
403!$OMP THREADPRIVATE(tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm)
404      REAL,SAVE,ALLOCATABLE :: tau_aero_lw_rrtm(:,:,:,:), piz_aero_lw_rrtm(:,:,:,:), cg_aero_lw_rrtm(:,:,:,:)
405!$OMP THREADPRIVATE(tau_aero_lw_rrtm, piz_aero_lw_rrtm, cg_aero_lw_rrtm)
406      REAL,SAVE,ALLOCATABLE :: ccm(:,:,:)
407!$OMP THREADPRIVATE(ccm)
408
409!!! nrlmd le 10/04/2012
410      REAL,SAVE,ALLOCATABLE :: ale_bl_trig(:)
411!$OMP THREADPRIVATE(ale_bl_trig)
412!!! fin nrlmd le 10/04/2012
413
414      REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), du_gwd_front(:, :)
415      !$OMP THREADPRIVATE(du_gwd_rando, du_gwd_front)
416      ! tendencies on wind due to gravity waves
417
418CONTAINS
419
420!======================================================================
421SUBROUTINE phys_state_var_init(read_climoz)
422USE dimphy
423USE aero_mod
424USE infotrac_phy, ONLY : nbtr
425USE indice_sol_mod
426IMPLICIT NONE
427
428integer, intent(in)::  read_climoz
429! read ozone climatology
430! Allowed values are 0, 1 and 2
431! 0: do not read an ozone climatology
432! 1: read a single ozone climatology that will be used day and night
433! 2: read two ozone climatologies, the average day and night
434! climatology and the daylight climatology
435
436include "clesphys.h"
437
438      ALLOCATE(pctsrf(klon,nbsrf))
439      ALLOCATE(ftsol(klon,nbsrf))
440      ALLOCATE(qsol(klon),fevap(klon,nbsrf))
441      ALLOCATE(z0m(klon,nbsrf+1),z0h(klon,nbsrf+1),agesno(klon,nbsrf))
442!FC
443      ALLOCATE(treedrg(klon,klev,nbsrf))
444      ALLOCATE(falb1(klon,nbsrf))
445      ALLOCATE(falb2(klon,nbsrf))
446!albedo SB >>>
447      ALLOCATE(falb_dir(klon,nsw,nbsrf),falb_dif(klon,nsw,nbsrf))
448      ALLOCATE(chl_con(klon))
449!albedo SB <<<
450      ALLOCATE(rain_fall(klon))
451      ALLOCATE(snow_fall(klon))
452      ALLOCATE(solsw(klon), sollw(klon))
453      ALLOCATE(radsol(klon))
454      ALLOCATE(swradcorr(klon))
455      ALLOCATE(zmea(klon), zstd(klon), zsig(klon), zgam(klon))
456      ALLOCATE(zthe(klon), zpic(klon), zval(klon))
457
458      ALLOCATE(rugoro(klon))
459      ALLOCATE(t_ancien(klon,klev), q_ancien(klon,klev))
460      ALLOCATE(ql_ancien(klon,klev), qs_ancien(klon,klev))
461      ALLOCATE(prw_ancien(klon), prlw_ancien(klon), prsw_ancien(klon))
462      ALLOCATE(u_ancien(klon,klev), v_ancien(klon,klev))
463!!! Rom P >>>
464      ALLOCATE(tr_ancien(klon,klev,nbtr))
465!!! Rom P <<<
466      ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev))
467      ALLOCATE(qtc_cv(klon,klev),sigt_cv(klon,klev))
468      ALLOCATE(ratqs(klon,klev))
469      ALLOCATE(pbl_tke(klon,klev+1,nbsrf+1))
470!nrlmd<
471      ALLOCATE(delta_tsurf(klon,nbsrf))
472!>nrlmd
473      ALLOCATE(coefh(klon,klev+1,nbsrf+1))
474      ALLOCATE(coefm(klon,klev+1,nbsrf+1))
475      ALLOCATE(zmax0(klon), f0(klon))
476      ALLOCATE(sig1(klon,klev), w01(klon,klev))
477      ALLOCATE(entr_therm(klon,klev), fm_therm(klon,klev+1))
478      ALLOCATE(detr_therm(klon,klev))
479!     pour phsystoke avec thermiques
480      ALLOCATE(clwcon0th(klon,klev),rnebcon0th(klon,klev))
481! radiation outputs
482      ALLOCATE(swdnc0(klon,klevp1), swdn0(klon,klevp1), swdn(klon,klevp1))
483      ALLOCATE(swupc0(klon,klevp1), swup0(klon,klevp1), swup(klon,klevp1))
484      ALLOCATE(lwdnc0(klon,klevp1), lwdn0(klon,klevp1), lwdn(klon,klevp1))
485      ALLOCATE(lwupc0(klon,klevp1), lwup0(klon,klevp1), lwup(klon,klevp1))
486      ALLOCATE(SWdn200clr(klon), SWdn200(klon))
487      ALLOCATE(SWup200clr(klon), SWup200(klon))
488      ALLOCATE(LWdn200clr(klon), LWdn200(klon))
489      ALLOCATE(LWup200clr(klon), LWup200(klon))
490      ALLOCATE(LWdnTOA(klon), LWdnTOAclr(klon))
491! pressure level
492      ALLOCATE(tsumSTD(klon,nlevSTD,nout))
493      ALLOCATE(usumSTD(klon,nlevSTD,nout), vsumSTD(klon,nlevSTD,nout))
494      ALLOCATE(wsumSTD(klon,nlevSTD,nout), phisumSTD(klon,nlevSTD,nout))
495      ALLOCATE(qsumSTD(klon,nlevSTD,nout), rhsumSTD(klon,nlevSTD,nout))
496      ALLOCATE(tnondef(klon,nlevSTD,nout))
497      ALLOCATE(uvsumSTD(klon,nlevSTD,nout))
498      ALLOCATE(vqsumSTD(klon,nlevSTD,nout))
499      ALLOCATE(vTsumSTD(klon,nlevSTD,nout))
500      ALLOCATE(wqsumSTD(klon,nlevSTD,nout))
501      ALLOCATE(vphisumSTD(klon,nlevSTD,nout))
502      ALLOCATE(wTsumSTD(klon,nlevSTD,nout))
503      ALLOCATE(u2sumSTD(klon,nlevSTD,nout))
504      ALLOCATE(v2sumSTD(klon,nlevSTD,nout))
505      ALLOCATE(T2sumSTD(klon,nlevSTD,nout))
506      ALLOCATE(O3sumSTD(klon,nlevSTD,nout))
507      ALLOCATE(O3daysumSTD(klon,nlevSTD,nout))
508!IM beg
509      ALLOCATE(wlevSTD(klon,nlevSTD), ulevSTD(klon,nlevSTD), vlevSTD(klon,nlevSTD))
510      ALLOCATE(tlevSTD(klon,nlevSTD), qlevSTD(klon,nlevSTD), rhlevSTD(klon,nlevSTD))
511      ALLOCATE(philevSTD(klon,nlevSTD))
512      ALLOCATE(uvSTD(klon,nlevSTD),vqSTD(klon,nlevSTD))
513      ALLOCATE(vTSTD(klon,nlevSTD),wqSTD(klon,nlevSTD))
514      ALLOCATE(vphiSTD(klon,nlevSTD),wTSTD(klon,nlevSTD))
515      ALLOCATE(u2STD(klon,nlevSTD),v2STD(klon,nlevSTD))
516      ALLOCATE(T2STD(klon,nlevSTD))
517      ALLOCATE(O3STD(klon,nlevSTD))
518      ALLOCATE(O3daySTD(klon,nlevSTD))
519!IM end
520      ALLOCATE(seed_old(klon,napisccp))
521      ALLOCATE(zuthe(klon),zvthe(klon))
522      ALLOCATE(alb_neig(klon))
523!cloud base mass flux
524      ALLOCATE(ema_cbmf(klon))
525!cloud base pressure & cloud top pressure
526      ALLOCATE(ema_pcb(klon), ema_pct(klon))
527!
528      ALLOCATE(Ma(klon,klev))
529      ALLOCATE(qcondc(klon,klev))
530      ALLOCATE(wd(klon))
531      ALLOCATE(sigd(klon))
532      ALLOCATE(cin(klon), ALE(klon), ALP(klon))
533      ALLOCATE(ftd(klon,klev), fqd(klon,klev))
534      ALLOCATE(Ale_bl(klon))
535      ALLOCATE(ale_wake(klon))
536      ALLOCATE(ale_bl_stat(klon))
537      ALLOCATE(Alp_bl(klon))
538      ALLOCATE(lalim_conv(klon))
539      ALLOCATE(wght_th(klon,klev))
540      ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev))
541      ALLOCATE(wake_s(klon), wake_dens(klon))
542      ALLOCATE(wake_Cstar(klon))
543      ALLOCATE(wake_pe(klon), wake_fip(klon))
544!jyg<
545      ALLOCATE(wake_delta_pbl_TKE(klon,klev+1,nbsrf+1))
546!>jyg
547      ALLOCATE(pfrac_impa(klon,klev), pfrac_nucl(klon,klev))
548      ALLOCATE(pfrac_1nucl(klon,klev))
549      ALLOCATE(total_rain(klon), nday_rain(klon))
550      ALLOCATE(albsol1(klon), albsol2(klon))
551!albedo SB >>>
552      ALLOCATE(albsol_dir(klon,nsw),albsol_dif(klon,nsw))
553!albedo SB <<<
554
555      if (read_climoz <= 1) then
556         ALLOCATE(wo(klon,klev, 1))
557      else
558         ! read_climoz == 2
559         ALLOCATE(wo(klon,klev, 2))
560      end if
561     
562      ALLOCATE(clwcon0(klon,klev),rnebcon0(klon,klev))
563      ALLOCATE(heat(klon,klev), heat0(klon,klev))
564      ALLOCATE(cool(klon,klev), cool0(klon,klev))
565      ALLOCATE(topsw(klon), toplw(klon))
566      ALLOCATE(sollwdown(klon), sollwdownclr(klon))
567      ALLOCATE(toplwdown(klon), toplwdownclr(klon))
568      ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon))
569      ALLOCATE(albpla(klon))
570!IM ajout variables CFMIP2/CMIP5
571      ALLOCATE(heatp(klon,klev), coolp(klon,klev))
572      ALLOCATE(heat0p(klon,klev), cool0p(klon,klev))
573      ALLOCATE(radsolp(klon), topswp(klon), toplwp(klon))
574      ALLOCATE(albplap(klon))
575      ALLOCATE(solswp(klon), sollwp(klon))
576      ALLOCATE(gustiness(klon))
577      ALLOCATE(sollwdownp(klon))
578      ALLOCATE(topsw0p(klon),toplw0p(klon))
579      ALLOCATE(solsw0p(klon),sollw0p(klon))
580      ALLOCATE(lwdnc0p(klon,klevp1), lwdn0p(klon,klevp1), lwdnp(klon,klevp1))
581      ALLOCATE(lwupc0p(klon,klevp1), lwup0p(klon,klevp1), lwupp(klon,klevp1))
582      ALLOCATE(swdnc0p(klon,klevp1), swdn0p(klon,klevp1), swdnp(klon,klevp1))
583      ALLOCATE(swupc0p(klon,klevp1), swup0p(klon,klevp1), swupp(klon,klevp1))
584
585      ALLOCATE(cape(klon))
586      ALLOCATE(pbase(klon),bbase(klon))
587      ALLOCATE(zqasc(klon,klev))
588      ALLOCATE(ibas_con(klon), itop_con(klon))
589      ALLOCATE(rain_con(klon), snow_con(klon))
590      ALLOCATE(rlonPOS(klon))
591      ALLOCATE(newsst(klon))
592      ALLOCATE(ustar(klon,nbsrf),u10m(klon,nbsrf), v10m(klon,nbsrf),wstar(klon,nbsrf+1))
593      ALLOCATE(topswad(klon), solswad(klon))
594      ALLOCATE(topswai(klon), solswai(klon))
595      ALLOCATE(tau_aero(klon,klev,naero_grp,nbands),piz_aero(klon,klev,naero_grp,nbands),cg_aero(klon,klev,naero_grp,nbands))
596      ALLOCATE(tau_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm),piz_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm))
597      ALLOCATE(cg_aero_sw_rrtm(klon,klev,2,nbands_sw_rrtm))
598      ALLOCATE(tau_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm),piz_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm))
599      ALLOCATE(cg_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm))
600      ALLOCATE(ccm(klon,klev,nbands))
601
602!!! nrlmd le 10/04/2012
603      ALLOCATE(ale_bl_trig(klon))
604!!! fin nrlmd le 10/04/2012
605      if (ok_gwd_rando) allocate(du_gwd_rando(klon, klev))
606      if (.not. ok_hines .and. ok_gwd_rando) allocate(du_gwd_front(klon, klev))
607
608END SUBROUTINE phys_state_var_init
609
610!======================================================================
611SUBROUTINE phys_state_var_end
612!USE dimphy
613USE indice_sol_mod
614IMPLICIT NONE
615include "clesphys.h"
616
617      deallocate(pctsrf, ftsol, falb1, falb2)
618      deallocate(qsol,fevap,z0m,z0h,agesno)
619!FC
620      deallocate(treedrg)
621      deallocate(rain_fall, snow_fall, solsw, sollw, radsol, swradcorr)
622      deallocate(zmea, zstd, zsig, zgam)
623      deallocate(zthe, zpic, zval)
624      deallocate(rugoro, t_ancien, q_ancien, clwcon, rnebcon)
625      deallocate(qs_ancien, ql_ancien)
626      deallocate(prw_ancien, prlw_ancien, prsw_ancien)
627      deallocate(qtc_cv,sigt_cv)
628      deallocate(u_ancien, v_ancien)
629      deallocate(tr_ancien)                           !RomP
630      deallocate(ratqs, pbl_tke,coefh,coefm)
631!nrlmd<
632      deallocate(delta_tsurf)
633!>nrlmd
634      deallocate(zmax0, f0)
635      deallocate(sig1, w01)
636      deallocate(entr_therm, fm_therm)
637      deallocate(detr_therm)
638      deallocate(clwcon0th, rnebcon0th)
639! radiation outputs
640      deallocate(swdnc0, swdn0, swdn)
641      deallocate(swupc0, swup0, swup)
642      deallocate(lwdnc0, lwdn0, lwdn)
643      deallocate(lwupc0, lwup0, lwup)
644      deallocate(SWdn200clr, SWdn200)
645      deallocate(SWup200clr, SWup200)
646      deallocate(LWdn200clr, LWdn200)
647      deallocate(LWup200clr, LWup200)
648      deallocate(LWdnTOA, LWdnTOAclr)
649! pressure level
650      deallocate(tsumSTD)
651      deallocate(usumSTD, vsumSTD)
652      deallocate(wsumSTD, phisumSTD)
653      deallocate(tnondef)
654      deallocate(qsumSTD, rhsumSTD)
655      deallocate(uvsumSTD)
656      deallocate(vqsumSTD)
657      deallocate(vTsumSTD)
658      deallocate(wqsumSTD)
659      deallocate(vphisumSTD)
660      deallocate(wTsumSTD)
661      deallocate(u2sumSTD)
662      deallocate(v2sumSTD)
663      deallocate(T2sumSTD)
664      deallocate(O3sumSTD)
665      deallocate(O3daysumSTD)
666!IM beg
667      deallocate(wlevSTD,ulevSTD,vlevSTD,tlevSTD,qlevSTD,rhlevSTD,philevSTD)
668      deallocate(uvSTD,vqSTD,vTSTD,wqSTD,vphiSTD,wTSTD,u2STD,v2STD,T2STD,O3STD,O3daySTD)
669!IM end
670      deallocate(seed_old)
671      deallocate(zuthe, zvthe)
672      deallocate(alb_neig)
673      deallocate(ema_cbmf)
674      deallocate(ema_pcb, ema_pct)
675      deallocate(Ma, qcondc)
676      deallocate(wd, sigd)
677      deallocate(cin, ALE, ALP)
678      deallocate(ftd, fqd)
679      deallocate(Ale_bl, Alp_bl)
680      deallocate(ale_wake)
681      deallocate(ale_bl_stat)
682      deallocate(lalim_conv, wght_th)
683      deallocate(wake_deltat, wake_deltaq)
684      deallocate(wake_s, wake_dens)
685      deallocate(wake_Cstar, wake_pe, wake_fip)
686!jyg<
687      deallocate(wake_delta_pbl_TKE)
688!>jyg
689      deallocate(pfrac_impa, pfrac_nucl)
690      deallocate(pfrac_1nucl)
691      deallocate(total_rain, nday_rain)
692      deallocate(albsol1, albsol2)
693!albedo SB >>>
694      deallocate(albsol_dir,albsol_dif,falb_dir,falb_dif,chl_con)
695!albedo SB <<<
696      deallocate(wo)
697      deallocate(clwcon0,rnebcon0)
698      deallocate(heat, heat0)
699      deallocate(cool, cool0)
700      deallocate(topsw, toplw)
701      deallocate(sollwdown, sollwdownclr)
702      deallocate(gustiness)
703      deallocate(toplwdown, toplwdownclr)
704      deallocate(topsw0,toplw0,solsw0,sollw0)
705      deallocate(albpla)
706!IM ajout variables CFMIP2/CMIP5
707      deallocate(heatp, coolp)
708      deallocate(heat0p, cool0p)
709      deallocate(radsolp, topswp, toplwp)
710      deallocate(albplap)
711      deallocate(solswp, sollwp)
712      deallocate(sollwdownp)
713      deallocate(topsw0p,toplw0p)
714      deallocate(solsw0p,sollw0p)
715      deallocate(lwdnc0p, lwdn0p, lwdnp)
716      deallocate(lwupc0p, lwup0p, lwupp)
717      deallocate(swdnc0p, swdn0p, swdnp)
718      deallocate(swupc0p, swup0p, swupp)
719      deallocate(cape)
720      deallocate(pbase,bbase)
721      deallocate(zqasc)
722      deallocate(ibas_con, itop_con)
723      deallocate(rain_con, snow_con)
724      deallocate(rlonPOS)
725      deallocate(newsst)
726      deallocate(ustar,u10m, v10m,wstar)
727      deallocate(topswad, solswad)
728      deallocate(topswai, solswai)
729      deallocate(tau_aero,piz_aero,cg_aero)
730      deallocate(tau_aero_sw_rrtm,piz_aero_sw_rrtm,cg_aero_sw_rrtm)
731      deallocate(tau_aero_lw_rrtm,piz_aero_lw_rrtm,cg_aero_lw_rrtm)
732      deallocate(ccm)
733      if (ok_gwd_rando) deallocate(du_gwd_rando)
734      if (.not. ok_hines .and. ok_gwd_rando) deallocate(du_gwd_front)
735       
736!!! nrlmd le 10/04/2012
737      deallocate(ale_bl_trig)
738!!! fin nrlmd le 10/04/2012
739
740END SUBROUTINE phys_state_var_end
741
742      END MODULE phys_state_var_mod
Note: See TracBrowser for help on using the repository browser.