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

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

Ajout cape_max, SWup, SWdn au TOA et SFC ciel clair IM
  • 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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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",itap,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.