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

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

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 19.5 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
210c
211c=================================================================
212c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
213c=================================================================
214      if ( lev_histday.gt.1 ) then
215
216c   En attendant un eventuel debugage.
217
218      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
219      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
220     .               iim*jjmp1,ndex2d)
221c
222      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
223      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
224     .               iim*jjmp1,ndex2d)
225c
226c=================================================================
227c=================================================================
228c=================================================================
229c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
230c=================================================================
231
232      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
234      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
235c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
236c    .              'Surface Temperature','K')
237C
238      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
239      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
240      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
241c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
242c    .              'Surface Temperature','K')
243C
244      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
245      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
246      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
247c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
248c    .              'Surface Temperature','K')
249C
250      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
251      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
252      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
253c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
254c    .              'Surface Temperature','K')
255C
256cccIM
257      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
258      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
259     .               ndex2d)
260c
261      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
262      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
263      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
264     .               iim*jjmp1,ndex2d)
265c
266      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
267      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
268      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
269     .               iim*jjmp1,ndex2d)
270c
271      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
272      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
273      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
274     .               ndex2d)
275c
276      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
277      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
278      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
279     .               ndex2d)
280c
281      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
282      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
283      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
284     .               ndex2d)
285c
286      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
287      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
288      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
289     .               ndex2d)
290c
291      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
292      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
293      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
294     .               ndex2d)
295c
296      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
297      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
298      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
299     .               ndex2d)
300c
301      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
302      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
303      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
304     .               ndex2d)
305c
306      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
307      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
308      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
309     .               ndex2d)
310c
311      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
312      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
313      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
314     .               ndex2d)
315c
316      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
317      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
318      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
319     .               ndex2d)
320C
321      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
322      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
323      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
324     .               ndex2d)
325C
326      DO nsrf = 1, nbsrf
327CXXX
328        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
329        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
330        CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
331     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
332c       call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf),
333c    .                pctsrf( 1 : klon, nsrf),
334c    .                'Fraction'//clnsurf(nsrf),'-')
335C
336        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
337        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
338        CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
339     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
340c       call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf),
341c    .                ftsol( 1 : klon, nsrf),
342c    .                'Surf. Temp'//clnsurf(nsrf),'K')
343C
344        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
345        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
346        CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
347     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
348c       call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf),
349c    .                fluxt( 1 : klon, 1, nsrf),
350c    .                'Sensible heat flux '//clnsurf(nsrf),'W/m2')
351C 
352        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
353        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
354        CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
355     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
356c       call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf),
357c    .                fluxlat( 1 : klon, nsrf),
358c    .                'Latent heat flux '//clnsurf(nsrf),'W/m2')
359C
360        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
361        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
362        CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
363     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
364c       call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf),
365c    .                fluxu( 1 : klon, 1, nsrf),
366c    .                'Zonal wind stress '//clnsurf(nsrf),'Pa')
367C     
368        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
369        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
370        CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
371     $      zx_tmp_2d,iim*jjmp1,ndex2d)
372c       call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf),
373c    .                fluxv( 1 : klon, 1, nsrf),
374c    .                'Meridional wind stress '//clnsurf(nsrf),'Pa')
375C
376        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
377        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
378        CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
379     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
380c       call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf),
381c    .                falbe( 1 : klon, nsrf),
382c    .                'Albedo surf. SW'//clnsurf(nsrf),'-')
383c       call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf),
384c    .                falblw( 1 : klon, nsrf),
385c    .                'Albedo surf. LW'//clnsurf(nsrf),'-')
386C
387        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
388        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
389        CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
390     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
391c       call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf),
392c    .                frugs( 1 : klon, nsrf),
393c    .                'Rugosity '//clnsurf(nsrf),' - ')
394C
395      END DO 
396c=================================================================
397c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
398c=================================================================
399
400c=================================================================
401c   ECRITURE DES CHAMPS 3D
402c=================================================================
403
404c
405c Champs 3D:
406c
407      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
408      CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
409     .                                   iim*jjmp1*klev,ndex3d)
410c Essai writephys
411      varname = 'temp'
412      vartitle= 'temperature 3D'
413      varunits= 'K'
414c     call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits)
415c
416      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
417      CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
418     .                                   iim*jjmp1*klev,ndex3d)
419c     call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap),
420c    .              'Specific humidity','Kg/Kg')
421c
422      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
423      CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
424     .                                   iim*jjmp1*klev,ndex3d)
425c     call writephy(fid_day,prof3d_av,'geop',zphi,
426c    .              'Geopotential height','m')
427c
428      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
429      CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
430     .                                   iim*jjmp1*klev,ndex3d)
431c     call writephy(fid_day,prof3d_av,'vitu',u_seri,
432c    .              'Zonal wind','m/s')
433c
434      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
435      CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
436     .                                   iim*jjmp1*klev,ndex3d)
437c     call writephy(fid_day,prof3d_av,'vitv',v_seri,
438c    .              'Meridional wind','m/s')
439c
440      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
441      CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
442     .                                   iim*jjmp1*klev,ndex3d)
443c     call writephy(fid_day,prof3d_av,'vitw',omega,
444c    .              'Vertical wind','m/s')
445c
446      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
447      CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
448     .                                   iim*jjmp1*klev,ndex3d)
449c     call writephy(fid_day,prof3d_av,'pres',pplay,
450c    .              'Air pressure','Pa')
451cccIM
452
453c=================================================================
454c   FIN ECRITURE DES CHAMPS 3D
455c=================================================================
456
457c=================================================================
458c=================================================================
459      endif !  lev_histday.gt.1
460c=================================================================
461
462
463      if (ok_sync) then
464c       call writephy_sync(fid_day)
465        call histsync(nid_day)
466      endif
467
468      ENDIF
Note: See TracBrowser for help on using the repository browser.