- Timestamp:
- Jun 4, 2007, 4:39:50 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ini_histhf3d.h
r760 r767 1 c2 1 c $Header$ 3 2 c … … 12 11 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 13 12 c 14 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 15 DO i = 1, iim 16 zx_lon(i,1) = rlon(i+1) 17 zx_lon(i,jjmp1) = rlon(i+1) 18 ENDDO 19 c 20 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 21 CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 22 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 13 14 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 15 cym DO i = 1, iim 16 cym zx_lon(i,1) = rlon(i+1) 17 cym zx_lon(i,jjmp1) = rlon(i+1) 18 cym ENDDO 19 20 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 21 22 cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat, 23 cym CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 24 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 25 cym . nhori, nid_hf3d) 26 CALL histbeg_phy("histhf3d", itau_phy, zjulian, dtime, 23 27 . nhori, nid_hf3d) 24 c 28 25 29 CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb", 26 30 . klev, presnivs/100., nvert) … … 29 33 c 30 34 CALL histdef(nid_hf3d, "temp", "Air temperature", "K", 31 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,32 . "inst(X)", zstohf, zstohf)35 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 36 . "inst(X)", zstohf,zout) 33 37 c 34 38 CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg", 35 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,36 . "inst(X)", zstohf, zstohf)39 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 40 . "inst(X)", zstohf,zout) 37 41 c 38 42 CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s", 39 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,40 . "inst(X)", zstohf, zstohf)43 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 44 . "inst(X)", zstohf,zout) 41 45 c 42 46 CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s", 43 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,44 . "inst(X)", zstohf, zstohf)47 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 48 . "inst(X)", zstohf,zout) 45 49 c 46 50 CALL histend(nid_hf3d)
Note: See TracChangeset
for help on using the changeset viewer.