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

Last change on this file since 3661 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
RevLine 
[416]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)
[506]14c
15      IF(lev_histday.GE.1) THEN
16c
[416]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= '-'
[506]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
[416]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
[506]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
[416]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
[506]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
[416]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
[506]115      DO i=1, klon
[416]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)
[506]122      CALL histwrite(nid_day,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[416]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
[506]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
[416]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
[506]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
[416]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
[506]226      ENDIF !lev_histday.GE.2
[416]227c
[506]228      IF(lev_histday.GE.3) THEN
229c=================================================================
230c   ECRITURE DES CHAMPS 3D
231c=================================================================
232
[433]233c
[506]234c Champs 3D:
[433]235c
[506]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)
[433]239c
[506]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)
[433]243c
[506]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)
[433]247c
[506]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)
[433]251c
[506]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 )
[433]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
[467]279      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
[433]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
[506]284      zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, klevp1 )
[433]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
[506]289      zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, klevp1 )
[433]290      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
[467]291      CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
292     .                               iim*jjmp1,ndex2d)
[506]293c
[467]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)
[506]298c
299      zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, klevp1 )
[467]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)
[506]303c
[467]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)
[433]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)
[506]312c
[433]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)
[506]316c
[433]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)
[506]320c
[433]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)
[506]324c
[433]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
[416]394      DO nsrf = 1, nbsrf
[506]395C
[416]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 
[433]437c=================================================================
438c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
439c=================================================================
[506]440      ENDIF !lev_histday.GE.4
[416]441c
[433]442c=================================================================
443c=================================================================
444c=================================================================
[506]445c
[416]446      if (ok_sync) then
447        call histsync(nid_day)
448      endif
[433]449
[416]450      ENDIF
Note: See TracBrowser for help on using the repository browser.