source: LMDZ4/trunk/libf/phytherm/write_paramLMDZ_phy.h @ 832

Last change on this file since 832 was 814, checked in by Laurent Fairhead, 17 years ago

Rajout de la physique utilisant les thermiques FH
LF
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 11.1 KB
RevLine 
[814]1c
2      IF (is_sequential) THEN
3     
4      ndex2d = 0
5      itau_w = itau_phy + itap
6c
7c Variables type caractere : plusieurs valeurs possibles
8c
9      IF(ocean.EQ.'force ') THEN
10       zx_tmp_2d(1:iim,1:jjmp1)=1.
11      ELSE IF(ocean.EQ.'slab  ') THEN
12       zx_tmp_2d(1:iim,1:jjmp1)=2.
13      ELSE IF(ocean.EQ.'couple') THEN
14       zx_tmp_2d(1:iim,1:jjmp1)=3.
15      ENDIF
16      CALL histwrite(nid_ctesGCM,"ocean",itau_w,
17     .               zx_tmp_2d,iim*jjmp1,ndex2d)
18c
19      IF(ok_slab_sicOBS) THEN
20       zx_tmp_2d(1:iim,1:jjmp1)=1.
21      ELSE
22       zx_tmp_2d(1:iim,1:jjmp1)=0.
23      ENDIF
24      CALL histwrite(nid_ctesGCM,"ok_slab_sicOBS",itau_w,
25     .               zx_tmp_2d,iim*jjmp1,ndex2d)
26c
27      IF(type_run.EQ.'CLIM'.OR.type_run.EQ.'ENSP') THEN
28       zx_tmp_2d(1:iim,1:jjmp1)=1.
29      ELSE IF(type_run.EQ.'AMIP'.OR.type_run.EQ.'CFMI') THEN
30       zx_tmp_2d(1:iim,1:jjmp1)=2.
31      ENDIF
32      CALL histwrite(nid_ctesGCM,"type_run",itau_w,
33     .               zx_tmp_2d,iim*jjmp1,ndex2d)
34c
35c Variables logiques (1=true, 2=false)
36c
37      IF(ok_veget) THEN
38       zx_tmp_2d(1:iim,1:jjmp1)=1.
39      ELSE
40       zx_tmp_2d(1:iim,1:jjmp1)=0.
41      ENDIF
42      CALL histwrite(nid_ctesGCM,"ok_veget",itau_w,
43     .               zx_tmp_2d,iim*jjmp1,ndex2d)
44c
45      IF(ok_journe) THEN
46       zx_tmp_2d(1:iim,1:jjmp1)=1.
47      ELSE
48       zx_tmp_2d(1:iim,1:jjmp1)=0.
49      ENDIF
50      CALL histwrite(nid_ctesGCM,"ok_journe",itau_w,
51     .               zx_tmp_2d,iim*jjmp1,ndex2d)
52c
53      IF(ok_mensuel) THEN
54       zx_tmp_2d(1:iim,1:jjmp1)=1.
55      ELSE
56       zx_tmp_2d(1:iim,1:jjmp1)=0.
57      ENDIF
58      CALL histwrite(nid_ctesGCM,"ok_mensuel",itau_w,
59     .               zx_tmp_2d,iim*jjmp1,ndex2d)
60c
61      IF(ok_instan) THEN
62       zx_tmp_2d(1:iim,1:jjmp1)=1.
63      ELSE
64       zx_tmp_2d(1:iim,1:jjmp1)=0.
65      ENDIF
66      CALL histwrite(nid_ctesGCM,"ok_instan",itau_w,
67     .               zx_tmp_2d,iim*jjmp1,ndex2d)
68c
69      IF(ok_ade) THEN
70       zx_tmp_2d(1:iim,1:jjmp1)=1.
71      ELSE
72       zx_tmp_2d(1:iim,1:jjmp1)=0.
73      ENDIF
74      CALL histwrite(nid_ctesGCM,"ok_ade",itau_w,
75     .               zx_tmp_2d,iim*jjmp1,ndex2d)
76c
77      IF(ok_aie) THEN
78       zx_tmp_2d(1:iim,1:jjmp1)=1.
79      ELSE
80       zx_tmp_2d(1:iim,1:jjmp1)=0.
81      ENDIF
82      CALL histwrite(nid_ctesGCM,"ok_aie",itau_w,
83     .               zx_tmp_2d,iim*jjmp1,ndex2d)
84c
85
86c
87c Champs 2D:
88c
89      zx_tmp_2d(1:iim,1:jjmp1)=bl95_b0
90      CALL histwrite(nid_ctesGCM,"bl95_b0",itau_w,
91     .               zx_tmp_2d,iim*jjmp1,ndex2d)
92c
93      zx_tmp_2d(1:iim,1:jjmp1)=bl95_b1
94      CALL histwrite(nid_ctesGCM,"bl95_b1",itau_w,
95     .               zx_tmp_2d,iim*jjmp1,ndex2d)
96c
97      zx_tmp_2d(1:iim,1:jjmp1)=ip_ebil_phy
98      CALL histwrite(nid_ctesGCM,"ip_ebil_phy",itau_w,
99     .               zx_tmp_2d,iim*jjmp1,ndex2d)
100c
101      zx_tmp_2d(1:iim,1:jjmp1)=R_ecc
102      CALL histwrite(nid_ctesGCM,"R_ecc",itau_w,
103     .               zx_tmp_2d,iim*jjmp1,ndex2d)
104c
105      zx_tmp_2d(1:iim,1:jjmp1)=R_peri
106      CALL histwrite(nid_ctesGCM,"R_peri",itau_w,
107     .               zx_tmp_2d,iim*jjmp1,ndex2d)
108c
109      zx_tmp_2d(1:iim,1:jjmp1)=R_incl
110      CALL histwrite(nid_ctesGCM,"R_incl",itau_w,
111     .               zx_tmp_2d,iim*jjmp1,ndex2d)
112c
113      zx_tmp_2d(1:iim,1:jjmp1)=solaire
114      CALL histwrite(nid_ctesGCM,"solaire",itau_w,
115     .               zx_tmp_2d,iim*jjmp1,ndex2d)
116c
117      zx_tmp_2d(1:iim,1:jjmp1)=co2_ppm
118      CALL histwrite(nid_ctesGCM,"co2_ppm",itau_w,
119     .               zx_tmp_2d,iim*jjmp1,ndex2d)
120c
121      zx_tmp_2d(1:iim,1:jjmp1)=CH4_ppb
122      CALL histwrite(nid_ctesGCM,"CH4_ppb",itau_w,
123     .               zx_tmp_2d,iim*jjmp1,ndex2d)
124c
125      zx_tmp_2d(1:iim,1:jjmp1)=N2O_ppb
126      CALL histwrite(nid_ctesGCM,"N2O_ppb",itau_w,
127     .               zx_tmp_2d,iim*jjmp1,ndex2d)
128c
129      zx_tmp_2d(1:iim,1:jjmp1)=CFC11_ppt
130      CALL histwrite(nid_ctesGCM,"CFC11_ppt",itau_w,
131     .               zx_tmp_2d,iim*jjmp1,ndex2d)
132c
133      zx_tmp_2d(1:iim,1:jjmp1)=CFC12_ppt
134      CALL histwrite(nid_ctesGCM,"CFC12_ppt",itau_w,
135     .               zx_tmp_2d,iim*jjmp1,ndex2d)
136c
137      zx_tmp_2d(1:iim,1:jjmp1)=epmax
138      CALL histwrite(nid_ctesGCM,"epmax",itau_w,
139     .               zx_tmp_2d,iim*jjmp1,ndex2d)
140c
141      IF(ok_adj_ema) THEN
142       zx_tmp_2d(1:iim,1:jjmp1)=1.
143      ELSE
144       zx_tmp_2d(1:iim,1:jjmp1)=0.
145      ENDIF
146      CALL histwrite(nid_ctesGCM,"ok_adj_ema",itau_w,
147     .               zx_tmp_2d,iim*jjmp1,ndex2d)
148c
149      zx_tmp_2d(1:iim,1:jjmp1)=iflag_clw
150      CALL histwrite(nid_ctesGCM,"iflag_clw",itau_w,
151     .               zx_tmp_2d,iim*jjmp1,ndex2d)
152c
153      zx_tmp_2d(1:iim,1:jjmp1)=cld_lc_lsc
154      CALL histwrite(nid_ctesGCM,"cld_lc_lsc",itau_w,
155     .               zx_tmp_2d,iim*jjmp1,ndex2d)
156c
157      zx_tmp_2d(1:iim,1:jjmp1)=cld_lc_con
158      CALL histwrite(nid_ctesGCM,"cld_lc_con",itau_w,
159     .               zx_tmp_2d,iim*jjmp1,ndex2d)
160c
161      zx_tmp_2d(1:iim,1:jjmp1)=cld_tau_lsc
162      CALL histwrite(nid_ctesGCM,"cld_tau_lsc",itau_w,
163     .               zx_tmp_2d,iim*jjmp1,ndex2d)
164c
165      zx_tmp_2d(1:iim,1:jjmp1)=cld_tau_con
166      CALL histwrite(nid_ctesGCM,"cld_tau_con",itau_w,
167     .               zx_tmp_2d,iim*jjmp1,ndex2d)
168c
169      zx_tmp_2d(1:iim,1:jjmp1)=ffallv_lsc
170      CALL histwrite(nid_ctesGCM,"ffallv_lsc",itau_w,
171     .               zx_tmp_2d,iim*jjmp1,ndex2d)
172c
173      zx_tmp_2d(1:iim,1:jjmp1)=ffallv_con
174      CALL histwrite(nid_ctesGCM,"ffallv_con",itau_w,
175     .               zx_tmp_2d,iim*jjmp1,ndex2d)
176c
177      zx_tmp_2d(1:iim,1:jjmp1)=coef_eva
178      CALL histwrite(nid_ctesGCM,"coef_eva",itau_w,
179     .               zx_tmp_2d,iim*jjmp1,ndex2d)
180c
181      IF(reevap_ice) THEN
182       zx_tmp_2d(1:iim,1:jjmp1)=1.
183      ELSE
184       zx_tmp_2d(1:iim,1:jjmp1)=0.
185      ENDIF
186      CALL histwrite(nid_ctesGCM,"reevap_ice",itau_w,
187     .               zx_tmp_2d,iim*jjmp1,ndex2d)
188c
189      zx_tmp_2d(1:iim,1:jjmp1)=iflag_cldcon
190      CALL histwrite(nid_ctesGCM,"iflag_cldcon",itau_w,
191     .               zx_tmp_2d,iim*jjmp1,ndex2d)
192c
193      zx_tmp_2d(1:iim,1:jjmp1)=iflag_pdf
194      CALL histwrite(nid_ctesGCM,"iflag_pdf",itau_w,
195     .               zx_tmp_2d,iim*jjmp1,ndex2d)
196c
197      zx_tmp_2d(1:iim,1:jjmp1)=fact_cldcon
198      CALL histwrite(nid_ctesGCM,"fact_cldcon",itau_w,
199     .               zx_tmp_2d,iim*jjmp1,ndex2d)
200c
201      zx_tmp_2d(1:iim,1:jjmp1)=facttemps
202      CALL histwrite(nid_ctesGCM,"facttemps",itau_w,
203     .               zx_tmp_2d,iim*jjmp1,ndex2d)
204c
205      IF(ok_newmicro) THEN
206       zx_tmp_2d(1:iim,1:jjmp1)=1.
207      ELSE
208       zx_tmp_2d(1:iim,1:jjmp1)=0.
209      ENDIF
210      CALL histwrite(nid_ctesGCM,"ok_newmicro",itau_w,
211     .               zx_tmp_2d,iim*jjmp1,ndex2d)
212c
213      zx_tmp_2d(1:iim,1:jjmp1)=ratqsbas
214      CALL histwrite(nid_ctesGCM,"ratqsbas",itau_w,
215     .               zx_tmp_2d,iim*jjmp1,ndex2d)
216c
217      zx_tmp_2d(1:iim,1:jjmp1)=ratqshaut
218      CALL histwrite(nid_ctesGCM,"ratqshaut",itau_w,
219     .               zx_tmp_2d,iim*jjmp1,ndex2d)
220c
221      zx_tmp_2d(1:iim,1:jjmp1)=rad_froid
222      CALL histwrite(nid_ctesGCM,"rad_froid",itau_w,
223     .               zx_tmp_2d,iim*jjmp1,ndex2d)
224c
225      zx_tmp_2d(1:iim,1:jjmp1)=rad_chau1
226      CALL histwrite(nid_ctesGCM,"rad_chau1",itau_w,
227     .               zx_tmp_2d,iim*jjmp1,ndex2d)
228c
229      zx_tmp_2d(1:iim,1:jjmp1)=rad_chau2
230      CALL histwrite(nid_ctesGCM,"rad_chau2",itau_w,
231     .               zx_tmp_2d,iim*jjmp1,ndex2d)
232c
233      zx_tmp_2d(1:iim,1:jjmp1)=top_height
234      CALL histwrite(nid_ctesGCM,"top_height",itau_w,
235     .               zx_tmp_2d,iim*jjmp1,ndex2d)
236c
237      zx_tmp_2d(1:iim,1:jjmp1)=overlap
238      CALL histwrite(nid_ctesGCM,"overlap",itau_w,
239     .               zx_tmp_2d,iim*jjmp1,ndex2d)
240c
241      zx_tmp_2d(1:iim,1:jjmp1)=cdmmax
242      CALL histwrite(nid_ctesGCM,"cdmmax",itau_w,
243     .               zx_tmp_2d,iim*jjmp1,ndex2d)
244c
245      zx_tmp_2d(1:iim,1:jjmp1)=cdhmax
246      CALL histwrite(nid_ctesGCM,"cdhmax",itau_w,
247     .               zx_tmp_2d,iim*jjmp1,ndex2d)
248c
249      zx_tmp_2d(1:iim,1:jjmp1)=ksta
250      CALL histwrite(nid_ctesGCM,"ksta",itau_w,
251     .               zx_tmp_2d,iim*jjmp1,ndex2d)
252c
253      zx_tmp_2d(1:iim,1:jjmp1)=ksta_ter
254      CALL histwrite(nid_ctesGCM,"ksta_ter",itau_w,
255     .               zx_tmp_2d,iim*jjmp1,ndex2d)
256c
257      IF(ok_kzmin) THEN
258       zx_tmp_2d(1:iim,1:jjmp1)=1.
259      ELSE
260       zx_tmp_2d(1:iim,1:jjmp1)=0.
261      ENDIF
262      CALL histwrite(nid_ctesGCM,"ok_kzmin",itau_w,
263     .               zx_tmp_2d,iim*jjmp1,ndex2d)
264c
265      zx_tmp_2d(1:iim,1:jjmp1)=iflag_pbl
266      CALL histwrite(nid_ctesGCM,"iflag_pbl",itau_w,
267     .               zx_tmp_2d,iim*jjmp1,ndex2d)
268c
269      zx_tmp_2d(1:iim,1:jjmp1)=lev_histhf
270      CALL histwrite(nid_ctesGCM,"lev_histhf",itau_w,
271     .               zx_tmp_2d,iim*jjmp1,ndex2d)
272c
273      zx_tmp_2d(1:iim,1:jjmp1)=lev_histday
274      CALL histwrite(nid_ctesGCM,"lev_histday",itau_w,
275     .               zx_tmp_2d,iim*jjmp1,ndex2d)
276c
277      zx_tmp_2d(1:iim,1:jjmp1)=lev_histmth
278      CALL histwrite(nid_ctesGCM,"lev_histmth",itau_w,
279     .               zx_tmp_2d,iim*jjmp1,ndex2d)
280c
281      IF(ok_isccp) THEN
282       zx_tmp_2d(1:iim,1:jjmp1)=1.
283      ELSE
284       zx_tmp_2d(1:iim,1:jjmp1)=0.
285      ENDIF
286      CALL histwrite(nid_ctesGCM,"ok_isccp",itau_w,
287     .               zx_tmp_2d,iim*jjmp1,ndex2d)
288c
289      zx_tmp_2d(1:iim,1:jjmp1)=lonmin_ins
290      CALL histwrite(nid_ctesGCM,"lonmin_ins",itau_w,
291     .               zx_tmp_2d,iim*jjmp1,ndex2d)
292c
293      zx_tmp_2d(1:iim,1:jjmp1)=lonmax_ins
294      CALL histwrite(nid_ctesGCM,"lonmax_ins",itau_w,
295     .               zx_tmp_2d,iim*jjmp1,ndex2d)
296c
297      zx_tmp_2d(1:iim,1:jjmp1)=latmin_ins
298      CALL histwrite(nid_ctesGCM,"latmin_ins",itau_w,
299     .               zx_tmp_2d,iim*jjmp1,ndex2d)
300c
301      zx_tmp_2d(1:iim,1:jjmp1)=latmax_ins
302      CALL histwrite(nid_ctesGCM,"latmax_ins",itau_w,
303     .               zx_tmp_2d,iim*jjmp1,ndex2d)
304c
305      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_ins
306      CALL histwrite(nid_ctesGCM,"ecrit_ins",itau_w,
307     .               zx_tmp_2d,iim*jjmp1,ndex2d)
308c
309      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_hf
310      CALL histwrite(nid_ctesGCM,"ecrit_hf",itau_w,
311     .               zx_tmp_2d,iim*jjmp1,ndex2d)
312c
313      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_day
314      CALL histwrite(nid_ctesGCM,"ecrit_day",itau_w,
315     .               zx_tmp_2d,iim*jjmp1,ndex2d)
316c
317      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_mth
318      CALL histwrite(nid_ctesGCM,"ecrit_mth",itau_w,
319     .               zx_tmp_2d,iim*jjmp1,ndex2d)
320c
321      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_tra
322      CALL histwrite(nid_ctesGCM,"ecrit_tra",itau_w,
323     .               zx_tmp_2d,iim*jjmp1,ndex2d)
324c
325      zx_tmp_2d(1:iim,1:jjmp1)=ecrit_reg
326      CALL histwrite(nid_ctesGCM,"ecrit_reg",itau_w,
327     .               zx_tmp_2d,iim*jjmp1,ndex2d)
328c
329      zx_tmp_2d(1:iim,1:jjmp1)=freqin_isccp
330      CALL histwrite(nid_ctesGCM,"freqin_isccp",itau_w,
331     .               zx_tmp_2d,iim*jjmp1,ndex2d)
332c
333      zx_tmp_2d(1:iim,1:jjmp1)=freqout_isccp
334      CALL histwrite(nid_ctesGCM,"freqout_isccp",itau_w,
335     .               zx_tmp_2d,iim*jjmp1,ndex2d)
336c
337c=================================================================
338c=================================================================
339c=================================================================
340c
341      if (ok_sync) then
342        call histsync(nid_ctesGCM)
343      endif
344c
345      ENDIF ! mono_cpu
Note: See TracBrowser for help on using the repository browser.