Context Navigation

ini_histhf3d.h @ 735

Last change on this file since 735 was 723, checked in by lmdzadmin, 19 years ago
On passe a des ecrit_ins, ecrit_day, etc en nombre de jours (REAL) On lit frequence ecriture traceurs ecrit_trac dans physiq.def Correction petits pbs ini_histrac.h IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 1.7 KB

Line
1	c
2	c $Header$
3	c
4	c sorties hf 3d
5	c
6	cIM 130904 zstohf = dtime * REAL(NINT(86400./dtime*ecrit_hf))
7	cIM 130904 zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
8	cIM zstohf = dtime * FLOAT(ecrit_hf)
9	cIM zout = dtime * FLOAT(ecrit_hf)
10	zstohf = dtime * ecrit_hf
11	zout = dtime * ecrit_hf
12	c
13	c PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf
14	c
15	idayref = day_ref
16	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
17	c
18	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
19	DO i = 1, iim
20	zx_lon(i,1) = rlon(i+1)
21	zx_lon(i,jjmp1) = rlon(i+1)
22	ENDDO
23	c
24	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
25	CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
26	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
27	. nhori, nid_hf3d)
28	c
29	CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb",
30	. klev, presnivs/100., nvert)
31	c
32	c Champs 3D:
33	c
34	CALL histdef(nid_hf3d, "temp", "Air temperature", "K",
35	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
36	. "ave(X)", zstohf,zout)
37	c
38	CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg",
39	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
40	. "ave(X)", zstohf,zout)
41	c
42	CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s",
43	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
44	. "ave(X)", zstohf,zout)
45	c
46	CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s",
47	. iim,jjmp1,nhori, klev,1,klev,nvert, 32,
48	. "ave(X)", zstohf,zout)
49	c
50	CALL histend(nid_hf3d)
51	c

Note: See TracBrowser for help on using the repository browser.