source: LMDZ5/branches/testing/libf/phylmd/phys_local_var_mod.F90 @ 1664

Last change on this file since 1664 was 1539, checked in by musat, 13 years ago

Ajouts CFMIP2/CMIP5

  • 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers (voir DefLists?): npCFMIP_param.data(_*) contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) pointlocations.txt(_*) contenant le numero, les coordonnees (lon,lat) et le nom de chaque station
  • flag LOGICAL dans tous les appels histwrite_phy pour pouvoir sortir le fichier histstn.nc

NB: 1) les flags de type phys_ que l'on met dans le physiq.def_L* pour ajouter plus de sorties

necessitent dorenavant 6 valeurs, la 6eme correspondant au fichier histstn.nc

2) par defaut le fichier histstn.nc ne sort pas; pour le sortir ajouter les lignes suivantes

dans physiq.def_L*

### Type de fichier : global (n) ou stations (y)
phys_out_filestations = n n n n n y

  • introduction de 2 jeux de flags pour les taux des GES; taux actuels avec suffixes _act, taux futurs avec "_per" avec 2 appels au rayonnement si taux "_per" different des taux "_act" (utiles pour diags. CFMIP 4CO2)
  • flags "betaCRF" pour calculs CRF pour experiences sensibilite proprietes optiques eau liquide nuageuse avec initialisations par defaut; sinon besoin de fichier beta_crf.data

Ajout flag LOGICAL lCOSP necessaire pour sortir un fichier stations
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 14.3 KB
Line 
1!
2! $Id: phys_local_var_mod.F90 1539 2011-06-08 22:13:33Z fairhead $
3!
4      MODULE phys_local_var_mod
5
6! Variables locales pour effectuer les appels en serie
7!======================================================================
8!
9!
10!======================================================================
11! Declaration des variables
12
13      REAL, SAVE, ALLOCATABLE :: t_seri(:,:), q_seri(:,:)
14      !$OMP THREADPRIVATE(t_seri, q_seri)
15      REAL, SAVE, ALLOCATABLE :: ql_seri(:,:),qs_seri(:,:)
16      !$OMP THREADPRIVATE(ql_seri,qs_seri)
17      REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:)
18      !$OMP THREADPRIVATE(u_seri, v_seri)
19
20      REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:)
21      !$OMP THREADPRIVATE(tr_seri)
22      REAL, SAVE, ALLOCATABLE :: d_t_dyn(:,:), d_q_dyn(:,:)
23      !$OMP THREADPRIVATE(d_t_dyn, d_q_dyn)
24      REAL, SAVE, ALLOCATABLE :: d_u_dyn(:,:), d_v_dyn(:,:)
25      !$OMP THREADPRIVATE(d_u_dyn, d_v_dyn)
26      REAL, SAVE, ALLOCATABLE :: d_t_con(:,:),d_q_con(:,:)
27      !$OMP THREADPRIVATE(d_t_con,d_q_con)
28      REAL, SAVE, ALLOCATABLE :: d_u_con(:,:),d_v_con(:,:)
29      !$OMP THREADPRIVATE(d_u_con,d_v_con)
30      REAL, SAVE, ALLOCATABLE :: d_t_wake(:,:),d_q_wake(:,:)
31      !$OMP THREADPRIVATE( d_t_wake,d_q_wake)
32      REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:)
33      !$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc)
34      REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:)
35      !$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb)
36      REAL, SAVE, ALLOCATABLE :: d_t_ajs(:,:), d_q_ajs(:,:)
37      !$OMP THREADPRIVATE(d_t_ajs, d_q_ajs)
38      REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:)
39      !$OMP THREADPRIVATE(d_u_ajs, d_v_ajs)
40      REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:)
41      !$OMP THREADPRIVATE(d_t_eva,d_q_eva)
42      REAL, SAVE, ALLOCATABLE :: d_t_lscst(:,:),d_q_lscst(:,:)
43      !$OMP THREADPRIVATE(d_t_lscst,d_q_lscst)
44      REAL, SAVE, ALLOCATABLE :: d_t_lscth(:,:),d_q_lscth(:,:)
45      !$OMP THREADPRIVATE(d_t_lscth,d_q_lscth)
46      REAL, SAVE, ALLOCATABLE :: plul_th(:),plul_st(:)
47      !$OMP THREADPRIVATE(plul_th,plul_st)
48!tendances dues a oro et lif
49      REAL, SAVE, ALLOCATABLE :: d_t_oli(:,:)
50      !$OMP THREADPRIVATE(d_t_oli)
51      REAL, SAVE, ALLOCATABLE :: d_u_oli(:,:), d_v_oli(:,:)
52      !$OMP THREADPRIVATE(d_u_oli, d_v_oli)
53      REAL, SAVE, ALLOCATABLE :: d_t_vdf(:,:), d_q_vdf(:,:)
54      !$OMP THREADPRIVATE( d_t_vdf, d_q_vdf)
55      REAL, SAVE, ALLOCATABLE :: d_u_vdf(:,:), d_v_vdf(:,:)
56      !$OMP THREADPRIVATE(d_u_vdf, d_v_vdf)
57      REAL, SAVE, ALLOCATABLE :: d_t_oro(:,:)
58      !$OMP THREADPRIVATE(d_t_oro)
59      REAL, SAVE, ALLOCATABLE :: d_u_oro(:,:), d_v_oro(:,:)
60      !$OMP THREADPRIVATE(d_u_oro, d_v_oro)
61      REAL, SAVE, ALLOCATABLE :: d_t_lif(:,:)
62      !$OMP THREADPRIVATE(d_t_lif)
63      REAL, SAVE, ALLOCATABLE :: d_u_lif(:,:), d_v_lif(:,:)
64      !$OMP THREADPRIVATE(d_u_lif, d_v_lif)
65! Tendances Ondes de G non oro (runs strato).
66      REAL, SAVE, ALLOCATABLE :: d_u_hin(:,:)
67      !$OMP THREADPRIVATE(d_u_hin)
68      REAL, SAVE, ALLOCATABLE :: d_v_hin(:,:)
69      !$OMP THREADPRIVATE(d_v_hin)
70      REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:)
71      !$OMP THREADPRIVATE(d_t_hin)
72
73! tendance du a la conersion Ec -> E thermique
74      REAL, SAVE, ALLOCATABLE :: d_t_ec(:,:)
75      !$OMP THREADPRIVATE(d_t_ec)
76      REAL, SAVE, ALLOCATABLE :: d_ts(:,:), d_tr(:,:,:)
77      !$OMP THREADPRIVATE(d_ts, d_tr)
78
79! diagnostique pour le rayonnement
80      REAL, SAVE, ALLOCATABLE :: topswad_aero(:),  solswad_aero(:)      ! diag
81      !$OMP THREADPRIVATE(topswad_aero,solswad_aero)
82      REAL, SAVE, ALLOCATABLE :: topswai_aero(:),  solswai_aero(:)      ! diag
83      !$OMP THREADPRIVATE(topswai_aero,solswai_aero)
84      REAL, SAVE, ALLOCATABLE :: topswad0_aero(:), solswad0_aero(:)     ! diag
85      !$OMP THREADPRIVATE(topswad0_aero,solswad0_aero)
86      REAL, SAVE, ALLOCATABLE :: topsw_aero(:,:),  solsw_aero(:,:)      ! diag
87      !$OMP THREADPRIVATE(topsw_aero,solsw_aero)
88      REAL, SAVE, ALLOCATABLE :: topsw0_aero(:,:), solsw0_aero(:,:)     ! diag
89      !$OMP THREADPRIVATE(topsw0_aero,solsw0_aero)
90      REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:),  solswcf_aero(:,:)  ! diag
91      !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero)
92      REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:)
93      !$OMP THREADPRIVATE(tausum_aero)
94      REAL, SAVE, ALLOCATABLE :: tau3d_aero(:,:,:,:)
95      !$OMP THREADPRIVATE(tau3d_aero)
96      REAL, SAVE, ALLOCATABLE :: scdnc(:,:)
97      !$OMP THREADPRIVATE(scdnc)
98      REAL, SAVE, ALLOCATABLE :: cldncl(:)
99      !$OMP THREADPRIVATE(cldncl)
100      REAL, SAVE, ALLOCATABLE :: reffclwtop(:)
101      !$OMP THREADPRIVATE(reffclwtop)
102      REAL, SAVE, ALLOCATABLE :: lcc(:)
103      !$OMP THREADPRIVATE(lcc)
104      REAL, SAVE, ALLOCATABLE :: reffclws(:,:)
105      !$OMP THREADPRIVATE(reffclws)
106      REAL, SAVE, ALLOCATABLE :: reffclwc(:,:)
107      !$OMP THREADPRIVATE(reffclwc)
108      REAL, SAVE, ALLOCATABLE :: cldnvi(:)
109      !$OMP THREADPRIVATE(cldnvi)
110      REAL, SAVE, ALLOCATABLE :: lcc3d(:,:)
111      !$OMP THREADPRIVATE(lcc3d)
112      REAL, SAVE, ALLOCATABLE :: lcc3dcon(:,:)
113      !$OMP THREADPRIVATE(lcc3dcon)
114      REAL, SAVE, ALLOCATABLE :: lcc3dstra(:,:)
115      !$OMP THREADPRIVATE(lcc3dstra)
116      REAL, SAVE, ALLOCATABLE :: od550aer(:)
117      !$OMP THREADPRIVATE(od550aer)
118      REAL, SAVE, ALLOCATABLE :: absvisaer(:)
119      !$OMP THREADPRIVATE(absvisaer)
120      REAL, SAVE, ALLOCATABLE :: od865aer(:)
121      !$OMP THREADPRIVATE(od865aer)
122      REAL, SAVE, ALLOCATABLE :: ec550aer(:,:)
123      !$OMP THREADPRIVATE(ec550aer)
124      REAL, SAVE, ALLOCATABLE :: od550lt1aer(:)
125      !$OMP THREADPRIVATE(od550lt1aer)
126      REAL, SAVE, ALLOCATABLE :: sconcso4(:)
127      !$OMP THREADPRIVATE(sconcso4)
128      REAL, SAVE, ALLOCATABLE :: sconcoa(:)
129      !$OMP THREADPRIVATE(sconcoa)
130      REAL, SAVE, ALLOCATABLE :: sconcbc(:)
131      !$OMP THREADPRIVATE(sconcbc)
132      REAL, SAVE, ALLOCATABLE :: sconcss(:)
133      !$OMP THREADPRIVATE(sconcss)
134      REAL, SAVE, ALLOCATABLE :: sconcdust(:)
135      !$OMP THREADPRIVATE(sconcdust)
136      REAL, SAVE, ALLOCATABLE :: concso4(:,:)
137      !$OMP THREADPRIVATE(concso4)
138      REAL, SAVE, ALLOCATABLE :: concoa(:,:)
139      !$OMP THREADPRIVATE(concoa)
140      REAL, SAVE, ALLOCATABLE :: concbc(:,:)
141      !$OMP THREADPRIVATE(concbc)
142      REAL, SAVE, ALLOCATABLE :: concss(:,:)
143      !$OMP THREADPRIVATE(concss)
144      REAL, SAVE, ALLOCATABLE :: concdust(:,:)
145      !$OMP THREADPRIVATE(concdust)
146      REAL, SAVE, ALLOCATABLE :: loadso4(:)
147      !$OMP THREADPRIVATE(loadso4)
148      REAL, SAVE, ALLOCATABLE :: loadoa(:)
149      !$OMP THREADPRIVATE(loadoa)
150      REAL, SAVE, ALLOCATABLE :: loadbc(:)
151      !$OMP THREADPRIVATE(loadbc)
152      REAL, SAVE, ALLOCATABLE :: loadss(:)
153      !$OMP THREADPRIVATE(loadss)
154      REAL, SAVE, ALLOCATABLE :: loaddust(:)
155      !$OMP THREADPRIVATE(loaddust)
156      REAL, SAVE, ALLOCATABLE :: load_tmp1(:)
157      !$OMP THREADPRIVATE(load_tmp1)
158      REAL, SAVE, ALLOCATABLE :: load_tmp2(:)
159      !$OMP THREADPRIVATE(load_tmp2)
160      REAL, SAVE, ALLOCATABLE :: load_tmp3(:)
161      !$OMP THREADPRIVATE(load_tmp3)
162      REAL, SAVE, ALLOCATABLE :: load_tmp4(:)
163      !$OMP THREADPRIVATE(load_tmp4)
164      REAL, SAVE, ALLOCATABLE :: load_tmp5(:)
165      !$OMP THREADPRIVATE(load_tmp5)
166      REAL, SAVE, ALLOCATABLE :: load_tmp6(:)
167      !$OMP THREADPRIVATE(load_tmp6)
168      REAL, SAVE, ALLOCATABLE :: load_tmp7(:)
169      !$OMP THREADPRIVATE(load_tmp7)
170
171!IM ajout variables CFMIP2/CMIP5
172      REAL,ALLOCATABLE,SAVE :: topswad_aerop(:), solswad_aerop(:)
173!$OMP THREADPRIVATE(topswad_aerop, solswad_aerop)
174      REAL,ALLOCATABLE,SAVE :: topswai_aerop(:), solswai_aerop(:)
175!$OMP THREADPRIVATE(topswai_aerop, solswai_aerop)
176      REAL,ALLOCATABLE,SAVE :: topswad0_aerop(:), solswad0_aerop(:)
177!$OMP THREADPRIVATE(topswad0_aerop, solswad0_aerop)
178      REAL,ALLOCATABLE,SAVE :: topsw_aerop(:,:), topsw0_aerop(:,:) 
179!$OMP THREADPRIVATE(topsw_aerop, topsw0_aerop)
180      REAL,ALLOCATABLE,SAVE :: solsw_aerop(:,:), solsw0_aerop(:,:)
181!$OMP THREADPRIVATE(solsw_aerop, solsw0_aerop)
182      REAL,ALLOCATABLE,SAVE :: topswcf_aerop(:,:), solswcf_aerop(:,:)
183!$OMP THREADPRIVATE(topswcf_aerop, solswcf_aerop)
184
185CONTAINS
186
187!======================================================================
188SUBROUTINE phys_local_var_init
189use dimphy
190use infotrac, ONLY : nbtr
191USE aero_mod
192
193IMPLICIT NONE
194#include "indicesol.h"
195      allocate(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
196      allocate(u_seri(klon,klev),v_seri(klon,klev))
197
198      allocate(tr_seri(klon,klev,nbtr))
199      allocate(d_t_dyn(klon,klev),d_q_dyn(klon,klev))
200      allocate(d_u_dyn(klon,klev),d_v_dyn(klon,klev))
201      allocate(d_t_con(klon,klev),d_q_con(klon,klev))
202      allocate(d_u_con(klon,klev),d_v_con(klon,klev))
203      allocate(d_t_wake(klon,klev),d_q_wake(klon,klev))
204      allocate(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
205      allocate(d_ql_lsc(klon,klev))
206      allocate(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
207      allocate(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
208      allocate(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
209      allocate(d_t_eva(klon,klev),d_q_eva(klon,klev))
210      allocate(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
211      allocate(d_t_lscth(klon,klev),d_q_lscth(klon,klev))
212      allocate(plul_st(klon),plul_th(klon))
213      allocate(d_t_vdf(klon,klev),d_q_vdf(klon,klev))
214      allocate(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
215      allocate(d_t_oli(klon,klev),d_t_oro(klon,klev))
216      allocate(d_u_oli(klon,klev),d_v_oli(klon,klev))
217      allocate(d_u_oro(klon,klev),d_v_oro(klon,klev))
218      allocate(d_t_lif(klon,klev),d_t_ec(klon,klev))
219      allocate(d_u_lif(klon,klev),d_v_lif(klon,klev))
220      allocate(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
221      allocate(topswad_aero(klon), solswad_aero(klon))
222      allocate(topswai_aero(klon), solswai_aero(klon))
223      allocate(topswad0_aero(klon), solswad0_aero(klon))
224      allocate(topsw_aero(klon,naero_grp), solsw_aero(klon,naero_grp))
225      allocate(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
226      allocate(topswcf_aero(klon,3), solswcf_aero(klon,3))
227      allocate(d_u_hin(klon,klev),d_v_hin(klon,klev),d_t_hin(klon,klev))
228      allocate(tausum_aero(klon,nwave,naero_spc))
229      allocate(tau3d_aero(klon,klev,nwave,naero_spc))
230      allocate(scdnc(klon, klev))
231      allocate(cldncl(klon))
232      allocate(reffclwtop(klon))
233      allocate(lcc(klon))
234      allocate(reffclws(klon, klev))
235      allocate(reffclwc(klon, klev))
236      allocate(cldnvi(klon))
237      allocate(lcc3d(klon, klev))
238      allocate(lcc3dcon(klon, klev))
239      allocate(lcc3dstra(klon, klev))
240      allocate(od550aer(klon))   
241      allocate(od865aer(klon))   
242      allocate(absvisaer(klon)) 
243      allocate(ec550aer(klon,klev))
244      allocate(od550lt1aer(klon))               
245      allocate(sconcso4(klon))
246      allocate(sconcoa(klon))
247      allocate(sconcbc(klon))
248      allocate(sconcss(klon))
249      allocate(sconcdust(klon))
250      allocate(concso4(klon,klev))
251      allocate(concoa(klon,klev))
252      allocate(concbc(klon,klev))
253      allocate(concss(klon,klev))
254      allocate(concdust(klon,klev))
255      allocate(loadso4(klon))
256      allocate(loadoa(klon))
257      allocate(loadbc(klon))
258      allocate(loadss(klon))
259      allocate(loaddust(klon))
260      allocate(load_tmp1(klon))
261      allocate(load_tmp2(klon))
262      allocate(load_tmp3(klon))
263      allocate(load_tmp4(klon))
264      allocate(load_tmp5(klon))
265      allocate(load_tmp6(klon))
266      allocate(load_tmp7(klon))
267
268!IM ajout variables CFMIP2/CMIP5
269      ALLOCATE(topswad_aerop(klon), solswad_aerop(klon))
270      ALLOCATE(topswai_aerop(klon), solswai_aerop(klon))
271      ALLOCATE(topswad0_aerop(klon), solswad0_aerop(klon))
272      ALLOCATE(topsw_aerop(klon,naero_grp), topsw0_aerop(klon,naero_grp))
273      ALLOCATE(solsw_aerop(klon,naero_grp), solsw0_aerop(klon,naero_grp))
274      ALLOCATE(topswcf_aerop(klon,naero_grp), solswcf_aerop(klon,naero_grp))
275
276END SUBROUTINE phys_local_var_init
277
278!======================================================================
279SUBROUTINE phys_local_var_end
280use dimphy
281IMPLICIT NONE
282#include "indicesol.h"
283      deallocate(t_seri,q_seri,ql_seri,qs_seri)
284      deallocate(u_seri,v_seri)
285
286      deallocate(tr_seri)
287      deallocate(d_t_dyn,d_q_dyn)
288      deallocate(d_u_dyn,d_v_dyn)
289      deallocate(d_t_con,d_q_con)
290      deallocate(d_u_con,d_v_con)
291      deallocate(d_t_wake,d_q_wake)
292      deallocate(d_t_lsc,d_q_lsc)
293      deallocate(d_ql_lsc)
294      deallocate(d_t_ajsb,d_q_ajsb)
295      deallocate(d_t_ajs,d_q_ajs)
296      deallocate(d_u_ajs,d_v_ajs)
297      deallocate(d_t_eva,d_q_eva)
298      deallocate(d_t_lscst,d_q_lscst)
299      deallocate(d_t_lscth,d_q_lscth)
300      deallocate(plul_st,plul_th)
301      deallocate(d_t_vdf,d_q_vdf)
302      deallocate(d_u_vdf,d_v_vdf)
303      deallocate(d_t_oli,d_t_oro)
304      deallocate(d_u_oli,d_v_oli)
305      deallocate(d_u_oro,d_v_oro)
306      deallocate(d_t_lif,d_t_ec)
307      deallocate(d_u_lif,d_v_lif)
308      deallocate(d_ts, d_tr)
309      deallocate(topswad_aero,solswad_aero)
310      deallocate(topswai_aero,solswai_aero)
311      deallocate(topswad0_aero,solswad0_aero)
312      deallocate(topsw_aero,solsw_aero)
313      deallocate(topsw0_aero,solsw0_aero)
314      deallocate(topswcf_aero,solswcf_aero)
315      deallocate(tausum_aero)
316      deallocate(tau3d_aero)
317      deallocate(scdnc)
318      deallocate(cldncl)
319      deallocate(reffclwtop)
320      deallocate(lcc)
321      deallocate(reffclws)
322      deallocate(reffclwc)
323      deallocate(cldnvi)
324      deallocate(lcc3d)
325      deallocate(lcc3dcon)
326      deallocate(lcc3dstra)
327      deallocate(od550aer)       
328      deallocate(od865aer)
329      deallocate(absvisaer)
330      deallocate(ec550aer)
331      deallocate(od550lt1aer)
332      deallocate(sconcso4)
333      deallocate(sconcoa)
334      deallocate(sconcbc)
335      deallocate(sconcss)
336      deallocate(sconcdust)
337      deallocate(concso4)
338      deallocate(concoa)
339      deallocate(concbc)
340      deallocate(concss)
341      deallocate(concdust)
342      deallocate(loadso4)
343      deallocate(loadoa)
344      deallocate(loadbc)
345      deallocate(loadss)
346      deallocate(loaddust)
347      deallocate(load_tmp1)
348      deallocate(load_tmp2)
349      deallocate(load_tmp3)
350      deallocate(load_tmp4)
351      deallocate(load_tmp5)
352      deallocate(load_tmp6)
353      deallocate(load_tmp7)
354      deallocate(d_u_hin,d_v_hin,d_t_hin)
355
356!IM ajout variables CFMIP2/CMIP5
357      deallocate(topswad_aerop, solswad_aerop)
358      deallocate(topswai_aerop, solswai_aerop)
359      deallocate(topswad0_aerop, solswad0_aerop)
360      deallocate(topsw_aerop, topsw0_aerop)
361      deallocate(solsw_aerop, solsw0_aerop)
362      deallocate(topswcf_aerop, solswcf_aerop)
363
364END SUBROUTINE phys_local_var_end
365
366END MODULE phys_local_var_mod
Note: See TracBrowser for help on using the repository browser.