[644] | 1 | c |
---|
| 2 | c $Header$ |
---|
| 3 | c |
---|
[524] | 4 | IF (ok_journe) THEN |
---|
[766] | 5 | |
---|
| 6 | c$OMP MASTER |
---|
[524] | 7 | c |
---|
[644] | 8 | zstophy = dtime |
---|
[684] | 9 | zstoday= ecrit_day |
---|
| 10 | zout = ecrit_day |
---|
[524] | 11 | c |
---|
| 12 | idayref = day_ref |
---|
| 13 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
| 14 | c |
---|
[766] | 15 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) |
---|
| 16 | cym DO i = 1, iim |
---|
| 17 | cym zx_lon(i,1) = rlon(i+1) |
---|
| 18 | cym zx_lon(i,jjmp1) = rlon(i+1) |
---|
| 19 | cym ENDDO |
---|
[524] | 20 | DO ll=1,klev |
---|
| 21 | znivsig(ll)=float(ll) |
---|
| 22 | ENDDO |
---|
[766] | 23 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) |
---|
| 24 | cym CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), |
---|
| 25 | cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, |
---|
| 26 | cym . nhori, nid_day) |
---|
| 27 | |
---|
| 28 | CALL histbeg_phy("histday",itau_phy, zjulian, dtime, |
---|
| 29 | . nhori, nid_day) |
---|
| 30 | |
---|
| 31 | write(*,*)'Journee ', itau_phy, zjulian |
---|
[524] | 32 | CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", |
---|
[644] | 33 | . klev, presnivs/100., nvert) |
---|
[524] | 34 | c |
---|
| 35 | IF(lev_histday.GE.1) THEN |
---|
| 36 | c |
---|
[644] | 37 | CALL histdef(nid_day, "phis", "Surface geop. height", "m2/s2", |
---|
[776] | 38 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 39 | . "once", zstophy,zout) |
---|
[524] | 40 | c |
---|
[918] | 41 | CALL histdef(nid_day, "aire", "Grid area", "m2", |
---|
[776] | 42 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 43 | . "once", zstophy,zout) |
---|
[524] | 44 | c |
---|
[918] | 45 | CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","1", |
---|
[776] | 46 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 47 | . "once", zstophy,zout) |
---|
[524] | 48 | c |
---|
[918] | 49 | CALL histdef(nid_day, "contfracOR","% sfce terre OR", "1", |
---|
[776] | 50 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 51 | . "once", zstophy,zout) |
---|
[524] | 52 | c |
---|
| 53 | c Champs 2D: |
---|
| 54 | c |
---|
| 55 | CALL histdef(nid_day, "tsol", "Surface Temperature", "K", |
---|
[776] | 56 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 57 | . "ave(X)", zstophy,zout) |
---|
[524] | 58 | c |
---|
| 59 | CALL histdef(nid_day, "t2m", "Temperature 2m", "K", |
---|
[776] | 60 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 61 | . "ave(X)", zstophy,zout) |
---|
[524] | 62 | c |
---|
[766] | 63 | c Champs retires momentannement en attendant un hypothetique |
---|
[524] | 64 | c debugage |
---|
| 65 | |
---|
| 66 | CALL histdef(nid_day, "t2m_min", "Temp. 2m min.", |
---|
| 67 | . "K", |
---|
[776] | 68 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 69 | . t2mincels, zstophy,zout) |
---|
[524] | 70 | c |
---|
| 71 | CALL histdef(nid_day, "t2m_max", "Temp. 2m max.", |
---|
| 72 | . "K", |
---|
[776] | 73 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 74 | . t2maxcels, zstophy,zout) |
---|
[524] | 75 | c |
---|
| 76 | CALL histdef(nid_day, "plul", "Large-scale Precip.", |
---|
| 77 | . "kg/(s*m2)", |
---|
[776] | 78 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 79 | . "ave(X)", zstophy,zout) |
---|
[524] | 80 | c |
---|
| 81 | CALL histdef(nid_day, "pluc", "Convective Precip.", |
---|
| 82 | . "kg/(s*m2)", |
---|
[776] | 83 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 84 | . "ave(X)", zstophy,zout) |
---|
[524] | 85 | c |
---|
| 86 | CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.", |
---|
[644] | 87 | . "kg/(m2*s)", |
---|
[776] | 88 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 89 | . "ave(X)", zstophy,zout) |
---|
[524] | 90 | c |
---|
| 91 | IF(1.EQ.0) THEN !snowc=0. |
---|
| 92 | CALL histdef(nid_day, "snowc", "Solid Convective Precip.", |
---|
[644] | 93 | . "kg/(m2*s)", |
---|
[776] | 94 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 95 | . "ave(X)", zstophy,zout) |
---|
[524] | 96 | ENDIF !snowc=0. |
---|
| 97 | c |
---|
| 98 | CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2", |
---|
[776] | 99 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 100 | . "ave(X)", zstophy,zout) |
---|
[524] | 101 | c |
---|
[918] | 102 | CALL histdef(nid_day, "sicf", "Sea-ice fraction", "1", |
---|
[776] | 103 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 104 | . "inst(X)", zstoday,zout) |
---|
[524] | 105 | c |
---|
| 106 | CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg", |
---|
[776] | 107 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 108 | . "ave(X)", zstophy,zout) |
---|
[524] | 109 | c |
---|
| 110 | CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s", |
---|
[776] | 111 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 112 | . "ave(X)", zstophy,zout) |
---|
[524] | 113 | c |
---|
| 114 | CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s", |
---|
[776] | 115 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 116 | . "ave(X)", zstophy,zout) |
---|
[524] | 117 | c |
---|
| 118 | CALL histdef(nid_day, "wind10m","10-m wind speed","m/s", |
---|
[776] | 119 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 120 | . "ave(X)", zstophy,zout) |
---|
[524] | 121 | c |
---|
| 122 | CALL histdef(nid_day, "wind10max", "10-m wind speed max.", |
---|
| 123 | . "m/s", |
---|
[776] | 124 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 125 | . t2maxcels, zstophy,zout) |
---|
[524] | 126 | c |
---|
| 127 | CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", |
---|
[776] | 128 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 129 | . "ave(X)", zstophy,zout) |
---|
[524] | 130 | c |
---|
| 131 | CALL histdef(nid_day, "precip","Precipitation Totale liq+sol" |
---|
| 132 | . , "kg/(s*m2)", |
---|
[776] | 133 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 134 | . "ave(X)", zstophy,zout) |
---|
[524] | 135 | c |
---|
[644] | 136 | CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)", |
---|
[776] | 137 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 138 | . "ave(X)", zstophy,zout) |
---|
[524] | 139 | c |
---|
[644] | 140 | c CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2", |
---|
[776] | 141 | c . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 142 | c . "ave(X)", zstophy,zout) |
---|
[524] | 143 | c |
---|
[644] | 144 | CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)", |
---|
[776] | 145 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 146 | . "ave(X)", zstophy,zout) |
---|
[524] | 147 | c |
---|
| 148 | CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", |
---|
[776] | 149 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 150 | . "ave(X)", zstophy,zout) |
---|
[524] | 151 | c |
---|
| 152 | CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", |
---|
[776] | 153 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 154 | . "ave(X)", zstophy,zout) |
---|
[524] | 155 | c |
---|
| 156 | CALL histdef(nid_day, "sols", "Net Solar rad. at surf.", |
---|
| 157 | . "W/m2", |
---|
[776] | 158 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 159 | . "ave(X)", zstophy,zout) |
---|
[524] | 160 | c |
---|
| 161 | CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2", |
---|
[776] | 162 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 163 | . "ave(X)", zstophy,zout) |
---|
[524] | 164 | c |
---|
[684] | 165 | CALL histdef(nid_day, "radsol", "Rayonnement au sol", "W/m2", |
---|
[776] | 166 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[684] | 167 | . "ave(X)", zstophy,zout) |
---|
| 168 | c |
---|
[918] | 169 | CALL histdef(nid_day, "cldl", "Low-level cloudiness", "1", |
---|
[776] | 170 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 171 | . "ave(X)", zstophy,zout) |
---|
[524] | 172 | c |
---|
[918] | 173 | CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "1", |
---|
[776] | 174 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 175 | . "ave(X)", zstophy,zout) |
---|
[524] | 176 | c |
---|
[918] | 177 | CALL histdef(nid_day, "cldh", "High-level cloudiness", "1", |
---|
[776] | 178 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 179 | . "ave(X)", zstophy,zout) |
---|
[524] | 180 | c |
---|
[918] | 181 | CALL histdef(nid_day, "cldt", "Total cloudiness", "1", |
---|
[776] | 182 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 183 | . "ave(X)", zstophy,zout) |
---|
[524] | 184 | c |
---|
| 185 | CALL histdef(nid_day, "cldq", "Cloud liquid water path", |
---|
| 186 | . "kg/m2", |
---|
[776] | 187 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 188 | . "ave(X)", zstophy,zout) |
---|
[524] | 189 | c |
---|
| 190 | CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2", |
---|
[776] | 191 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 192 | . "ave(X)", zstophy,zout) |
---|
[524] | 193 | c |
---|
[897] | 194 | CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa", |
---|
| 195 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 196 | . "ave(X)", zstophy,zout) |
---|
| 197 | c |
---|
| 198 | CALL histdef(nid_day, "cape_max", "CAPE max.", |
---|
| 199 | . "J/kg", |
---|
| 200 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 201 | . capemaxcels, zstophy,zout) |
---|
| 202 | c |
---|
| 203 | CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", |
---|
| 204 | . "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 205 | . "ave(X)", zstophy,zout) |
---|
| 206 | c |
---|
| 207 | CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", |
---|
| 208 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 209 | . "ave(X)", zstophy,zout) |
---|
| 210 | c |
---|
| 211 | CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2", |
---|
| 212 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
| 213 | . 32, "ave(X)", zstophy,zout) |
---|
| 214 | c |
---|
[524] | 215 | c Champs dynamiques sur niveaux de pression |
---|
| 216 | c |
---|
[644] | 217 | c DO k=1, nlevENS |
---|
| 218 | DO k=1, nlevSTD |
---|
| 219 | c |
---|
| 220 | c bb=clevSTD(k) |
---|
| 221 | c |
---|
| 222 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 223 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 224 | c aa=clevSTD(k) |
---|
| 225 | c bb2=clevSTD(k) |
---|
| 226 | c ENDIF |
---|
| 227 | c |
---|
| 228 | IF(bb2.EQ."850".OR.bb2.EQ."700" |
---|
| 229 | . .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
| 230 | c |
---|
| 231 | CALL histdef(nid_day, "u"//bb2, |
---|
| 232 | . "Zonal wind "//bb2//"mb","m/s", |
---|
[776] | 233 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 234 | . "inst(X)", zout,zout) |
---|
[524] | 235 | c |
---|
[644] | 236 | CALL histdef(nid_day, "v"//bb2, |
---|
| 237 | . "Meridional wind "//bb2//"mb","m/s", |
---|
[776] | 238 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 239 | . "inst(X)", zout,zout) |
---|
[524] | 240 | c |
---|
[644] | 241 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
[524] | 242 | c |
---|
[644] | 243 | IF(bb2.EQ."500") THEN |
---|
| 244 | c |
---|
| 245 | CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb", |
---|
| 246 | . "m/s", |
---|
[776] | 247 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 248 | . "inst(X)", zout,zout) |
---|
[524] | 249 | c |
---|
[644] | 250 | CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb", |
---|
| 251 | . "m", |
---|
[776] | 252 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 253 | . "inst(X)", zout,zout) |
---|
[524] | 254 | c |
---|
[644] | 255 | ENDIF !(bb2.EQ."500") |
---|
| 256 | c |
---|
| 257 | ENDDO !nlevSTD |
---|
| 258 | c |
---|
[897] | 259 | |
---|
[524] | 260 | c |
---|
[644] | 261 | cIM diagnostiques flux ocean-atm ou ocean-glace de mer |
---|
| 262 | cIM pour utilisation dans un modele de "slab" ocean |
---|
| 263 | c |
---|
| 264 | IF (OCEAN.EQ.'force ') THEN |
---|
| 265 | c |
---|
| 266 | CALL histdef(nid_day, "lmt_bils", |
---|
| 267 | . "Bilan au sol atmosphere forcee", "W/m2", |
---|
[776] | 268 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 269 | . "ave(X)", zstophy,zout) |
---|
| 270 | c |
---|
| 271 | ELSE IF (OCEAN.EQ.'slab ') THEN |
---|
| 272 | c |
---|
| 273 | CALL histdef(nid_day, "slab_bils", |
---|
| 274 | . "Bilan au sol slab", "W/m2", |
---|
[776] | 275 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 276 | . "ave(X)", zstophy,zout) |
---|
| 277 | c |
---|
| 278 | CALL histdef(nid_day, "tslab", "Slab SST ", "K", |
---|
[776] | 279 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 280 | . "ave(X)", zstophy,zout) |
---|
| 281 | c |
---|
| 282 | CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2", |
---|
[776] | 283 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 284 | . "ave(X)", zstophy,zout) |
---|
| 285 | c |
---|
| 286 | CALL histdef(nid_day, "siceh", "Slab seaice height", "m", |
---|
[776] | 287 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 288 | . "ave(X)", zstophy,zout) |
---|
| 289 | c |
---|
| 290 | ENDIF !OCEAN.EQ.force/slab |
---|
| 291 | c |
---|
| 292 | CALL histdef(nid_day, "fluxo", |
---|
| 293 | . "Flux turbulents ocean-atmosphere", "W/m2", |
---|
[776] | 294 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 295 | . "ave(X)", zstophy,zout) |
---|
| 296 | c |
---|
| 297 | CALL histdef(nid_day, "fluxg", |
---|
| 298 | . "Flux ocean-glace de mer", "W/m2", |
---|
[776] | 299 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 300 | . "ave(X)", zstophy,zout) |
---|
| 301 | c |
---|
[524] | 302 | ENDIF !lev_histday.GE.1 |
---|
| 303 | c |
---|
| 304 | IF (lev_histday.GE.2) THEN |
---|
| 305 | c |
---|
| 306 | CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2", |
---|
[776] | 307 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 308 | . "ave(X)", zstophy,zout) |
---|
[524] | 309 | c |
---|
| 310 | CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K", |
---|
[776] | 311 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 312 | . "ave(X)", zstophy,zout) |
---|
[524] | 313 | c |
---|
[897] | 314 | CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%", |
---|
| 315 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 316 | . "ave(X)", zstophy,zout) |
---|
[524] | 317 | c |
---|
| 318 | CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2", |
---|
[776] | 319 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 320 | . 32, "ave(X)", zstophy,zout) |
---|
[524] | 321 | c |
---|
| 322 | CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2", |
---|
[776] | 323 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 324 | . 32, "ave(X)", zstophy,zout) |
---|
[524] | 325 | c |
---|
| 326 | CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2", |
---|
[776] | 327 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 328 | . 32, "ave(X)", zstophy,zout) |
---|
[524] | 329 | c |
---|
| 330 | CALL histdef(nid_day, "SWupTOAclr", |
---|
| 331 | . "SWup clear sky at TOA","W/m2", |
---|
[776] | 332 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 333 | . 32, "ave(X)", zstophy,zout) |
---|
[897] | 334 | c |
---|
[524] | 335 | CALL histdef(nid_day, "SWupSFCclr", |
---|
| 336 | . "SWup clear sky at surface","W/m2", |
---|
[776] | 337 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 338 | . 32, "ave(X)", zstophy,zout) |
---|
[897] | 339 | c |
---|
[524] | 340 | CALL histdef(nid_day, "SWdnTOAclr", |
---|
| 341 | . "SWdn clear sky at TOA","W/m2", |
---|
[776] | 342 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 343 | . 32, "ave(X)", zstophy,zout) |
---|
[897] | 344 | c |
---|
[524] | 345 | CALL histdef(nid_day, "SWdnSFCclr", |
---|
| 346 | . "SWdn clear sky at surface","W/m2", |
---|
[776] | 347 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[644] | 348 | . 32, "ave(X)", zstophy,zout) |
---|
[524] | 349 | c |
---|
[684] | 350 | cIM AMMA-MIP |
---|
| 351 | CALL histdef(nid_day, "LWdnSFC", |
---|
| 352 | . "LW down at surface","W/m2", |
---|
[776] | 353 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[684] | 354 | . 32, "ave(X)", zstophy,zout) |
---|
| 355 | c |
---|
| 356 | CALL histdef(nid_day, "LWupSFC", |
---|
[897] | 357 | . "LW up at surface","W/m2", |
---|
[776] | 358 | . iim,jj_nb,nhori, 1,1,1,-99, |
---|
[684] | 359 | . 32, "ave(X)", zstophy,zout) |
---|
| 360 | c |
---|
[897] | 361 | CALL histdef(nid_day,"lwp","Cloud water path","kg/m2", |
---|
| 362 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 363 | . "ave(X)", zstophy,zout) |
---|
| 364 | c |
---|
| 365 | CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2", |
---|
| 366 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 367 | . "ave(X)", zstophy,zout) |
---|
| 368 | c |
---|
| 369 | CALL histdef(nid_day, "taux","Zonal wind stress","Pa", |
---|
| 370 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 371 | . "ave(X)", zstophy,zout) |
---|
| 372 | c |
---|
| 373 | CALL histdef(nid_day, "tauy","Meridional wind stress ", "Pa", |
---|
| 374 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
| 375 | . "ave(X)", zstophy,zout) |
---|
| 376 | c |
---|
| 377 | ENDIF !lev_histday.GE.2 |
---|
[524] | 378 | c================================================================= |
---|
[897] | 379 | IF(lev_histday.GE.3) THEN |
---|
[524] | 380 | c================================================================= |
---|
[897] | 381 | c INITIALISATION DES CHAMPS 3D |
---|
[524] | 382 | c================================================================= |
---|
[897] | 383 | c Champs 3D: |
---|
| 384 | c |
---|
| 385 | CALL histdef(nid_day, "temp", "Air temperature", "K", |
---|
| 386 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 387 | . "ave(X)", zstophy,zout) |
---|
| 388 | c |
---|
| 389 | CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg", |
---|
| 390 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 391 | . "ave(X)", zstophy,zout) |
---|
| 392 | c |
---|
| 393 | CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2", |
---|
| 394 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 395 | . "ave(X)", zstophy,zout) |
---|
| 396 | c |
---|
| 397 | CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", |
---|
| 398 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 399 | . "ave(X)", zstophy,zout) |
---|
| 400 | c |
---|
| 401 | CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", |
---|
| 402 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 403 | . "ave(X)", zstophy,zout) |
---|
| 404 | c |
---|
| 405 | CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", |
---|
| 406 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 407 | . "ave(X)", zstophy,zout) |
---|
| 408 | c |
---|
| 409 | CALL histdef(nid_day, "pres", "Air pressure", "Pa", |
---|
| 410 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
| 411 | . "ave(X)", zstophy,zout) |
---|
| 412 | c |
---|
| 413 | c================================================================= |
---|
| 414 | c FIN INITIALISATION DES CHAMPS 3D |
---|
| 415 | c================================================================= |
---|
| 416 | ENDIF ! lev_histday.GE.3 |
---|
| 417 | c |
---|
| 418 | IF (lev_histday.GE.4) THEN |
---|
| 419 | c================================================================= |
---|
[524] | 420 | c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES |
---|
| 421 | c================================================================= |
---|
| 422 | c |
---|
[644] | 423 | CALL histdef(nid_day, "tter", "Surface Temperature ter", "K", |
---|
[776] | 424 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 425 | . "ave(X)", zstophy,zout) |
---|
[524] | 426 | c |
---|
[644] | 427 | CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K", |
---|
[776] | 428 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 429 | . "ave(X)", zstophy,zout) |
---|
[524] | 430 | c |
---|
[644] | 431 | CALL histdef(nid_day, "toce", "Surface Temperature oce", "K", |
---|
[776] | 432 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 433 | . "ave(X)", zstophy,zout) |
---|
[524] | 434 | c |
---|
[644] | 435 | CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K", |
---|
[776] | 436 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 437 | . "ave(X)", zstophy,zout) |
---|
[524] | 438 | c |
---|
[644] | 439 | CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K", |
---|
[776] | 440 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 441 | . "ave(X)", zstophy,zout) |
---|
[524] | 442 | c |
---|
| 443 | CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K", |
---|
[776] | 444 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 445 | . "ave(X)", zstophy,zout) |
---|
[524] | 446 | c |
---|
| 447 | CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K", |
---|
[776] | 448 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 449 | . "ave(X)", zstophy,zout) |
---|
[524] | 450 | c |
---|
| 451 | CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K", |
---|
[776] | 452 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 453 | . "ave(X)", zstophy,zout) |
---|
[524] | 454 | c |
---|
| 455 | CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.", |
---|
| 456 | . "K", |
---|
[776] | 457 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 458 | . t2mincels, zstophy,zout) |
---|
[524] | 459 | c |
---|
| 460 | CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.", |
---|
| 461 | . "K", |
---|
[776] | 462 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 463 | . t2maxcels, zstophy,zout) |
---|
[524] | 464 | c |
---|
| 465 | CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s", |
---|
[776] | 466 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 467 | . "ave(X)", zstophy,zout) |
---|
[524] | 468 | c |
---|
| 469 | CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s", |
---|
[776] | 470 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 471 | . "ave(X)", zstophy,zout) |
---|
[524] | 472 | c |
---|
| 473 | CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s", |
---|
[776] | 474 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 475 | . "ave(X)", zstophy,zout) |
---|
[524] | 476 | c |
---|
| 477 | CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m", |
---|
[776] | 478 | . "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 479 | . "ave(X)", zstophy,zout) |
---|
[524] | 480 | c |
---|
| 481 | CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m", |
---|
[776] | 482 | . "m/s", iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 483 | . "ave(X)", zstophy,zout) |
---|
[524] | 484 | c |
---|
| 485 | CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m", |
---|
[776] | 486 | . "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 487 | . "ave(X)", zstophy,zout) |
---|
[524] | 488 | c |
---|
| 489 | CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m", |
---|
[776] | 490 | . "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 491 | . "ave(X)", zstophy,zout) |
---|
[524] | 492 | c |
---|
| 493 | CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m", |
---|
[776] | 494 | . "m/s",iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 495 | . "ave(X)", zstophy,zout) |
---|
[524] | 496 | c |
---|
| 497 | DO nsrf = 1, nbsrf |
---|
| 498 | C |
---|
| 499 | call histdef(nid_day, "pourc_"//clnsurf(nsrf), |
---|
[644] | 500 | $ "% "//clnsurf(nsrf), "%", |
---|
[776] | 501 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 502 | $ "ave(X)", zstophy,zout) |
---|
[524] | 503 | C |
---|
[644] | 504 | call histdef(nid_day, "fract_"//clnsurf(nsrf), |
---|
| 505 | $ "Fraction "//clnsurf(nsrf), "1", |
---|
[776] | 506 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 507 | $ "ave(X)", zstophy,zout) |
---|
| 508 | C |
---|
[524] | 509 | call histdef(nid_day, "tsol_"//clnsurf(nsrf), |
---|
[644] | 510 | $ "Temperature "//clnsurf(nsrf), "K", |
---|
[776] | 511 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 512 | $ "ave(X)", zstophy,zout) |
---|
[524] | 513 | C |
---|
| 514 | call histdef(nid_day, "sens_"//clnsurf(nsrf), |
---|
| 515 | $ "Sensible heat flux "//clnsurf(nsrf), "W/m2", |
---|
[776] | 516 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 517 | $ "ave(X)", zstophy,zout) |
---|
[524] | 518 | c |
---|
| 519 | call histdef(nid_day, "lat_"//clnsurf(nsrf), |
---|
| 520 | $ "Latent heat flux "//clnsurf(nsrf), "W/m2", |
---|
[776] | 521 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 522 | $ "ave(X)", zstophy,zout) |
---|
[524] | 523 | C |
---|
| 524 | call histdef(nid_day, "taux_"//clnsurf(nsrf), |
---|
[644] | 525 | $ "Zonal wind stress "//clnsurf(nsrf),"Pa", |
---|
[776] | 526 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 527 | $ "ave(X)", zstophy,zout) |
---|
[524] | 528 | |
---|
| 529 | call histdef(nid_day, "tauy_"//clnsurf(nsrf), |
---|
[644] | 530 | $ "Meridional wind stress "//clnsurf(nsrf), "Pa", |
---|
[776] | 531 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 532 | $ "ave(X)", zstophy,zout) |
---|
[524] | 533 | C |
---|
| 534 | call histdef(nid_day, "albe_"//clnsurf(nsrf), |
---|
[644] | 535 | $ "Albedo surf. "//clnsurf(nsrf), "-", |
---|
[776] | 536 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 537 | $ "ave(X)", zstophy,zout) |
---|
[524] | 538 | C |
---|
| 539 | call histdef(nid_day, "rugs_"//clnsurf(nsrf), |
---|
[644] | 540 | $ "Rugosite "//clnsurf(nsrf), "m", |
---|
[776] | 541 | $ iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 542 | $ "ave(X)", zstophy,zout) |
---|
[524] | 543 | C |
---|
| 544 | END DO |
---|
| 545 | C |
---|
| 546 | c================================================================= |
---|
| 547 | c |
---|
| 548 | c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES |
---|
| 549 | c |
---|
| 550 | c================================================================= |
---|
| 551 | ENDIF !lev_histday.GE.4 |
---|
[644] | 552 | c |
---|
| 553 | IF(lev_histday.GE.5) THEN !lev_histday.GE.5 |
---|
| 554 | c |
---|
| 555 | c rajout sorties F. Aires |
---|
| 556 | CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA", |
---|
| 557 | . "W/m2", |
---|
[776] | 558 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 559 | . "ave(X)", zstophy,zout) |
---|
| 560 | c |
---|
| 561 | CALL histdef(nid_day, "topl0", "CS IR rad. at TOA", |
---|
| 562 | . "W/m2", |
---|
[776] | 563 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 564 | . "ave(X)", zstophy,zout) |
---|
| 565 | c |
---|
| 566 | CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.", |
---|
| 567 | . "W/m2", |
---|
[776] | 568 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 569 | . "ave(X)", zstophy,zout) |
---|
| 570 | c |
---|
| 571 | CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface", |
---|
| 572 | . "W/m2", |
---|
[776] | 573 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[644] | 574 | . "ave(X)", zstophy,zout) |
---|
| 575 | c |
---|
[918] | 576 | CALL histdef(nid_day, "rneb", "Cloud fraction", "1", |
---|
[776] | 577 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 578 | . "ave(X)", zstophy,zout) |
---|
| 579 | c |
---|
| 580 | CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction" |
---|
[918] | 581 | . , "1", |
---|
[776] | 582 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
[644] | 583 | . "ave(X)", zstophy,zout) |
---|
| 584 | c |
---|
| 585 | CALL histdef(nid_day,"lwcon","Cloud liquid water content", |
---|
| 586 | . "kg/kg", |
---|
[776] | 587 | . iim,jj_nb,nhori, klev,1,klev, nvert, 32, |
---|
[644] | 588 | . "ave(X)", zstophy,zout) |
---|
| 589 | c |
---|
| 590 | CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg", |
---|
[776] | 591 | . iim,jj_nb,nhori, klev,1,klev, nvert, 32, |
---|
[644] | 592 | . "ave(X)", zstophy,zout) |
---|
| 593 | c |
---|
| 594 | CALL histdef(nid_day,"cldtau","Cloud optical thickness","1", |
---|
[776] | 595 | . iim,jj_nb,nhori, klev,1,klev, nvert, 32, |
---|
[644] | 596 | . "ave(X)", zstophy,zout) |
---|
| 597 | c |
---|
| 598 | CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1", |
---|
[776] | 599 | . iim,jj_nb,nhori, klev,1,klev, nvert, 32, |
---|
[644] | 600 | . "ave(X)", zstophy,zout) |
---|
| 601 | c |
---|
[684] | 602 | cIM rajout AMMA-MIP |
---|
| 603 | c |
---|
| 604 | CALL histdef(nid_day, "dtcon", "Convection dT", "K/s", |
---|
[776] | 605 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
[684] | 606 | . "ave(X)", zstophy,zout) |
---|
| 607 | c |
---|
| 608 | CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s", |
---|
[776] | 609 | . iim,jj_nb,nhori, klev,1,klev,nvert, 32, |
---|
[684] | 610 | . "ave(X)", zstophy,zout) |
---|
| 611 | c |
---|
| 612 | CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%", |
---|
[776] | 613 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[684] | 614 | . "ave(X)", zstophy,zout) |
---|
| 615 | c |
---|
| 616 | CALL histdef(nid_day, "tpot", |
---|
| 617 | . "Surface air potential temperature", "K", |
---|
[776] | 618 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[684] | 619 | . "ave(X)", zstophy,zout) |
---|
| 620 | c |
---|
| 621 | CALL histdef(nid_day, "tpote", |
---|
| 622 | . "Surface air equivalent potential temperature", "K", |
---|
[776] | 623 | . iim,jj_nb,nhori, 1,1,1, -99, 32, |
---|
[684] | 624 | . "ave(X)", zstophy,zout) |
---|
| 625 | c |
---|
[644] | 626 | ENDIF !lev_histday.GE.5 |
---|
[524] | 627 | c================================================================= |
---|
| 628 | c |
---|
| 629 | CALL histend(nid_day) |
---|
| 630 | c |
---|
| 631 | ndex2d = 0 |
---|
| 632 | ndex3d = 0 |
---|
| 633 | c |
---|
| 634 | c================================================================= |
---|
[766] | 635 | c$OMP END MASTER |
---|
[524] | 636 | ENDIF ! fin de test sur ok_journe |
---|