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

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

Modifs sur les seuils (cdrag etc...), inclusion des diagnostics ISCCP par Ionela
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 21.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
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
162c     zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
163      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 2 )
164
165      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
166      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
167     .                               iim*jjmp1,ndex2d)
168c
169c      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
170      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
171
172      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
173      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
174     .                               iim*jjmp1,ndex2d)
175c
176c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
177      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, 2 )
178
179      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
180      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
181     .                               iim*jjmp1,ndex2d)
182c
183c      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
184      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, 1 )
185
186      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
187      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
188     .                               iim*jjmp1,ndex2d)
189
190c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, klevp1)
191      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 2 )
192
193      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
194      CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
195     .                               iim*jjmp1,ndex2d)
196
197c      zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1)
198      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 1 )
199
200      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
201      CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
202     .                               iim*jjmp1,ndex2d)
203 
204c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, klevp1)
205      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 2 )
206
207      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
208      CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
209     .                               iim*jjmp1,ndex2d)
210 
211c      zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1)
212      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 1 )
213
214      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
215      CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
216     .                               iim*jjmp1,ndex2d)
217cIM
218      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
219      CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
220     .               iim*jjmp1,ndex2d)
221
222c   Ecriture de champs dynamiques sur des niveaux de pression
223
224      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
225      CALL histwrite(nid_day,"u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
226
227      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
228      CALL histwrite(nid_day,"v850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
229
230      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
231      CALL histwrite(nid_day,"u500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
232
233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
234      CALL histwrite(nid_day,"v500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
235
236      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
237      CALL histwrite(nid_day,"u200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
238
239      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
240      CALL histwrite(nid_day,"v200",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
241
242      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
243      CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
244
245      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
246      CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
247
248      CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
249      CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
250     .               iim*jjmp1,ndex2d)
251
252
253c
254c=================================================================
255c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
256c=================================================================
257      if ( lev_histday.gt.1 ) then
258
259c   En attendant un eventuel debugage.
260
261      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
262      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
263     .               iim*jjmp1,ndex2d)
264c
265      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
266      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
267     .               iim*jjmp1,ndex2d)
268c
269c=================================================================
270c=================================================================
271c=================================================================
272c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
273c=================================================================
274
275      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
276      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
277      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
278c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
279c    .              'Surface Temperature','K')
280C
281      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
282      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
283      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
284c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
285c    .              'Surface Temperature','K')
286C
287      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
288      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
289      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
290c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
291c    .              'Surface Temperature','K')
292C
293      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
294      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
295      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
296c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
297c    .              'Surface Temperature','K')
298C
299cccIM
300      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
301      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
302     .               ndex2d)
303c
304      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
305      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
306      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
307     .               iim*jjmp1,ndex2d)
308c
309      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
310      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
311      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
312     .               iim*jjmp1,ndex2d)
313c
314      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
315      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
316      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
317     .               ndex2d)
318c
319      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
320      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
321      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
322     .               ndex2d)
323c
324      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
325      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
326      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
327     .               ndex2d)
328c
329      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
330      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
331      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
332     .               ndex2d)
333c
334      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
335      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
336      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
337     .               ndex2d)
338c
339      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
340      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
341      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
342     .               ndex2d)
343c
344      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
345      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
346      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
347     .               ndex2d)
348c
349      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
350      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
351      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
352     .               ndex2d)
353c
354      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
355      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
356      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
357     .               ndex2d)
358c
359      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
360      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
361      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
362     .               ndex2d)
363C
364      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
365      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
366      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
367     .               ndex2d)
368C
369      DO nsrf = 1, nbsrf
370CXXX
371        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
372        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
373        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
374     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
375c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
376c    .                pctsrf( 1 : klon, nsrf),
377c    .                'Fraction'//clnsurf(nsrf),'-')
378C
379        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
380        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
381        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
382     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
383c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
384c    .                ftsol( 1 : klon, nsrf),
385c    .                'Surf. Temp'//clnsurf(nsrf),'K')
386C
387        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
388        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
389        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
390     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
391c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
392c    .                fluxt( 1 : klon, 1, nsrf),
393c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
394C 
395        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
396        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
397        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
398     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
399c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
400c    .                fluxlat( 1 : klon, nsrf),
401c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
402C
403        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
404        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
405        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
406     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
407c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
408c    .                fluxu( 1 : klon, 1, nsrf),
409c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
410C     
411        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
412        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
413        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
414     $      zx_tmp_2d,iim*jjmp1,ndex2d)
415c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
416c    .                fluxv( 1 : klon, 1, nsrf),
417c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
418C
419        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
420        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
421        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
422     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
423c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
424c    .                falbe( 1 : klon, nsrf),
425c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
426c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
427c    .                falblw( 1 : klon, nsrf),
428c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
429C
430        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
431        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
432        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
433     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
434c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
435c    .                frugs( 1 : klon, nsrf),
436c    .                'Rugosity '//clnsurf(nsrf),' - ')
437C
438      END DO 
439c=================================================================
440c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
441c=================================================================
442
443c=================================================================
444c   ECRITURE DES CHAMPS 3D
445c=================================================================
446
447c
448c Champs 3D:
449c
450      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
451      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
452     .                                   iim*jjmp1*klev,ndex3d)
453c Essai writephys
454      varname = 'temp'
455      vartitle= 'temperature 3D'
456      varunits= 'K'
457c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
458c
459      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
460      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
461     .                                   iim*jjmp1*klev,ndex3d)
462c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
463c    .              'Specific humidity','Kg/Kg')
464c
465      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
466      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
467     .                                   iim*jjmp1*klev,ndex3d)
468c     call writephy(fid_day,prof3d_av,'geop',zphi,
469c    .              'Geopotential height','m')
470c
471      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
472      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
473     .                                   iim*jjmp1*klev,ndex3d)
474c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
475c    .              'Zonal wind','m/s')
476c
477      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
478      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
479     .                                   iim*jjmp1*klev,ndex3d)
480c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
481c    .              'Meridional wind','m/s')
482c
483      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
484      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
485     .                                   iim*jjmp1*klev,ndex3d)
486c     call writephy(fid_day,prof3d_av,'vitw',omega,
487c    .              'Vertical wind','m/s')
488c
489      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
490      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
491     .                                   iim*jjmp1*klev,ndex3d)
492c     call writephy(fid_day,prof3d_av,'pres',pplay,
493c    .              'Air pressure','Pa')
494cccIM
495
496c=================================================================
497c   FIN ECRITURE DES CHAMPS 3D
498c=================================================================
499
500c=================================================================
501c=================================================================
502      endif !  lev_histday.gt.1
503c=================================================================
504
505
506      if (ok_sync) then
507c       call writephy_sync(fid_day)
508        call histsync(nid_day)
509      endif
510
511      ENDIF
Note: See TracBrowser for help on using the repository browser.