[1] | 1 | c |
---|
| 2 | c $Id: ini_bilKP_ave.h 1403 2010-07-01 09:02:53Z fairhead $ |
---|
| 3 | c |
---|
| 4 | IF (ok_journe) THEN |
---|
| 5 | c |
---|
| 6 | zsto = dtime |
---|
| 7 | zout = ecrit_day |
---|
| 8 | typeval=tave |
---|
| 9 | c |
---|
| 10 | idayref = day_ref |
---|
| 11 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
| 12 | c |
---|
| 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 |
---|
| 18 | DO ll=1,klev |
---|
| 19 | znivsig(ll)=REAL(ll) |
---|
| 20 | ENDDO |
---|
| 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 | cym CALL histbeg("histbilKP_ave", iim,zx_lon(:,1), jjmp1, |
---|
| 25 | cym . zx_lat(1,:), |
---|
| 26 | cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, |
---|
| 27 | cym . nhori, nid_bilKPave) |
---|
| 28 | CALL histbeg_phy("histbilKP_ave", itau_phy, zjulian, dtime, |
---|
| 29 | . nhori, nid_bilKPave) |
---|
| 30 | |
---|
| 31 | write(*,*)'Journee ', itau_phy, zjulian |
---|
| 32 | CALL histvert(nid_bilKPave, "presnivs", |
---|
| 33 | . "Vertical levels","mb", |
---|
| 34 | . klev, presnivs/100., nvert) |
---|
| 35 | c |
---|
| 36 | c |
---|
| 37 | c Champs 3D: |
---|
| 38 | c |
---|
| 39 | CALL histdef(nid_bilKPave,"ue", |
---|
| 40 | . "Zonal energy transport","-", |
---|
| 41 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 42 | . typeval, zsto,zout) |
---|
| 43 | c |
---|
| 44 | CALL histdef(nid_bilKPave,"ve", |
---|
| 45 | . "Merid energy transport","-", |
---|
| 46 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 47 | . typeval, zsto,zout) |
---|
| 48 | c |
---|
| 49 | CALL histdef(nid_bilKPave,"uq", |
---|
| 50 | . "Zonal humidity transport","-", |
---|
| 51 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 52 | . typeval, zsto,zout) |
---|
| 53 | c |
---|
| 54 | CALL histdef(nid_bilKPave,"vq", |
---|
| 55 | . "Merid humidity transport","-", |
---|
| 56 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 57 | . typeval, zsto,zout) |
---|
| 58 | c |
---|
| 59 | c Champs 3D: |
---|
| 60 | c |
---|
| 61 | CALL histdef(nid_bilKPave,"temp", |
---|
| 62 | . "Air temperature","K", |
---|
| 63 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 64 | . typeval, zsto,zout) |
---|
| 65 | c |
---|
| 66 | CALL histdef(nid_bilKPave,"ovap", |
---|
| 67 | . "Specific humidity","Kg/Kg", |
---|
| 68 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 69 | . typeval, zsto,zout) |
---|
| 70 | c |
---|
| 71 | CALL histdef(nid_bilKPave,"geop", |
---|
| 72 | . "Geopotential height","m", |
---|
| 73 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 74 | . typeval, zsto,zout) |
---|
| 75 | c |
---|
| 76 | CALL histdef(nid_bilKPave,"vitu", |
---|
| 77 | . "Zonal wind","m/s", |
---|
| 78 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 79 | . typeval, zsto,zout) |
---|
| 80 | c |
---|
| 81 | CALL histdef(nid_bilKPave,"vitv", |
---|
| 82 | . "Meridional wind","m/s", |
---|
| 83 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 84 | . typeval, zsto,zout) |
---|
| 85 | c |
---|
| 86 | CALL histdef(nid_bilKPave, "vitw", |
---|
| 87 | . "Vertical wind", "m/s", |
---|
| 88 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 89 | . typeval, zsto,zout) |
---|
| 90 | c |
---|
| 91 | CALL histdef(nid_bilKPave, "pres", |
---|
| 92 | . "Inter-Layer Air pressure", |
---|
| 93 | . "Pa", |
---|
| 94 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 95 | . typeval, zsto,zout) |
---|
| 96 | c |
---|
| 97 | CALL histdef(nid_bilKPave, "play", |
---|
| 98 | . "Mean-Layer Air pressure", |
---|
| 99 | . "Pa", |
---|
| 100 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 101 | . typeval, zsto,zout) |
---|
| 102 | c |
---|
| 103 | CALL histdef(nid_bilKPave, "oliq", |
---|
| 104 | . "Liquid water content", |
---|
| 105 | . "kg/kg", |
---|
| 106 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 107 | . typeval, zsto,zout) |
---|
| 108 | c |
---|
| 109 | CALL histdef(nid_bilKPave, "dtdyn", |
---|
| 110 | . "Dynamics dT", "K/s", |
---|
| 111 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 112 | . typeval, zsto,zout) |
---|
| 113 | c |
---|
| 114 | CALL histdef(nid_bilKPave, "dqdyn", |
---|
| 115 | . "Dynamics dQ", "Kg/Kg/s", |
---|
| 116 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 117 | . typeval, zsto,zout) |
---|
| 118 | c |
---|
| 119 | CALL histdef(nid_bilKPave, "dtcon", |
---|
| 120 | . "Convection dT", "K/s", |
---|
| 121 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 122 | . typeval, zsto,zout) |
---|
| 123 | c |
---|
| 124 | CALL histdef(nid_bilKPave, "ducon", |
---|
| 125 | . "Convection du", "m/s2", |
---|
| 126 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 127 | . typeval, zsto,zout) |
---|
| 128 | c |
---|
| 129 | CALL histdef(nid_bilKPave, "dvcon", |
---|
| 130 | . "Convection dv", "m/s2", |
---|
| 131 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 132 | . typeval, zsto,zout) |
---|
| 133 | c |
---|
| 134 | CALL histdef(nid_bilKPave,"dqcon", |
---|
| 135 | . "Convection dQ","Kg/Kg/s", |
---|
| 136 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 137 | . typeval, zsto,zout) |
---|
| 138 | c |
---|
| 139 | CALL histdef(nid_bilKPave, "dtlsc", |
---|
| 140 | . "Condensation dT", "K/s", |
---|
| 141 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 142 | . typeval, zsto,zout) |
---|
| 143 | c |
---|
| 144 | CALL histdef(nid_bilKPave,"dqlsc", |
---|
| 145 | . "Condensation dQ","Kg/Kg/s", |
---|
| 146 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 147 | . typeval, zsto,zout) |
---|
| 148 | c |
---|
| 149 | CALL histdef(nid_bilKPave,"dtvdf", |
---|
| 150 | . "Boundary-layer dT","K/s", |
---|
| 151 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 152 | . typeval, zsto,zout) |
---|
| 153 | c |
---|
| 154 | CALL histdef(nid_bilKPave, "dqvdf", |
---|
| 155 | . "Boundary-layer dQ", |
---|
| 156 | . "Kg/Kg/s", |
---|
| 157 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 158 | . typeval, zsto,zout) |
---|
| 159 | c |
---|
| 160 | CALL histdef(nid_bilKPave,"dtajs", |
---|
| 161 | . "Ajustement sec dT","K/s", |
---|
| 162 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 163 | . typeval, zsto,zout) |
---|
| 164 | c |
---|
| 165 | CALL histdef(nid_bilKPave, "dqajs", |
---|
| 166 | . "Ajustement sec dQ", |
---|
| 167 | . "Kg/Kg/s", |
---|
| 168 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 169 | . typeval, zsto,zout) |
---|
| 170 | c |
---|
| 171 | CALL histdef(nid_bilKPave, "dteva", |
---|
| 172 | . "Reevaporation dT", "K/s", |
---|
| 173 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 174 | . typeval, zsto,zout) |
---|
| 175 | c |
---|
| 176 | CALL histdef(nid_bilKPave,"dqeva", |
---|
| 177 | . "Reevaporation dQ", |
---|
| 178 | . "Kg/Kg/s", |
---|
| 179 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 180 | . typeval, zsto,zout) |
---|
| 181 | |
---|
| 182 | c |
---|
| 183 | CALL histdef(nid_bilKPave, "dtswr", |
---|
| 184 | . "SW radiation dT", "K/s", |
---|
| 185 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 186 | . typeval, zsto,zout) |
---|
| 187 | c |
---|
| 188 | CALL histdef(nid_bilKPave, "dtsw0", |
---|
| 189 | . "SW radiation dT", "K/s", |
---|
| 190 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 191 | . typeval, zsto,zout) |
---|
| 192 | c |
---|
| 193 | CALL histdef(nid_bilKPave, "dtlwr", |
---|
| 194 | . "LW radiation dT", "K/s", |
---|
| 195 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 196 | . typeval, zsto,zout) |
---|
| 197 | c |
---|
| 198 | CALL histdef(nid_bilKPave, "dtlw0", |
---|
| 199 | . "LW radiation dT", "K/s", |
---|
| 200 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 201 | . typeval, zsto,zout) |
---|
| 202 | c |
---|
| 203 | CALL histdef(nid_bilKPave,"duvdf", |
---|
| 204 | . "Boundary-layer dU","m/s2", |
---|
| 205 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 206 | . typeval, zsto,zout) |
---|
| 207 | c |
---|
| 208 | CALL histdef(nid_bilKPave,"dvvdf", |
---|
| 209 | . "Boundary-layer dV","m/s2", |
---|
| 210 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 211 | . typeval, zsto,zout) |
---|
| 212 | c |
---|
| 213 | IF (ok_orodr) THEN |
---|
| 214 | IF (ok_orolf) THEN |
---|
| 215 | CALL histdef(nid_bilKPave, "duoli", |
---|
| 216 | . "Orography dU", "m/s2", |
---|
| 217 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 218 | . typeval, zsto,zout) |
---|
| 219 | c |
---|
| 220 | CALL histdef(nid_bilKPave, "dvoli", |
---|
| 221 | . "Orography dV", "m/s2", |
---|
| 222 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 223 | . typeval, zsto,zout) |
---|
| 224 | c |
---|
| 225 | ENDIF |
---|
| 226 | ENDIF |
---|
| 227 | C |
---|
| 228 | CALL histdef(nid_bilKPave, "duphy", |
---|
| 229 | . "Physiq dU","-", |
---|
| 230 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 231 | . typeval, zsto,zout) |
---|
| 232 | C |
---|
| 233 | CALL histdef(nid_bilKPave, "dvphy", |
---|
| 234 | . "Physiq dV","-", |
---|
| 235 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 236 | . typeval, zsto,zout) |
---|
| 237 | C |
---|
| 238 | CALL histdef(nid_bilKPave, "dtphy", |
---|
| 239 | . "Physiq dT","-", |
---|
| 240 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 241 | . typeval, zsto,zout) |
---|
| 242 | C |
---|
| 243 | CALL histdef(nid_bilKPave, "dqphy", |
---|
| 244 | . "Physiq dQ","-", |
---|
| 245 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 246 | . typeval, zsto,zout) |
---|
| 247 | C |
---|
| 248 | CALL histdef(nid_bilKPave, "dqlphy", |
---|
| 249 | . "Physiq dQl","-", |
---|
| 250 | . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 251 | . typeval, zsto,zout) |
---|
| 252 | C |
---|
| 253 | C |
---|
| 254 | CALL histend(nid_bilKPave) |
---|
| 255 | c |
---|
| 256 | ndex2d = 0 |
---|
| 257 | ndex3d = 0 |
---|
| 258 | c |
---|
| 259 | ENDIF ! fin de test sur ok_journe |
---|