Context Navigation

source: LMDZ4/trunk/libf/phylmd/write_paramLMDZ_phy.h @ 719

Last change on this file since 719 was 719, checked in by lmdzadmin, 18 years ago
Correction typage ecrit_day, ecrit_mth (*.h) et definition rh2m (physiq.F) FH/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
1	c
2	ndex2d = 0
3	itau_w = itau_phy + itap
4	c
5	c Variables type caractere : plusieurs valeurs possibles
6	c
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)
16	c
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)
24	c
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)
32	c
33	c Variables logiques (1=true, 2=false)
34	c
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)
42	c
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)
50	c
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)
58	c
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)
66	c
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)
74	c
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)
82	c
83
84	c
85	c Champs 2D:
86	c
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)
90	c
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)
94	c
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)
98	c
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)
102	c
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)
106	c
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)
110	c
111	zx_tmp_2d(1:iim,1:jjmp1)=solaire
112	CALL histwrite(nid_ctesGCM,"solaire",itau_w,
113	. zx_tmp_2d,iim*jjmp1,ndex2d)
114	c
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)
118	c
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)
122	c
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)
126	c
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)
130	c
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)
134	c
135	zx_tmp_2d(1:iim,1:jjmp1)=epmax
136	CALL histwrite(nid_ctesGCM,"epmax",itau_w,
137	. zx_tmp_2d,iim*jjmp1,ndex2d)
138	c
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)
146	c
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)
150	c
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)
154	c
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)
158	c
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)
162	c
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)
166	c
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)
170	c
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)
174	c
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)
178	c
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)
186	c
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)
190	c
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)
194	c
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)
198	c
199	zx_tmp_2d(1:iim,1:jjmp1)=facttemps
200	CALL histwrite(nid_ctesGCM,"facttemps",itau_w,
201	. zx_tmp_2d,iim*jjmp1,ndex2d)
202	c
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)
210	c
211	zx_tmp_2d(1:iim,1:jjmp1)=ratqsbas
212	CALL histwrite(nid_ctesGCM,"ratqsbas",itau_w,
213	. zx_tmp_2d,iim*jjmp1,ndex2d)
214	c
215	zx_tmp_2d(1:iim,1:jjmp1)=ratqshaut
216	CALL histwrite(nid_ctesGCM,"ratqshaut",itau_w,
217	. zx_tmp_2d,iim*jjmp1,ndex2d)
218	c
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)
222	c
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)
226	c
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)
230	c
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)
234	c
235	zx_tmp_2d(1:iim,1:jjmp1)=overlap
236	CALL histwrite(nid_ctesGCM,"overlap",itau_w,
237	. zx_tmp_2d,iim*jjmp1,ndex2d)
238	c
239	zx_tmp_2d(1:iim,1:jjmp1)=cdmmax
240	CALL histwrite(nid_ctesGCM,"cdmmax",itau_w,
241	. zx_tmp_2d,iim*jjmp1,ndex2d)
242	c
243	zx_tmp_2d(1:iim,1:jjmp1)=cdhmax
244	CALL histwrite(nid_ctesGCM,"cdhmax",itau_w,
245	. zx_tmp_2d,iim*jjmp1,ndex2d)
246	c
247	zx_tmp_2d(1:iim,1:jjmp1)=ksta
248	CALL histwrite(nid_ctesGCM,"ksta",itau_w,
249	. zx_tmp_2d,iim*jjmp1,ndex2d)
250	c
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)
254	c
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)
262	c
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)
266	c
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)
270	c
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)
274	c
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)
278	c
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)
286	c
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)
290	c
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)
294	c
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)
298	c
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)
302	c
303	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_ins
304	CALL histwrite(nid_ctesGCM,"ecrit_ins",itau_w,
305	. zx_tmp_2d,iim*jjmp1,ndex2d)
306	c
307	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_hf
308	CALL histwrite(nid_ctesGCM,"ecrit_hf",itau_w,
309	. zx_tmp_2d,iim*jjmp1,ndex2d)
310	c
311	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_day
312	CALL histwrite(nid_ctesGCM,"ecrit_day",itau_w,
313	. zx_tmp_2d,iim*jjmp1,ndex2d)
314	c
315	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_mth
316	CALL histwrite(nid_ctesGCM,"ecrit_mth",itau_w,
317	. zx_tmp_2d,iim*jjmp1,ndex2d)
318	c
319	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_tra
320	CALL histwrite(nid_ctesGCM,"ecrit_tra",itau_w,
321	. zx_tmp_2d,iim*jjmp1,ndex2d)
322	c
323	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_reg
324	CALL histwrite(nid_ctesGCM,"ecrit_reg",itau_w,
325	. zx_tmp_2d,iim*jjmp1,ndex2d)
326	c
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)
330	c
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)
334	c
335	c=================================================================
336	c=================================================================
337	c=================================================================
338	c
339	if (ok_sync) then
340	call histsync(nid_ctesGCM)
341	endif
342	c

Note: See TracBrowser for help on using the repository browser.

Download in other formats: