| 1 | ! |
|---|
| 2 | ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histrac.h,v 1.2 2004/10/18 15:42:20 lmdzadmin Exp $ |
|---|
| 3 | ! |
|---|
| 4 | zsto = pdtphys |
|---|
| 5 | z4julian = zjulian |
|---|
| 6 | CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian) |
|---|
| 7 | c |
|---|
| 8 | CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon) |
|---|
| 9 | DO i = 1, iim |
|---|
| 10 | zx_lon(i,1) = xlon(i+1) |
|---|
| 11 | zx_lon(i,jjm+1) = xlon(i+1) |
|---|
| 12 | ENDDO |
|---|
| 13 | CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat) |
|---|
| 14 | z4_lon = zx_lon |
|---|
| 15 | z4_lat = zx_lat |
|---|
| 16 | zpresnivs=presnivs |
|---|
| 17 | CALL histbeg("histrac", iim,z4_lon(:,1), jjm+1,z4_lat(1,:), |
|---|
| 18 | . 1,iim,1,jjm+1, itau_phy, z4julian, zsto, |
|---|
| 19 | . nhori, nid_tra) |
|---|
| 20 | CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", |
|---|
| 21 | . klev, zpresnivs, nvert) |
|---|
| 22 | |
|---|
| 23 | zout = pdtphys * FLOAT(ecrit_tra) |
|---|
| 24 | c |
|---|
| 25 | CALL histdef(nid_tra, "phis", "Surface geop. height", "-", |
|---|
| 26 | . iim,jjm+1,nhori, 1,1,1, -99, 32, |
|---|
| 27 | . "once", zsto,zout) |
|---|
| 28 | c |
|---|
| 29 | CALL histdef(nid_tra, "aire", "Grid area", "-", |
|---|
| 30 | . iim,jjm+1,nhori, 1,1,1, -99, 32, |
|---|
| 31 | . "once", zsto,zout) |
|---|
| 32 | c |
|---|
| 33 | CALL histend(nid_tra) |
|---|
| 34 | ndex2d = 0 |
|---|
| 35 | ndex3d = 0 |
|---|
| 36 | ndex = 0 |
|---|
