source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histday.h @ 463

Last change on this file since 463 was 463, checked in by lmdzadmin, 21 years ago

En rajoutant des diagnostiques dans les fichiers, quelqu'un s'etait goure dans
le pas de temps d'ecriture d'ou probleme dans le calendrier interne des
fichiers.
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 20.6 KB
Line 
1      IF (ok_journe) THEN
2c
3      ndex2d = 0
4      ndex3d = 0
5c
6c Champs 2D:
7c
8         zsto = dtime
9         zout = dtime * FLOAT(ecrit_day)
10         itau_w = itau_phy + itap
11
12         i = NINT(zout/zsto)
13         CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
14       CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
15         varname = 'phis'
16         vartitle= 'Surface geop. height'
17         varunits= '-'
18c        call writephy(fid_day,prof2d_on,varname,pphis,vartitle,
19c    .                                                    varunits)
20c
21         i = NINT(zout/zsto)
22         CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
23       CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
24         varname = 'aire'
25         vartitle= 'Grid area'
26         varunits= '-'
27c        call writephy(fid_day,prof2d_on,varname,paire,vartitle,
28c    .                                                    varunits)
29C
30      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
31      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
32c     call writephy(fid_day,prof2d_av,'tsol',zxtsol,
33c    .              'Surface Temperature','K')
34c
35C
36C
37      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
38      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
39c
40      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
41      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
42c
43      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
44      CALL histwrite(nid_day,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,
45     .               ndex2d)
46c
47      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
48      CALL histwrite(nid_day,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,
49     .               ndex2d)
50c
51      DO i = 1, klon
52         zx_tmp_fi2d(i) = paprs(i,1)
53      ENDDO
54      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
55      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
56c Essai writephys
57      varname = 'psol'
58      vartitle= 'pression au sol'
59      varunits= 'hPa'
60c     call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle,
61c    .                                                    varunits)
62c
63      DO i = 1, klon
64         zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
65      ENDDO
66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
67      CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
68c     call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
69c    .              'Precipitation','mm/day')
70
71
72c
73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
74      CALL histwrite(nid_day,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
75c     call writephy(fid_day,prof2d_av,'snow',snow_fall,
76c    .              'Snow','mm/day')
77c
78      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
79      CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
80     .               ndex2d)
81c     call writephy(fid_day,prof2d_av,'snow_mass',zxsnow,
82c    .              'Snow cover','mm')
83c
84      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
85      CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
86c     call writephy(fid_day,prof2d_av,'evap',evap,
87c    .              'Evaporation','mm/day')
88c
89      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
90      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
91c     call writephy(fid_day,prof2d_av,'tops',topsw,
92c    .              'Solar rad. at TOA','W/m2')
93c
94      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
95      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
96c     call writephy(fid_day,prof2d_av,'topl',toplw,
97c    .              'IR rad. at TOA','W/m2')
98c
99      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
100      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
101c     call writephy(fid_day,prof2d_av,'sols',solsw,
102c    .              'Solar rad. at surf.','W/m2')
103c
104      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
105      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
106c     call writephy(fid_day,prof2d_av,'soll',sollw,
107c    .              'IR rad. at surface','W/m2')
108c
109      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
110      CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
111     .               ndex2d)
112c     call writephy(fid_day,prof2d_av,'solldown',sollwdown,
113c    .              'Down. IR rad. at surface','W/m2')
114c
115      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
116      CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
117c     call writephy(fid_day,prof2d_av,'bils',bils,
118c    .              'Surf. total heat flux','W/m2')
119c
120      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
121      CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
122c     call writephy(fid_day,prof2d_av,'sens',sens,
123c    .              'Sensible heat flux','W/m2')
124c
125      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
126      CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
127c     call writephy(fid_day,prof2d_av,'fder',fder,
128c    .              'Heat flux derivation','W/m2')
129c
130c
131cXXX      DO i = 1, klon
132cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
133cXXX      ENDDO
134cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
135cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
136c
137      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
138      CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
139c     call writephy(fid_day,prof2d_av,'cldl',cldl,
140c    .              'Low-level cloudiness','-')
141c
142      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
143      CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
144c     call writephy(fid_day,prof2d_av,'cldm',cldm,
145c    .              'Mid-level cloudiness','-')
146c
147      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
148      CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
149c     call writephy(fid_day,prof2d_av,'cldh',cldh,
150c    .              'High-level cloudiness','-')
151c
152      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
153      CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
154c     call writephy(fid_day,prof2d_av,'cldt',cldt,
155c    .              'Total cloudiness','-')
156c
157      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
158      CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
159c     call writephy(fid_day,prof2d_av,'cldq',cldq,
160c    .              'Cloud liquid water path','-')
161c
162      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
163      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
164      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
165     .                               iim*jjmp1,ndex2d)
166c
167      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
168      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
169      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
170     .                               iim*jjmp1,ndex2d)
171c
172      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
173      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
174      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
175     .                               iim*jjmp1,ndex2d)
176c
177      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
178      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
179      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
180     .                               iim*jjmp1,ndex2d)
181
182
183c   Ecriture de champs dynamiques sur des niveaux de pression
184
185      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
186      CALL histwrite(nid_day,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
187
188      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
189      CALL histwrite(nid_day,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
190
191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
192      CALL histwrite(nid_day,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
193
194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
195      CALL histwrite(nid_day,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
196
197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
198      CALL histwrite(nid_day,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
199
200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
201      CALL histwrite(nid_day,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
202
203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
204      CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
205
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
207      CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
208
209      CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
210      CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
211     .               iim*jjmp1,ndex2d)
212
213      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
214      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
215      CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
216     .                               iim*jjmp1,ndex2d)
217
218      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
219      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
220      CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
221     .                               iim*jjmp1,ndex2d)
222 
223      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
224      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
225      CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
226     .                               iim*jjmp1,ndex2d)
227 
228      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
229      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
230      CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
231     .                               iim*jjmp1,ndex2d)
232
233c
234c=================================================================
235c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
236c=================================================================
237      if ( lev_histday.gt.1 ) then
238
239c   En attendant un eventuel debugage.
240
241      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
242      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
243     .               iim*jjmp1,ndex2d)
244c
245      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
246      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
247     .               iim*jjmp1,ndex2d)
248c
249c=================================================================
250c=================================================================
251c=================================================================
252c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
253c=================================================================
254
255      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
256      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
257      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
258c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
259c    .              'Surface Temperature','K')
260C
261      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
262      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
263      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
264c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
265c    .              'Surface Temperature','K')
266C
267      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
268      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
269      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
270c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
271c    .              'Surface Temperature','K')
272C
273      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
274      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
275      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
276c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
277c    .              'Surface Temperature','K')
278C
279cccIM
280      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
281      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
282     .               ndex2d)
283c
284      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
285      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
286      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
287     .               iim*jjmp1,ndex2d)
288c
289      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
290      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
291      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
292     .               iim*jjmp1,ndex2d)
293c
294      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
295      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
296      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
297     .               ndex2d)
298c
299      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
300      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
301      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
302     .               ndex2d)
303c
304      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
305      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
306      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
307     .               ndex2d)
308c
309      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
310      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
311      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
312     .               ndex2d)
313c
314      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
315      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
316      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
317     .               ndex2d)
318c
319      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
320      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
321      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
322     .               ndex2d)
323c
324      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
325      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
326      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
327     .               ndex2d)
328c
329      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
330      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
331      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
332     .               ndex2d)
333c
334      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
335      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
336      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
337     .               ndex2d)
338c
339      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
340      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
341      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
342     .               ndex2d)
343C
344      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
345      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
346      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
347     .               ndex2d)
348C
349      DO nsrf = 1, nbsrf
350CXXX
351        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
352        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
353        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
354     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
355c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
356c    .                pctsrf( 1 : klon, nsrf),
357c    .                'Fraction'//clnsurf(nsrf),'-')
358C
359        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
360        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
361        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
362     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
363c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
364c    .                ftsol( 1 : klon, nsrf),
365c    .                'Surf. Temp'//clnsurf(nsrf),'K')
366C
367        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
368        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
369        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
370     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
371c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
372c    .                fluxt( 1 : klon, 1, nsrf),
373c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
374C 
375        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
376        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
377        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
378     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
379c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
380c    .                fluxlat( 1 : klon, nsrf),
381c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
382C
383        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
384        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
385        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
386     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
387c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
388c    .                fluxu( 1 : klon, 1, nsrf),
389c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
390C     
391        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
392        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
393        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
394     $      zx_tmp_2d,iim*jjmp1,ndex2d)
395c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
396c    .                fluxv( 1 : klon, 1, nsrf),
397c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
398C
399        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
400        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
401        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
402     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
403c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
404c    .                falbe( 1 : klon, nsrf),
405c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
406c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
407c    .                falblw( 1 : klon, nsrf),
408c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
409C
410        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
411        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
412        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
413     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
414c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
415c    .                frugs( 1 : klon, nsrf),
416c    .                'Rugosity '//clnsurf(nsrf),' - ')
417C
418      END DO 
419c=================================================================
420c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
421c=================================================================
422
423c=================================================================
424c   ECRITURE DES CHAMPS 3D
425c=================================================================
426
427c
428c Champs 3D:
429c
430      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
431      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
432     .                                   iim*jjmp1*klev,ndex3d)
433c Essai writephys
434      varname = 'temp'
435      vartitle= 'temperature 3D'
436      varunits= 'K'
437c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
438c
439      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
440      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
441     .                                   iim*jjmp1*klev,ndex3d)
442c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
443c    .              'Specific humidity','Kg/Kg')
444c
445      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
446      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
447     .                                   iim*jjmp1*klev,ndex3d)
448c     call writephy(fid_day,prof3d_av,'geop',zphi,
449c    .              'Geopotential height','m')
450c
451      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
452      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
453     .                                   iim*jjmp1*klev,ndex3d)
454c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
455c    .              'Zonal wind','m/s')
456c
457      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
458      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
459     .                                   iim*jjmp1*klev,ndex3d)
460c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
461c    .              'Meridional wind','m/s')
462c
463      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
464      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
465     .                                   iim*jjmp1*klev,ndex3d)
466c     call writephy(fid_day,prof3d_av,'vitw',omega,
467c    .              'Vertical wind','m/s')
468c
469      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
470      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
471     .                                   iim*jjmp1*klev,ndex3d)
472c     call writephy(fid_day,prof3d_av,'pres',pplay,
473c    .              'Air pressure','Pa')
474cccIM
475
476c=================================================================
477c   FIN ECRITURE DES CHAMPS 3D
478c=================================================================
479
480c=================================================================
481c=================================================================
482      endif !  lev_histday.gt.1
483c=================================================================
484
485
486      if (ok_sync) then
487c       call writephy_sync(fid_day)
488        call histsync(nid_day)
489      endif
490
491      ENDIF
Note: See TracBrowser for help on using the repository browser.