source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_local_var_mod.F90 @ 3536

Last change on this file since 3536 was 1534, checked in by musat, 14 years ago

Ajouts CFMIP2/CMIP5

  • 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers: PARAM/npCFMIP_param.data contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) PARAM/pointlocations.txt contenat 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

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