Here we attempt to provide useful code to generate figures from WRF outputs based on known galleries. For instance, NCL and WRF-Python provides extensive examples for plotting WRF outputs. Therefore, we aim to replicate some of these. Our approach to read wrfout files is based on eixport which relies r packages with GDAL bindings such as raster and stars. We do not try to provide a full gallery, instead, some basics and necessary plots to inspire other R used and receive more examples so share with the community.
library(eixport)
library(raster)
#> Loading required package: sp
library(stars)
#> Loading required package: abind
#> Loading required package: sf
#> Linking to GEOS 3.9.0, GDAL 3.2.2, PROJ 7.2.1
library(cptcity)
library(sf)Based on NCL:
First lets get a summary of a WRF output file
wrfo <- "/media/sergio/ext41/wrfo_sebrazil/wrfout_d02_2014-11-03_000000"
(dt <- wrf_meta(wrfo)$vars)
#> Warning in rbind(c(FieldType = "104", MemoryOrder = "XYZ", description =
#> "Aerosol Optical Depth", : number of columns of result is not a multiple of
#> vector length (arg 158)
#> vars FieldType MemoryOrder
#> 110 ac0 104 XYZ
#> 148 ac0cw 104 XYZ
#> 227 ACGRDFLX 104 XY
#> 299 ACHFX 104 XY
#> 300 ACLHF 104 XY
#> 60 aco3 104 XYZ
#> 228 ACSNOM 104 XY
#> 289 ALBBCK 104 XY
#> 287 ALBEDO 104 XY
#> 39 ald 104 XYZ
#> 106 antha 104 XYZ
#> 144 anthcw 104 XYZ
#> 1 AOD_OUT 104 XYZ
#> 2 AOD2D_OUT 104 XY
#> 24 ASOA 104 XYZ
#> 3 ATOP2D_OUT 104 XY
#> 23 BSOA 104 XYZ
#> 312 C1F 104 Z
#> 310 C1H 104 Z
#> 313 C2F 104 Z
#> 311 C2H 104 Z
#> 316 C3F 104 Z
#> 314 C3H 104 Z
#> 317 C4F 104 Z
#> 315 C4H 104 Z
#> 231 CANWAT 104 XY
#> 201 CF1 104 0
#> 202 CF2 104 0
#> 203 CF3 104 0
#> 187 CFN 104 0
#> 188 CFN1 104 0
#> 83 clai 104 XYZ
#> 82 claj 104 XYZ
#> 288 CLAT 104 XY
#> 121 clcwi 104 XYZ
#> 120 clcwj 104 XYZ
#> 278 CLDFRA 104 XYZ
#> 6 CLDFRA2 104 XYZ
#> 152 CN2O5 104 XYZ
#> 54 co 104 XYZ
#> 111 corn 104 XYZ
#> 149 corncw 104 XYZ
#> 260 COSALPHA 104 XY
#> 233 COSZEN 104 XY
#> 69 csl 104 XYZ
#> 67 dcb 104 XYZ
#> 25 DMS_0 104 XY
#> 186 DN 104 Z
#> 185 DNW 104 Z
#> 14 DRY_DEP_FLUX 104 XYZ
#> 13 DRY_DEP_LEN 104 XYZ
#> 12 DRYDEPVEL 104 XY
#> 19 dvel_o3 104 XYZ
#> 161 DZS 104 Z
#> 258 E 104 XY
#> 17 EBIO_API 104 XY
#> 16 EBIO_ISO 104 XY
#> 141 eccwi 104 XYZ
#> 140 eccwj 104 XYZ
#> 103 eci 104 XYZ
#> 102 ecj 104 XYZ
#> 236 EDMF_A 104 XYZ
#> 240 EDMF_ENT 104 XYZ
#> 241 EDMF_QC 104 XYZ
#> 239 EDMF_QT 104 XYZ
#> 238 EDMF_THL 104 XYZ
#> 237 EDMF_W 104 XYZ
#> 246 EL_PBL 104 XYZ
#> 290 EMISS 104 XY
#> 53 eth 104 XYZ
#> 8 EVAPPROD 104 XYZ
#> 257 F 104 XY
#> 181 FNM 104 Z
#> 182 FNP 104 Z
#> 151 GAMN2O5 104 XYZ
#> 280 GLW 104 XY
#> 65 gly 104 XYZ
#> 276 GRAUPELNC 104 XY
#> 226 GRDFLX 104 XY
#> 38 h2o2 104 XYZ
#> 277 HAILNC 104 XY
#> 50 hc3 104 XYZ
#> 51 hc5 104 XYZ
#> 52 hc8 104 XYZ
#> 40 hcho 104 XYZ
#> 71 hcl 104 XYZ
#> 296 HFX 104 XY
#> 170 HFX_FORCE 104 0
#> 173 HFX_FORCE_TEND 104 0
#> 261 HGT 104 XY
#> 37 hno3 104 XYZ
#> 63 hno4 104 XYZ
#> 72 ho 104 XYZ
#> 73 ho2 104 XYZ
#> 62 hono 104 XYZ
#> 4 ICN_DIAG 104 XYZ
#> 306 ISEEDARR_RAND_PERTURB 106 Z
#> 305 ISEEDARR_SKEBS 106 Z
#> 304 ISEEDARR_SPPT 106 Z
#> 307 ISEEDARRAY_SPP_CONV 106 Z
#> 309 ISEEDARRAY_SPP_LSM 106 Z
#> 308 ISEEDARRAY_SPP_PBL 106 Z
#> 224 ISLTYP 106 XY
#> 70 iso 104 XYZ
#> 204 ITIMESTEP 106 0
#> 223 IVGTYP 106 XY
#> 64 ket 104 XYZ
#> 153 KN2O5 104 XYZ
#> 234 LAI 104 XY
#> 18 LAI_VEGMASK 104 XY
#> 321 LAKEMASK 104 XY
#> 320 LANDMASK 104 XY
#> 298 LH 104 XY
#> 171 LH_FORCE 104 0
#> 174 LH_FORCE_TEND 104 0
#> 157 LU_INDEX 104 XY
#> 247 MAPFAC_M 104 XY
#> 250 MAPFAC_MX 104 XY
#> 251 MAPFAC_MY 104 XY
#> 248 MAPFAC_U 104 XY
#> 252 MAPFAC_UX 104 XY
#> 253 MAPFAC_UY 104 XY
#> 249 MAPFAC_V 104 XY
#> 254 MAPFAC_VX 104 XY
#> 256 MAPFAC_VY 104 XY
#> 270 MAX_MSTFX 104 0
#> 271 MAX_MSTFY 104 0
#> 243 MAXMF 104 XY
#> 255 MF_VX_INV 104 XY
#> 66 mgly 104 XYZ
#> 176 MU 104 XY
#> 177 MUB 104 XY
#> 47 n2o5 104 XYZ
#> 81 naai 104 XYZ
#> 80 naaj 104 XYZ
#> 119 nacwi 104 XYZ
#> 118 nacwj 104 XYZ
#> 5 NC_DIAG 104 XYZ
#> 178 NEST_POS 104 XY
#> 46 nh3 104 XYZ
#> 77 nh4ai 104 XYZ
#> 76 nh4aj 104 XYZ
#> 115 nh4cwi 104 XYZ
#> 114 nh4cwj 104 XYZ
#> 35 no 104 XYZ
#> 34 no2 104 XYZ
#> 48 no3 104 XYZ
#> 79 no3ai 104 XYZ
#> 78 no3aj 104 XYZ
#> 117 no3cwi 104 XYZ
#> 116 no3cwj 104 XYZ
#> 291 NOAHRES 104 XY
#> 109 nu0 104 XYZ
#> 147 nu0cw 104 XYZ
#> 242 NUPDRAFT 106 XY
#> 36 o3 104 XYZ
#> 55 ol2 104 XYZ
#> 57 oli 104 XYZ
#> 282 OLR 104 XY
#> 56 olt 104 XYZ
#> 68 onit 104 XYZ
#> 41 op1 104 XYZ
#> 42 op2 104 XYZ
#> 44 ora1 104 XYZ
#> 45 ora2 104 XYZ
#> 127 orgalk1cwi 104 XYZ
#> 126 orgalk1cwj 104 XYZ
#> 89 orgalk1i 104 XYZ
#> 88 orgalk1j 104 XYZ
#> 123 orgaro1cwi 104 XYZ
#> 122 orgaro1cwj 104 XYZ
#> 85 orgaro1i 104 XYZ
#> 84 orgaro1j 104 XYZ
#> 125 orgaro2cwi 104 XYZ
#> 124 orgaro2cwj 104 XYZ
#> 87 orgaro2i 104 XYZ
#> 86 orgaro2j 104 XYZ
#> 131 orgba1cwi 104 XYZ
#> 130 orgba1cwj 104 XYZ
#> 93 orgba1i 104 XYZ
#> 92 orgba1j 104 XYZ
#> 133 orgba2cwi 104 XYZ
#> 132 orgba2cwj 104 XYZ
#> 95 orgba2i 104 XYZ
#> 94 orgba2j 104 XYZ
#> 135 orgba3cwi 104 XYZ
#> 134 orgba3cwj 104 XYZ
#> 97 orgba3i 104 XYZ
#> 96 orgba3j 104 XYZ
#> 137 orgba4cwi 104 XYZ
#> 136 orgba4cwj 104 XYZ
#> 99 orgba4i 104 XYZ
#> 98 orgba4j 104 XYZ
#> 129 orgole1cwi 104 XYZ
#> 128 orgole1cwj 104 XYZ
#> 91 orgole1i 104 XYZ
#> 90 orgole1j 104 XYZ
#> 139 orgpacwi 104 XYZ
#> 138 orgpacwj 104 XYZ
#> 101 orgpai 104 XYZ
#> 100 orgpaj 104 XYZ
#> 179 P 104 XYZ
#> 190 P_HYD 104 XYZ
#> 269 P_STRAT 104 0
#> 263 P_TOP 104 0
#> 265 P00 104 0
#> 143 p25cwi 104 XYZ
#> 142 p25cwj 104 XYZ
#> 105 p25i 104 XYZ
#> 104 p25j 104 XYZ
#> 43 paa 104 XYZ
#> 49 pan 104 XYZ
#> 180 PB 104 XYZ
#> 295 PBLH 104 XY
#> 319 PC 104 XY
#> 318 PCB 104 XY
#> 166 PH 104 XYZ
#> 167 PHB 104 XYZ
#> 27 PHOTR201 104 XYZ
#> 28 PHOTR202 104 XYZ
#> 29 PHOTR203 104 XYZ
#> 26 PHOTR204 104 XYZ
#> 21 PM10 104 XYZ
#> 20 PM2_5_DRY 104 XYZ
#> 194 PSFC 104 XY
#> 150 PV 104 XYZ
#> 191 Q2 104 XY
#> 207 QCLOUD 104 XYZ
#> 297 QFX 104 XY
#> 211 QGRAUP 104 XYZ
#> 209 QICE 104 XYZ
#> 235 QKE 104 XYZ
#> 212 QNDROP 104 XYZ
#> 208 QRAIN 104 XYZ
#> 210 QSNOW 104 XYZ
#> 206 QVAPOR 104 XYZ
#> 272 RAINC 104 XY
#> 274 RAINNC 104 XY
#> 7 RAINPROD 104 XYZ
#> 273 RAINSH 104 XY
#> 184 RDN 104 Z
#> 183 RDNW 104 Z
#> 197 RDX 104 0
#> 198 RDY 104 0
#> 199 RESM 104 0
#> 10 ROUGH_COR 104 XY
#> 31 SAC 104 XYZ
#> 303 SAVE_TOPO_FROM_REAL 106 0
#> 219 SEAICE 104 XY
#> 107 seas 104 XYZ
#> 145 seascw 104 XYZ
#> 221 SFROFF 104 XY
#> 218 SH2O 104 XYZ
#> 213 SHDMAX 104 XY
#> 214 SHDMIN 104 XY
#> 259 SINALPHA 104 XY
#> 217 SMOIS 104 XYZ
#> 11 SMOIS_COR 104 XY
#> 215 SNOALB 104 XY
#> 229 SNOW 104 XY
#> 301 SNOWC 104 XY
#> 230 SNOWH 104 XY
#> 275 SNOWNC 104 XY
#> 30 SNU 104 XYZ
#> 32 so2 104 XYZ
#> 75 so4ai 104 XYZ
#> 74 so4aj 104 XYZ
#> 113 so4cwi 104 XYZ
#> 112 so4cwj 104 XYZ
#> 108 soila 104 XYZ
#> 146 soilcw 104 XYZ
#> 302 SR 104 XY
#> 322 SST 104 XY
#> 323 SST_INPUT 104 XY
#> 232 SSTSK 104 XY
#> 33 sulf 104 XYZ
#> 279 SWDOWN 104 XY
#> 281 SWNORM 104 XY
#> 168 T 104 XYZ
#> 264 T00 104 0
#> 192 T2 104 XY
#> 193 TH2 104 XY
#> 189 THIS_IS_AN_IDEAL_RUN 106 0
#> 169 THM 104 XYZ
#> 267 TISO 104 0
#> 245 TKE_PBL 104 XYZ
#> 266 TLP 104 0
#> 268 TLP_STRAT 104 0
#> 292 TMN 104 XY
#> 58 tol 104 XYZ
#> 61 tpan 104 XYZ
#> 262 TSK 104 XY
#> 172 TSK_FORCE 104 0
#> 175 TSK_FORCE_TEND 104 0
#> 216 TSLB 104 XYZ
#> 22 TSOA 104 XYZ
#> 163 U 104 XYZ
#> 195 U10 104 XY
#> 222 UDROFF 104 XY
#> 294 UST 104 XY
#> 9 UST_T 104 XY
#> 164 V 104 XYZ
#> 196 V10 104 XY
#> 244 VAR 104 XY
#> 162 VAR_SSO 104 XY
#> 225 VEGFRA 104 XY
#> 165 W 104 XYZ
#> 15 WET_DEP_FLUX 104 XYZ
#> 220 XICEM 104 XY
#> 293 XLAND 104 XY
#> 155 XLAT 104 XY
#> 283 XLAT_U 104 XY
#> 285 XLAT_V 104 XY
#> 156 XLONG 104 XY
#> 284 XLONG_U 104 XY
#> 286 XLONG_V 104 XY
#> 205 XTIME 104 0
#> 59 xyl 104 XYZ
#> 154 YCLNO2 104 XYZ
#> 200 ZETATOP 104 0
#> 158 ZNU 104 Z
#> 159 ZNW 104 Z
#> 160 ZS 104 Z
#> description
#> 110 Accumulation mode number
#> 148 Accumulation mode number in cloud
#> 227 ACCUMULATED GROUND HEAT FLUX
#> 299 ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE
#> 300 ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE
#> 60 ACO3 mixing ratio
#> 228 ACCUMULATED MELTED SNOW
#> 289 BACKGROUND ALBEDO
#> 287 ALBEDO
#> 39 ALD mixing ratio
#> 106 Coarse anthropogenic aerosols
#> 144 Coarse anthropogenic aerosols in cloud
#> 1 Aerosol Optical Depth
#> 2 Aerosol Optical Depth, 2d
#> 24 asoa dry mass
#> 3 Aerosol Optical Depth, top
#> 23 bsoa dry mass
#> 312 full levels, c1f = d bf / d eta, using znu
#> 310 half levels, c1h = d bf / d eta, using znw
#> 313 full levels, c2f = (1-c1f)*(p0-pt)
#> 311 half levels, c2h = (1-c1h)*(p0-pt)
#> 316 full levels, c3f = bf
#> 314 half levels, c3h = bh
#> 317 full levels, c4f = (eta-bf)*(p0-pt)+pt, using znw
#> 315 half levels, c4h = (eta-bh)*(p0-pt)+pt, using znu
#> 231 CANOPY WATER
#> 201 2nd order extrapolation constant
#> 202 2nd order extrapolation constant
#> 203 2nd order extrapolation constant
#> 187 extrapolation constant
#> 188 extrapolation constant
#> 83 Chloride conc. Aitken mode
#> 82 Chloride conc. Acc. mode
#> 288 COMPUTATIONAL GRID LATITUDE, SOUTH IS NEGATIVE
#> 121 Chloride conc. Aitken mode in cloud
#> 120 Chloride conc. Acc. mode in cloud
#> 278 CLOUD FRACTION
#> 6 CLOUD FRACTION
#> 152 n2o5 velocity
#> 54 CO mixing ratio
#> 111 Coarse mode number
#> 149 Coarse mode number in cloud
#> 260 Local cosine of map rotation
#> 233 COS of SOLAR ZENITH ANGLE
#> 69 CSL mixing ratio
#> 67 DCB mixing ratio
#> 25 dms oceanic concentrations
#> 186 d(eta) values between half (mass) levels
#> 185 d(eta) values between full (w) levels
#> 14 dry deposition flux
#> 13 dry deposition velocity
#> 12 dust dry deposition velocity
#> 19 O3 deposition velocity
#> 161 THICKNESSES OF SOIL LAYERS
#> 258 Coriolis cosine latitude term
#> 17 Actual biog emiss
#> 16 Actual biog emiss
#> 141 Elemental carbon Aitken mode in cloud
#> 140 Elemental carbon Acc. mode in cloud
#> 103 Elemental carbon Aitken mode
#> 102 Elemental carbon Acc. mode
#> 236 EDMF relative updraft area - moist updrafts
#> 240 EDMF entrainment - mean moist updrafts
#> 241 EDMF qc - mean moist updrafts
#> 239 EDMF qt - mean moist updrafts
#> 238 EDMF thetaL - mean moist updrafts
#> 237 EDMF vertical velocity - mean moist updrafts
#> 246 Length scale from PBL
#> 290 SURFACE EMISSIVITY
#> 53 ETH mixing ratio
#> 8 RAIN EVAPORATION RATE
#> 257 Coriolis sine latitude term
#> 181 upper weight for vertical stretching
#> 182 lower weight for vertical stretching
#> 151 n2o5 uptake by aerosol
#> 280 DOWNWARD LONG WAVE FLUX AT GROUND SURFACE
#> 65 GLY mixing ratio
#> 276 ACCUMULATED TOTAL GRID SCALE GRAUPEL
#> 226 GROUND HEAT FLUX
#> 38 H2O2 mixing ratio
#> 277 ACCUMULATED TOTAL GRID SCALE HAIL
#> 50 HC3 mixing ratio
#> 51 HC5 mixing ratio
#> 52 HC8 mixing ratio
#> 40 HCHO mixing ratio
#> 71 HCL mixing ratio
#> 296 UPWARD HEAT FLUX AT THE SURFACE
#> 170 SCM ideal surface sensible heat flux
#> 173 SCM ideal surface sensible heat flux tendency
#> 261 Terrain Height
#> 37 HNO3 mixing ratio
#> 63 HNO4 mixing ratio
#> 72 HO mixing ratio
#> 73 HO2 mixing ratio
#> 62 HONO mixing ratio
#> 4 Some chem thing, ICN_DIAG
#> 306 Array to hold seed for restart, RAND_PERT
#> 305 Array to hold seed for restart, SKEBS
#> 304 Array to hold seed for restart, SPPT
#> 307 Array to hold seed for restart, RAND_PERT2
#> 309 Array to hold seed for restart, RAND_PERT4
#> 308 Array to hold seed for restart, RAND_PERT3
#> 224 DOMINANT SOIL CATEGORY
#> 70 ISO mixing ratio
#> 204
#> 223 DOMINANT VEGETATION CATEGORY
#> 64 KET mixing ratio
#> 153 n2o5 het reaction rate
#> 234 LEAF AREA INDEX
#> 18 MODIS LAI vegetation mask for this date; 0=no dust produced (vegetation)
#> 321 LAKE MASK (1 FOR LAND, 0 FOR WATER)
#> 320 LAND MASK (1 FOR LAND, 0 FOR WATER)
#> 298 LATENT HEAT FLUX AT THE SURFACE
#> 171 SCM ideal surface latent heat flux
#> 174 SCM ideal surface latent heat flux tendency
#> 157 LAND USE CATEGORY
#> 247 Map scale factor on mass grid
#> 250 Map scale factor on mass grid, x direction
#> 251 Map scale factor on mass grid, y direction
#> 248 Map scale factor on u-grid
#> 252 Map scale factor on u-grid, x direction
#> 253 Map scale factor on u-grid, y direction
#> 249 Map scale factor on v-grid
#> 254 Map scale factor on v-grid, x direction
#> 256 Map scale factor on v-grid, y direction
#> 270 Max map factor in domain
#> 271 Max map factor in domain
#> 243 Maximum mass-flux (neg: all dry, pos: moist)
#> 255 Inverse map scale factor on v-grid, x direction
#> 66 MGLY mixing ratio
#> 176 perturbation dry air mass in column
#> 177 base state dry air mass in column
#> 47 N2O5 mixing ratio
#> 81 Sodium conc. Aitken mode
#> 80 Sodium conc. Acc. mode
#> 119 Sodium conc. Aitken mode in cloud
#> 118 Sodium conc. Acc. mode in cloud
#> 5 Some chem thing, NC_DIAG
#> 178 -
#> 46 NH3 mixing ratio
#> 77 Ammonium conc. Aitken mode
#> 76 Ammonium conc. Acc. mode
#> 115 Ammonium conc. Aitken mode in cloud
#> 114 Ammonium conc. Acc. mode in cloud
#> 35 NO mixing ratio
#> 34 NO2 mixing ratio
#> 48 NO3 mixing ratio
#> 79 Nitrate conc. Aitken mode
#> 78 Nitrate conc. Acc. mode
#> 117 Nitrate conc. Aitken mode in cloud
#> 116 Nitrate conc. Acc. mode in cloud
#> 291 RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET
#> 109 Aitken mode number
#> 147 Aitken mode number in cloud
#> 242 Number of updrafts per grid cell
#> 36 O3 mixing ratio
#> 55 OL2 mixing ratio
#> 57 OLI mixing ratio
#> 282 TOA OUTGOING LONG WAVE
#> 56 OLT mixing ratio
#> 68 ONIT mixing ratio
#> 41 OP1 mixing ratio
#> 42 OP2 mixing ratio
#> 44 ORA1 mixing ratio
#> 45 ORA2 mixing ratio
#> 127 SOA Anth. org. conc. from alkanes and others except aromatics Aitken mode in cloud
#> 126 SOA Anth. org. conc. from alkanes and others except aromatics Acc. mode in cloud
#> 89 SOA Anth. org. conc. from alkanes and others except aromatics Aitken mode
#> 88 SOA Anth. org. conc. from alkanes and others except aromatics Acc. mode
#> 123 SOA Anth. org. conc. from aromatics Aitken mode in cloud
#> 122 SOA Anth. org. conc. from aromatics Acc. mode in cloud
#> 85 SOA Anth. org. conc. from aromatics Aitken mode
#> 84 SOA Anth. org. conc. from aromatics Acc. mode
#> 125 SOA Anth. org. conc. from aromatics Aitken mode in cloud
#> 124 SOA Anth. org. conc. from aromatics Acc. mode in cloud
#> 87 SOA Anth. org. conc. from aromatics Aitken mode
#> 86 SOA Anth. org. conc. from aromatics Acc. mode
#> 131 SOA Biog. org. conc. from aromatics Aitken mode in cloud
#> 130 SOA Biog. org. conc. from aromatics Acc. mode in cloud
#> 93 SOA Biog. org. conc. from aromatics Aitken mode
#> 92 SOA Biog. org. conc. from aromatics Acc. mode
#> 133 SOA Biog. org. conc. from aromatics Aitken mode in cloud
#> 132 SOA Biog. org. conc. from aromatics Acc. mode in cloud
#> 95 SOA Biog. org. conc. from aromatics Aitken mode
#> 94 SOA Biog. org. conc. from aromatics Acc. mode
#> 135 SOA Biog. org. conc. from aromatics Aitken mode in cloud
#> 134 SOA Biog. org. conc. from aromatics Acc. mode in cloud
#> 97 SOA Biog. org. conc. from aromatics Aitken mode
#> 96 SOA Biog. org. conc. from aromatics Acc. mode
#> 137 SOA Biog. org. conc. from aromatics Aitken mode in cloud
#> 136 SOA Biog. org. conc. from aromatics Acc. mode in cloud
#> 99 SOA Biog. org. conc. from aromatics Aitken mode
#> 98 SOA Biog. org. conc. from aromatics Acc. mode
#> 129 SOA Anth. org. conc. from alkenes and others except aromatics Aitken mode in cloud
#> 128 SOA Anth. org. conc. from alkenes and others except aromatics Acc. mode in cloud
#> 91 SOA Anth. org. conc. from alkenes and others except aromatics Aitken mode
#> 90 SOA Anth. org. conc. from alkenes and others except aromatics Acc. mode
#> 139 Prim. anth. org. conc. from aromatics Aitken mode in cloud
#> 138 Prim. anth. org. conc. from aromatics Acc. mode in cloud
#> 101 Prim. anth. org. conc. from aromatics Aitken mode
#> 100 Prim. anth. org. conc. from aromatics Acc. mode
#> 179 perturbation pressure
#> 190 hydrostatic pressure
#> 269 BASE STATE PRESSURE AT BOTTOM OF STRATOSPHERE
#> 263 PRESSURE TOP OF THE MODEL
#> 265 BASE STATE PRESURE
#> 143 Primary PM2.5 Aitken mode in cloud
#> 142 Primary PM2.5 Acc. mode in cloud
#> 105 Primary PM2.5 Aitken mode
#> 104 Primary PM2.5 Acc. mode
#> 43 PAA mixing ratio
#> 49 PAN mixing ratio
#> 180 BASE STATE PRESSURE
#> 295 PBL HEIGHT
#> 319 perturbation dry air mass in column
#> 318 base state dry air mass in column
#> 166 perturbation geopotential
#> 167 base-state geopotential
#> 27 cl2 photolysis rate
#> 28 hocl photolysis rate
#> 29 fmcl photolysis rate
#> 26 CLNO2 Photolysis Rate
#> 21 pm10 dry mass
#> 20 pm2.5 aerosol dry mass
#> 194 SFC PRESSURE
#> 150 Potential Vorticity
#> 191 QV at 2 M
#> 207 Cloud water mixing ratio
#> 297 UPWARD MOISTURE FLUX AT THE SURFACE
#> 211 Graupel mixing ratio
#> 209 Ice mixing ratio
#> 235 twice TKE from MYNN
#> 212 Droplet number mixing ratio
#> 208 Rain water mixing ratio
#> 210 Snow mixing ratio
#> 206 Water vapor mixing ratio
#> 272 ACCUMULATED TOTAL CUMULUS PRECIPITATION
#> 274 ACCUMULATED TOTAL GRID SCALE PRECIPITATION
#> 7 TOTAL RAIN PRODUCTION RATE
#> 273 ACCUMULATED SHALLOW CUMULUS PRECIPITATION
#> 184 inverse d(eta) values between half (mass) levels
#> 183 inverse d(eta) values between full (w) levels
#> 197 INVERSE X GRID LENGTH
#> 198 INVERSE Y GRID LENGTH
#> 199 TIME WEIGHT CONSTANT FOR SMALL STEPS
#> 10 roughness elements correction
#> 31 2nd moment Aitken mode
#> 303 1=original topo from real/0=topo modified by WRF
#> 219 SEA ICE FLAG
#> 107 Coarse marine aerosols
#> 145 Coarse marine aerosols in cloud
#> 221 SURFACE RUNOFF
#> 218 SOIL LIQUID WATER
#> 213 ANNUAL MAX VEG FRACTION
#> 214 ANNUAL MIN VEG FRACTION
#> 259 Local sine of map rotation
#> 217 SOIL MOISTURE
#> 11 soil moisture correction
#> 215 ANNUAL MAX SNOW ALBEDO IN FRACTION
#> 229 SNOW WATER EQUIVALENT
#> 301 FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)
#> 230 PHYSICAL SNOW DEPTH
#> 275 ACCUMULATED TOTAL GRID SCALE SNOW AND ICE
#> 30 2nd moment Aitken mode
#> 32 SO2 mixing ratio
#> 75 Sulfate conc. Aitken mode
#> 74 Sulfate conc. Acc. mode
#> 113 Sulfate conc. Aitken mode in cloud
#> 112 Sulfate conc. Acc. mode in cloud
#> 108 Coarse soil-derived aerosols
#> 146 Coarse soil-derived aerosols in cloud
#> 302 fraction of frozen precipitation
#> 322 SEA SURFACE TEMPERATURE
#> 323 SEA SURFACE TEMPERATURE FROM WRFLOWINPUT FILE
#> 232 SKIN SEA SURFACE TEMPERATURE
#> 33 SULF mixing ratio
#> 279 DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE
#> 281 NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT)
#> 168 perturbation potential temperature theta-t0
#> 264 BASE STATE TEMPERATURE
#> 192 TEMP at 2 M
#> 193 POT TEMP at 2 M
#> 189 T/F flag: this is an ARW ideal simulation
#> 169 either 1) pert moist pot temp=(1+Rv/Rd Qv)*(theta)-T0, or 2) pert dry pot temp=t
#> 267 TEMP AT WHICH THE BASE T TURNS CONST
#> 245 TKE from PBL
#> 266 BASE STATE LAPSE RATE
#> 268 BASE STATE LAPSE RATE (DT/D(LN(P)) IN STRATOSPHERE
#> 292 SOIL TEMPERATURE AT LOWER BOUNDARY
#> 58 TOL mixing ratio
#> 61 TPAN mixing ratio
#> 262 SURFACE SKIN TEMPERATURE
#> 172 SCM ideal surface skin temperature
#> 175 SCM ideal surface skin temperature tendency
#> 216 SOIL TEMPERATURE
#> 22 soa dry mass
#> 163 x-wind component
#> 195 U at 10 M
#> 222 UNDERGROUND RUNOFF
#> 294 U* IN SIMILARITY THEORY
#> 9 Threshold Friction Velocity
#> 164 y-wind component
#> 196 V at 10 M
#> 244 OROGRAPHIC VARIANCE
#> 162 variance of subgrid-scale orography
#> 225 VEGETATION FRACTION
#> 165 z-wind component
#> 15 column wet scavening flux
#> 220 SEA ICE FLAG (PREVIOUS STEP)
#> 293 LAND MASK (1 FOR LAND, 2 FOR WATER)
#> 155 LATITUDE, SOUTH IS NEGATIVE
#> 283 LATITUDE, SOUTH IS NEGATIVE
#> 285 LATITUDE, SOUTH IS NEGATIVE
#> 156 LONGITUDE, WEST IS NEGATIVE
#> 284 LONGITUDE, WEST IS NEGATIVE
#> 286 LONGITUDE, WEST IS NEGATIVE
#> 205 minutes since 2014-10-25 00:00:00
#> 59 XYL mixing ratio
#> 154 clno2 yield from n2o5 het
#> 200 ZETA AT MODEL TOP
#> 158 eta values on half (mass) levels
#> 159 eta values on full (w) levels
#> 160 DEPTHS OF CENTERS OF SOIL LAYERS
#> units stagger coordinates
#> 110 /kg-dryair XLONG XLAT XTIME
#> 148 /kg-dryair XLONG XLAT XTIME
#> 227 J m-2 XLONG XLAT XTIME
#> 299 J m-2 XLONG XLAT XTIME
#> 300 J m-2 XLONG XLAT XTIME
#> 60 ppmv XLONG XLAT XTIME
#> 228 kg m-2 XLONG XLAT XTIME
#> 289 XLONG XLAT XTIME
#> 287 - XLONG XLAT XTIME
#> 39 ppmv XLONG XLAT XTIME
#> 106 ug/kg-dryair XLONG XLAT XTIME
#> 144 ug/kg-dryair XLONG XLAT XTIME
#> 1 XLONG XLAT XTIME
#> 2 XLONG XLAT XTIME
#> 24 ug m^-3 XLONG XLAT XTIME
#> 3 XLONG XLAT XTIME
#> 23 ug m^-3 XLONG XLAT XTIME
#> 312 Dimensionless Z 104
#> 310 Dimensionless 104
#> 313 Pa Z 104
#> 311 Pa 104
#> 316 Dimensionless Z 104
#> 314 Dimensionless 104
#> 317 Pa Z 104
#> 315 Pa 104
#> 231 kg m-2 XLONG XLAT XTIME
#> 201 104
#> 202 104
#> 203 104
#> 187 104
#> 188 104
#> 83 ug/kg-dryair XLONG XLAT XTIME
#> 82 ug/kg-dryair XLONG XLAT XTIME
#> 288 degree_north XLONG XLAT XTIME
#> 121 ug/kg-dryair XLONG XLAT XTIME
#> 120 ug/kg-dryair XLONG XLAT XTIME
#> 278 XLONG XLAT XTIME
#> 6 - XLONG XLAT XTIME
#> 152 m/s XLONG XLAT XTIME
#> 54 ppmv XLONG XLAT XTIME
#> 111 /kg-dryair XLONG XLAT XTIME
#> 149 /kg-dryair XLONG XLAT XTIME
#> 260 XLONG XLAT XTIME
#> 233 dimensionless XLONG XLAT XTIME
#> 69 ppmv XLONG XLAT XTIME
#> 67 ppmv XLONG XLAT XTIME
#> 25 nM/L XLONG XLAT XTIME
#> 186 104
#> 185 104
#> 14 mol or ug m^-2 Z XLONG XLAT XTIME
#> 13 cm/s Z XLONG XLAT XTIME
#> 12 m/s XLONG XLAT XTIME
#> 19 cm/s XLONG XLAT XTIME
#> 161 m Z 104
#> 258 s-1 XLONG XLAT XTIME
#> 17 mol km^-2 hr^-1 XLONG XLAT XTIME
#> 16 mol km^-2 hr^-1 XLONG XLAT XTIME
#> 141 ug/kg-dryair XLONG XLAT XTIME
#> 140 ug/kg-dryair XLONG XLAT XTIME
#> 103 ug/kg-dryair XLONG XLAT XTIME
#> 102 ug/kg-dryair XLONG XLAT XTIME
#> 236 - XLONG XLAT XTIME
#> 240 m-1 XLONG XLAT XTIME
#> 241 kg kg-1 XLONG XLAT XTIME
#> 239 kg kg-1 XLONG XLAT XTIME
#> 238 K XLONG XLAT XTIME
#> 237 m s-1 XLONG XLAT XTIME
#> 246 m Z XLONG XLAT XTIME
#> 290 XLONG XLAT XTIME
#> 53 ppmv XLONG XLAT XTIME
#> 8 s-1 XLONG XLAT XTIME
#> 257 s-1 XLONG XLAT XTIME
#> 181 104
#> 182 104
#> 151 numerical value XLONG XLAT XTIME
#> 280 W m-2 XLONG XLAT XTIME
#> 65 ppmv XLONG XLAT XTIME
#> 276 mm XLONG XLAT XTIME
#> 226 W m-2 XLONG XLAT XTIME
#> 38 ppmv XLONG XLAT XTIME
#> 277 mm XLONG XLAT XTIME
#> 50 ppmv XLONG XLAT XTIME
#> 51 ppmv XLONG XLAT XTIME
#> 52 ppmv XLONG XLAT XTIME
#> 40 ppmv XLONG XLAT XTIME
#> 71 ppmv XLONG XLAT XTIME
#> 296 W m-2 XLONG XLAT XTIME
#> 170 W m-2 104
#> 173 W m-2 s-1 104
#> 261 m XLONG XLAT XTIME
#> 37 ppmv XLONG XLAT XTIME
#> 63 ppmv XLONG XLAT XTIME
#> 72 ppmv XLONG XLAT XTIME
#> 73 ppmv XLONG XLAT XTIME
#> 62 ppmv XLONG XLAT XTIME
#> 4 XLONG XLAT XTIME
#> 306 Z 106
#> 305 Z 106
#> 304 Z 106
#> 307 Z 106
#> 309 Z 106
#> 308 Z 106
#> 224 XLONG XLAT XTIME
#> 70 ppmv XLONG XLAT XTIME
#> 204 106
#> 223 XLONG XLAT XTIME
#> 64 ppmv XLONG XLAT XTIME
#> 153 s-1 XLONG XLAT XTIME
#> 234 m-2/m-2 XLONG XLAT XTIME
#> 18 none XLONG XLAT XTIME
#> 321 XLONG XLAT XTIME
#> 320 XLONG XLAT XTIME
#> 298 W m-2 XLONG XLAT XTIME
#> 171 W m-2 104
#> 174 W m-2 s-1 104
#> 157 XLONG XLAT XTIME
#> 247 XLONG XLAT XTIME
#> 250 XLONG XLAT XTIME
#> 251 XLONG XLAT XTIME
#> 248 X XLONG_U XLAT_U XTIME
#> 252 X XLONG_U XLAT_U XTIME
#> 253 X XLONG_U XLAT_U XTIME
#> 249 Y XLONG_V XLAT_V XTIME
#> 254 Y XLONG_V XLAT_V XTIME
#> 256 Y XLONG_V XLAT_V XTIME
#> 270 104
#> 271 104
#> 243 m/s * area XLONG XLAT XTIME
#> 255 Y XLONG_V XLAT_V XTIME
#> 66 ppmv XLONG XLAT XTIME
#> 176 Pa XLONG XLAT XTIME
#> 177 Pa XLONG XLAT XTIME
#> 47 ppmv XLONG XLAT XTIME
#> 81 ug/kg-dryair XLONG XLAT XTIME
#> 80 ug/kg-dryair XLONG XLAT XTIME
#> 119 ug/kg-dryair XLONG XLAT XTIME
#> 118 ug/kg-dryair XLONG XLAT XTIME
#> 5 XLONG XLAT XTIME
#> 178 - XLONG XLAT XTIME
#> 46 ppmv XLONG XLAT XTIME
#> 77 ug/kg-dryair XLONG XLAT XTIME
#> 76 ug/kg-dryair XLONG XLAT XTIME
#> 115 ug/kg-dryair XLONG XLAT XTIME
#> 114 ug/kg-dryair XLONG XLAT XTIME
#> 35 ppmv XLONG XLAT XTIME
#> 34 ppmv XLONG XLAT XTIME
#> 48 ppmv XLONG XLAT XTIME
#> 79 ug/kg-dryair XLONG XLAT XTIME
#> 78 ug/kg-dryair XLONG XLAT XTIME
#> 117 ug/kg-dryair XLONG XLAT XTIME
#> 116 ug/kg-dryair XLONG XLAT XTIME
#> 291 W m{-2} XLONG XLAT XTIME
#> 109 /kg-dryair XLONG XLAT XTIME
#> 147 /kg-dryair XLONG XLAT XTIME
#> 242 XLONG XLAT XTIME
#> 36 ppmv XLONG XLAT XTIME
#> 55 ppmv XLONG XLAT XTIME
#> 57 ppmv XLONG XLAT XTIME
#> 282 W m-2 XLONG XLAT XTIME
#> 56 ppmv XLONG XLAT XTIME
#> 68 ppmv XLONG XLAT XTIME
#> 41 ppmv XLONG XLAT XTIME
#> 42 ppmv XLONG XLAT XTIME
#> 44 ppmv XLONG XLAT XTIME
#> 45 ppmv XLONG XLAT XTIME
#> 127 ug/kg-dryair XLONG XLAT XTIME
#> 126 ug/kg-dryair XLONG XLAT XTIME
#> 89 ug/kg-dryair XLONG XLAT XTIME
#> 88 ug/kg-dryair XLONG XLAT XTIME
#> 123 ug/kg-dryair XLONG XLAT XTIME
#> 122 ug/kg-dryair XLONG XLAT XTIME
#> 85 ug/kg-dryair XLONG XLAT XTIME
#> 84 ug/kg-dryair XLONG XLAT XTIME
#> 125 ug/kg-dryair XLONG XLAT XTIME
#> 124 ug/kg-dryair XLONG XLAT XTIME
#> 87 ug/kg-dryair XLONG XLAT XTIME
#> 86 ug/kg-dryair XLONG XLAT XTIME
#> 131 ug/kg-dryair XLONG XLAT XTIME
#> 130 ug/kg-dryair XLONG XLAT XTIME
#> 93 ug/kg-dryair XLONG XLAT XTIME
#> 92 ug/kg-dryair XLONG XLAT XTIME
#> 133 ug/kg-dryair XLONG XLAT XTIME
#> 132 ug/kg-dryair XLONG XLAT XTIME
#> 95 ug/kg-dryair XLONG XLAT XTIME
#> 94 ug/kg-dryair XLONG XLAT XTIME
#> 135 ug/kg-dryair XLONG XLAT XTIME
#> 134 ug/kg-dryair XLONG XLAT XTIME
#> 97 ug/kg-dryair XLONG XLAT XTIME
#> 96 ug/kg-dryair XLONG XLAT XTIME
#> 137 ug/kg-dryair XLONG XLAT XTIME
#> 136 ug/kg-dryair XLONG XLAT XTIME
#> 99 ug/kg-dryair XLONG XLAT XTIME
#> 98 ug/kg-dryair XLONG XLAT XTIME
#> 129 ug/kg-dryair XLONG XLAT XTIME
#> 128 ug/kg-dryair XLONG XLAT XTIME
#> 91 ug/kg-dryair XLONG XLAT XTIME
#> 90 ug/kg-dryair XLONG XLAT XTIME
#> 139 ug/kg-dryair XLONG XLAT XTIME
#> 138 ug/kg-dryair XLONG XLAT XTIME
#> 101 ug/kg-dryair XLONG XLAT XTIME
#> 100 ug/kg-dryair XLONG XLAT XTIME
#> 179 Pa XLONG XLAT XTIME
#> 190 Pa XLONG XLAT XTIME
#> 269 Pa 104
#> 263 Pa 104
#> 265 Pa 104
#> 143 ug/kg-dryair XLONG XLAT XTIME
#> 142 ug/kg-dryair XLONG XLAT XTIME
#> 105 ug/kg-dryair XLONG XLAT XTIME
#> 104 ug/kg-dryair XLONG XLAT XTIME
#> 43 ppmv XLONG XLAT XTIME
#> 49 ppmv XLONG XLAT XTIME
#> 180 Pa XLONG XLAT XTIME
#> 295 m XLONG XLAT XTIME
#> 319 Pa XLONG XLAT XTIME
#> 318 Pa XLONG XLAT XTIME
#> 166 m2 s-2 Z XLONG XLAT XTIME
#> 167 m2 s-2 Z XLONG XLAT XTIME
#> 27 min{-1} XLONG XLAT XTIME
#> 28 min{-1} XLONG XLAT XTIME
#> 29 min{-1} XLONG XLAT XTIME
#> 26 min{-1} XLONG XLAT XTIME
#> 21 ug m^-3 XLONG XLAT XTIME
#> 20 ug m^-3 XLONG XLAT XTIME
#> 194 Pa XLONG XLAT XTIME
#> 150 pvu XLONG XLAT XTIME
#> 191 kg kg-1 XLONG XLAT XTIME
#> 207 kg kg-1 XLONG XLAT XTIME
#> 297 kg m-2 s-1 XLONG XLAT XTIME
#> 211 kg kg-1 XLONG XLAT XTIME
#> 209 kg kg-1 XLONG XLAT XTIME
#> 235 m2 s-2 XLONG XLAT XTIME
#> 212 kg-1 XLONG XLAT XTIME
#> 208 kg kg-1 XLONG XLAT XTIME
#> 210 kg kg-1 XLONG XLAT XTIME
#> 206 kg kg-1 XLONG XLAT XTIME
#> 272 mm XLONG XLAT XTIME
#> 274 mm XLONG XLAT XTIME
#> 7 s-1 XLONG XLAT XTIME
#> 273 mm XLONG XLAT XTIME
#> 184 104
#> 183 104
#> 197 104
#> 198 104
#> 199 104
#> 10 XLONG XLAT XTIME
#> 31 m2 m-3 XLONG XLAT XTIME
#> 303 flag 106
#> 219 XLONG XLAT XTIME
#> 107 ug/kg-dryair XLONG XLAT XTIME
#> 145 ug/kg-dryair XLONG XLAT XTIME
#> 221 mm XLONG XLAT XTIME
#> 218 m3 m-3 Z XLONG XLAT XTIME
#> 213 XLONG XLAT XTIME
#> 214 XLONG XLAT XTIME
#> 259 XLONG XLAT XTIME
#> 217 m3 m-3 Z XLONG XLAT XTIME
#> 11 XLONG XLAT XTIME
#> 215 XLONG XLAT XTIME
#> 229 kg m-2 XLONG XLAT XTIME
#> 301 XLONG XLAT XTIME
#> 230 m XLONG XLAT XTIME
#> 275 mm XLONG XLAT XTIME
#> 30 m2 m-3 XLONG XLAT XTIME
#> 32 ppmv XLONG XLAT XTIME
#> 75 ug/kg-dryair XLONG XLAT XTIME
#> 74 ug/kg-dryair XLONG XLAT XTIME
#> 113 ug/kg-dryair XLONG XLAT XTIME
#> 112 ug/kg-dryair XLONG XLAT XTIME
#> 108 ug/kg-dryair XLONG XLAT XTIME
#> 146 ug/kg-dryair XLONG XLAT XTIME
#> 302 - XLONG XLAT XTIME
#> 322 K XLONG XLAT XTIME
#> 323 K XLONG XLAT XTIME
#> 232 K XLONG XLAT XTIME
#> 33 ppmv XLONG XLAT XTIME
#> 279 W m-2 XLONG XLAT XTIME
#> 281 W m-2 XLONG XLAT XTIME
#> 168 K XLONG XLAT XTIME
#> 264 K 104
#> 192 K XLONG XLAT XTIME
#> 193 K XLONG XLAT XTIME
#> 189 - 106
#> 169 K XLONG XLAT XTIME
#> 267 K 104
#> 245 m2 s-2 Z XLONG XLAT XTIME
#> 266 104
#> 268 K 104
#> 292 K XLONG XLAT XTIME
#> 58 ppmv XLONG XLAT XTIME
#> 61 ppmv XLONG XLAT XTIME
#> 262 K XLONG XLAT XTIME
#> 172 W m-2 104
#> 175 W m-2 s-1 104
#> 216 K Z XLONG XLAT XTIME
#> 22 ug m^-3 XLONG XLAT XTIME
#> 163 m s-1 X XLONG_U XLAT_U XTIME
#> 195 m s-1 XLONG XLAT XTIME
#> 222 mm XLONG XLAT XTIME
#> 294 m s-1 XLONG XLAT XTIME
#> 9 m s-1 XLONG XLAT XTIME
#> 164 m s-1 Y XLONG_V XLAT_V XTIME
#> 196 m s-1 XLONG XLAT XTIME
#> 244 XLONG XLAT XTIME
#> 162 m2 XLONG XLAT XTIME
#> 225 XLONG XLAT XTIME
#> 165 m s-1 Z XLONG XLAT XTIME
#> 15 mmol or ug m^-2 Z XLONG XLAT XTIME
#> 220 XLONG XLAT XTIME
#> 293 XLONG XLAT XTIME
#> 155 degree_north XLONG XLAT
#> 283 degree_north X XLONG_U XLAT_U
#> 285 degree_north Y XLONG_V XLAT_V
#> 156 degree_east XLONG XLAT
#> 284 degree_east X XLONG_U XLAT_U
#> 286 degree_east Y XLONG_V XLAT_V
#> 205 minutes since 2014-10-25 00:00:00 104
#> 59 ppmv XLONG XLAT XTIME
#> 154 numerical value XLONG XLAT XTIME
#> 200 104
#> 158 104
#> 159 Z 104
#> 160 m Z 104Now we can select some variables
(vars <- dt[dt$MemoryOrder == "XYZ", ]$vars)
#> [1] "ac0" "ac0cw" "aco3" "ald" "antha"
#> [6] "anthcw" "AOD_OUT" "ASOA" "BSOA" "clai"
#> [11] "claj" "clcwi" "clcwj" "CLDFRA" "CLDFRA2"
#> [16] "CN2O5" "co" "corn" "corncw" "csl"
#> [21] "dcb" "DRY_DEP_FLUX" "DRY_DEP_LEN" "dvel_o3" "eccwi"
#> [26] "eccwj" "eci" "ecj" "EDMF_A" "EDMF_ENT"
#> [31] "EDMF_QC" "EDMF_QT" "EDMF_THL" "EDMF_W" "EL_PBL"
#> [36] "eth" "EVAPPROD" "GAMN2O5" "gly" "h2o2"
#> [41] "hc3" "hc5" "hc8" "hcho" "hcl"
#> [46] "hno3" "hno4" "ho" "ho2" "hono"
#> [51] "ICN_DIAG" "iso" "ket" "KN2O5" "mgly"
#> [56] "n2o5" "naai" "naaj" "nacwi" "nacwj"
#> [61] "NC_DIAG" "nh3" "nh4ai" "nh4aj" "nh4cwi"
#> [66] "nh4cwj" "no" "no2" "no3" "no3ai"
#> [71] "no3aj" "no3cwi" "no3cwj" "nu0" "nu0cw"
#> [76] "o3" "ol2" "oli" "olt" "onit"
#> [81] "op1" "op2" "ora1" "ora2" "orgalk1cwi"
#> [86] "orgalk1cwj" "orgalk1i" "orgalk1j" "orgaro1cwi" "orgaro1cwj"
#> [91] "orgaro1i" "orgaro1j" "orgaro2cwi" "orgaro2cwj" "orgaro2i"
#> [96] "orgaro2j" "orgba1cwi" "orgba1cwj" "orgba1i" "orgba1j"
#> [101] "orgba2cwi" "orgba2cwj" "orgba2i" "orgba2j" "orgba3cwi"
#> [106] "orgba3cwj" "orgba3i" "orgba3j" "orgba4cwi" "orgba4cwj"
#> [111] "orgba4i" "orgba4j" "orgole1cwi" "orgole1cwj" "orgole1i"
#> [116] "orgole1j" "orgpacwi" "orgpacwj" "orgpai" "orgpaj"
#> [121] "P" "P_HYD" "p25cwi" "p25cwj" "p25i"
#> [126] "p25j" "paa" "pan" "PB" "PH"
#> [131] "PHB" "PHOTR201" "PHOTR202" "PHOTR203" "PHOTR204"
#> [136] "PM10" "PM2_5_DRY" "PV" "QCLOUD" "QGRAUP"
#> [141] "QICE" "QKE" "QNDROP" "QRAIN" "QSNOW"
#> [146] "QVAPOR" "RAINPROD" "SAC" "seas" "seascw"
#> [151] "SH2O" "SMOIS" "SNU" "so2" "so4ai"
#> [156] "so4aj" "so4cwi" "so4cwj" "soila" "soilcw"
#> [161] "sulf" "T" "THM" "TKE_PBL" "tol"
#> [166] "tpan" "TSLB" "TSOA" "U" "V"
#> [171] "W" "WET_DEP_FLUX" "xyl" "YCLNO2"Now get some statistics
wrf_summary(wrfo, c("o3", "PM2_5_DRY", "T", "U", "V"))
#>
|
| | 0%
|
|============== | 20%
|
|============================ | 40%
|
|========================================== | 60%
|
|======================================================== | 80%
|
|======================================================================| 100%
#> Min. 1st Qu. Median Mean 3rd Qu.
#> o3 1.000000e-16 0.0001494586 0.003905759 0.00645222 0.01011842
#> PM2_5_DRY 1.000149e-04 0.0001599413 0.012459201 0.13730251 0.05450001
#> T -7.659851e+00 6.2148132324 22.712966919 34.82713373 49.47671509
#> U -2.072658e+01 -2.2874763012 1.582949162 2.71344012 6.05402517
#> V -1.938167e+01 -3.7864468098 -0.764581144 0.10207729 2.48909324
#> Max. sum
#> o3 0.08323991 37102.53
#> PM2_5_DRY 72.09506989 789537.79
#> T 178.68499756 200268278.11
#> U 32.16123199 15782583.34
#> V 27.93202972 594227.03