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

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

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 17.4 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)
14c
15      IF(lev_histday.GE.1) THEN
16c
17       CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
18         varname = 'phis'
19         vartitle= 'Surface geop. height'
20         varunits= '-'
21c
22         i = NINT(zout/zsto)
23         CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
24       CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
25         varname = 'aire'
26         vartitle= 'Grid area'
27         varunits= '-'
28c
29      DO i=1, klon
30       zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
31      ENDDO
32c
33      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
34      CALL histwrite(nid_day,"contfracATM",itau_w,
35     $      zx_tmp_2d,iim*jjmp1,ndex2d)
36c
37      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
38      CALL histwrite(nid_day,"contfracOR",itau_w,
39     $      zx_tmp_2d,iim*jjmp1,ndex2d)
40c
41      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
42      CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
43C
44      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
45      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
46c
47c   En attendant un eventuel debugage.
48
49      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
50      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
51     .               iim*jjmp1,ndex2d)
52c
53      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
54      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
55     .               iim*jjmp1,ndex2d)
56c
57      DO i = 1, klon
58         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
59      ENDDO
60      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
61      CALL histwrite(nid_day,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
62c
63      DO i = 1, klon
64         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
65      ENDDO
66      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
67      CALL histwrite(nid_day,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
68c
69      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_lsc,zx_tmp_2d)
70      CALL histwrite(nid_day,"snowl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
71c
72      IF(1.EQ.0) THEN
73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_con,zx_tmp_2d)
74      CALL histwrite(nid_day,"snowc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
75      ENDIF
76c
77      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
78      CALL histwrite(nid_day,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
79c
80      DO i = 1, klon
81         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
82      ENDDO
83c
84      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
85      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
86c
87      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
88      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
89c
90      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
91      CALL histwrite(nid_day,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,
92     .               ndex2d)
93c
94      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
95      CALL histwrite(nid_day,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,
96     .               ndex2d)
97c
98      DO i=1, klon
99       zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
100      ENDDO
101      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
102      CALL histwrite(nid_day,"wind10m",itau_w,zx_tmp_2d,
103     $     iim*jjmp1,ndex2d)
104c
105      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
106      CALL histwrite(nid_day,"wind10max",itau_w,zx_tmp_2d,
107     .               iim*jjmp1,ndex2d)
108c
109      DO i=1, klon
110         zx_tmp_fi2d(i) = paprs(i,1)
111      ENDDO
112      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
113      CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
114c
115      DO i=1, klon
116         zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
117      ENDDO
118      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
119      CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
120c
121      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
122      CALL histwrite(nid_day,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
123c
124      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
125      CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
126     .               ndex2d)
127c
128      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
129      CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
130c
131      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
132      CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
133c
134      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
135      CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
136c
137      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
138      CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
139c
140      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
141      CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
142c
143      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
144      CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
145c
146      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
147      CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
148c
149      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
150      CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
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
155      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
156      CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
157c
158      CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
159      CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
160     .               iim*jjmp1,ndex2d)
161
162c   Ecriture de champs dynamiques sur des niveaux de pression
163      DO k=1, nlevENS
164c
165        CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
166        CALL histwrite(nid_day,"u"//clev(k),itau_w,zx_tmp_2d,
167     $       iim*jjmp1,ndex2d)
168c
169        CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
170        CALL histwrite(nid_day,"v"//clev(k),itau_w,zx_tmp_2d,
171     $       iim*jjmp1,ndex2d)
172c
173      ENDDO !nlevENS
174c
175      DO k=1, nlevENS
176      IF(clev(k).EQ."500") THEN
177c
178      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,indENS(k)),zx_tmp_2d)
179      CALL histwrite(nid_day,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
180c
181c     DO i=1, klon
182c     print*,' ind, phi 500hPa',i,indENS(k),philev(:,indENS(k))
183c     print*,' ind, phi 500hPa',i,indENS(k)
184c     ENDDO
185c
186cBAD  CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,indENS(k)),zx_tmp_2d)
187      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
188      CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
189c
190      ENDIF !clev(k).EQ.500
191      ENDDO !k=1, nlevENS
192c
193      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
194      CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
195
196      CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
197      CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
198     .               iim*jjmp1,ndex2d)
199c
200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
201      CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
202     .               ndex2d)
203c
204      DO i=1, klon
205        zx_tmp_fi2d(i)=-1*sens(i)
206      ENDDO
207      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
208      CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
209c
210      zx_tmp_fi2d(1 : klon) = swdn(1 : klon, 1)
211
212      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
213      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
214     .                               iim*jjmp1,ndex2d)
215c
216      ENDIF !lev_histday.GE.1
217c
218      IF(lev_histday.GE.2) THEN
219c
220      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
221      CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
222c
223      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
224      CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
225c
226      ENDIF !lev_histday.GE.2
227c
228      IF(lev_histday.GE.3) THEN
229c=================================================================
230c   ECRITURE DES CHAMPS 3D
231c=================================================================
232
233c
234c Champs 3D:
235c
236      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
237      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
238     .                                   iim*jjmp1*klev,ndex3d)
239c
240      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
241      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
242     .                                   iim*jjmp1*klev,ndex3d)
243c
244      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
245      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
246     .                                   iim*jjmp1*klev,ndex3d)
247c
248      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
249      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
250     .                                   iim*jjmp1*klev,ndex3d)
251c
252      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
253      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
254     .                                   iim*jjmp1*klev,ndex3d)
255c
256      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
257      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
258     .                                   iim*jjmp1*klev,ndex3d)
259c
260      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
261      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
262     .                                   iim*jjmp1*klev,ndex3d)
263c
264      ENDIF !lev_histday.GE.3
265c=================================================================
266c   FIN ECRITURE DES CHAMPS 3D
267c=================================================================
268      IF(lev_histday.ge.4) THEN
269c=================================================================
270c
271c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
272c
273c=================================================================
274      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 )
275      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
276      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
277     .                               iim*jjmp1,ndex2d)
278c
279      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
280      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
281      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
282     .                               iim*jjmp1,ndex2d)
283c
284      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, klevp1 )
285      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
286      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
287     .                               iim*jjmp1,ndex2d)
288c
289      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, klevp1 )
290      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
291      CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
292     .                               iim*jjmp1,ndex2d)
293c
294      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 1 )
295      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
296      CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
297     .                               iim*jjmp1,ndex2d)
298c
299      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, klevp1 )
300      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
301      CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
302     .                               iim*jjmp1,ndex2d)
303c
304      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 1 )
305      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
306      CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
307     .                               iim*jjmp1,ndex2d)
308c
309      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
310      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
311      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
312c
313      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
314      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
315      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
316c
317      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
318      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
319      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
320c
321      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
322      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
323      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
324c
325      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
326      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
327     .               ndex2d)
328c
329      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
330      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
331      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
332     .               iim*jjmp1,ndex2d)
333c
334      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
335      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
336      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
337     .               iim*jjmp1,ndex2d)
338c
339      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
340      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
341      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
342     .               ndex2d)
343c
344      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
345      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
346      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
347     .               ndex2d)
348c
349      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
350      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
351      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
352     .               ndex2d)
353c
354      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
355      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
356      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
357     .               ndex2d)
358c
359      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
360      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
361      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
362     .               ndex2d)
363c
364      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
365      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
366      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
367     .               ndex2d)
368c
369      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
370      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
371      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
372     .               ndex2d)
373c
374      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
375      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
376      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
377     .               ndex2d)
378c
379      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
380      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
381      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
382     .               ndex2d)
383c
384      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
385      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
386      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
387     .               ndex2d)
388C
389      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
390      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
391      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
392     .               ndex2d)
393C
394      DO nsrf = 1, nbsrf
395C
396        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
397        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
398        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
399     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
400C
401        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
402        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
403        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
404     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
405C
406        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
407        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
408        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
409     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
410C 
411        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
412        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
413        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
414     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
415C
416        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
417        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
418        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
419     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
420C     
421        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
422        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
423        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
424     $      zx_tmp_2d,iim*jjmp1,ndex2d)
425C
426        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
427        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
428        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
429     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
430C
431        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
432        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
433        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
434     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
435C
436      END DO 
437c=================================================================
438c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
439c=================================================================
440      ENDIF !lev_histday.GE.4
441c
442c=================================================================
443c=================================================================
444c=================================================================
445c
446      if (ok_sync) then
447        call histsync(nid_day)
448      endif
449
450      ENDIF
Note: See TracBrowser for help on using the repository browser.