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)
#> The legacy packages maptools, rgdal, and rgeos, underpinning the sp package,
#> which was just loaded, will retire in October 2023.
#> Please refer to R-spatial evolution reports for details, especially
#> https://r-spatial.org/r/2023/05/15/evolution4.html.
#> It may be desirable to make the sf package available;
#> package maintainers should consider adding sf to Suggests:.
#> The sp package is now running under evolution status 2
#> (status 2 uses the sf package in place of rgdal)
library(raster)
#> Loading required package: sp
library(stars)
#> Loading required package: abind
#> Loading required package: sf
#> Linking to GEOS 3.11.1, GDAL 3.6.2, PROJ 9.1.1; sf_use_s2() is TRUE
library(cptcity)
library(sf)Based on NCL:
First lets get a summary of a WRF output file
wrfo <- "/media/sergio/ext5/WRF4/WRF/test/em_real/wrfout_d01_2014-10-03_00:00:00"
(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 87)
#> vars FieldType MemoryOrder
#> 156 ACGRDFLX 104 XY
#> 228 ACHFX 104 XY
#> 229 ACLHF 104 XY
#> 64 aco3 104 XYZ
#> 157 ACSNOM 104 XY
#> 218 ALBBCK 104 XY
#> 216 ALBEDO 104 XY
#> 43 ald 104 XYZ
#> 1 AOD_OUT 104 XYZ
#> 2 AOD2D_OUT 104 XY
#> 3 ATOP2D_OUT 104 XY
#> 241 C1F 104 Z
#> 239 C1H 104 Z
#> 242 C2F 104 Z
#> 240 C2H 104 Z
#> 245 C3F 104 Z
#> 243 C3H 104 Z
#> 246 C4F 104 Z
#> 244 C4H 104 Z
#> 160 CANWAT 104 XY
#> 130 CF1 104 0
#> 131 CF2 104 0
#> 132 CF3 104 0
#> 116 CFN 104 0
#> 117 CFN1 104 0
#> 76 ch4 104 XYZ
#> 217 CLAT 104 XY
#> 207 CLDFRA 104 XYZ
#> 6 CLDFRA2 104 XYZ
#> 81 CN2O5 104 XYZ
#> 58 co 104 XYZ
#> 75 co2 104 XYZ
#> 189 COSALPHA 104 XY
#> 162 COSZEN 104 XY
#> 73 csl 104 XYZ
#> 71 dcb 104 XYZ
#> 29 DMS_0 104 XY
#> 115 DN 104 Z
#> 114 DNW 104 Z
#> 22 DRY_DEP_LEN 104 XYZ
#> 12 DRYDEPVEL 104 XY
#> 26 dvel_o3 104 XYZ
#> 90 DZS 104 Z
#> 187 E 104 XY
#> 24 EBIO_API 104 XY
#> 23 EBIO_ISO 104 XY
#> 13 ebu_oc 104 XYZ
#> 165 EDMF_A 104 XYZ
#> 169 EDMF_ENT 104 XYZ
#> 170 EDMF_QC 104 XYZ
#> 168 EDMF_QT 104 XYZ
#> 167 EDMF_THL 104 XYZ
#> 166 EDMF_W 104 XYZ
#> 175 EL_PBL 104 XYZ
#> 219 EMISS 104 XY
#> 57 eth 104 XYZ
#> 8 EVAPPROD 104 XYZ
#> 186 F 104 XY
#> 19 FIRESIZE_AGEF 104 XY
#> 21 FIRESIZE_AGGR 104 XY
#> 20 FIRESIZE_AGSV 104 XY
#> 18 FIRESIZE_AGTF 104 XY
#> 110 FNM 104 Z
#> 111 FNP 104 Z
#> 80 GAMN2O5 104 XYZ
#> 209 GLW 104 XY
#> 69 gly 104 XYZ
#> 205 GRAUPELNC 104 XY
#> 155 GRDFLX 104 XY
#> 42 h2o2 104 XYZ
#> 206 HAILNC 104 XY
#> 54 hc3 104 XYZ
#> 55 hc5 104 XYZ
#> 56 hc8 104 XYZ
#> 44 hcho 104 XYZ
#> 225 HFX 104 XY
#> 99 HFX_FORCE 104 0
#> 102 HFX_FORCE_TEND 104 0
#> 190 HGT 104 XY
#> 41 hno3 104 XYZ
#> 67 hno4 104 XYZ
#> 77 ho 104 XYZ
#> 78 ho2 104 XYZ
#> 66 hono 104 XYZ
#> 4 ICN_DIAG 104 XYZ
#> 235 ISEEDARR_RAND_PERTURB 106 Z
#> 234 ISEEDARR_SKEBS 106 Z
#> 233 ISEEDARR_SPPT 106 Z
#> 236 ISEEDARRAY_SPP_CONV 106 Z
#> 238 ISEEDARRAY_SPP_LSM 106 Z
#> 237 ISEEDARRAY_SPP_PBL 106 Z
#> 153 ISLTYP 106 XY
#> 74 iso 104 XYZ
#> 133 ITIMESTEP 106 0
#> 152 IVGTYP 106 XY
#> 68 ket 104 XYZ
#> 82 KN2O5 104 XYZ
#> 163 LAI 104 XY
#> 25 LAI_VEGMASK 104 XY
#> 250 LAKEMASK 104 XY
#> 249 LANDMASK 104 XY
#> 227 LH 104 XY
#> 100 LH_FORCE 104 0
#> 103 LH_FORCE_TEND 104 0
#> 86 LU_INDEX 104 XY
#> 176 MAPFAC_M 104 XY
#> 179 MAPFAC_MX 104 XY
#> 180 MAPFAC_MY 104 XY
#> 177 MAPFAC_U 104 XY
#> 181 MAPFAC_UX 104 XY
#> 182 MAPFAC_UY 104 XY
#> 178 MAPFAC_V 104 XY
#> 183 MAPFAC_VX 104 XY
#> 185 MAPFAC_VY 104 XY
#> 199 MAX_MSTFX 104 0
#> 200 MAX_MSTFY 104 0
#> 172 MAXMF 104 XY
#> 15 MEAN_FCT_AGEF 104 XY
#> 17 MEAN_FCT_AGGR 104 XY
#> 16 MEAN_FCT_AGSV 104 XY
#> 14 MEAN_FCT_AGTF 104 XY
#> 184 MF_VX_INV 104 XY
#> 70 mgly 104 XYZ
#> 105 MU 104 XY
#> 106 MUB 104 XY
#> 51 n2o5 104 XYZ
#> 5 NC_DIAG 104 XYZ
#> 107 NEST_POS 104 XY
#> 50 nh3 104 XYZ
#> 39 no 104 XYZ
#> 38 no2 104 XYZ
#> 52 no3 104 XYZ
#> 220 NOAHRES 104 XY
#> 171 NUPDRAFT 106 XY
#> 40 o3 104 XYZ
#> 59 ol2 104 XYZ
#> 61 oli 104 XYZ
#> 211 OLR 104 XY
#> 60 olt 104 XYZ
#> 72 onit 104 XYZ
#> 45 op1 104 XYZ
#> 46 op2 104 XYZ
#> 48 ora1 104 XYZ
#> 49 ora2 104 XYZ
#> 108 P 104 XYZ
#> 119 P_HYD 104 XYZ
#> 198 P_STRAT 104 0
#> 192 P_TOP 104 0
#> 194 P00 104 0
#> 47 paa 104 XYZ
#> 53 pan 104 XYZ
#> 109 PB 104 XYZ
#> 224 PBLH 104 XY
#> 248 PC 104 XY
#> 247 PCB 104 XY
#> 95 PH 104 XYZ
#> 96 PHB 104 XYZ
#> 31 PHOTR201 104 XYZ
#> 32 PHOTR202 104 XYZ
#> 33 PHOTR203 104 XYZ
#> 30 PHOTR204 104 XYZ
#> 28 PM10 104 XYZ
#> 27 PM2_5_DRY 104 XYZ
#> 123 PSFC 104 XY
#> 79 PV 104 XYZ
#> 120 Q2 104 XY
#> 136 QCLOUD 104 XYZ
#> 226 QFX 104 XY
#> 140 QGRAUP 104 XYZ
#> 138 QICE 104 XYZ
#> 164 QKE 104 XYZ
#> 141 QNDROP 104 XYZ
#> 137 QRAIN 104 XYZ
#> 139 QSNOW 104 XYZ
#> 135 QVAPOR 104 XYZ
#> 201 RAINC 104 XY
#> 203 RAINNC 104 XY
#> 7 RAINPROD 104 XYZ
#> 202 RAINSH 104 XY
#> 113 RDN 104 Z
#> 112 RDNW 104 Z
#> 126 RDX 104 0
#> 127 RDY 104 0
#> 128 RESM 104 0
#> 10 ROUGH_COR 104 XY
#> 35 SAC 104 XYZ
#> 232 SAVE_TOPO_FROM_REAL 106 0
#> 148 SEAICE 104 XY
#> 150 SFROFF 104 XY
#> 147 SH2O 104 XYZ
#> 142 SHDMAX 104 XY
#> 143 SHDMIN 104 XY
#> 188 SINALPHA 104 XY
#> 146 SMOIS 104 XYZ
#> 11 SMOIS_COR 104 XY
#> 144 SNOALB 104 XY
#> 158 SNOW 104 XY
#> 230 SNOWC 104 XY
#> 159 SNOWH 104 XY
#> 204 SNOWNC 104 XY
#> 34 SNU 104 XYZ
#> 36 so2 104 XYZ
#> 231 SR 104 XY
#> 251 SST 104 XY
#> 252 SST_INPUT 104 XY
#> 161 SSTSK 104 XY
#> 37 sulf 104 XYZ
#> 208 SWDOWN 104 XY
#> 210 SWNORM 104 XY
#> 97 T 104 XYZ
#> 193 T00 104 0
#> 121 T2 104 XY
#> 122 TH2 104 XY
#> 118 THIS_IS_AN_IDEAL_RUN 106 0
#> 98 THM 104 XYZ
#> 196 TISO 104 0
#> 174 TKE_PBL 104 XYZ
#> 195 TLP 104 0
#> 197 TLP_STRAT 104 0
#> 221 TMN 104 XY
#> 62 tol 104 XYZ
#> 65 tpan 104 XYZ
#> 191 TSK 104 XY
#> 101 TSK_FORCE 104 0
#> 104 TSK_FORCE_TEND 104 0
#> 145 TSLB 104 XYZ
#> 92 U 104 XYZ
#> 124 U10 104 XY
#> 151 UDROFF 104 XY
#> 223 UST 104 XY
#> 9 UST_T 104 XY
#> 93 V 104 XYZ
#> 125 V10 104 XY
#> 173 VAR 104 XY
#> 91 VAR_SSO 104 XY
#> 154 VEGFRA 104 XY
#> 94 W 104 XYZ
#> 149 XICEM 104 XY
#> 222 XLAND 104 XY
#> 84 XLAT 104 XY
#> 212 XLAT_U 104 XY
#> 214 XLAT_V 104 XY
#> 85 XLONG 104 XY
#> 213 XLONG_U 104 XY
#> 215 XLONG_V 104 XY
#> 134 XTIME 104 0
#> 63 xyl 104 XYZ
#> 83 YCLNO2 104 XYZ
#> 129 ZETATOP 104 0
#> 87 ZNU 104 Z
#> 88 ZNW 104 Z
#> 89 ZS 104 Z
#> description
#> 156 ACCUMULATED GROUND HEAT FLUX
#> 228 ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE
#> 229 ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE
#> 64 ACO3 mixing ratio
#> 157 ACCUMULATED MELTED SNOW
#> 218 BACKGROUND ALBEDO
#> 216 ALBEDO
#> 43 ALD mixing ratio
#> 1 Aerosol Optical Depth
#> 2 Aerosol Optical Depth, 2d
#> 3 Aerosol Optical Depth, top
#> 241 full levels, c1f = d bf / d eta, using znu
#> 239 half levels, c1h = d bf / d eta, using znw
#> 242 full levels, c2f = (1-c1f)*(p0-pt)
#> 240 half levels, c2h = (1-c1h)*(p0-pt)
#> 245 full levels, c3f = bf
#> 243 half levels, c3h = bh
#> 246 full levels, c4f = (eta-bf)*(p0-pt)+pt, using znw
#> 244 half levels, c4h = (eta-bh)*(p0-pt)+pt, using znu
#> 160 CANOPY WATER
#> 130 2nd order extrapolation constant
#> 131 2nd order extrapolation constant
#> 132 2nd order extrapolation constant
#> 116 extrapolation constant
#> 117 extrapolation constant
#> 76 CH4 mixing ratio
#> 217 COMPUTATIONAL GRID LATITUDE, SOUTH IS NEGATIVE
#> 207 CLOUD FRACTION
#> 6 CLOUD FRACTION
#> 81 n2o5 velocity
#> 58 CO mixing ratio
#> 75 CO2 mixing ratio
#> 189 Local cosine of map rotation
#> 162 COS of SOLAR ZENITH ANGLE
#> 73 CSL mixing ratio
#> 71 DCB mixing ratio
#> 29 dms oceanic concentrations
#> 115 d(eta) values between half (mass) levels
#> 114 d(eta) values between full (w) levels
#> 22 dry deposition velocity
#> 12 dust dry deposition velocity
#> 26 O3 deposition velocity
#> 90 THICKNESSES OF SOIL LAYERS
#> 187 Coriolis cosine latitude term
#> 24 Actual biog emiss
#> 23 Actual biog emiss
#> 13 biomass burning emiss
#> 165 EDMF relative updraft area - moist updrafts
#> 169 EDMF entrainment - mean moist updrafts
#> 170 EDMF qc - mean moist updrafts
#> 168 EDMF qt - mean moist updrafts
#> 167 EDMF thetaL - mean moist updrafts
#> 166 EDMF vertical velocity - mean moist updrafts
#> 175 Length scale from PBL
#> 219 SURFACE EMISSIVITY
#> 57 ETH mixing ratio
#> 8 RAIN EVAPORATION RATE
#> 186 Coriolis sine latitude term
#> 19 mean firesize for extratropical forest
#> 21 mean firesize for grassland
#> 20 mean firesize for savanna
#> 18 mean firesize for tropical forest
#> 110 upper weight for vertical stretching
#> 111 lower weight for vertical stretching
#> 80 n2o5 uptake by aerosol
#> 209 DOWNWARD LONG WAVE FLUX AT GROUND SURFACE
#> 69 GLY mixing ratio
#> 205 ACCUMULATED TOTAL GRID SCALE GRAUPEL
#> 155 GROUND HEAT FLUX
#> 42 H2O2 mixing ratio
#> 206 ACCUMULATED TOTAL GRID SCALE HAIL
#> 54 HC3 mixing ratio
#> 55 HC5 mixing ratio
#> 56 HC8 mixing ratio
#> 44 HCHO mixing ratio
#> 225 UPWARD HEAT FLUX AT THE SURFACE
#> 99 SCM ideal surface sensible heat flux
#> 102 SCM ideal surface sensible heat flux tendency
#> 190 Terrain Height
#> 41 HNO3 mixing ratio
#> 67 HNO4 mixing ratio
#> 77 HO mixing ratio
#> 78 HO2 mixing ratio
#> 66 HONO mixing ratio
#> 4 Some chem thing, ICN_DIAG
#> 235 Array to hold seed for restart, RAND_PERT
#> 234 Array to hold seed for restart, SKEBS
#> 233 Array to hold seed for restart, SPPT
#> 236 Array to hold seed for restart, RAND_PERT2
#> 238 Array to hold seed for restart, RAND_PERT4
#> 237 Array to hold seed for restart, RAND_PERT3
#> 153 DOMINANT SOIL CATEGORY
#> 74 ISO mixing ratio
#> 133
#> 152 DOMINANT VEGETATION CATEGORY
#> 68 KET mixing ratio
#> 82 n2o5 het reaction rate
#> 163 LEAF AREA INDEX
#> 25 MODIS LAI vegetation mask for this date; 0=no dust produced (vegetation)
#> 250 LAKE MASK (1 FOR LAND, 0 FOR WATER)
#> 249 LAND MASK (1 FOR LAND, 0 FOR WATER)
#> 227 LATENT HEAT FLUX AT THE SURFACE
#> 100 SCM ideal surface latent heat flux
#> 103 SCM ideal surface latent heat flux tendency
#> 86 LAND USE CATEGORY
#> 176 Map scale factor on mass grid
#> 179 Map scale factor on mass grid, x direction
#> 180 Map scale factor on mass grid, y direction
#> 177 Map scale factor on u-grid
#> 181 Map scale factor on u-grid, x direction
#> 182 Map scale factor on u-grid, y direction
#> 178 Map scale factor on v-grid
#> 183 Map scale factor on v-grid, x direction
#> 185 Map scale factor on v-grid, y direction
#> 199 Max map factor in domain
#> 200 Max map factor in domain
#> 172 Maximum mass-flux (neg: all dry, pos: moist)
#> 15 mean fraction of extra tropical forest
#> 17 mean fraction of grassland
#> 16 mean fraction of savanna
#> 14 mean fraction of tropical forest
#> 184 Inverse map scale factor on v-grid, x direction
#> 70 MGLY mixing ratio
#> 105 perturbation dry air mass in column
#> 106 base state dry air mass in column
#> 51 N2O5 mixing ratio
#> 5 Some chem thing, NC_DIAG
#> 107 -
#> 50 NH3 mixing ratio
#> 39 NO mixing ratio
#> 38 NO2 mixing ratio
#> 52 NO3 mixing ratio
#> 220 RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET
#> 171 Number of updrafts per grid cell
#> 40 O3 mixing ratio
#> 59 OL2 mixing ratio
#> 61 OLI mixing ratio
#> 211 TOA OUTGOING LONG WAVE
#> 60 OLT mixing ratio
#> 72 ONIT mixing ratio
#> 45 OP1 mixing ratio
#> 46 OP2 mixing ratio
#> 48 ORA1 mixing ratio
#> 49 ORA2 mixing ratio
#> 108 perturbation pressure
#> 119 hydrostatic pressure
#> 198 BASE STATE PRESSURE AT BOTTOM OF STRATOSPHERE
#> 192 PRESSURE TOP OF THE MODEL
#> 194 BASE STATE PRESURE
#> 47 PAA mixing ratio
#> 53 PAN mixing ratio
#> 109 BASE STATE PRESSURE
#> 224 PBL HEIGHT
#> 248 perturbation dry air mass in column
#> 247 base state dry air mass in column
#> 95 perturbation geopotential
#> 96 base-state geopotential
#> 31 cl2 photolysis rate
#> 32 hocl photolysis rate
#> 33 fmcl photolysis rate
#> 30 CLNO2 Photolysis Rate
#> 28 pm10 dry mass
#> 27 pm2.5 aerosol dry mass
#> 123 SFC PRESSURE
#> 79 Potential Vorticity
#> 120 QV at 2 M
#> 136 Cloud water mixing ratio
#> 226 UPWARD MOISTURE FLUX AT THE SURFACE
#> 140 Graupel mixing ratio
#> 138 Ice mixing ratio
#> 164 twice TKE from MYNN
#> 141 Droplet number mixing ratio
#> 137 Rain water mixing ratio
#> 139 Snow mixing ratio
#> 135 Water vapor mixing ratio
#> 201 ACCUMULATED TOTAL CUMULUS PRECIPITATION
#> 203 ACCUMULATED TOTAL GRID SCALE PRECIPITATION
#> 7 TOTAL RAIN PRODUCTION RATE
#> 202 ACCUMULATED SHALLOW CUMULUS PRECIPITATION
#> 113 inverse d(eta) values between half (mass) levels
#> 112 inverse d(eta) values between full (w) levels
#> 126 INVERSE X GRID LENGTH
#> 127 INVERSE Y GRID LENGTH
#> 128 TIME WEIGHT CONSTANT FOR SMALL STEPS
#> 10 roughness elements correction
#> 35 2nd moment Aitken mode
#> 232 1=original topo from real/0=topo modified by WRF
#> 148 SEA ICE FLAG
#> 150 SURFACE RUNOFF
#> 147 SOIL LIQUID WATER
#> 142 ANNUAL MAX VEG FRACTION
#> 143 ANNUAL MIN VEG FRACTION
#> 188 Local sine of map rotation
#> 146 SOIL MOISTURE
#> 11 soil moisture correction
#> 144 ANNUAL MAX SNOW ALBEDO IN FRACTION
#> 158 SNOW WATER EQUIVALENT
#> 230 FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)
#> 159 PHYSICAL SNOW DEPTH
#> 204 ACCUMULATED TOTAL GRID SCALE SNOW AND ICE
#> 34 2nd moment Aitken mode
#> 36 SO2 mixing ratio
#> 231 fraction of frozen precipitation
#> 251 SEA SURFACE TEMPERATURE
#> 252 SEA SURFACE TEMPERATURE FROM WRFLOWINPUT FILE
#> 161 SKIN SEA SURFACE TEMPERATURE
#> 37 SULF mixing ratio
#> 208 DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE
#> 210 NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT)
#> 97 perturbation potential temperature theta-t0
#> 193 BASE STATE TEMPERATURE
#> 121 TEMP at 2 M
#> 122 POT TEMP at 2 M
#> 118 T/F flag: this is an ARW ideal simulation
#> 98 either 1) pert moist pot temp=(1+Rv/Rd Qv)*(theta)-T0, or 2) pert dry pot temp=t
#> 196 TEMP AT WHICH THE BASE T TURNS CONST
#> 174 TKE from PBL
#> 195 BASE STATE LAPSE RATE
#> 197 BASE STATE LAPSE RATE (DT/D(LN(P)) IN STRATOSPHERE
#> 221 SOIL TEMPERATURE AT LOWER BOUNDARY
#> 62 TOL mixing ratio
#> 65 TPAN mixing ratio
#> 191 SURFACE SKIN TEMPERATURE
#> 101 SCM ideal surface skin temperature
#> 104 SCM ideal surface skin temperature tendency
#> 145 SOIL TEMPERATURE
#> 92 x-wind component
#> 124 U at 10 M
#> 151 UNDERGROUND RUNOFF
#> 223 U* IN SIMILARITY THEORY
#> 9 Threshold Friction Velocity
#> 93 y-wind component
#> 125 V at 10 M
#> 173 OROGRAPHIC VARIANCE
#> 91 variance of subgrid-scale orography
#> 154 VEGETATION FRACTION
#> 94 z-wind component
#> 149 SEA ICE FLAG (PREVIOUS STEP)
#> 222 LAND MASK (1 FOR LAND, 2 FOR WATER)
#> 84 LATITUDE, SOUTH IS NEGATIVE
#> 212 LATITUDE, SOUTH IS NEGATIVE
#> 214 LATITUDE, SOUTH IS NEGATIVE
#> 85 LONGITUDE, WEST IS NEGATIVE
#> 213 LONGITUDE, WEST IS NEGATIVE
#> 215 LONGITUDE, WEST IS NEGATIVE
#> 134 minutes since 2014-10-03 00:00:00
#> 63 XYL mixing ratio
#> 83 clno2 yield from n2o5 het
#> 129 ZETA AT MODEL TOP
#> 87 eta values on half (mass) levels
#> 88 eta values on full (w) levels
#> 89 DEPTHS OF CENTERS OF SOIL LAYERS
#> units stagger coordinates
#> 156 J m-2 XLONG XLAT XTIME
#> 228 J m-2 XLONG XLAT XTIME
#> 229 J m-2 XLONG XLAT XTIME
#> 64 ppmv XLONG XLAT XTIME
#> 157 kg m-2 XLONG XLAT XTIME
#> 218 XLONG XLAT XTIME
#> 216 - XLONG XLAT XTIME
#> 43 ppmv XLONG XLAT XTIME
#> 1 XLONG XLAT XTIME
#> 2 XLONG XLAT XTIME
#> 3 XLONG XLAT XTIME
#> 241 Dimensionless Z 104
#> 239 Dimensionless 104
#> 242 Pa Z 104
#> 240 Pa 104
#> 245 Dimensionless Z 104
#> 243 Dimensionless 104
#> 246 Pa Z 104
#> 244 Pa 104
#> 160 kg m-2 XLONG XLAT XTIME
#> 130 104
#> 131 104
#> 132 104
#> 116 104
#> 117 104
#> 76 ppmv XLONG XLAT XTIME
#> 217 degree_north XLONG XLAT XTIME
#> 207 XLONG XLAT XTIME
#> 6 - XLONG XLAT XTIME
#> 81 m/s XLONG XLAT XTIME
#> 58 ppmv XLONG XLAT XTIME
#> 75 ppmv XLONG XLAT XTIME
#> 189 XLONG XLAT XTIME
#> 162 dimensionless XLONG XLAT XTIME
#> 73 ppmv XLONG XLAT XTIME
#> 71 ppmv XLONG XLAT XTIME
#> 29 nM/L XLONG XLAT XTIME
#> 115 104
#> 114 104
#> 22 cm/s Z XLONG XLAT XTIME
#> 12 m/s XLONG XLAT XTIME
#> 26 cm/s XLONG XLAT XTIME
#> 90 m Z 104
#> 187 s-1 XLONG XLAT XTIME
#> 24 mol km^-2 hr^-1 XLONG XLAT XTIME
#> 23 mol km^-2 hr^-1 XLONG XLAT XTIME
#> 13 ug/m2/s XLONG XLAT XTIME
#> 165 - XLONG XLAT XTIME
#> 169 m-1 XLONG XLAT XTIME
#> 170 kg kg-1 XLONG XLAT XTIME
#> 168 kg kg-1 XLONG XLAT XTIME
#> 167 K XLONG XLAT XTIME
#> 166 m s-1 XLONG XLAT XTIME
#> 175 m Z XLONG XLAT XTIME
#> 219 XLONG XLAT XTIME
#> 57 ppmv XLONG XLAT XTIME
#> 8 s-1 XLONG XLAT XTIME
#> 186 s-1 XLONG XLAT XTIME
#> 19 ? XLONG XLAT XTIME
#> 21 ? XLONG XLAT XTIME
#> 20 ? XLONG XLAT XTIME
#> 18 ? XLONG XLAT XTIME
#> 110 104
#> 111 104
#> 80 numerical value XLONG XLAT XTIME
#> 209 W m-2 XLONG XLAT XTIME
#> 69 ppmv XLONG XLAT XTIME
#> 205 mm XLONG XLAT XTIME
#> 155 W m-2 XLONG XLAT XTIME
#> 42 ppmv XLONG XLAT XTIME
#> 206 mm XLONG XLAT XTIME
#> 54 ppmv XLONG XLAT XTIME
#> 55 ppmv XLONG XLAT XTIME
#> 56 ppmv XLONG XLAT XTIME
#> 44 ppmv XLONG XLAT XTIME
#> 225 W m-2 XLONG XLAT XTIME
#> 99 W m-2 104
#> 102 W m-2 s-1 104
#> 190 m XLONG XLAT XTIME
#> 41 ppmv XLONG XLAT XTIME
#> 67 ppmv XLONG XLAT XTIME
#> 77 ppmv XLONG XLAT XTIME
#> 78 ppmv XLONG XLAT XTIME
#> 66 ppmv XLONG XLAT XTIME
#> 4 XLONG XLAT XTIME
#> 235 Z 106
#> 234 Z 106
#> 233 Z 106
#> 236 Z 106
#> 238 Z 106
#> 237 Z 106
#> 153 XLONG XLAT XTIME
#> 74 ppmv XLONG XLAT XTIME
#> 133 106
#> 152 XLONG XLAT XTIME
#> 68 ppmv XLONG XLAT XTIME
#> 82 s-1 XLONG XLAT XTIME
#> 163 m-2/m-2 XLONG XLAT XTIME
#> 25 none XLONG XLAT XTIME
#> 250 XLONG XLAT XTIME
#> 249 XLONG XLAT XTIME
#> 227 W m-2 XLONG XLAT XTIME
#> 100 W m-2 104
#> 103 W m-2 s-1 104
#> 86 XLONG XLAT XTIME
#> 176 XLONG XLAT XTIME
#> 179 XLONG XLAT XTIME
#> 180 XLONG XLAT XTIME
#> 177 X XLONG_U XLAT_U XTIME
#> 181 X XLONG_U XLAT_U XTIME
#> 182 X XLONG_U XLAT_U XTIME
#> 178 Y XLONG_V XLAT_V XTIME
#> 183 Y XLONG_V XLAT_V XTIME
#> 185 Y XLONG_V XLAT_V XTIME
#> 199 104
#> 200 104
#> 172 m/s * area XLONG XLAT XTIME
#> 15 ? XLONG XLAT XTIME
#> 17 ? XLONG XLAT XTIME
#> 16 ? XLONG XLAT XTIME
#> 14 ? XLONG XLAT XTIME
#> 184 Y XLONG_V XLAT_V XTIME
#> 70 ppmv XLONG XLAT XTIME
#> 105 Pa XLONG XLAT XTIME
#> 106 Pa XLONG XLAT XTIME
#> 51 ppmv XLONG XLAT XTIME
#> 5 XLONG XLAT XTIME
#> 107 - XLONG XLAT XTIME
#> 50 ppmv XLONG XLAT XTIME
#> 39 ppmv XLONG XLAT XTIME
#> 38 ppmv XLONG XLAT XTIME
#> 52 ppmv XLONG XLAT XTIME
#> 220 W m{-2} XLONG XLAT XTIME
#> 171 XLONG XLAT XTIME
#> 40 ppmv XLONG XLAT XTIME
#> 59 ppmv XLONG XLAT XTIME
#> 61 ppmv XLONG XLAT XTIME
#> 211 W m-2 XLONG XLAT XTIME
#> 60 ppmv XLONG XLAT XTIME
#> 72 ppmv XLONG XLAT XTIME
#> 45 ppmv XLONG XLAT XTIME
#> 46 ppmv XLONG XLAT XTIME
#> 48 ppmv XLONG XLAT XTIME
#> 49 ppmv XLONG XLAT XTIME
#> 108 Pa XLONG XLAT XTIME
#> 119 Pa XLONG XLAT XTIME
#> 198 Pa 104
#> 192 Pa 104
#> 194 Pa 104
#> 47 ppmv XLONG XLAT XTIME
#> 53 ppmv XLONG XLAT XTIME
#> 109 Pa XLONG XLAT XTIME
#> 224 m XLONG XLAT XTIME
#> 248 Pa XLONG XLAT XTIME
#> 247 Pa XLONG XLAT XTIME
#> 95 m2 s-2 Z XLONG XLAT XTIME
#> 96 m2 s-2 Z XLONG XLAT XTIME
#> 31 min{-1} XLONG XLAT XTIME
#> 32 min{-1} XLONG XLAT XTIME
#> 33 min{-1} XLONG XLAT XTIME
#> 30 min{-1} XLONG XLAT XTIME
#> 28 ug m^-3 XLONG XLAT XTIME
#> 27 ug m^-3 XLONG XLAT XTIME
#> 123 Pa XLONG XLAT XTIME
#> 79 pvu XLONG XLAT XTIME
#> 120 kg kg-1 XLONG XLAT XTIME
#> 136 kg kg-1 XLONG XLAT XTIME
#> 226 kg m-2 s-1 XLONG XLAT XTIME
#> 140 kg kg-1 XLONG XLAT XTIME
#> 138 kg kg-1 XLONG XLAT XTIME
#> 164 m2 s-2 XLONG XLAT XTIME
#> 141 kg-1 XLONG XLAT XTIME
#> 137 kg kg-1 XLONG XLAT XTIME
#> 139 kg kg-1 XLONG XLAT XTIME
#> 135 kg kg-1 XLONG XLAT XTIME
#> 201 mm XLONG XLAT XTIME
#> 203 mm XLONG XLAT XTIME
#> 7 s-1 XLONG XLAT XTIME
#> 202 mm XLONG XLAT XTIME
#> 113 104
#> 112 104
#> 126 104
#> 127 104
#> 128 104
#> 10 XLONG XLAT XTIME
#> 35 m2 m-3 XLONG XLAT XTIME
#> 232 flag 106
#> 148 XLONG XLAT XTIME
#> 150 mm XLONG XLAT XTIME
#> 147 m3 m-3 Z XLONG XLAT XTIME
#> 142 XLONG XLAT XTIME
#> 143 XLONG XLAT XTIME
#> 188 XLONG XLAT XTIME
#> 146 m3 m-3 Z XLONG XLAT XTIME
#> 11 XLONG XLAT XTIME
#> 144 XLONG XLAT XTIME
#> 158 kg m-2 XLONG XLAT XTIME
#> 230 XLONG XLAT XTIME
#> 159 m XLONG XLAT XTIME
#> 204 mm XLONG XLAT XTIME
#> 34 m2 m-3 XLONG XLAT XTIME
#> 36 ppmv XLONG XLAT XTIME
#> 231 - XLONG XLAT XTIME
#> 251 K XLONG XLAT XTIME
#> 252 K XLONG XLAT XTIME
#> 161 K XLONG XLAT XTIME
#> 37 ppmv XLONG XLAT XTIME
#> 208 W m-2 XLONG XLAT XTIME
#> 210 W m-2 XLONG XLAT XTIME
#> 97 K XLONG XLAT XTIME
#> 193 K 104
#> 121 K XLONG XLAT XTIME
#> 122 K XLONG XLAT XTIME
#> 118 - 106
#> 98 K XLONG XLAT XTIME
#> 196 K 104
#> 174 m2 s-2 Z XLONG XLAT XTIME
#> 195 104
#> 197 K 104
#> 221 K XLONG XLAT XTIME
#> 62 ppmv XLONG XLAT XTIME
#> 65 ppmv XLONG XLAT XTIME
#> 191 K XLONG XLAT XTIME
#> 101 W m-2 104
#> 104 W m-2 s-1 104
#> 145 K Z XLONG XLAT XTIME
#> 92 m s-1 X XLONG_U XLAT_U XTIME
#> 124 m s-1 XLONG XLAT XTIME
#> 151 mm XLONG XLAT XTIME
#> 223 m s-1 XLONG XLAT XTIME
#> 9 m s-1 XLONG XLAT XTIME
#> 93 m s-1 Y XLONG_V XLAT_V XTIME
#> 125 m s-1 XLONG XLAT XTIME
#> 173 XLONG XLAT XTIME
#> 91 m2 XLONG XLAT XTIME
#> 154 XLONG XLAT XTIME
#> 94 m s-1 Z XLONG XLAT XTIME
#> 149 XLONG XLAT XTIME
#> 222 XLONG XLAT XTIME
#> 84 degree_north XLONG XLAT
#> 212 degree_north X XLONG_U XLAT_U
#> 214 degree_north Y XLONG_V XLAT_V
#> 85 degree_east XLONG XLAT
#> 213 degree_east X XLONG_U XLAT_U
#> 215 degree_east Y XLONG_V XLAT_V
#> 134 minutes since 2014-10-03 00:00:00 104
#> 63 ppmv XLONG XLAT XTIME
#> 83 numerical value XLONG XLAT XTIME
#> 129 104
#> 87 104
#> 88 Z 104
#> 89 m Z 104Now we can select some variables
(vars <- dt[dt$MemoryOrder == "XYZ", ]$vars)
#> [1] "aco3" "ald" "AOD_OUT" "ch4" "CLDFRA"
#> [6] "CLDFRA2" "CN2O5" "co" "co2" "csl"
#> [11] "dcb" "DRY_DEP_LEN" "dvel_o3" "ebu_oc" "EDMF_A"
#> [16] "EDMF_ENT" "EDMF_QC" "EDMF_QT" "EDMF_THL" "EDMF_W"
#> [21] "EL_PBL" "eth" "EVAPPROD" "GAMN2O5" "gly"
#> [26] "h2o2" "hc3" "hc5" "hc8" "hcho"
#> [31] "hno3" "hno4" "ho" "ho2" "hono"
#> [36] "ICN_DIAG" "iso" "ket" "KN2O5" "mgly"
#> [41] "n2o5" "NC_DIAG" "nh3" "no" "no2"
#> [46] "no3" "o3" "ol2" "oli" "olt"
#> [51] "onit" "op1" "op2" "ora1" "ora2"
#> [56] "P" "P_HYD" "paa" "pan" "PB"
#> [61] "PH" "PHB" "PHOTR201" "PHOTR202" "PHOTR203"
#> [66] "PHOTR204" "PM10" "PM2_5_DRY" "PV" "QCLOUD"
#> [71] "QGRAUP" "QICE" "QKE" "QNDROP" "QRAIN"
#> [76] "QSNOW" "QVAPOR" "RAINPROD" "SAC" "SH2O"
#> [81] "SMOIS" "SNU" "so2" "sulf" "T"
#> [86] "THM" "TKE_PBL" "tol" "tpan" "TSLB"
#> [91] "U" "V" "W" "xyl" "YCLNO2"Now get some statistics