Context Navigation

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

Last change on this file since 900 was 828, checked in by lmdzadmin, 17 years ago
Ajout parametres "fmagic" reglage albedo, "ok_hf", frequence calcul "freq_ISCCP" et sortie "ecrit_ISCCP" ISCCP 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	IF (is_sequential) THEN
3
4	ndex2d = 0
5	itau_w = itau_phy + itap
6	c
7	c Variables type caractere : plusieurs valeurs possibles
8	c
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)
18	c
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)
26	c
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)
34	c
35	c Variables logiques (1=true, 2=false)
36	c
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)
44	c
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)
52	c
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)
60	c
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)
68	c
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)
76	c
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)
84	c
85
86	c
87	c Champs 2D:
88	c
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)
92	c
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)
96	c
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)
100	c
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)
104	c
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)
108	c
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)
112	c
113	zx_tmp_2d(1:iim,1:jjmp1)=solaire
114	CALL histwrite(nid_ctesGCM,"solaire",itau_w,
115	. zx_tmp_2d,iim*jjmp1,ndex2d)
116	c
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)
120	c
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)
124	c
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)
128	c
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)
132	c
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)
136	c
137	zx_tmp_2d(1:iim,1:jjmp1)=epmax
138	CALL histwrite(nid_ctesGCM,"epmax",itau_w,
139	. zx_tmp_2d,iim*jjmp1,ndex2d)
140	c
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)
148	c
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)
152	c
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)
156	c
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)
160	c
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)
164	c
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)
168	c
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)
172	c
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)
176	c
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)
180	c
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)
188	c
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)
192	c
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)
196	c
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)
200	c
201	zx_tmp_2d(1:iim,1:jjmp1)=facttemps
202	CALL histwrite(nid_ctesGCM,"facttemps",itau_w,
203	. zx_tmp_2d,iim*jjmp1,ndex2d)
204	c
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)
212	c
213	zx_tmp_2d(1:iim,1:jjmp1)=ratqsbas
214	CALL histwrite(nid_ctesGCM,"ratqsbas",itau_w,
215	. zx_tmp_2d,iim*jjmp1,ndex2d)
216	c
217	zx_tmp_2d(1:iim,1:jjmp1)=ratqshaut
218	CALL histwrite(nid_ctesGCM,"ratqshaut",itau_w,
219	. zx_tmp_2d,iim*jjmp1,ndex2d)
220	c
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)
224	c
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)
228	c
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)
232	c
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)
236	c
237	zx_tmp_2d(1:iim,1:jjmp1)=overlap
238	CALL histwrite(nid_ctesGCM,"overlap",itau_w,
239	. zx_tmp_2d,iim*jjmp1,ndex2d)
240	c
241	zx_tmp_2d(1:iim,1:jjmp1)=cdmmax
242	CALL histwrite(nid_ctesGCM,"cdmmax",itau_w,
243	. zx_tmp_2d,iim*jjmp1,ndex2d)
244	c
245	zx_tmp_2d(1:iim,1:jjmp1)=cdhmax
246	CALL histwrite(nid_ctesGCM,"cdhmax",itau_w,
247	. zx_tmp_2d,iim*jjmp1,ndex2d)
248	c
249	zx_tmp_2d(1:iim,1:jjmp1)=ksta
250	CALL histwrite(nid_ctesGCM,"ksta",itau_w,
251	. zx_tmp_2d,iim*jjmp1,ndex2d)
252	c
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)
256	c
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)
264	c
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)
268	c
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)
272	c
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)
276	c
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)
280	c
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)
288	c
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)
292	c
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)
296	c
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)
300	c
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)
304	c
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)
308	c
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)
312	c
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)
316	c
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)
320	c
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)
324	c
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)
328	c
329	zx_tmp_2d(1:iim,1:jjmp1)=freq_ISCCP
330	CALL histwrite(nid_ctesGCM,"freq_ISCCP",itau_w,
331	. zx_tmp_2d,iim*jjmp1,ndex2d)
332	c
333	zx_tmp_2d(1:iim,1:jjmp1)=ecrit_ISCCP
334	CALL histwrite(nid_ctesGCM,"ecrit_ISCCP",itau_w,
335	. zx_tmp_2d,iim*jjmp1,ndex2d)
336	c
337	c=================================================================
338	c=================================================================
339	c=================================================================
340	c
341	if (ok_sync) then
342	call histsync(nid_ctesGCM)
343	endif
344	c
345	ENDIF ! mono_cpu

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

Download in other formats: