source: LMDZ4/branches/unlabeled-1.1.1/libf/phylmd/ini_histhf3d.h @ 4667

Last change on this file since 4667 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
Line 
1!
2! $Header$
3!
4
5c     IF (ok_hf) THEN
6c
7        zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
8        zsto = zout
9        zsto1 = dtime
10        PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf
11c
12cIM cf LF
13         idayref = day_ref
14         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15
16cccIM    CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
17c        CALL ymds2ju(annee_ref, 1, 1, 0.0, zjulian)
18c        zjulian = zjulian + day_ini
19cIM cf LF
20
21c
22         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
23         DO i = 1, iim
24            zx_lon(i,1) = rlon(i+1)
25            zx_lon(i,jjmp1) = rlon(i+1)
26         ENDDO
27
28         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
29
30cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
31         CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
32     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
33     .                 nhori, nid_hf3d)
34
35         CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb",
36     .                 klev, presnivs, nvert)
37c
38
39c     IF(lev_histhf.GE.4) THEN
40c
41c Champs 3D:
42c
43         CALL histdef(nid_hf3d, "temp", "Air temperature", "K",
44     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
45     .                "ave(X)", zsto,zout)
46c
47         CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg",
48     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
49     .                "ave(X)", zsto,zout)
50c
51         CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s",
52     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
53     .                "ave(X)", zsto,zout)
54c
55         CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s",
56     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
57     .                "ave(X)", zsto,zout)
58c
59c     ENDIF
60c
61         CALL histend(nid_hf3d)
62c
63c     endif ! ok_hf
Note: See TracBrowser for help on using the repository browser.