source: LMDZ.3.3/trunk/libf/phylmd/write_histday.h @ 877

Last change on this file since 877 was 416, checked in by lmdzadmin, 22 years ago

Inclusion initiale
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 17.2 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
36      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
37      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
38      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
39c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
40c    .              'Surface Temperature','K')
41C
42      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
43      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
44      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
45c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
46c    .              'Surface Temperature','K')
47C
48      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
49      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
50      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
51c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
52c    .              'Surface Temperature','K')
53C
54      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
55      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
56      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
57c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
58c    .              'Surface Temperature','K')
59C
60cccIM
61      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
62      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
63c
64      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
65      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
66     .               iim*jjmp1,ndex2d)
67c
68      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
69      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
70     .               iim*jjmp1,ndex2d)
71c
72      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
74      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
75     .               ndex2d)
76c
77      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
78      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
79      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
80     .               iim*jjmp1,ndex2d)
81c
82      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
83      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
84      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
85     .               iim*jjmp1,ndex2d)
86c
87      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
88      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
89      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
90     .               ndex2d)
91c
92      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
93      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
94      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
95     .               ndex2d)
96c
97      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
98      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
99      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
100     .               ndex2d)
101c
102      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
103      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
104c
105      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
106      CALL histwrite(nid_day,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,
107     .               ndex2d)
108c
109      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
110      CALL histwrite(nid_day,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,
111     .               ndex2d)
112c
113      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
114      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
115      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
116     .               ndex2d)
117c
118      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
119      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
120      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
121     .               ndex2d)
122c
123      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
124      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
125      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
126     .               ndex2d)
127c
128      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
129      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
130      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
131     .               ndex2d)
132c
133      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
134      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
135      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
136     .               ndex2d)
137c
138      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
139      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
140      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
141     .               ndex2d)
142c
143      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
144      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
145      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
146     .               ndex2d)
147C
148      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
149      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
150      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
151     .               ndex2d)
152C
153      DO i = 1, klon
154         zx_tmp_fi2d(i) = paprs(i,1)
155      ENDDO
156      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
157      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
158c Essai writephys
159      varname = 'psol'
160      vartitle= 'pression au sol'
161      varunits= 'hPa'
162c     call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle,
163c    .                                                    varunits)
164c
165      DO i = 1, klon
166         zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
167      ENDDO
168      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
169      CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
170c     call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
171c    .              'Precipitation','mm/day')
172
173
174c
175      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
176      CALL histwrite(nid_day,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
177c     call writephy(fid_day,prof2d_av,'snow',snow_fall,
178c    .              'Snow','mm/day')
179c
180      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
181      CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
182     .               ndex2d)
183c     call writephy(fid_day,prof2d_av,'snow_mass',zxsnow,
184c    .              'Snow cover','mm')
185c
186      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
187      CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
188c     call writephy(fid_day,prof2d_av,'evap',evap,
189c    .              'Evaporation','mm/day')
190c
191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
192      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
193c     call writephy(fid_day,prof2d_av,'tops',topsw,
194c    .              'Solar rad. at TOA','W/m2')
195c
196      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
197      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
198c     call writephy(fid_day,prof2d_av,'topl',toplw,
199c    .              'IR rad. at TOA','W/m2')
200c
201      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
202      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
203c     call writephy(fid_day,prof2d_av,'sols',solsw,
204c    .              'Solar rad. at surf.','W/m2')
205c
206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
207      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
208c     call writephy(fid_day,prof2d_av,'soll',sollw,
209c    .              'IR rad. at surface','W/m2')
210c
211      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
212      CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
213     .               ndex2d)
214c     call writephy(fid_day,prof2d_av,'solldown',sollwdown,
215c    .              'Down. IR rad. at surface','W/m2')
216c
217      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
218      CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
219c     call writephy(fid_day,prof2d_av,'bils',bils,
220c    .              'Surf. total heat flux','W/m2')
221c
222      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
223      CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
224c     call writephy(fid_day,prof2d_av,'sens',sens,
225c    .              'Sensible heat flux','W/m2')
226c
227      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
228      CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
229c     call writephy(fid_day,prof2d_av,'fder',fder,
230c    .              'Heat flux derivation','W/m2')
231c
232c
233      DO nsrf = 1, nbsrf
234CXXX
235        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
236        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
237        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
238     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
239c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
240c    .                pctsrf( 1 : klon, nsrf),
241c    .                'Fraction'//clnsurf(nsrf),'-')
242C
243        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
244        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
245        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
246     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
247c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
248c    .                ftsol( 1 : klon, nsrf),
249c    .                'Surf. Temp'//clnsurf(nsrf),'K')
250C
251        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
252        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
253        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
254     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
255c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
256c    .                fluxt( 1 : klon, 1, nsrf),
257c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
258C 
259        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
260        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
261        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
262     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
263c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
264c    .                fluxlat( 1 : klon, nsrf),
265c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
266C
267        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
268        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
269        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
270     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
271c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
272c    .                fluxu( 1 : klon, 1, nsrf),
273c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
274C     
275        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
276        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
277        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
278     $      zx_tmp_2d,iim*jjmp1,ndex2d)
279c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
280c    .                fluxv( 1 : klon, 1, nsrf),
281c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
282C
283        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
284        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
285        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
286     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
287c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
288c    .                falbe( 1 : klon, nsrf),
289c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
290c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
291c    .                falblw( 1 : klon, nsrf),
292c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
293C
294        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
295        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
296        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
297     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
298c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
299c    .                frugs( 1 : klon, nsrf),
300c    .                'Rugosity '//clnsurf(nsrf),' - ')
301C
302      END DO 
303C
304cXXX      DO i = 1, klon
305cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
306cXXX      ENDDO
307cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
308cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
309c
310      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
311      CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
312c     call writephy(fid_day,prof2d_av,'cldl',cldl,
313c    .              'Low-level cloudiness','-')
314c
315      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
316      CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
317c     call writephy(fid_day,prof2d_av,'cldm',cldm,
318c    .              'Mid-level cloudiness','-')
319c
320      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
321      CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
322c     call writephy(fid_day,prof2d_av,'cldh',cldh,
323c    .              'High-level cloudiness','-')
324c
325      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
326      CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
327c     call writephy(fid_day,prof2d_av,'cldt',cldt,
328c    .              'Total cloudiness','-')
329c
330      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
331      CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
332c     call writephy(fid_day,prof2d_av,'cldq',cldq,
333c    .              'Cloud liquid water path','-')
334c
335c Champs 3D:
336c
337      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
338      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
339     .                                   iim*jjmp1*klev,ndex3d)
340c Essai writephys
341      varname = 'temp'
342      vartitle= 'temperature 3D'
343      varunits= 'K'
344c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
345c
346      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
347      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
348     .                                   iim*jjmp1*klev,ndex3d)
349c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
350c    .              'Specific humidity','Kg/Kg')
351c
352      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
353      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
354     .                                   iim*jjmp1*klev,ndex3d)
355c     call writephy(fid_day,prof3d_av,'geop',zphi,
356c    .              'Geopotential height','m')
357c
358      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
359      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
360     .                                   iim*jjmp1*klev,ndex3d)
361c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
362c    .              'Zonal wind','m/s')
363c
364      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
365      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
366     .                                   iim*jjmp1*klev,ndex3d)
367c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
368c    .              'Meridional wind','m/s')
369c
370      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
371      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
372     .                                   iim*jjmp1*klev,ndex3d)
373c     call writephy(fid_day,prof3d_av,'vitw',omega,
374c    .              'Vertical wind','m/s')
375c
376      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
377      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
378     .                                   iim*jjmp1*klev,ndex3d)
379c     call writephy(fid_day,prof3d_av,'pres',pplay,
380c    .              'Air pressure','Pa')
381cccIM
382      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
383      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
384      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
385     .                               iim*jjmp1,ndex2d)
386c
387      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
388      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
389      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
390     .                               iim*jjmp1,ndex2d)
391c
392      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
393      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
394      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
395     .                               iim*jjmp1,ndex2d)
396c
397      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
398      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
399      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
400     .                               iim*jjmp1,ndex2d)
401c
402      if (ok_sync) then
403c       call writephy_sync(fid_day)
404        call histsync(nid_day)
405      endif
406      ENDIF
Note: See TracBrowser for help on using the repository browser.