source: LMDZ4/trunk/libf/dyn3d/write_paramLMDZ_dyn.h @ 923

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

Probleme de calcul des pas de temps d'ecriture dans le fichier paramLMDZ_dyn.nc
LF

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 7.6 KB
Line 
1c
2      ndex2d = 0
3      itau_w=itau_dyn+itau
4c
5      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(prt_level) 
6      CALL histwrite(nid_ctesGCM, "prt_level", itau_w,
7     .               zx_tmp_2d,iip1*jjp1,ndex2d)
8c
9      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(dayref)
10      CALL histwrite(nid_ctesGCM, "dayref", itau_w,
11     .               zx_tmp_2d,iip1*jjp1,ndex2d)
12c
13      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(anneeref)
14      CALL histwrite(nid_ctesGCM, "anneeref", itau_w,
15     .               zx_tmp_2d,iip1*jjp1,ndex2d)
16c
17      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(raz_date)
18      CALL histwrite(nid_ctesGCM, "raz_date", itau_w,
19     .               zx_tmp_2d,iip1*jjp1,ndex2d)
20c
21      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nday)
22      CALL histwrite(nid_ctesGCM, "nday", itau_w,
23     .               zx_tmp_2d,iip1*jjp1,ndex2d)
24c
25      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(day_step)
26      CALL histwrite(nid_ctesGCM, "day_step", itau_w,
27     .               zx_tmp_2d,iip1*jjp1,ndex2d)
28c
29      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iperiod)
30      CALL histwrite(nid_ctesGCM, "iperiod", itau_w,
31     .               zx_tmp_2d,iip1*jjp1,ndex2d)
32c
33      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iapp_tracvl)
34      CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_w,
35     .               zx_tmp_2d,iip1*jjp1,ndex2d)
36c
37      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iconser)
38      CALL histwrite(nid_ctesGCM, "iconser", itau_w,
39     .               zx_tmp_2d,iip1*jjp1,ndex2d)
40c
41      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iecri)
42      CALL histwrite(nid_ctesGCM, "iecri", itau_w,
43     .               zx_tmp_2d,iip1*jjp1,ndex2d)
44c
45      zx_tmp_2d(1:iip1,1:jjp1)=periodav
46      CALL histwrite(nid_ctesGCM, "periodav", itau_w,
47     .               zx_tmp_2d,iip1*jjp1,ndex2d)
48c
49      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(idissip)
50      CALL histwrite(nid_ctesGCM, "idissip", itau_w,
51     .               zx_tmp_2d,iip1*jjp1,ndex2d)
52c
53      IF(lstardis) THEN
54       zx_tmp_2d(1:iip1,1:jjp1)=1.
55      ELSE
56       zx_tmp_2d(1:iip1,1:jjp1)=0.
57      ENDIF
58      CALL histwrite(nid_ctesGCM, "lstardis", itau_w,
59     .               zx_tmp_2d,iip1*jjp1,ndex2d)
60c
61      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergdiv)
62      CALL histwrite(nid_ctesGCM, "nitergdiv", itau_w,
63     .               zx_tmp_2d,iip1*jjp1,ndex2d)
64c
65      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergrot)
66      CALL histwrite(nid_ctesGCM, "nitergrot", itau_w,
67     .               zx_tmp_2d,iip1*jjp1,ndex2d)
68c
69      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(niterh) 
70      CALL histwrite(nid_ctesGCM, "niterh", itau_w,
71     .               zx_tmp_2d,iip1*jjp1,ndex2d)
72c
73      zx_tmp_2d(1:iip1,1:jjp1)=tetagdiv
74      CALL histwrite(nid_ctesGCM, "tetagdiv", itau_w,
75     .               zx_tmp_2d,iip1*jjp1,ndex2d)
76c
77      zx_tmp_2d(1:iip1,1:jjp1)=tetagrot
78      CALL histwrite(nid_ctesGCM, "tetagrot", itau_w,
79     .               zx_tmp_2d,iip1*jjp1,ndex2d)
80c
81      zx_tmp_2d(1:iip1,1:jjp1)=tetatemp
82      CALL histwrite(nid_ctesGCM, "tetatemp", itau_w,
83     .               zx_tmp_2d,iip1*jjp1,ndex2d)
84c
85      zx_tmp_2d(1:iip1,1:jjp1)=coefdis
86      CALL histwrite(nid_ctesGCM, "coefdis", itau_w,
87     .               zx_tmp_2d,iip1*jjp1,ndex2d)
88c
89      IF(purmats) THEN
90       zx_tmp_2d(1:iip1,1:jjp1)=1.
91      ELSE
92       zx_tmp_2d(1:iip1,1:jjp1)=0.
93      ENDIF
94      CALL histwrite(nid_ctesGCM, "purmats", itau_w,
95     .               zx_tmp_2d,iip1*jjp1,ndex2d)
96c
97      IF(ok_guide) THEN
98       zx_tmp_2d(1:iip1,1:jjp1)=1.
99      ELSE
100       zx_tmp_2d(1:iip1,1:jjp1)=0.
101      ENDIF
102      CALL histwrite(nid_ctesGCM, "ok_guide", itau_w,
103     .               zx_tmp_2d,iip1*jjp1,ndex2d)
104c
105      IF(true_calendar) THEN
106       zx_tmp_2d(1:iip1,1:jjp1)=1.
107      ELSE
108       zx_tmp_2d(1:iip1,1:jjp1)=0.
109      ENDIF
110      CALL histwrite(nid_ctesGCM, "true_calendar", itau_w,
111     .               zx_tmp_2d,iip1*jjp1,ndex2d)
112c
113      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iflag_phys)
114      CALL histwrite(nid_ctesGCM, "iflag_phys", itau_w,
115     .               zx_tmp_2d,iip1*jjp1,ndex2d)
116c
117      zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iphysiq)
118      CALL histwrite(nid_ctesGCM, "iphysiq", itau_w,
119     .               zx_tmp_2d,iip1*jjp1,ndex2d)
120c
121      zx_tmp_2d(1:iip1,1:jjp1)=ecritphy
122      CALL histwrite(nid_ctesGCM, "ecritphy", itau_w,
123     .               zx_tmp_2d,iip1*jjp1,ndex2d)
124c
125      IF(cycle_diurne) THEN
126       zx_tmp_2d(1:iip1,1:jjp1)=1.
127      ELSE
128       zx_tmp_2d(1:iip1,1:jjp1)=0.
129      ENDIF
130      CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_w,
131     .               zx_tmp_2d,iip1*jjp1,ndex2d)
132c
133      IF(soil_model) THEN
134       zx_tmp_2d(1:iip1,1:jjp1)=1.
135      ELSE
136       zx_tmp_2d(1:iip1,1:jjp1)=0.
137      ENDIF
138      CALL histwrite(nid_ctesGCM, "soil_model", itau_w,
139     .               zx_tmp_2d,iip1*jjp1,ndex2d)
140c
141      IF(new_oliq) THEN
142       zx_tmp_2d(1:iip1,1:jjp1)=1.
143      ELSE
144       zx_tmp_2d(1:iip1,1:jjp1)=0.
145      ENDIF
146      CALL histwrite(nid_ctesGCM, "new_oliq", itau_w,
147     .               zx_tmp_2d,iip1*jjp1,ndex2d)
148c
149      IF(ok_orodr) THEN
150       zx_tmp_2d(1:iip1,1:jjp1)=1.
151      ELSE
152       zx_tmp_2d(1:iip1,1:jjp1)=0.
153      ENDIF
154      CALL histwrite(nid_ctesGCM, "ok_orodr", itau_w,
155     .               zx_tmp_2d,iip1*jjp1,ndex2d)
156c
157      IF(ok_orolf) THEN
158       zx_tmp_2d(1:iip1,1:jjp1)=1.
159      ELSE
160       zx_tmp_2d(1:iip1,1:jjp1)=0.
161      ENDIF
162      CALL histwrite(nid_ctesGCM, "ok_orolf", itau_w,
163     .               zx_tmp_2d,iip1*jjp1,ndex2d)
164c
165      IF(ok_limitvrai) THEN
166       zx_tmp_2d(1:iip1,1:jjp1)=1.
167      ELSE
168       zx_tmp_2d(1:iip1,1:jjp1)=0.
169      ENDIF
170      CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_w,
171     .               zx_tmp_2d,iip1*jjp1,ndex2d)
172c
173      zx_tmp_2d(1:iip1,1:jjp1)=nbapp_rad
174      CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_w,
175     .               zx_tmp_2d,iip1*jjp1,ndex2d)
176c
177      zx_tmp_2d(1:iip1,1:jjp1)=iflag_con
178      CALL histwrite(nid_ctesGCM, "iflag_con", itau_w,
179     .               zx_tmp_2d,iip1*jjp1,ndex2d)
180c
181      zx_tmp_2d(1:iip1,1:jjp1)=clon
182      CALL histwrite(nid_ctesGCM, "clon", itau_w,
183     .               zx_tmp_2d,iip1*jjp1,ndex2d)
184c
185      zx_tmp_2d(1:iip1,1:jjp1)=clat
186      CALL histwrite(nid_ctesGCM, "clat", itau_w,
187     .               zx_tmp_2d,iip1*jjp1,ndex2d)
188c
189      zx_tmp_2d(1:iip1,1:jjp1)=grossismx
190      CALL histwrite(nid_ctesGCM, "grossismx", itau_w,
191     .               zx_tmp_2d,iip1*jjp1,ndex2d)
192c
193      zx_tmp_2d(1:iip1,1:jjp1)=grossismy
194      CALL histwrite(nid_ctesGCM, "grossismy", itau_w,
195     .               zx_tmp_2d,iip1*jjp1,ndex2d)
196c
197      IF(fxyhypb) THEN
198       zx_tmp_2d(1:iip1,1:jjp1)=1.
199      ELSE
200       zx_tmp_2d(1:iip1,1:jjp1)=0.
201      ENDIF
202      CALL histwrite(nid_ctesGCM, "fxyhypb", itau_w,
203     .               zx_tmp_2d,iip1*jjp1,ndex2d)
204c
205      zx_tmp_2d(1:iip1,1:jjp1)=dzoomx
206      CALL histwrite(nid_ctesGCM, "dzoomx", itau_w,
207     .               zx_tmp_2d,iip1*jjp1,ndex2d)
208c
209      zx_tmp_2d(1:iip1,1:jjp1)=dzoomy
210      CALL histwrite(nid_ctesGCM, "dzoomy", itau_w,
211     .               zx_tmp_2d,iip1*jjp1,ndex2d)
212c
213      zx_tmp_2d(1:iip1,1:jjp1)=taux
214      CALL histwrite(nid_ctesGCM, "taux", itau_w,
215     .               zx_tmp_2d,iip1*jjp1,ndex2d)
216c
217      zx_tmp_2d(1:iip1,1:jjp1)=tauy
218      CALL histwrite(nid_ctesGCM, "tauy", itau_w,
219     .               zx_tmp_2d,iip1*jjp1,ndex2d)
220c
221      IF(ysinus) THEN
222       zx_tmp_2d(1:iip1,1:jjp1)=1.
223      ELSE
224       zx_tmp_2d(1:iip1,1:jjp1)=0.
225      ENDIF
226      CALL histwrite(nid_ctesGCM, "ysinus", itau_w,
227     .               zx_tmp_2d,iip1*jjp1,ndex2d)
228c
229      zx_tmp_2d(1:iip1,1:jjp1)=ip_ebil_dyn
230      CALL histwrite(nid_ctesGCM, "ip_ebil_dyn", itau_w,
231     .               zx_tmp_2d,iip1*jjp1,ndex2d)
232c
233c=================================================================
234c
235      if (ok_sync) then
236        call histsync(nid_ctesGCM)
237      endif
238c
239c=================================================================
Note: See TracBrowser for help on using the repository browser.