source: LMDZ6/trunk/libf/phylmdiso/phys_state_var_mod.F90 @ 3927

Last change on this file since 3927 was 3927, checked in by Laurent Fairhead, 3 years ago

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