Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/ini_bilKP_ins.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_bilKP_ins.h
r684 r704 11 11 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 12 12 c 13 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)14 DO i = 1, iim15 zx_lon(i,1) = rlon(i+1)16 zx_lon(i,jjmp1) = rlon(i+1)17 ENDDO13 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 14 cym DO i = 1, iim 15 cym zx_lon(i,1) = rlon(i+1) 16 cym zx_lon(i,jjmp1) = rlon(i+1) 17 cym ENDDO 18 18 DO ll=1,klev 19 19 znivsig(ll)=float(ll) 20 20 ENDDO 21 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)22 write(*,*)'zx_lon = ',zx_lon(:,1)23 write(*,*)'zx_lat = ',zx_lat(1,:)21 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 22 cym write(*,*)'zx_lon = ',zx_lon(:,1) 23 cym write(*,*)'zx_lat = ',zx_lat(1,:) 24 24 c 25 25 cIM 280405 BEG … … 30 30 jmin_ins=1 31 31 jmax_ins=jjmp1 32 do i=1,iim-133 if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i34 if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+135 enddo36 do j=1,jjmp137 if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j38 if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j39 enddo32 cym do i=1,iim-1 33 cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i 34 cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1 35 cym enddo 36 cym do j=1,jjmp1 37 cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j 38 cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j 39 cym enddo 40 40 c 41 41 print*,'On stoke le fichier bilKP instantanne sur ', … … 48 48 cIM 280405 END 49 49 c 50 IF(1.EQ.0) THEN51 CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,52 . zx_lat(1,:),53 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,54 . nhori, nid_bilKPins)50 cym IF(1.EQ.0) THEN 51 cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1, 52 cym . zx_lat(1,:), 53 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 54 cym . nhori, nid_bilKPins) 55 55 ENDIF 56 56 c … … 58 58 c 59 59 cIM cf. AM 081204 BEG region 60 CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), 61 . jjmp1,zx_lat(1,:), 62 . imin_ins,imax_ins-imin_ins+1, 63 . jmin_ins,jmax_ins-jmin_ins+1, 64 . itau_phy, zjulian, dtime, 60 cym CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), 61 cym . jjmp1,zx_lat(1,:), 62 cym . imin_ins,imax_ins-imin_ins+1, 63 cym . jmin_ins,jmax_ins-jmin_ins+1, 64 cym . itau_phy, zjulian, dtime, 65 cym . nhori, nid_bilKPins) 66 CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime, 65 67 . nhori, nid_bilKPins) 66 68 cIM 081204 END … … 77 79 CALL histdef(nid_bilKPins,"ue", 78 80 . "Zonal energy transport","-", 79 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,81 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 80 82 . typeval, zsto,zout) 81 83 c 82 84 CALL histdef(nid_bilKPins,"ve", 83 85 . "Merid energy transport","-", 84 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,86 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 85 87 . typeval, zsto,zout) 86 88 c 87 89 CALL histdef(nid_bilKPins,"uq", 88 90 . "Zonal humidity transport","-", 89 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,91 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 90 92 . typeval, zsto,zout) 91 93 c 92 94 CALL histdef(nid_bilKPins,"vq", 93 95 . "Merid humidity transport","-", 94 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,96 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 95 97 . typeval, zsto,zout) 96 98 c … … 99 101 CALL histdef(nid_bilKPins, "temp", 100 102 . "Air temperature", "K", 101 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,103 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 102 104 . typeval, zsto,zout) 103 105 c 104 106 CALL histdef(nid_bilKPins,"ovap", 105 107 . "Specific humidity","Kg/Kg", 106 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,108 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 107 109 . typeval, zsto,zout) 108 110 c 109 111 CALL histdef(nid_bilKPins,"geop", 110 112 . "Geopotential height", "m", 111 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,113 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 112 114 . typeval, zsto,zout) 113 115 c 114 116 CALL histdef(nid_bilKPins,"vitu", 115 117 . "Zonal wind", "m/s", 116 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,118 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 117 119 . typeval, zsto,zout) 118 120 c 119 121 CALL histdef(nid_bilKPins,"vitv", 120 122 . "Meridional wind", "m/s", 121 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,123 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 122 124 . typeval, zsto,zout) 123 125 c 124 126 CALL histdef(nid_bilKPins, "vitw", 125 127 . "Vertical wind", "m/s", 126 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,128 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 127 129 . typeval, zsto,zout) 128 130 c … … 130 132 . "Inter-Layer Air pressure", 131 133 . "Pa", 132 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,134 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 133 135 . typeval, zsto,zout) 134 136 c … … 136 138 . "Mean-Layer Air pressure", 137 139 . "Pa", 138 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,140 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 139 141 . typeval, zsto,zout) 140 142 c … … 142 144 . "Liquid water content", 143 145 . "kg/kg", 144 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,146 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 145 147 . typeval, zsto,zout) 146 148 c 147 149 CALL histdef(nid_bilKPins, "dtdyn", 148 150 . "Dynamics dT", "K/s", 149 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,151 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 150 152 . typeval, zsto,zout) 151 153 c 152 154 CALL histdef(nid_bilKPins, "dqdyn", 153 155 . "Dynamics dQ", "Kg/Kg/s", 154 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,156 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 155 157 . typeval, zsto,zout) 156 158 c 157 159 CALL histdef(nid_bilKPins, "dtcon", 158 160 . "Convection dT", "K/s", 159 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,161 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 160 162 . typeval, zsto,zout) 161 163 c 162 164 CALL histdef(nid_bilKPins, "ducon", 163 165 . "Convection du", "m/s2", 164 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,166 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 165 167 . typeval, zsto,zout) 166 168 c 167 169 CALL histdef(nid_bilKPins, "dvcon", 168 170 . "Convection dv", "m/s2", 169 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,171 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 170 172 . typeval, zsto,zout) 171 173 c 172 174 CALL histdef(nid_bilKPins,"dqcon", 173 175 . "Convection dQ","Kg/Kg/s", 174 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,176 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 175 177 . typeval, zsto,zout) 176 178 c 177 179 CALL histdef(nid_bilKPins, "dtlsc", 178 180 . "Condensation dT", "K/s", 179 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,181 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 180 182 . typeval, zsto,zout) 181 183 c 182 184 CALL histdef(nid_bilKPins,"dqlsc", 183 185 . "Condensation dQ","Kg/Kg/s", 184 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,186 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 185 187 . typeval, zsto,zout) 186 188 c 187 189 CALL histdef(nid_bilKPins,"dtvdf", 188 190 . "Boundary-layer dT","K/s", 189 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,191 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 190 192 . typeval, zsto,zout) 191 193 c … … 193 195 . "Boundary-layer dQ", 194 196 . "Kg/Kg/s", 195 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,197 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 196 198 . typeval, zsto,zout) 197 199 c 198 200 CALL histdef(nid_bilKPins,"dtajs", 199 201 . "Ajustement sec dT","K/s", 200 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,202 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 201 203 . typeval, zsto,zout) 202 204 c … … 204 206 . "Ajustement sec dQ", 205 207 . "Kg/Kg/s", 206 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,208 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 207 209 . typeval, zsto,zout) 208 210 c 209 211 CALL histdef(nid_bilKPins,"dteva", 210 212 . "Reevaporation dT","K/s", 211 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,213 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 212 214 . typeval, zsto,zout) 213 215 c … … 215 217 . "Reevaporation dQ", 216 218 . "Kg/Kg/s", 217 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,219 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 218 220 . typeval, zsto,zout) 219 221 … … 221 223 CALL histdef(nid_bilKPins, "dtswr", 222 224 . "SW radiation dT", "K/s", 223 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,225 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 224 226 . typeval, zsto,zout) 225 227 c 226 228 CALL histdef(nid_bilKPins, "dtsw0", 227 229 . "SW radiation dT", "K/s", 228 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,230 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 229 231 . typeval, zsto,zout) 230 232 c 231 233 CALL histdef(nid_bilKPins, "dtlwr", 232 234 . "LW radiation dT", "K/s", 233 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,235 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 234 236 . typeval, zsto,zout) 235 237 c 236 238 CALL histdef(nid_bilKPins, "dtlw0", 237 239 . "LW radiation dT", "K/s", 238 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,240 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 239 241 . typeval, zsto,zout) 240 242 c 241 243 CALL histdef(nid_bilKPins,"duvdf", 242 244 . "Boundary-layer dU","m/s2", 243 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,245 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 244 246 . typeval, zsto,zout) 245 247 c 246 248 CALL histdef(nid_bilKPins,"dvvdf", 247 249 . "Boundary-layer dV","m/s2", 248 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,250 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 249 251 . typeval, zsto,zout) 250 252 c … … 253 255 CALL histdef(nid_bilKPins, "duoli", 254 256 . "Orography dU", "m/s2", 255 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,257 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 256 258 . typeval, zsto,zout) 257 259 c 258 260 CALL histdef(nid_bilKPins, "dvoli", 259 261 . "Orography dV", "m/s2", 260 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,262 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 261 263 . typeval, zsto,zout) 262 264 c … … 266 268 CALL histdef(nid_bilKPins, "duphy", 267 269 . "Physiq dU","-", 268 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,270 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 269 271 . typeval, zsto,zout) 270 272 C 271 273 CALL histdef(nid_bilKPins, "dvphy", 272 274 . "Physiq dV","-", 273 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,275 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 274 276 . typeval, zsto,zout) 275 277 C 276 278 CALL histdef(nid_bilKPins, "dtphy", 277 279 . "Physiq dT","-", 278 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,280 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 279 281 . typeval, zsto,zout) 280 282 C 281 283 CALL histdef(nid_bilKPins, "dqphy", 282 284 . "Physiq dQ","-", 283 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,285 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 284 286 . typeval, zsto,zout) 285 287 C 286 288 CALL histdef(nid_bilKPins, "dqlphy", 287 289 . "Physiq dQl","-", 288 . iim,jj mp1,nhori, klev,1,klev,nvert, 32,290 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 289 291 . typeval, zsto,zout) 290 292 cIM 280405 BEG … … 303 305 CALL histdef(nid_bilKPins, "u"//bb2, 304 306 . "Zonal wind "//bb2//"mb","m/s", 305 . iim,jj mp1,nhori, 1,1,1, -99, 32,307 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 306 308 . typeval, zsto,zout) 307 309 c 308 310 CALL histdef(nid_bilKPins, "v"//bb2, 309 311 . "Meridional wind "//bb2//"mb","m/s", 310 . iim,jj mp1,nhori, 1,1,1, -99, 32,312 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 311 313 . typeval, zsto,zout) 312 314 c
Note: See TracChangeset
for help on using the changeset viewer.