Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/ini_histmthNMC.h
- Timestamp:
- Aug 17, 2006, 5:41:51 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/V3_test/libf/phylmd/ini_histmthNMC.h
r684 r704 12 12 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 13 13 c 14 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)15 DO i = 1, iim16 zx_lon(i,1) = rlon(i+1)17 zx_lon(i,jjmp1) = rlon(i+1)18 ENDDO14 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 19 DO ll=1,klev 20 20 znivsig(ll)=float(ll) 21 21 ENDDO 22 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 23 CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 24 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 22 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 23 cym CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 24 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 25 cym . nhori, nid_nmc) 26 27 CALL histbeg_phy("histNMC",itau_phy, zjulian, dtime, 25 28 . nhori, nid_nmc) 26 29 c … … 34 37 CALL histdef(nid_nmc, "temp", 35 38 . "Temperature","K", 36 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,39 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 37 40 . "inst(X)", zout,zout) 38 41 c 39 42 CALL histdef(nid_nmc, "phi", 40 43 . "Geopotential", "m", 41 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,44 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 42 45 . "inst(X)", zout,zout) 43 46 c 44 47 CALL histdef(nid_nmc, "q", 45 48 . "Specific humidity","kg/kg", 46 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,49 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 47 50 . "inst(X)", zout,zout) 48 51 c 49 52 CALL histdef(nid_nmc, "rh", 50 53 . "Relative humidity", "%", 51 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,54 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 52 55 . "inst(X)", zout,zout) 53 56 c 54 57 CALL histdef(nid_nmc, "u", 55 58 . "Zonal wind","m/s", 56 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,59 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 57 60 . "inst(X)", zout,zout) 58 61 c 59 62 CALL histdef(nid_nmc, "v", 60 63 . "Meridional wind","m/s", 61 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,64 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 62 65 . "inst(X)", zout,zout) 63 66 c … … 73 76 CALL histdef(nid_nmc, "temp", 74 77 . "Temperature","K", 75 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,78 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 76 79 . "inst(X)", zout,zout) 77 80 c 78 81 CALL histdef(nid_nmc, "phi", 79 82 . "Geopotential ", "m", 80 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,83 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 81 84 . "inst(X)", zout,zout) 82 85 c 83 86 CALL histdef(nid_nmc, "q", 84 87 . "Specific humidity","kg/kg", 85 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,88 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 86 89 . "inst(X)", zout,zout) 87 90 c 88 91 CALL histdef(nid_nmc, "rh", 89 92 . "Relative humidity", "%", 90 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,93 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 91 94 . "inst(X)", zout,zout) 92 95 c 93 96 CALL histdef(nid_nmc, "u", 94 97 . "Zonal wind","m/s", 95 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,98 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 96 99 . "inst(X)", zout,zout) 97 100 c 98 101 CALL histdef(nid_nmc, "v", 99 102 . "Meridional wind","m/s", 100 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,103 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 101 104 . "inst(X)", zout,zout) 102 105 c 103 106 CALL histdef(nid_nmc, "w", 104 107 . "Vertical motion","Pa/s", 105 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,108 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 106 109 . "inst(X)", zout,zout) 107 110 c … … 116 119 CALL histdef(nid_nmc, "psbg", 117 120 . "Pressure sfce below ground","%", 118 . iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,121 . iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 119 122 . "inst(X)", zout,zout) 120 123 c 121 124 CALL histdef(nid_nmc, "uv", 122 125 . "uv ", 123 . "m2/s2",iim,jj mp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,126 . "m2/s2",iim,jjphy_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32, 124 127 . "inst(X)", zout,zout) 125 128 c 126 129 CALL histdef(nid_nmc, "vq", 127 130 . "vq ", 128 . "m/s * (kg/kg)",iim,jj mp1,nhori,131 . "m/s * (kg/kg)",iim,jjphy_nb,nhori, 129 132 . nlevSTD,1,nlevSTD, nvert, 32, 130 133 . "inst(X)", zout,zout) … … 132 135 CALL histdef(nid_nmc, "vT", 133 136 . "vT ", 134 . "mK/s",iim,jj mp1,nhori,137 . "mK/s",iim,jjphy_nb,nhori, 135 138 . nlevSTD,1,nlevSTD, nvert, 32, 136 139 . "inst(X)", zout,zout) … … 138 141 CALL histdef(nid_nmc, "wq", 139 142 . "wq ", 140 . "(Pa/s)*(kg/kg)",iim,jj mp1,nhori,143 . "(Pa/s)*(kg/kg)",iim,jjphy_nb,nhori, 141 144 . nlevSTD,1,nlevSTD, nvert, 32, 142 145 . "inst(X)", zout,zout) … … 144 147 CALL histdef(nid_nmc, "vphi", 145 148 . "vphi ", 146 . "m2/s",iim,jj mp1,nhori,149 . "m2/s",iim,jjphy_nb,nhori, 147 150 . nlevSTD,1,nlevSTD, nvert, 32, 148 151 . "inst(X)", zout,zout) … … 150 153 CALL histdef(nid_nmc, "wT", 151 154 . "wT ", 152 . "K*Pa/s",iim,jj mp1,nhori,155 . "K*Pa/s",iim,jjphy_nb,nhori, 153 156 . nlevSTD,1,nlevSTD, nvert, 32, 154 157 . "inst(X)", zout,zout) … … 156 159 CALL histdef(nid_nmc, "uxu", 157 160 . "u2 ", 158 . "m2/s2",iim,jj mp1,nhori,161 . "m2/s2",iim,jjphy_nb,nhori, 159 162 . nlevSTD,1,nlevSTD, nvert, 32, 160 163 . "inst(X)", zout,zout) … … 162 165 CALL histdef(nid_nmc, "vxv", 163 166 . "v2 ", 164 . "m2/s2",iim,jj mp1,nhori,167 . "m2/s2",iim,jjphy_nb,nhori, 165 168 . nlevSTD,1,nlevSTD, nvert, 32, 166 169 . "inst(X)", zout,zout) … … 168 171 CALL histdef(nid_nmc, "TxT", 169 172 . "T2 ", 170 . "K2",iim,jj mp1,nhori,173 . "K2",iim,jjphy_nb,nhori, 171 174 . nlevSTD,1,nlevSTD, nvert, 32, 172 175 . "inst(X)", zout,zout)
Note: See TracChangeset
for help on using the changeset viewer.