[644] | 1 | c |
---|
[1403] | 2 | c $Id: ini_bilKP_ins.h 1907 2013-11-26 13:10:46Z acaubel $ |
---|
[644] | 3 | c |
---|
| 4 | IF (ok_journe) THEN |
---|
| 5 | c |
---|
| 6 | zsto = dtime |
---|
| 7 | zout = dtime |
---|
| 8 | typeval=tinst |
---|
| 9 | c |
---|
| 10 | idayref = day_ref |
---|
| 11 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
| 12 | c |
---|
[766] | 13 | 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 |
---|
[644] | 18 | DO ll=1,klev |
---|
[1403] | 19 | znivsig(ll)=REAL(ll) |
---|
[644] | 20 | ENDDO |
---|
[766] | 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,:) |
---|
[644] | 24 | c |
---|
| 25 | cIM 280405 BEG |
---|
| 26 | c |
---|
| 27 | cIM cf. AM 081204 BEG region |
---|
| 28 | imin_ins=1 |
---|
| 29 | imax_ins=iim |
---|
| 30 | jmin_ins=1 |
---|
| 31 | jmax_ins=jjmp1 |
---|
[766] | 32 | 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 |
---|
[644] | 40 | c |
---|
| 41 | print*,'On stoke le fichier bilKP instantanne sur ', |
---|
| 42 | s imin_ins,imax_ins,jmin_ins,jmax_ins |
---|
| 43 | print*,'On stoke le fichier bilKP instantanne sur ', |
---|
| 44 | s zx_lon(imin_ins,1),zx_lon(imax_ins,1), |
---|
| 45 | s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins) |
---|
| 46 | cIM cf. AM 081204 END region |
---|
| 47 | c |
---|
| 48 | cIM 280405 END |
---|
| 49 | c |
---|
[766] | 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) |
---|
[644] | 55 | ENDIF |
---|
| 56 | c |
---|
| 57 | cIM 280405 BEG |
---|
| 58 | c |
---|
| 59 | cIM cf. AM 081204 BEG region |
---|
[766] | 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, |
---|
[644] | 67 | . nhori, nid_bilKPins) |
---|
| 68 | cIM 081204 END |
---|
| 69 | c |
---|
| 70 | cIM 280405 END |
---|
| 71 | c |
---|
| 72 | write(*,*)'Journee ', itau_phy, zjulian |
---|
| 73 | CALL histvert(nid_bilKPins, "presnivs", |
---|
| 74 | . "Vertical levels","mb", |
---|
| 75 | . klev, presnivs/100., nvert) |
---|
| 76 | c |
---|
| 77 | c Champs 3D: |
---|
| 78 | c |
---|
| 79 | CALL histdef(nid_bilKPins,"ue", |
---|
| 80 | . "Zonal energy transport","-", |
---|
[766] | 81 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 82 | . typeval, zsto,zout) |
---|
| 83 | c |
---|
| 84 | CALL histdef(nid_bilKPins,"ve", |
---|
| 85 | . "Merid energy transport","-", |
---|
[766] | 86 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 87 | . typeval, zsto,zout) |
---|
| 88 | c |
---|
| 89 | CALL histdef(nid_bilKPins,"uq", |
---|
| 90 | . "Zonal humidity transport","-", |
---|
[766] | 91 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 92 | . typeval, zsto,zout) |
---|
| 93 | c |
---|
| 94 | CALL histdef(nid_bilKPins,"vq", |
---|
| 95 | . "Merid humidity transport","-", |
---|
[766] | 96 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 97 | . typeval, zsto,zout) |
---|
| 98 | c |
---|
| 99 | c Champs 3D: |
---|
| 100 | c |
---|
| 101 | CALL histdef(nid_bilKPins, "temp", |
---|
| 102 | . "Air temperature", "K", |
---|
[766] | 103 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 104 | . typeval, zsto,zout) |
---|
| 105 | c |
---|
| 106 | CALL histdef(nid_bilKPins,"ovap", |
---|
| 107 | . "Specific humidity","Kg/Kg", |
---|
[766] | 108 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 109 | . typeval, zsto,zout) |
---|
| 110 | c |
---|
| 111 | CALL histdef(nid_bilKPins,"geop", |
---|
| 112 | . "Geopotential height", "m", |
---|
[766] | 113 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 114 | . typeval, zsto,zout) |
---|
| 115 | c |
---|
| 116 | CALL histdef(nid_bilKPins,"vitu", |
---|
| 117 | . "Zonal wind", "m/s", |
---|
[766] | 118 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 119 | . typeval, zsto,zout) |
---|
| 120 | c |
---|
| 121 | CALL histdef(nid_bilKPins,"vitv", |
---|
| 122 | . "Meridional wind", "m/s", |
---|
[766] | 123 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 124 | . typeval, zsto,zout) |
---|
| 125 | c |
---|
| 126 | CALL histdef(nid_bilKPins, "vitw", |
---|
| 127 | . "Vertical wind", "m/s", |
---|
[766] | 128 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 129 | . typeval, zsto,zout) |
---|
| 130 | c |
---|
| 131 | CALL histdef(nid_bilKPins, "pres", |
---|
| 132 | . "Inter-Layer Air pressure", |
---|
| 133 | . "Pa", |
---|
[766] | 134 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 135 | . typeval, zsto,zout) |
---|
| 136 | c |
---|
| 137 | CALL histdef(nid_bilKPins, "play", |
---|
| 138 | . "Mean-Layer Air pressure", |
---|
| 139 | . "Pa", |
---|
[766] | 140 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 141 | . typeval, zsto,zout) |
---|
| 142 | c |
---|
| 143 | CALL histdef(nid_bilKPins, "oliq", |
---|
| 144 | . "Liquid water content", |
---|
| 145 | . "kg/kg", |
---|
[766] | 146 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 147 | . typeval, zsto,zout) |
---|
| 148 | c |
---|
| 149 | CALL histdef(nid_bilKPins, "dtdyn", |
---|
| 150 | . "Dynamics dT", "K/s", |
---|
[766] | 151 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 152 | . typeval, zsto,zout) |
---|
| 153 | c |
---|
| 154 | CALL histdef(nid_bilKPins, "dqdyn", |
---|
| 155 | . "Dynamics dQ", "Kg/Kg/s", |
---|
[766] | 156 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 157 | . typeval, zsto,zout) |
---|
| 158 | c |
---|
| 159 | CALL histdef(nid_bilKPins, "dtcon", |
---|
| 160 | . "Convection dT", "K/s", |
---|
[766] | 161 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 162 | . typeval, zsto,zout) |
---|
| 163 | c |
---|
| 164 | CALL histdef(nid_bilKPins, "ducon", |
---|
| 165 | . "Convection du", "m/s2", |
---|
[766] | 166 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 167 | . typeval, zsto,zout) |
---|
| 168 | c |
---|
| 169 | CALL histdef(nid_bilKPins, "dvcon", |
---|
| 170 | . "Convection dv", "m/s2", |
---|
[766] | 171 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 172 | . typeval, zsto,zout) |
---|
| 173 | c |
---|
| 174 | CALL histdef(nid_bilKPins,"dqcon", |
---|
| 175 | . "Convection dQ","Kg/Kg/s", |
---|
[766] | 176 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 177 | . typeval, zsto,zout) |
---|
| 178 | c |
---|
| 179 | CALL histdef(nid_bilKPins, "dtlsc", |
---|
| 180 | . "Condensation dT", "K/s", |
---|
[766] | 181 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 182 | . typeval, zsto,zout) |
---|
| 183 | c |
---|
| 184 | CALL histdef(nid_bilKPins,"dqlsc", |
---|
| 185 | . "Condensation dQ","Kg/Kg/s", |
---|
[766] | 186 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 187 | . typeval, zsto,zout) |
---|
| 188 | c |
---|
| 189 | CALL histdef(nid_bilKPins,"dtvdf", |
---|
| 190 | . "Boundary-layer dT","K/s", |
---|
[766] | 191 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 192 | . typeval, zsto,zout) |
---|
| 193 | c |
---|
| 194 | CALL histdef(nid_bilKPins, "dqvdf", |
---|
| 195 | . "Boundary-layer dQ", |
---|
| 196 | . "Kg/Kg/s", |
---|
[766] | 197 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 198 | . typeval, zsto,zout) |
---|
| 199 | c |
---|
| 200 | CALL histdef(nid_bilKPins,"dtajs", |
---|
| 201 | . "Ajustement sec dT","K/s", |
---|
[766] | 202 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 203 | . typeval, zsto,zout) |
---|
| 204 | c |
---|
| 205 | CALL histdef(nid_bilKPins,"dqajs", |
---|
| 206 | . "Ajustement sec dQ", |
---|
| 207 | . "Kg/Kg/s", |
---|
[766] | 208 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 209 | . typeval, zsto,zout) |
---|
| 210 | c |
---|
| 211 | CALL histdef(nid_bilKPins,"dteva", |
---|
| 212 | . "Reevaporation dT","K/s", |
---|
[766] | 213 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 214 | . typeval, zsto,zout) |
---|
| 215 | c |
---|
| 216 | CALL histdef(nid_bilKPins,"dqeva", |
---|
| 217 | . "Reevaporation dQ", |
---|
| 218 | . "Kg/Kg/s", |
---|
[766] | 219 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 220 | . typeval, zsto,zout) |
---|
| 221 | |
---|
| 222 | c |
---|
| 223 | CALL histdef(nid_bilKPins, "dtswr", |
---|
| 224 | . "SW radiation dT", "K/s", |
---|
[766] | 225 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 226 | . typeval, zsto,zout) |
---|
| 227 | c |
---|
| 228 | CALL histdef(nid_bilKPins, "dtsw0", |
---|
| 229 | . "SW radiation dT", "K/s", |
---|
[766] | 230 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 231 | . typeval, zsto,zout) |
---|
| 232 | c |
---|
| 233 | CALL histdef(nid_bilKPins, "dtlwr", |
---|
| 234 | . "LW radiation dT", "K/s", |
---|
[766] | 235 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 236 | . typeval, zsto,zout) |
---|
| 237 | c |
---|
| 238 | CALL histdef(nid_bilKPins, "dtlw0", |
---|
| 239 | . "LW radiation dT", "K/s", |
---|
[766] | 240 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 241 | . typeval, zsto,zout) |
---|
| 242 | c |
---|
| 243 | CALL histdef(nid_bilKPins,"duvdf", |
---|
| 244 | . "Boundary-layer dU","m/s2", |
---|
[766] | 245 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 246 | . typeval, zsto,zout) |
---|
| 247 | c |
---|
| 248 | CALL histdef(nid_bilKPins,"dvvdf", |
---|
| 249 | . "Boundary-layer dV","m/s2", |
---|
[766] | 250 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 251 | . typeval, zsto,zout) |
---|
| 252 | c |
---|
| 253 | IF (ok_orodr) THEN |
---|
| 254 | IF (ok_orolf) THEN |
---|
| 255 | CALL histdef(nid_bilKPins, "duoli", |
---|
| 256 | . "Orography dU", "m/s2", |
---|
[766] | 257 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 258 | . typeval, zsto,zout) |
---|
| 259 | c |
---|
| 260 | CALL histdef(nid_bilKPins, "dvoli", |
---|
| 261 | . "Orography dV", "m/s2", |
---|
[766] | 262 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 263 | . typeval, zsto,zout) |
---|
| 264 | c |
---|
| 265 | ENDIF |
---|
| 266 | ENDIF |
---|
| 267 | C |
---|
| 268 | CALL histdef(nid_bilKPins, "duphy", |
---|
| 269 | . "Physiq dU","-", |
---|
[766] | 270 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 271 | . typeval, zsto,zout) |
---|
| 272 | C |
---|
| 273 | CALL histdef(nid_bilKPins, "dvphy", |
---|
| 274 | . "Physiq dV","-", |
---|
[766] | 275 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 276 | . typeval, zsto,zout) |
---|
| 277 | C |
---|
| 278 | CALL histdef(nid_bilKPins, "dtphy", |
---|
| 279 | . "Physiq dT","-", |
---|
[766] | 280 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 281 | . typeval, zsto,zout) |
---|
| 282 | C |
---|
| 283 | CALL histdef(nid_bilKPins, "dqphy", |
---|
| 284 | . "Physiq dQ","-", |
---|
[766] | 285 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 286 | . typeval, zsto,zout) |
---|
| 287 | C |
---|
| 288 | CALL histdef(nid_bilKPins, "dqlphy", |
---|
| 289 | . "Physiq dQl","-", |
---|
[766] | 290 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 291 | . typeval, zsto,zout) |
---|
| 292 | cIM 280405 BEG |
---|
| 293 | c |
---|
| 294 | c Champs 2D: |
---|
| 295 | c |
---|
| 296 | c u850, v850 |
---|
| 297 | c DO k=1, nlevSTD |
---|
| 298 | DO k=1, 12 |
---|
| 299 | c |
---|
| 300 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 301 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 302 | c |
---|
| 303 | IF(bb2.EQ."850") THEN |
---|
| 304 | c |
---|
| 305 | CALL histdef(nid_bilKPins, "u"//bb2, |
---|
| 306 | . "Zonal wind "//bb2//"mb","m/s", |
---|
[766] | 307 | . iim,jjphy_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 308 | . typeval, zsto,zout) |
---|
| 309 | c |
---|
| 310 | CALL histdef(nid_bilKPins, "v"//bb2, |
---|
| 311 | . "Meridional wind "//bb2//"mb","m/s", |
---|
[766] | 312 | . iim,jjphy_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 313 | . typeval, zsto,zout) |
---|
| 314 | c |
---|
| 315 | ENDIF !(bb2.EQ."850") |
---|
| 316 | c |
---|
| 317 | ENDDO !k=1, 12 |
---|
| 318 | c |
---|
| 319 | cIM 280405 END |
---|
| 320 | c |
---|
| 321 | CALL histend(nid_bilKPins) |
---|
| 322 | c |
---|
| 323 | ndex2d = 0 |
---|
| 324 | ndex3d = 0 |
---|
| 325 | c |
---|
| 326 | ENDIF ! fin de test sur ok_journe |
---|