source: LMDZ4/branches/pre_V3/libf/phylmd/write_paramLMDZ_phy.h @ 2629

Last change on this file since 2629 was 683, checked in by lmdzadmin, 19 years ago

New routines contenant les flags de run.def
IM

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