source: LMDZ4/trunk/libf/phylmd/ini_histhf3d.h @ 1940

Last change on this file since 1940 was 776, checked in by Laurent Fairhead, 18 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.7 KB
Line 
1c $Header$
2c
3c sorties hf 3d
4c
5        zstohf = ecrit_hf
6        zout = ecrit_hf
7c
8c       PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf
9c
10        idayref = day_ref
11        CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12c
13
14cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
15cym         DO i = 1, iim
16cym            zx_lon(i,1) = rlon(i+1)
17cym            zx_lon(i,jjmp1) = rlon(i+1)
18cym         ENDDO
19
20cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
21
22cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
23cym         CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
24cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
25cym     .                 nhori, nid_hf3d)
26         CALL histbeg_phy("histhf3d", itau_phy, zjulian, dtime, 
27     .                 nhori, nid_hf3d)
28
29        CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb",
30     .                 klev, presnivs/100., nvert)
31c
32c Champs 3D:
33c
34        CALL histdef(nid_hf3d, "temp", "Air temperature", "K",
35     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
36     .                "ave(X)", zstohf,zout)
37c
38        CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg",
39     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
40     .                "ave(X)", zstohf,zout)
41c
42        CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s",
43     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
44     .                "ave(X)", zstohf,zout)
45c
46        CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s",
47     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
48     .                "ave(X)", zstohf,zout)
49c
50        CALL histend(nid_hf3d)
51c
Note: See TracBrowser for help on using the repository browser.