IVT calculation app for HighResMIP

[taobrienlbl]

We need an app for calculating integrated vapor transport from the HighResMIP models.

It will need to accept files containing both 4D (time, lev, lat ,lon) and 3D (time, lat, lon) variables:

• 4D variables
  • humidity Q
  • zonal wind U
  • meridional wind V
• 3D variables
  • surface humidity QSFC
  • surface zonal wind USFC
  • surface meridional wind VSFC
  • surface pressure PS, or
    • surface sea-level pressure
    • surface elevation (a 2D [lat, lon] field)

It will output netCDF files containing IVT, a 3D field that contains the vertical integral of horizontal moisture transport.

A typical run command should be something like:

teca_ivt_integral \
    --q_file_regex QFILEREGEX \
    --u_file_regex UFILEREGEX \
    --v_file_regex VFILEREGEX \
    --qsfc_file_regex QSFCFILEREGEX \
    --usfc_file_regex USFCFILEREGEX \
    --vsfc_file_regex VSFCFILEREGEX \
    --ps_file_regex PSFILEREGEX
    --q_var_name QVARNAME \
    --u_var_name UVARNAME \
    --v_var_name VVARNAME \
    --ps_var_name PSVARNAME \
    --output_ivt_var_name IVTVARNAME \
    --output_file_name OUTPUTFILETEMPLATE \

The pipeline will be something like:

CF Reader --> Normalize Coordinates --> Moisture Transport vector (Q*U, V*Q) --> Vertical Integral (uIVT, vIVT) --> l2_norm (IVT) --> CF writer

Depends on:

• support cf_reader reading multiple variables stored across multiple files #353
• Vertical integral algorithm #390

[taobrienlbl]

Sketch of pipeline from 10/6/20

This supersedes the pipeline described in the original issue text.

[taobrienlbl]

An update to the current task; we are only focusing on integration of 4D atmospheric data (ignoring surface 3D variables for now).

It will need to accept files containing both 4D (time, lev, lat ,lon):

4D variables
    humidity Q
    zonal wind U
    meridional wind V

It will output netCDF files containing IVT, a 3D field that contains the vertical integral of horizontal moisture transport.

A typical run command should be something like:

teca_ivt_integral \
    --q_file_regex QFILEREGEX \
    --u_file_regex UFILEREGEX \
    --v_file_regex VFILEREGEX \
    --q_var_name QVARNAME \
    --u_var_name UVARNAME \
    --v_var_name VVARNAME \
    --output_ivt_var_name IVTVARNAME \
    --output_file_name OUTPUTFILETEMPLATE