Changeset 506 for LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h
- Timestamp:
- Mar 24, 2004, 9:43:54 AM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h
r486 r506 1 1 IF (ok_journe) THEN 2 c 3 IF(region.EQ.'3d') THEN 4 zsto = dtime 5 c zout = dtime * FLOAT(ecrit_day) 6 c zsto1= dtime * FLOAT(ecrit_day) 7 zout = dtime * REAL(ecrit_day) 8 zsto1= dtime * REAL(ecrit_day) 9 ELSEIF(region.EQ.'1d') THEN 10 zsto = dtime 11 c zout = dtime 12 zout = dtime *4. 13 ENDIF 2 14 c 3 15 idayref = day_ref … … 13 25 ENDDO 14 26 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 27 write(*,*)'zx_lon = ',zx_lon(:,1) 28 write(*,*)'zx_lat = ',zx_lat(1,:) 15 29 CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 16 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 30 c071103 31 c . 1,iim,1,jjmp1, 32 . imin_debut,nbpti, 33 . jmin_debut,nbptj, 34 c071103 35 . itau_phy, zjulian, dtime, 17 36 . nhori, nid_day) 18 37 write(*,*)'Journee ', itau_phy, zjulian 19 38 CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", 20 39 . klev, presnivs, nvert) 21 c call histvert(nid_day, 'sig_s', 'Niveaux sigma','-', 22 c . klev, znivsig, nvert) 23 c 24 zsto = dtime 25 zout = dtime * FLOAT(ecrit_day) 26 C Essai writephys 27 c nom_fichier = 'histday1' 28 c call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev, 29 c . rlon,rlat, presnivs, 30 c . zjulian, dtime) 31 c call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0) 32 c call writephy_def(prof3d_on, fid_day, "once", zsto, zout, 33 c . klev) 34 c call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0) 35 c call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout, 36 c . klev) 37 40 c 41 IF(lev_histday.GE.1) THEN 38 42 c 39 43 CALL histdef(nid_day, "phis", "Surface geop. height", "-", … … 45 49 . "once", zsto,zout) 46 50 c 51 CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-", 52 . iim,jjmp1,nhori, 1,1,1, -99, 32, 53 . "once", zsto,zout) 54 c 55 CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-", 56 . iim,jjmp1,nhori, 1,1,1, -99, 32, 57 . "once", zsto,zout) 58 c 47 59 c Champs 2D: 48 60 c … … 51 63 . "ave(X)", zsto,zout) 52 64 c 53 c54 65 CALL histdef(nid_day, "t2m", "Temperature 2m", "K", 55 66 . iim,jjmp1,nhori, 1,1,1, -99, 32, 56 67 . "ave(X)", zsto,zout) 57 68 c 58 CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",59 . iim,jjmp1,nhori, 1,1,1, -99, 32,60 . "ave(X)", zsto,zout)61 c62 CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",63 . iim,jjmp1,nhori, 1,1,1, -99, 32,64 . "ave(X)", zsto,zout)65 c66 CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",67 . iim,jjmp1,nhori, 1,1,1, -99, 32,68 . "ave(X)", zsto,zout)69 c70 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",71 . iim,jjmp1,nhori, 1,1,1, -99, 32,72 . "ave(X)", zsto,zout)73 c74 CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"75 . , "kg/(s*m2)",76 . iim,jjmp1,nhori, 1,1,1, -99, 32,77 . "ave(X)", zsto,zout)78 c79 CALL histdef(nid_day, "snow", "Snow fall", "kg/(s*m2)",80 . iim,jjmp1,nhori, 1,1,1, -99, 32,81 . "ave(X)", zsto,zout)82 c83 CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",84 . iim,jjmp1,nhori, 1,1,1, -99, 32,85 . "ave(X)", zsto,zout)86 c87 CALL histdef(nid_day, "evap", "Evaporation", "kg/(s*m2)",88 . iim,jjmp1,nhori, 1,1,1, -99, 32,89 . "ave(X)", zsto,zout)90 c91 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",92 . iim,jjmp1,nhori, 1,1,1, -99, 32,93 . "ave(X)", zsto,zout)94 c95 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",96 . iim,jjmp1,nhori, 1,1,1, -99, 32,97 . "ave(X)", zsto,zout)98 c99 CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",100 . "W/m2",101 . iim,jjmp1,nhori, 1,1,1, -99, 32,102 . "ave(X)", zsto,zout)103 c104 CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",105 . iim,jjmp1,nhori, 1,1,1, -99, 32,106 . "ave(X)", zsto,zout)107 c108 CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",109 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,110 . "ave(X)", zsto,zout)111 c112 CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",113 . iim,jjmp1,nhori, 1,1,1, -99, 32,114 . "ave(X)", zsto,zout)115 c116 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",117 . iim,jjmp1,nhori, 1,1,1, -99, 32,118 . "ave(X)", zsto,zout)119 c120 CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",121 . iim,jjmp1,nhori, 1,1,1, -99, 32,122 . "ave(X)", zsto,zout)123 c124 c CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",125 c . iim,jjmp1,nhori, 1,1,1, -99, 32,126 c . "ave(X)", zsto,zout)127 c128 c CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",129 c . iim,jjmp1,nhori, 1,1,1, -99, 32,130 c . "ave(X)", zsto,zout)131 c132 c CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",133 c . iim,jjmp1,nhori, 1,1,1, -99, 32,134 c . "ave(X)", zsto,zout)135 c136 CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",137 . iim,jjmp1,nhori, 1,1,1, -99, 32,138 . "ave(X)", zsto,zout)139 c140 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",141 . iim,jjmp1,nhori, 1,1,1, -99, 32,142 . "ave(X)", zsto,zout)143 c144 CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",145 . iim,jjmp1,nhori, 1,1,1, -99, 32,146 . "ave(X)", zsto,zout)147 c148 CALL histdef(nid_day, "cldt", "Total cloudiness", "-",149 . iim,jjmp1,nhori, 1,1,1, -99, 32,150 . "ave(X)", zsto,zout)151 c152 CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",153 . iim,jjmp1,nhori, 1,1,1, -99, 32,154 . "ave(X)", zsto,zout)155 c156 CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",157 . iim,jjmp1,nhori, 1,1,1,-99,158 . 32, "ave(X)", zsto,zout)159 c160 CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",161 . iim,jjmp1,nhori, 1,1,1,-99,162 . 32, "ave(X)", zsto,zout)163 c164 CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",165 . iim,jjmp1,nhori, 1,1,1,-99,166 . 32, "ave(X)", zsto,zout)167 c168 CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",169 . iim,jjmp1,nhori, 1,1,1,-99,170 . 32, "ave(X)", zsto,zout)171 c172 CALL histdef(nid_day, "SWupTOAclr",173 . "SWup clear sky at TOA","W/m2",174 . iim,jjmp1,nhori, 1,1,1,-99,175 . 32, "ave(X)", zsto,zout)176 177 CALL histdef(nid_day, "SWupSFCclr",178 . "SWup clear sky at surface","W/m2",179 . iim,jjmp1,nhori, 1,1,1,-99,180 . 32, "ave(X)", zsto,zout)181 182 CALL histdef(nid_day, "SWdnTOAclr",183 . "SWdn clear sky at TOA","W/m2",184 . iim,jjmp1,nhori, 1,1,1,-99,185 . 32, "ave(X)", zsto,zout)186 187 CALL histdef(nid_day, "SWdnSFCclr",188 . "SWdn clear sky at surface","W/m2",189 . iim,jjmp1,nhori, 1,1,1,-99,190 . 32, "ave(X)", zsto,zout)191 cccIM192 CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",193 . iim,jjmp1,nhori, 1,1,1, -99, 32,194 . "ave(X)", zsto,zout)195 c196 c Champs dynamiques sur niveaux de pression197 198 CALL histdef(nid_day, "u850", "Zonal wind 850mb", "m/s",199 . iim,jjmp1,nhori, 1,1,1, -99, 32,200 . "ave(X)", zsto,zout)201 202 CALL histdef(nid_day, "v850", "Meridional wind 850mb", "m/s",203 . iim,jjmp1,nhori, 1,1,1, -99, 32,204 . "ave(X)", zsto,zout)205 c206 CALL histdef(nid_day, "u500", "Zonal wind 500mb", "m/s",207 . iim,jjmp1,nhori, 1,1,1, -99, 32,208 . "ave(X)", zsto,zout)209 210 CALL histdef(nid_day, "v500", "Meridional wind 500mb", "m/s",211 . iim,jjmp1,nhori, 1,1,1, -99, 32,212 . "ave(X)", zsto,zout)213 214 CALL histdef(nid_day, "u200", "Zonal wind 200mb", "m/s",215 . iim,jjmp1,nhori, 1,1,1, -99, 32,216 . "ave(X)", zsto,zout)217 218 CALL histdef(nid_day, "v200", "Meridional wind 200mb", "m/s",219 . iim,jjmp1,nhori, 1,1,1, -99, 32,220 . "ave(X)", zsto,zout)221 222 CALL histdef(nid_day,"phi500", "Geopotentiel à 500mb", "m2/s2",223 . iim,jjmp1,nhori, 1,1,1, -99, 32,224 . "ave(X)", zsto,zout)225 226 CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",227 . iim,jjmp1,nhori, 1,1,1, -99, 32,228 . "ave(X)", zsto,zout)229 230 CALL histdef(nid_day, "cape_max", "CAPE max.",231 . "J/Kg",232 . iim,jjmp1,nhori, 1,1,1, -99, 32,233 . capemaxcels, zsto,zout)234 235 236 c=================================================================237 c CI DESSOUS, SORTIES A UTILISER EN DEBUG238 c=================================================================239 if ( lev_histday.gt.1 ) then240 241 69 c Champs retires momentannéement en attendant un hypothetique 242 70 c debugage … … 252 80 . t2maxcels, zsto,zout) 253 81 c 254 c================================================================= 255 c================================================================= 256 c================================================================= 82 CALL histdef(nid_day, "plul", "Large-scale Precip.", 83 . "kg/(s*m2)", 84 . iim,jjmp1,nhori, 1,1,1, -99, 32, 85 . "ave(X)", zsto,zout) 86 c 87 CALL histdef(nid_day, "pluc", "Convective Precip.", 88 . "kg/(s*m2)", 89 . iim,jjmp1,nhori, 1,1,1, -99, 32, 90 . "ave(X)", zsto,zout) 91 c 92 CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.", 93 . "kg/(s*m2)", 94 . iim,jjmp1,nhori, 1,1,1, -99, 32, 95 . "ave(X)", zsto,zout) 96 c 97 IF(1.EQ.0) THEN !snowc=0. 98 CALL histdef(nid_day, "snowc", "Solid Convective Precip.", 99 . "kg/(s*m2)", 100 . iim,jjmp1,nhori, 1,1,1, -99, 32, 101 . "ave(X)", zsto,zout) 102 ENDIF !snowc=0. 103 c 104 CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2", 105 . iim,jjmp1,nhori, 1,1,1, -99, 32, 106 . "ave(X)", zsto,zout) 107 c 108 CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-", 109 . iim,jjmp1,nhori, 1,1,1, -99, 32, 110 . "inst(X)", zsto1,zout) 111 c 112 CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg", 113 . iim,jjmp1,nhori, 1,1,1, -99, 32, 114 . "ave(X)", zsto,zout) 115 c 116 CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s", 117 . iim,jjmp1,nhori, 1,1,1, -99, 32, 118 . "ave(X)", zsto,zout) 119 c 120 CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s", 121 . iim,jjmp1,nhori, 1,1,1, -99, 32, 122 . "ave(X)", zsto,zout) 123 c 124 CALL histdef(nid_day, "wind10m","10-m wind speed","m/s", 125 . iim,jjmp1,nhori, 1,1,1, -99, 32, 126 . "ave(X)", zsto,zout) 127 c 128 CALL histdef(nid_day, "wind10max", "10-m wind speed max.", 129 . "m/s", 130 . iim,jjmp1,nhori, 1,1,1, -99, 32, 131 . t2maxcels, zsto,zout) 132 c 133 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", 134 . iim,jjmp1,nhori, 1,1,1, -99, 32, 135 . "ave(X)", zsto,zout) 136 c 137 CALL histdef(nid_day, "precip","Precipitation Totale liq+sol" 138 . , "kg/(s*m2)", 139 . iim,jjmp1,nhori, 1,1,1, -99, 32, 140 . "ave(X)", zsto,zout) 141 c 142 CALL histdef(nid_day, "snowf", "Snow fall", "kg/(s*m2)", 143 . iim,jjmp1,nhori, 1,1,1, -99, 32, 144 . "ave(X)", zsto,zout) 145 c 146 CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2", 147 . iim,jjmp1,nhori, 1,1,1, -99, 32, 148 . "ave(X)", zsto,zout) 149 c 150 CALL histdef(nid_day, "evap", "Evaporation", "kg/(s*m2)", 151 . iim,jjmp1,nhori, 1,1,1, -99, 32, 152 . "ave(X)", zsto,zout) 153 c 154 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 155 . iim,jjmp1,nhori, 1,1,1, -99, 32, 156 . "ave(X)", zsto,zout) 157 c 158 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 159 . iim,jjmp1,nhori, 1,1,1, -99, 32, 160 . "ave(X)", zsto,zout) 161 c 162 CALL histdef(nid_day, "sols", "Net Solar rad. at surf.", 163 . "W/m2", 164 . iim,jjmp1,nhori, 1,1,1, -99, 32, 165 . "ave(X)", zsto,zout) 166 c 167 CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2", 168 . iim,jjmp1,nhori, 1,1,1, -99, 32, 169 . "ave(X)", zsto,zout) 170 c 171 CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-", 172 . iim,jjmp1,nhori, 1,1,1, -99, 32, 173 . "ave(X)", zsto,zout) 174 c 175 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-", 176 . iim,jjmp1,nhori, 1,1,1, -99, 32, 177 . "ave(X)", zsto,zout) 178 c 179 CALL histdef(nid_day, "cldh", "High-level cloudiness", "-", 180 . iim,jjmp1,nhori, 1,1,1, -99, 32, 181 . "ave(X)", zsto,zout) 182 c 183 CALL histdef(nid_day, "cldt", "Total cloudiness", "-", 184 . iim,jjmp1,nhori, 1,1,1, -99, 32, 185 . "ave(X)", zsto,zout) 186 c 187 CALL histdef(nid_day, "cldq", "Cloud liquid water path", 188 . "kg/m2", 189 . iim,jjmp1,nhori, 1,1,1, -99, 32, 190 . "ave(X)", zsto,zout) 191 c 192 CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2", 193 . iim,jjmp1,nhori, 1,1,1, -99, 32, 194 . "ave(X)", zsto,zout) 195 c 196 c Champs dynamiques sur niveaux de pression 197 DO k=1, nlevENS 198 c 199 CALL histdef(nid_day, "u"//clev(k), 200 . "Zonal wind"//clev(k)//"mb","m/s", 201 . iim,jjmp1,nhori, 1,1,1, -99, 32, 202 . "ave(X)", zsto,zout) 203 c 204 CALL histdef(nid_day, "v"//clev(k), 205 . "Meridional wind"//clev(k)//"mb","m/s", 206 . iim,jjmp1,nhori, 1,1,1, -99, 32, 207 . "ave(X)", zsto,zout) 208 c 209 ENDDO !nlevENS 210 c 211 CALL histdef(nid_day, "w500", "Verical wind 500mb", "m/s", 212 . iim,jjmp1,nhori, 1,1,1, -99, 32, 213 . "ave(X)", zsto,zout) 214 c 215 CALL histdef(nid_day,"phi500", "Geopotentiel a 500mb","m2/s2", 216 . iim,jjmp1,nhori, 1,1,1, -99, 32, 217 . "ave(X)", zsto,zout) 218 c 219 CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa", 220 . iim,jjmp1,nhori, 1,1,1, -99, 32, 221 . "ave(X)", zsto,zout) 222 c 223 CALL histdef(nid_day, "cape_max", "CAPE max.", 224 . "J/kg", 225 . iim,jjmp1,nhori, 1,1,1, -99, 32, 226 . capemaxcels, zsto,zout) 227 c 228 CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", 229 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 230 . "ave(X)", zsto,zout) 231 c 232 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", 233 . iim,jjmp1,nhori, 1,1,1, -99, 32, 234 . "ave(X)", zsto,zout) 235 c 236 CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2", 237 . iim,jjmp1,nhori, 1,1,1,-99, 238 . 32, "ave(X)", zsto,zout) 239 c 240 ENDIF !lev_histday.GE.1 241 c 242 IF (lev_histday.GE.2) THEN 243 c 244 CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2", 245 . iim,jjmp1,nhori, 1,1,1, -99, 32, 246 . "ave(X)", zsto,zout) 247 c 248 CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K", 249 . iim,jjmp1,nhori, 1,1,1, -99, 32, 250 . "ave(X)", zsto,zout) 251 c 252 ENDIF !lev_histday.GE.2 253 c================================================================= 254 IF(lev_histday.GE.3) THEN 255 c================================================================= 256 c INITIALISATION DES CHAMPS 3D 257 c================================================================= 258 c================================================================= 259 c================================================================= 260 c Champs 3D: 261 c 262 CALL histdef(nid_day, "temp", "Air temperature", "K", 263 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 264 . "ave(X)", zsto,zout) 265 c 266 CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg", 267 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 268 . "ave(X)", zsto,zout) 269 c 270 CALL histdef(nid_day, "geop", "Geopotential height", "m", 271 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 272 . "ave(X)", zsto,zout) 273 c 274 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 275 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 276 . "ave(X)", zsto,zout) 277 c 278 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 279 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 280 . "ave(X)", zsto,zout) 281 c 282 CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", 283 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 284 . "ave(X)", zsto,zout) 285 c 286 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 287 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 288 . "ave(X)", zsto,zout) 289 c 290 c================================================================= 291 c FIN INITIALISATION DES CHAMPS 3D 292 c================================================================= 293 ENDIF ! lev_histday.GE.3 294 c 295 IF (lev_histday.GE.4) THEN 296 c================================================================= 297 c 257 298 c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES 258 c================================================================= 299 c 300 c================================================================= 301 c 302 CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2", 303 . iim,jjmp1,nhori, 1,1,1,-99, 304 . 32, "ave(X)", zsto,zout) 305 c 306 CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2", 307 . iim,jjmp1,nhori, 1,1,1,-99, 308 . 32, "ave(X)", zsto,zout) 309 c 310 CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2", 311 . iim,jjmp1,nhori, 1,1,1,-99, 312 . 32, "ave(X)", zsto,zout) 313 c 314 CALL histdef(nid_day, "SWupTOAclr", 315 . "SWup clear sky at TOA","W/m2", 316 . iim,jjmp1,nhori, 1,1,1,-99, 317 . 32, "ave(X)", zsto,zout) 259 318 319 CALL histdef(nid_day, "SWupSFCclr", 320 . "SWup clear sky at surface","W/m2", 321 . iim,jjmp1,nhori, 1,1,1,-99, 322 . 32, "ave(X)", zsto,zout) 323 324 CALL histdef(nid_day, "SWdnTOAclr", 325 . "SWdn clear sky at TOA","W/m2", 326 . iim,jjmp1,nhori, 1,1,1,-99, 327 . 32, "ave(X)", zsto,zout) 328 329 CALL histdef(nid_day, "SWdnSFCclr", 330 . "SWdn clear sky at surface","W/m2", 331 . iim,jjmp1,nhori, 1,1,1,-99, 332 . 32, "ave(X)", zsto,zout) 333 c 334 c================================================================= 335 c================================================================= 336 c================================================================= 337 c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES 338 c================================================================= 339 c 260 340 CALL histdef(nid_day, "tter", "Surface Temperature", "K", 261 341 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 273 353 . iim,jjmp1,nhori, 1,1,1, -99, 32, 274 354 . "ave(X)", zsto,zout) 275 c276 cccIM277 355 c 278 356 CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K", … … 334 412 . "ave(X)", zsto,zout) 335 413 c 336 c337 CXXX PB flux pour chaque sous surface338 C339 414 DO nsrf = 1, nbsrf 340 415 C … … 378 453 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 379 454 $ "ave(X)", zsto,zout) 380 381 CXXX 455 C 382 456 END DO 383 384 c================================================================= 385 c FIN DES CHAMPS SUR LES SOUS SURFACES 386 c================================================================= 387 c================================================================= 388 c FIN INITIALISATION DES CHAMPS 3D 389 c================================================================= 390 391 c Champs 3D: 392 c 393 CALL histdef(nid_day, "temp", "Air temperature", "K", 394 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 395 . "ave(X)", zsto,zout) 396 c 397 CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg", 398 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 399 . "ave(X)", zsto,zout) 400 c 401 CALL histdef(nid_day, "geop", "Geopotential height", "m", 402 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 403 . "ave(X)", zsto,zout) 404 c 405 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 406 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 407 . "ave(X)", zsto,zout) 408 c 409 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 410 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 411 . "ave(X)", zsto,zout) 412 c 413 CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", 414 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 415 . "ave(X)", zsto,zout) 416 c 417 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 418 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 419 . "ave(X)", zsto,zout) 420 c 421 c 422 c================================================================= 423 c FIN INITIALISATION DES CHAMPS 3D 424 c================================================================= 425 426 427 c================================================================= 428 c================================================================= 429 endif ! lev_histday.gt.1 430 c================================================================= 431 457 C 458 c================================================================= 459 c 460 c FIN INITIALISATION DES CHAMPS SUR LES SOUS SURFACES 461 c 462 c================================================================= 463 ENDIF !lev_histday.GE.4 464 c================================================================= 465 c 432 466 CALL histend(nid_day) 433 467 c
Note: See TracChangeset
for help on using the changeset viewer.