Skip to contents

Read output from WRF model and calculate the column of trace gases.

Usage

calculate_column(
  file = file.choose(),
  name,
  met,
  DU,
  flip_v = FALSE,
  flip_h = FALSE,
  verbose = TRUE,
  ...
)

Arguments

file

WRF output file, see notes

name

trace gas name to be integrated

met

(optional) WRF output file for meteorological variables, see notes

DU

true to change the output units from 'molecules cm-1' to 'DU'

flip_v

passed to wrf_rast, see notes

flip_h

passed to wrf_rast, see notes

verbose

display additional information

...

extra arguments passed to eva3dm::wrf_rast or eva3dm::wrf_sds

Value

SpatRaster object (from terra package).

Note

files in file should contain Times, XLAT, XLONG variables in addition to the concentration to be integrated, the variable should include at least 3 dimensions including vertical.

met is a optional file, it should have containing PHB, PH, PB, P, and T variables with the same dimension of the concentration integrated.

post processing can affect the orientation of the variables, the arguments flip_v and flip_h and other arguments from eva3dm::wrf_rast and eva3dm::wrf_sds can be used to take effect into account.