source: LMDZ4/branches/LMDZ4_par_0/libf/phylmd/write_histday.h @ 634

Last change on this file since 634 was 634, checked in by Laurent Fairhead, 19 years ago

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