Registry.xml

<?xml version="1.0"?>
<registry model="mpas" core="init_atmosphere" core_abbrev="init_atm" version="8.3.1">

<!-- **************************************************************************************** -->
<!-- ************************************** Dimensions ************************************** -->
<!-- **************************************************************************************** -->

        <dims>
                <dim name="nCells"
                     description="The number of Voronoi cells in the primal mesh"/>
                <dim name="nEdges"
                     description="The number of cell faces (edges) in the primal and dual meshes"/>
                <dim name="maxEdges"
                     description="The largest number of neighbors that a primal mesh cell may have"/>
                <dim name="maxEdges2"
                     description="Twice the value of maxEdges"/>
                <dim name="nVertices"
                     description="The number of Delaunay triangles in the dual mesh, equivalent to the number of cell corners in the primal mesh"/>
                <dim name="TWO" definition="2"
                     description="A constant value of 2"/>
                <dim name="THREE" definition="3"
                     description="A constant value of 3"/>
                <dim name="vertexDegree"
                     description="The number of edges or cells that meet at vertices"/>
                <dim name="FIFTEEN" definition="15"
                     description="A constant value of 15"/>
                <dim name="TWENTYONE" definition="21"
                     description="A constant value of 21"/>
                <dim name="R3" definition="3"
                     description="The dimensionality of vectors in $\Re^3$"/>
                <dim name="nVertLevels" definition="namelist:config_nvertlevels"
                     description="The number of atmospheric layers"/>
                <dim name="nSoilLevels" definition="namelist:config_nsoillevels"
                     description="The number of soil layers used by the land-surface scheme"/>
                <dim name="nFGLevels" definition="namelist:config_nfglevels"
                     description="The number of first-guess atmospheric layers"/>
                <dim name="nFGSoilLevels" definition="namelist:config_nfgsoillevels"
                     description="The number of first-guess soil layers"/>
                <dim name="nVertLevelsP1" definition="nVertLevels+1"
                     description="The number of atmospheric levels, always one more than the number of layers"/>
                <dim name="nSoilComps" definition="8"
                     description="The number of soil textures as needed as input in the NOAH-MP land surface model"/>
                <dim name="nGocartLevels" definition="namelist:config_gocartlevels"
                     description="The number of prescribed pressure levels in the climatological GOCART data set"/>
                <dim name="nMonths" definition="namelist:config_months"
                     description="The number of months in a year"/>
        </dims>


<!-- **************************************************************************************** -->
<!-- ************************************** Namelists *************************************** -->
<!-- **************************************************************************************** -->

        <nml_record name="nhyd_model" in_defaults="true">

                <nml_option name="config_init_case"             type="integer"       default_value="7"
                     units="-"
                     description="Type of initial conditions to create: \newline
                                        1 = Jablonowski &amp; Williamson barolinic wave (no initial perturbation), \newline
                                        2 = Jablonowski &amp; Williamson barolinic wave (with initial perturbation), \newline
                                        3 = Jablonowski &amp; Williamson barolinic wave (with normal-mode perturbation), \newline
                                        4 = squall line, \newline
                                        5 = super-cell, \newline
                                        6 = mountain wave, \newline
                                        7 = real-data initial conditions from, e.g., GFS, \newline
                                        8 = surface field (SST, sea-ice) update file for use with real-data simulations \newline
                                        9 = lateral boundary conditions update file for use with real-data simulations \newline
                                        13 = CAM-MPAS 3-d grid with specified topography and zeta levels"
                     possible_values="1 -- 9, or 13"/>

                <nml_option name="config_calendar_type" type="character" default_value="gregorian" in_defaults="false"
                     units="-"
                     description="Simulation calendar type"
                     possible_values="`gregorian',`gregorian_noleap'"/>

                <nml_option name="config_start_time"            type="character"     default_value="2010-10-23_00:00:00"
                     units="-"
                     description="Time to begin processing first-guess data (cases 7, 8, and 9 only)"
                     possible_values="`YYYY-MM-DD_hh:mm:ss'"/>

                <nml_option name="config_stop_time"             type="character"     default_value="2010-10-23_00:00:00"
                     units="-"
                     description="Time to end processing first-guess data (cases 8 and 9 only)"
                     possible_values="`YYYY-MM-DD_hh:mm:ss'"/>

                <nml_option name="config_theta_adv_order" type="integer" default_value="3"
                     units="-"
                     description="Horizontal advection order for theta"
                     possible_values="2, 3, or 4"/>

                <nml_option name="config_coef_3rd_order" type="real" default_value="0.25"
                     units="-"
                     description="Upwinding coefficient in the 3rd order advection scheme"
                     possible_values="0 $\leq$ config_coef_3rd_order $\leq$ 1"/>

                <nml_option name="config_num_halos" type="integer" default_value="2" in_defaults="false"
                     units="-"
                     description="Number of halo layers for fields"
                     possible_values="Integer values, typically 2 or 3; DO NOT CHANGE"/>

                <nml_option name="config_interface_projection" type="character"     default_value="linear_interpolation"
                     units="-"
                     description="projecting layer values to the interface, linear vertical interpolation or integral" 
                     possible_values="'linear_interpolation' or 'layer_integral' (layer average value)"/>

        </nml_record>

        <nml_record name="dimensions" in_defaults="true">

                <nml_option name="config_nvertlevels"           type="integer"       default_value="55"
                     units="-"
                     description="The number of vertical levels to be used in the model"
                     possible_values="Positive integer values"/>

                <nml_option name="config_nsoillevels"           type="integer"       default_value="4"
                     units="-"
                     description="The number of vertical soil levels needed by LSM in the model (case 7 only)"
                     possible_values="Positive integer values"/>

                <nml_option name="config_nfglevels"             type="integer"       default_value="38"
                     units="-"
                     description="The number of atmospheric levels (including surface and sea-level) in the first-guess dataset (cases 7 and 9 only)"
                     possible_values="Positive integer values"/>

                <nml_option name="config_nfgsoillevels"         type="integer"       default_value="4"
                     units="-"
                     description="The number of vertical soil levels in the first-guess dataset (case 7 only)"
                     possible_values="Positive integer values"/>

                <nml_option name="config_gocartlevels"          type="integer"       default_value="30"
                     units="-"
                     description="The number of vertical gocart levels in the climatological gocart data set"
                     possible_values="Positive integer values"/>

                <nml_option name="config_months"                type="integer"       default_value="12"      in_defaults="false"
                     units="-"
                     description="The number of months in a year"
                     possible_values="Positive integer values"/>

        </nml_record>

        <nml_record name="data_sources" in_defaults="true">

                <nml_option name="config_geog_data_path"        type="character"     default_value="/glade/work/wrfhelp/WPS_GEOG/"
                     units="-"
                     description="Path to the WPS static data files (case 7 only)"
                     possible_values="Any valid path"/>

                <nml_option name="config_met_prefix"            type="character"     default_value="CFSR"
                     units="-"
                     description="Filename prefix of ungrib intermediate file to use for initial conditions (cases 7 and 9 only)"
                     possible_values="Any alpha-numeric string"/>

                <nml_option name="config_sfc_prefix"            type="character"     default_value="SST"
                     units="-"
                     description="Filename prefix of ungrib intermediate file to use for SST and sea-ice (cases 7 and 8 only)"
                     possible_values="Any alpha-numeric string"/>

                <nml_option name="config_fg_interval"           type="integer"       default_value="86400"
                     units="-"
                     description="Interval between intermediate files (cases 8 and 9 only)"
                     possible_values="[DDD_]hh:mm:ss"/>

                <nml_option name="config_landuse_data"          type="character"     default_value="MODIFIED_IGBP_MODIS_NOAH"
                     units="-"
                     description="The land use classification to use (case 7 only)"
                     possible_values="`USGS', `MODIFIED_IGBP_MODIS_NOAH', or `MODIFIED_IGBP_MODIS_NOAH_15s'"/>

                <nml_option name="config_soilcat_data"          type="character"     default_value="STATSGO"
                     units="-"
                     description="The soil category classification to use"
                     possible_values="`STATSGO' or `BNU'"/>

                <nml_option name="config_topo_data"             type="character"     default_value="GMTED2010"
                     units="-"
                     description="The topography dataset to use for the model terrain and for GWDO static fields (case 7 only)"
                     possible_values="`GTOPO30' or `GMTED2010'"/>

                <nml_option name="config_vegfrac_data"          type="character"     default_value="MODIS"
                     units="-"
                     description="The climatological monthly vegetation fraction dataset to use (case 7 only)"
                     possible_values="`MODIS' or `NCEP'"/>

                <nml_option name="config_albedo_data"           type="character"     default_value="MODIS"
                     units="-"
                     description="The climatological monthly albedo dataset to use (case 7 only)"
                     possible_values="`MODIS' or `NCEP'"/>

                <nml_option name="config_maxsnowalbedo_data"    type="character"     default_value="MODIS"
                     units="-"
                     description="The maximum snow albedo dataset to use (case 7 only)"
                     possible_values="`MODIS' or `NCEP'"/>

                <nml_option name="config_supersample_factor"    type="integer"       default_value="3"
                     units="-"
                     description="The supersampling factor to be used for MODIS maximum snow albedo and monthly albedo datasets (case 7 only)"
                     possible_values="Positive integer values"/>

	        <nml_option name="config_lu_supersample_factor" type="integer"       default_value="1"
                     units="-"
                     description="The supersampling factor to be used for 30s or 15s MODIS land use, or for 30s USGS land use, as selected by `config_landuse_data' (case 7 only)"
                     possible_values="Positive integer values"/>

                <nml_option name="config_30s_supersample_factor"    type="integer"       default_value="1"
                     units="-"
                     description="The supersampling factor to be used for 30s terrain, soil category, and MODIS FPAR monthly vegetation fraction (case 7 only)"
                     possible_values="Positive integer values"/>

                <nml_option name="config_noahmp_static" type="logical" default_value="true" in_defaults="false"
                     units="-"
                     description="Whether to process, read, and write static fields only required by Noah-MP"
                     possible_values="true or false"/>

                <nml_option name="config_use_spechumd"          type="logical"       default_value="false"
                     units="-"
                     description="Whether to use specific-humidity as the first-guess moisture variable. If this option is False, relative humidity will be used."
                     possible_values="true or false"/>

        </nml_record>

        <nml_record name="vertical_grid" in_defaults="true">

                <nml_option name="config_ztop"                  type="real"          default_value="30000.0"
                     units="m"
                     description="Model top height"
                     possible_values="Positive real values"/>

                <nml_option name="config_hybrid_coordinate"       type="logical"       default_value="true"
                     units="-"
                     description="Whether to employ a hybrid vertical coordinate"
                     possible_values="true or false"/>

                <nml_option name="config_hybrid_top_z"       type="real"       default_value="30000.0"
                     units="-"
                     description="Height at which coordinate surfaces become constant height surfaces"
                     possible_values="Positive real values"/>

                <nml_option name="config_nsmterrain"            type="integer"       default_value="1"
                     units="-"
                     description="Number of smoothing passes to apply to the interpolated terrain field"
                     possible_values="Non-negative integer values"/>

                <nml_option name="config_smooth_surfaces"       type="logical"       default_value="true"
                     units="-"
                     description="Whether to smooth zeta surfaces"
                     possible_values="true or false"/>

                <nml_option name="config_dzmin"                 type="real"          default_value="0.3"
                     units="-"
                     description="Minimum thickness of layers as a fraction of nominal thickness"
                     possible_values="Real values in the interval (0,1)"/>

                <nml_option name="config_nsm"                   type="integer"       default_value="30"
                     units="-"
                     description="Maximum number of smoothing passes for coordinate surfaces"
                     possible_values="Positive integer values"/>

                <nml_option name="config_tc_vertical_grid"      type="logical"       default_value="true"
                     units="-"
                     description="Whether to use the vertical layer profile that was developed for use in real-time TC experiments"
                     possible_values="true or false"/>

                <nml_option name="config_specified_zeta_levels" type="character"     default_value=""    in_defaults="false"
                     units="-"
                     description="The name of a text file containing a list of vertical coordinate (zeta) values in increasing order at layer interfaces."
                     possible_values="Any valid filename"/>

                <nml_option name="config_blend_bdy_terrain"     type="logical"       default_value="false"
                     units="-"
                     description="Whether to blend terrain along domain boundaries with first-guess terrain. Only useful for limited-area domains."
                     possible_values="true or false"/>

        </nml_record>

        <nml_record name="interpolation_control" in_defaults="true">

                <nml_option name="config_extrap_airtemp"        type="character"     default_value="lapse-rate"
                     units="-"
                     description="Method of extrapolation of air temperature above/below first-guess levels."
                     possible_values="`constant' (last valid value), `linear' (linear extrapolation based on last two values), `lapse-rate' (0.0065 K/m from last valid value)"/>

        </nml_record>

        <nml_record name="preproc_stages" in_defaults="true">

                <nml_option name="config_static_interp"         type="logical"       default_value="true"
                     units="-"
                     description="Whether to interpolate WPS static data (case 7 only)"
                     possible_values="true or false"/>

                <nml_option name="config_native_gwd_static"     type="logical"       default_value="true"
                     units="-"
                     description="Whether to recompute subgrid-scale orography statistics directly on the native MPAS mesh (case 7 only)"
                     possible_values="true or false"/>

                <nml_option name="config_native_gwd_gsl_static"     type="logical"       default_value="false"
                     units="-"
                     description="Whether to recompute subgrid-scale orography statistics for the GSL drag directly on the native MPAS mesh (case 7 only)"
                     possible_values="true or false"/>

                <nml_option name="config_gwd_cell_scaling"      type="real"          default_value="1.0"      in_defaults="false"
                     units="-"
                     description="Scaling factor for the effective grid cell diameter used in computation of GWD static fields"
                     possible_values="Positive real values"/>

                <nml_option name="config_vertical_grid"         type="logical"       default_value="true"
                     units="-"
                     description="Whether to generate vertical grid"
                     possible_values="true or false"/>

                <nml_option name="config_met_interp"            type="logical"       default_value="true"
                     units="-"
                     description="Whether to interpolate first-guess fields from intermediate file"
                     possible_values="true or false"/>

                <nml_option name="config_input_sst"             type="logical"       default_value="false"
                     units="-"
                     description="Whether to re-compute SST and sea-ice fields from surface input data set; should be set to .true. when running case 8"
                     possible_values="true or false"/>

                <nml_option name="config_frac_seaice"           type="logical"       default_value="true"
                     units="-"
                     description="Whether to switch sea-ice threshold from 0.5 to 0.02"
                     possible_values="true or false"/>
        </nml_record>

        <nml_record name="physics" in_defaults="false">
                <nml_option name="config_tsk_seaice_threshold" type="real" default_value="100." units="K"
                     description="surface temperature threshold below which water points are set to seaice points"
                     possible_values="Positive real values"/>
        </nml_record>

        <nml_record name="io" in_defaults="true">

                <nml_option name="config_pio_num_iotasks" type="integer" default_value="0"
                     units="-"
                     description="Number of tasks to perform file I/O"
                     possible_values="Integer valued, 0 $\leq$ config_pio_num_iotasks $\leq$ \# MPI tasks, 0 indicates all tasks perform I/O"/>

                <nml_option name="config_pio_stride" type="integer" default_value="1"
                     units="-"
                     description="Stride between file I/O tasks"
                     possible_values="Integer valued, $\leq$ (\# MPI tasks) / config_pio_num_iotasks"/>

        </nml_record>

        <nml_record name="decomposition" in_defaults="true">

                <nml_option name="config_block_decomp_file_prefix" type="character" default_value="x1.40962.graph.info.part."
                     units="-"
                     description="Prefix of graph decomposition file, to be suffixed with the MPI task count"
                     possible_values="Any valid filename"/>

                <nml_option name="config_number_of_blocks" type="integer" default_value="0" in_defaults="false"
                     units="-"
                     description="Number of blocks to assign to each MPI task"
                     possible_values="Positive integer values"/>

                <nml_option name="config_explicit_proc_decomp" type="logical" default_value="false" in_defaults="false"
                     units="-"
                     description="Whether to use an explicit mapping of blocks to MPI tasks"
                     possible_values=".true. or .false."/>

                <nml_option name="config_proc_decomp_file_prefix" type="character" default_value="graph.info.part." in_defaults="false"
                     units="-"
                     description="Prefix of block mapping file"
                     possible_values="Any valid filename"/>

        </nml_record>


<!-- **************************************************************************************** -->
<!-- **************************************  Packages  ************************************** -->
<!-- **************************************************************************************** -->
        <packages>
                <package name="initial_conds" description="Used by test cases that produce initial conditions."/>
                <package name="sfc_update" description="Used by test cases that produce surface updates."/>
                <package name="lbcs" description="Active when lateral boundary conditions are interpolated and written."/>
                <package name="gwd_stage_in" description="Active only if GWD static fields are be generated and other static fields are not being generated"/>
                <package name="gwd_stage_out" description="Active only if GWD static fields are be generated and other static fields are not being generated"/>
                <package name="gwd_gsl_stage_out" description="Active only if GSL GWD static fields are be generated"/>
                <package name="vertical_stage_in" description="Active only if a vertical grid is being generated"/>
                <package name="vertical_stage_out" description="Active only if a vertical grid is being generated"/>
                <package name="met_stage_in" description="Active only if meteorological fields are being interpolated"/>
                <package name="met_stage_out" description="Active only if meteorological fields are being interpolated"/>
                <package name="first_guess_field" description="3-d atmospheric or land-surface fields on first-guess levels"/>
                <package name="mp_thompson_aers_in" description="initialization of GOCART-based Thompson water- and ice-friendly aerosols"/>
                <package name="noahmp" description="Whether to process, read, and write static fields only required by Noah-MP"/>
        </packages>


<!-- **************************************************************************************** -->
<!-- *************************************** Streams **************************************** -->
<!-- **************************************************************************************** -->

	<streams>
		<stream name="input" 
                        type="input" 
                        filename_template="x1.40962.grid.nc"
                        input_interval="initial_only" 
                        immutable="true">

                        <var name="latCell"/>
                        <var name="lonCell"/>
                        <var name="xCell"/>
                        <var name="yCell"/>
                        <var name="zCell"/>
                        <var name="indexToCellID"/>
                        <var name="latEdge"/>
                        <var name="lonEdge"/>
                        <var name="xEdge"/>
                        <var name="yEdge"/>
                        <var name="zEdge"/>
                        <var name="indexToEdgeID"/>
                        <var name="latVertex"/>
                        <var name="lonVertex"/>
                        <var name="xVertex"/>
                        <var name="yVertex"/>
                        <var name="zVertex"/>
                        <var name="indexToVertexID"/>
                        <var name="cellsOnEdge"/>
                        <var name="nEdgesOnCell"/>
                        <var name="nEdgesOnEdge"/>
                        <var name="edgesOnCell"/>
                        <var name="edgesOnEdge"/>
                        <var name="weightsOnEdge"/>
                        <var name="dvEdge"/>
                        <var name="dcEdge"/>
                        <var name="angleEdge"/>
                        <var name="areaCell"/>
                        <var name="areaTriangle"/>
                        <var name="cellsOnCell"/>
                        <var name="verticesOnCell"/>
                        <var name="verticesOnEdge"/>
                        <var name="edgesOnVertex"/>
                        <var name="cellsOnVertex"/>
                        <var name="kiteAreasOnVertex"/>
                        <var name="meshDensity"/>
                        <var name="nominalMinDc"/>
                        <var name="bdyMaskCell"/>
                        <var name="bdyMaskEdge"/>
                        <var name="bdyMaskVertex"/>
                        <var name="edgeNormalVectors" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="localVerticalUnitVectors" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="cellTangentPlane" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="fEdge" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="fVertex" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="ter" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="landmask" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="mminlu" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="isice_lu" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="iswater_lu" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="ivgtyp" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="isltyp" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="snoalb" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soiltemp" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="greenfrac" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="shdmin" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="shdmax" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="albedo12m" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soilcomp" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soilcl1" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soilcl2" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soilcl3" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="soilcl4" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="var2d" packages="vertical_stage_in;met_stage_in"/>
                        <var name="con" packages="vertical_stage_in;met_stage_in"/>
                        <var name="oa1" packages="vertical_stage_in;met_stage_in"/>
                        <var name="oa2" packages="vertical_stage_in;met_stage_in"/>
                        <var name="oa3" packages="vertical_stage_in;met_stage_in"/>
                        <var name="oa4" packages="vertical_stage_in;met_stage_in"/>
                        <var name="ol1" packages="vertical_stage_in;met_stage_in"/>
                        <var name="ol2" packages="vertical_stage_in;met_stage_in"/>
                        <var name="ol3" packages="vertical_stage_in;met_stage_in"/>
                        <var name="ol4" packages="vertical_stage_in;met_stage_in"/>
                        <var name="deriv_two" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="cell_gradient_coef_x" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="cell_gradient_coef_y" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="coeffs_reconstruct" packages="gwd_stage_in;vertical_stage_in;met_stage_in"/>
                        <var name="cf1" packages="met_stage_in"/>
                        <var name="cf2" packages="met_stage_in"/>
                        <var name="cf3" packages="met_stage_in"/>
                        <var name="zgrid" packages="met_stage_in"/>
                        <var name="rdzw" packages="met_stage_in"/>
                        <var name="dzu" packages="met_stage_in"/>
                        <var name="rdzu" packages="met_stage_in"/>
                        <var name="fzm" packages="met_stage_in"/>
                        <var name="fzp" packages="met_stage_in"/>
                        <var name="zxu" packages="met_stage_in"/>
                        <var name="zz" packages="met_stage_in"/>
                        <var name="zb" packages="met_stage_in"/>
                        <var name="zb3" packages="met_stage_in"/>
                        <var name="dss" packages="met_stage_in"/>
		</stream>

		<stream name="output" 
                        type="output" 
                        filename_template="x1.40962.init.nc"
                        output_interval="initial_only"
                        packages="initial_conds"
                        immutable="true">

                        <var name="initial_time" packages="met_stage_out"/>
                        <var name="xtime"/>
                        <var name="Time"/>
                        <var name="latCell"/>
                        <var name="lonCell"/>
                        <var name="xCell"/>
                        <var name="yCell"/>
                        <var name="zCell"/>
                        <var name="indexToCellID"/>
                        <var name="latEdge"/>
                        <var name="lonEdge"/>
                        <var name="xEdge"/>
                        <var name="yEdge"/>
                        <var name="zEdge"/>
                        <var name="indexToEdgeID"/>
                        <var name="latVertex"/>
                        <var name="lonVertex"/>
                        <var name="xVertex"/>
                        <var name="yVertex"/>
                        <var name="zVertex"/>
                        <var name="indexToVertexID"/>
                        <var name="cellsOnEdge"/>
                        <var name="nEdgesOnCell"/>
                        <var name="nEdgesOnEdge"/>
                        <var name="edgesOnCell"/>
                        <var name="edgesOnEdge"/>
                        <var name="weightsOnEdge"/>
                        <var name="dvEdge"/>
                        <var name="dcEdge"/>
                        <var name="angleEdge"/>
                        <var name="areaCell"/>
                        <var name="areaTriangle"/>
                        <var name="cellsOnCell"/>
                        <var name="verticesOnCell"/>
                        <var name="verticesOnEdge"/>
                        <var name="edgesOnVertex"/>
                        <var name="cellsOnVertex"/>
                        <var name="kiteAreasOnVertex"/>
                        <var name="meshDensity"/>
                        <var name="nominalMinDc"/>
                        <var name="bdyMaskCell"/>
                        <var name="bdyMaskEdge"/>
                        <var name="bdyMaskVertex"/>
                        <var name="edgeNormalVectors"/>
                        <var name="localVerticalUnitVectors"/>
                        <var name="cellTangentPlane"/>
                        <var name="fEdge"/>
                        <var name="fVertex"/>
                        <var name="ter"/>
                        <var name="landmask"/>
                        <var name="mminlu"/>
                        <var name="isice_lu"/>
                        <var name="iswater_lu"/>
                        <var name="ivgtyp"/>
                        <var name="isltyp"/>
                        <var name="snoalb"/>
                        <var name="soiltemp"/>
                        <var name="greenfrac"/>
                        <var name="shdmin"/>
                        <var name="shdmax"/>
                        <var name="albedo12m"/>
                        <var name="soilcomp"/>
                        <var name="soilcl1"/>
                        <var name="soilcl2"/>
                        <var name="soilcl3"/>
                        <var name="soilcl4"/>
                        <var name="var2d" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="con" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="oa1" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="oa2" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="oa3" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="oa4" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="ol1" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="ol2" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="ol3" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="ol4" packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>
                        <var name="deriv_two"/>
                        <var name="cell_gradient_coef_x"/>
                        <var name="cell_gradient_coef_y"/>
                        <var name="coeffs_reconstruct"/>
                        <var name="cf1" packages="vertical_stage_out;met_stage_out"/>
                        <var name="cf2" packages="vertical_stage_out;met_stage_out"/>
                        <var name="cf3" packages="vertical_stage_out;met_stage_out"/>
                        <var name="zgrid" packages="vertical_stage_out;met_stage_out"/>
                        <var name="rdzw" packages="vertical_stage_out;met_stage_out"/>
                        <var name="dzu" packages="vertical_stage_out;met_stage_out"/>
                        <var name="rdzu" packages="vertical_stage_out;met_stage_out"/>
                        <var name="fzm" packages="vertical_stage_out;met_stage_out"/>
                        <var name="fzp" packages="vertical_stage_out;met_stage_out"/>
                        <var name="zxu" packages="vertical_stage_out;met_stage_out"/>
                        <var name="zz" packages="vertical_stage_out;met_stage_out"/>
                        <var name="zb" packages="vertical_stage_out;met_stage_out"/>
                        <var name="zb3" packages="vertical_stage_out;met_stage_out"/>
                        <var name="dss" packages="vertical_stage_out;met_stage_out"/>
                        <var name="u_init" packages="met_stage_out"/>
                        <var name="v_init" packages="met_stage_out"/>
                        <var name="t_init" packages="met_stage_out"/>
                        <var name="qv_init" packages="met_stage_out"/>
                        <var_array name="scalars" packages="met_stage_out;mp_thompson_aers_in"/>
                        <var name="u" packages="met_stage_out"/>
                        <var name="w" packages="met_stage_out"/>
                        <var name="dz" packages="met_stage_out"/>
                        <var name="dzs" packages="met_stage_out"/>
                        <var name="zs" packages="met_stage_out"/>
                        <var name="sh2o" packages="met_stage_out"/>
                        <var name="smois" packages="met_stage_out"/>
                        <var name="tslb" packages="met_stage_out"/>
                        <var name="tmn" packages="met_stage_out"/>
                        <var name="skintemp" packages="met_stage_out"/>
                        <var name="sst" packages="met_stage_out"/>
                        <var name="snow" packages="met_stage_out"/>
                        <var name="snowc" packages="met_stage_out"/>
                        <var name="snowh" packages="met_stage_out"/>
                        <var name="xice" packages="met_stage_out"/>
                        <var name="seaice" packages="met_stage_out"/>
                        <var name="vegfra" packages="met_stage_out"/>
                        <var name="sfc_albbck" packages="met_stage_out"/>
                        <var name="xland" packages="met_stage_out"/>
                        <var name="rho" packages="met_stage_out"/>
                        <var name="theta" packages="met_stage_out"/>
                        <var name="relhum" packages="met_stage_out"/>
                        <var name="rho_base" packages="met_stage_out"/>
                        <var name="theta_base" packages="met_stage_out"/>
                        <var name="surface_pressure" packages="met_stage_out"/>
                        <var name="precipw" packages="met_stage_out"/>
                        <var name="h_oml_initial" packages="met_stage_out"/>
                        <var name="u10" packages="met_stage_out"/>
                        <var name="v10" packages="met_stage_out"/>
                        <var name="q2" packages="met_stage_out"/>
                        <var name="rh2" packages="met_stage_out"/>
                        <var name="t2m" packages="met_stage_out"/>
		</stream>

                <stream name="ugwp_oro_data"
                        type="output"
                        filename_template="x1.40962.ugwp_oro_data.nc"
                        output_interval="initial_only"
                        immutable="true"
                        packages="gwd_gsl_stage_out">

                        <var name="var2dls"/>
                        <var name="conls"/>
                        <var name="oa1ls"/>
                        <var name="oa2ls"/>
                        <var name="oa3ls"/>
                        <var name="oa4ls"/>
                        <var name="ol1ls"/>
                        <var name="ol2ls"/>
                        <var name="ol3ls"/>
                        <var name="ol4ls"/>
                        <var name="var2dss"/>
                        <var name="conss"/>
                        <var name="oa1ss"/>
                        <var name="oa2ss"/>
                        <var name="oa3ss"/>
                        <var name="oa4ss"/>
                        <var name="ol1ss"/>
                        <var name="ol2ss"/>
                        <var name="ol3ss"/>
                        <var name="ol4ss"/>
                </stream>

		<stream name="surface" 
                        type="output" 
                        filename_template="x1.40962.sfc_update.nc"
                        output_interval="86400" 
                        filename_interval="none"
                        packages="sfc_update"
                        immutable="true">

			<var name="xtime"/>
			<var name="Time"/>
			<var name="sst"/>
			<var name="xice"/>
		</stream>

		<stream name="lbc"
                        type="output"
                        filename_template="lbc.$Y-$M-$D_$h.$m.$s.nc"
                        output_interval="3:00:00"
                        filename_interval="output_interval"
                        packages="lbcs"
                        immutable="true">

			<var name="xtime"/>
			<var name="Time"/>
			<var_struct name="lbc_state"/>
		</stream>

	</streams>


<!-- **************************************************************************************** -->
<!-- ************************************** Variables *************************************** -->
<!-- **************************************************************************************** -->

        <var_struct name="mesh" time_levs="1">

                <var name="latCell" type="real" dimensions="nCells" units="rad"
                     description="Latitude of cells"/>

                <var name="lonCell" type="real" dimensions="nCells" units="rad"
                     description="Longitude of cells"/>

                <var name="xCell" type="real" dimensions="nCells" units="m"
                     description="Cartesian x-coordinate of cells"/>

                <var name="yCell" type="real" dimensions="nCells" units="m"
                     description="Cartesian y-coordinate of cells"/>

                <var name="zCell" type="real" dimensions="nCells" units="m"
                     description="Cartesian z-coordinate of cells"/>

                <var name="indexToCellID" type="integer" dimensions="nCells" units="-"
                     description="Mapping from local array index to global cell ID"/>

                <var name="latEdge" type="real" dimensions="nEdges" units="rad"
                     description="Latitude of edges"/>

                <var name="lonEdge" type="real" dimensions="nEdges" units="rad"
                     description="Longitude of edges"/>

                <var name="xEdge" type="real" dimensions="nEdges" units="m"
                     description="Cartesian x-coordinate of edges"/>

                <var name="yEdge" type="real" dimensions="nEdges" units="m"
                     description="Cartesian y-coordinate of edges"/>

                <var name="zEdge" type="real" dimensions="nEdges" units="m"
                     description="Cartesian z-coordinate of edges"/>

                <var name="indexToEdgeID" type="integer" dimensions="nEdges" units="-"
                     description="Mapping from local array index to global edge ID"/>

                <var name="latVertex" type="real" dimensions="nVertices" units="rad"
                     description="Latitude of vertices"/>

                <var name="lonVertex" type="real" dimensions="nVertices" units="rad"
                     description="Longitude of vertices"/>

                <var name="xVertex" type="real" dimensions="nVertices" units="m"
                     description="Cartesian x-coordinate of vertices"/>

                <var name="yVertex" type="real" dimensions="nVertices" units="m"
                     description="Cartesian y-coordinate of vertices"/>

                <var name="zVertex" type="real" dimensions="nVertices" units="m"
                     description="Cartesian z-coordinate of vertices"/>

                <var name="indexToVertexID" type="integer" dimensions="nVertices" units="-"
                     description="Mapping from local array index to global vertex ID"/>

                <var name="cellsOnEdge" type="integer" dimensions="TWO nEdges" units="-"
                     description="IDs of cells divided by an edge"/>

                <var name="nEdgesOnCell" type="integer" dimensions="nCells" units="-"
                     description="Number of edges forming the boundary of a cell"/>

                <var name="nEdgesOnEdge" type="integer" dimensions="nEdges" units="-"
                     description="Number of edges involved in reconstruction of tangential velocity for an edge"/>

                <var name="edgesOnCell" type="integer" dimensions="maxEdges nCells" units="-"
                     description="IDs of edges forming the boundary of a cell"/>

                <var name="edgesOnEdge" type="integer" dimensions="maxEdges2 nEdges" units="-"
                     description="IDs of edges involved in reconstruction of tangential velocity for an edge"/>

                <var name="weightsOnEdge" type="real" dimensions="maxEdges2 nEdges" units="-"
                     description="Weights used in reconstruction of tangential velocity for an edge"/>

                <var name="dvEdge" type="real" dimensions="nEdges" units="-"
                     description="Spherical distance between vertex endpoints of an edge"/>

                <var name="dcEdge" type="real" dimensions="nEdges" units="m"
                     description="Spherical distance between cells separated by an edge"/>

                <var name="angleEdge" type="real" dimensions="nEdges" units="rad"
                     description="Angle between local north and the positive tangential direction of an edge"/>

                <var name="areaCell" type="real" dimensions="nCells" units="m^2"
                     description="Spherical area of a Voronoi cell"/>

                <var name="areaTriangle" type="real" dimensions="nVertices" units="m^2"
                     description="Spherical area of a Delaunay triangle"/>

                <var name="edgeNormalVectors" type="real" dimensions="R3 nEdges" units="unitless"
                     description="Cartesian components of the vector normal to an edge and tangential to the surface of the sphere"/>

                <var name="localVerticalUnitVectors" type="real" dimensions="R3 nCells" units="unitless"
                     description="Cartesian components of the vector pointing in the local vertical direction for a cell"/>

                <var name="cellTangentPlane" type="real" dimensions="R3 TWO nCells" units="unitless"
                     description="Components of a pair of vectors defining the tangent plane at a cell"/>

                <var name="cellsOnCell" type="integer" dimensions="maxEdges nCells" units="-"
                     description="IDs of cells neighboring a cell"/>

                <var name="verticesOnCell" type="integer" dimensions="maxEdges nCells" units="-"
                     description="IDs of vertices (corner points) of a cell"/>

                <var name="verticesOnEdge" type="integer" dimensions="TWO nEdges" units="-"
                     description="IDs of the two vertex endpoints of an edge"/>

                <var name="edgesOnVertex" type="integer" dimensions="vertexDegree nVertices" units="-"
                     description="IDs of the edges that meet at a vertex"/>

                <var name="cellsOnVertex" type="integer" dimensions="vertexDegree nVertices" units="-"
                     description="IDs of the cells that meet at a vertex"/>

                <var name="kiteAreasOnVertex" type="real" dimensions="vertexDegree nVertices" units="m^2"
                     description="Intersection areas between primal (Voronoi) and dual (triangular) mesh cells"/>

                <var name="fEdge" type="real" dimensions="nEdges" units="unitless"
                     description="Coriolis parameter at an edge"/>

                <var name="fVertex" type="real" dimensions="nVertices" units="unitless"
                     description="Coriolis parameter at a vertex"/>

                <var name="meshDensity" type="real" dimensions="nCells" units="unitless"
                     description="Mesh density function (used when generating the mesh) evaluated at a cell"/>

                <var name="nominalMinDc" type="real" dimensions="" units="m"
                     description="Nominal minimum dcEdge value where meshDensity == 1.0"/>

                <!-- coefficients for vertical extrapolation to the surface -->
                <var name="cf1" type="real" dimensions="" units="unitless"
                     description="Surface interpolation weight for level k=1 value"/>

                <var name="cf2" type="real" dimensions="" units="unitless"
                     description="Surface interpolation weight for level k=2 value"/>

                <var name="cf3" type="real" dimensions="" units="unitless"
                     description="Surface interpolation weight for level k=3 value"/>

                <!-- static terrestrial fields -->
                <var name="ter" type="real" dimensions="nCells" units="m"
                     description="terrain height"/>

                <var name="PHIS" type="real" dimensions="nCells" units="m^2 s^{-2}"
                     description="CAM Surface geopotential"/>

                <var name="landmask" type="integer"  dimensions="nCells" units="unitless"
                     description="land-ocean mask (1=land ; 0=ocean)"/>

                <var name="ivgtyp" name_in_code="lu_index" type="integer"  dimensions="nCells" units="unitless"
                     description="dominant vegetation category"/>

                <var name="mminlu" type="text" dimensions="" units="unitless"
                     description="land use classification"/>

                <var name="isltyp" name_in_code="soilcat_top"  type="integer"  dimensions="nCells" units="unitless"
                     description="dominate soil category"/>

                <var name="snoalb" type="real" dimensions="nCells" units="unitless"
                     description="annual maximum snow albedo"/>

                <var name="soiltemp" type="real" dimensions="nCells" units="K"
                     description="annual mean deep soil temperature"/>

                <var name="greenfrac" type="real" dimensions="nMonths nCells" units="percent"
                     description="monthly-mean climatological greenness fraction"/>

                <var name="shdmin" type="real" dimensions="nCells" units="percent"
                     description="minimum fractional coverage of annual green vegetation fraction"/>

                <var name="shdmax" type="real" dimensions="nCells" units="percent"
                     description="maximum fractional coverage of annual green vegetation fraction"/>

                <var name="albedo12m" type="real" dimensions="nMonths nCells" units="percent"
                     description="monthly-mean climatological surface albedo"/>

                <var name="isice_lu" type="integer" dimensions="" units="unitless" default_value="24"
                     description="Index category for snow/ice"/>

                <var name="iswater_lu" type="integer" dimensions="" units="unitless" default_value="16"
                     description="Index category for water"/>

                <!-- SOIL COMPOSITION fields needed for the NOAH-MP land surface scheme -->
                <var name="soilcomp" type="real" dimensions="nSoilComps nCells" units="percent"
                     description="soil composition needed as input in the NOAH-MP land surface model"
                     packages="noahmp"/>

                <var name="soilcl1" type="real" dimensions="nCells" units="unitless"
                     description="soil texture class level 1 needed as input for the NOAH-MP land surface model"
                     packages="noahmp"/>

                <var name="soilcl2" type="real" dimensions="nCells" units="unitless"
                     description="soil texture class level 2 needed as input for the NOAH-MP land surface model"
                     packages="noahmp"/>

                <var name="soilcl3" type="real" dimensions="nCells" units="unitless"
                     description="soil texture class level 3 needed as input for the NOAH-MP land surface model"
                     packages="noahmp"/>

                <var name="soilcl4" type="real" dimensions="nCells" units="unitless"
                     description="soil texture class level 4 needed as input for the NOAH-MP land surface model"
                     packages="noahmp"/>

                <!-- GWDO fields -->
                <var name="var2d" type="real" dimensions="nCells" units="m"
                     description="standard deviation of subgrid-scale orography"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="con" type="real" dimensions="nCells" units="unitless"
                     description="orographic convexity"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="oa1" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="oa2" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for southerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="oa3" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for south-westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="oa4" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for north-westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="ol1" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="ol2" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for southerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="ol3" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for south-westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <var name="ol4" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for north-westerly flow"
                     packages="gwd_stage_out;vertical_stage_out;met_stage_out"/>

                <!-- GSL GWDO fields - meso-scale -->
                <var name="var2dls" type="real" dimensions="nCells" units="m"
                     description="standard deviation of subgrid-scale orography (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="conls" type="real" dimensions="nCells" units="unitless"
                     description="orographic convexity (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa1ls" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa2ls" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for southerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa3ls" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for south-westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa4ls" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for north-westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol1ls" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol2ls" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for southerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol3ls" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for south-westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol4ls" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for north-westerly flow (meso-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <!-- GSL GWDO fields - small-scale -->
                <var name="var2dss" type="real" dimensions="nCells" units="m"
                     description="standard deviation of subgrid-scale orography (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="conss" type="real" dimensions="nCells" units="unitless"
                     description="orographic convexity (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa1ss" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa2ss" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for southerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa3ss" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for south-westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="oa4ss" type="real" dimensions="nCells" units="unitless"
                     description="asymmetry of subgrid-scale orography for north-westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol1ss" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol2ss" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for southerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol3ss" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for south-westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <var name="ol4ss" type="real" dimensions="nCells" units="unitless"
                     description="effective orographic length for north-westerly flow (small-scale GWD)"
                     packages="gwd_gsl_stage_out"/>

                <!-- GOCART aerosol fields -->
                <var name="nifa_gocart_clim" type="real" dimensions="nMonths nGocartLevels nCells" units="nb kg^{-1}"
                     description="climatological Gocart ice-friendly number concentration"
                     packages="mp_thompson_aers_in"/>

                <var name="nwfa_gocart_clim" type="real" dimensions="nMonths nGocartLevels nCells" units="nb kg^{-1}"
                     description="climatological Gocart water-friendly number concentration"
                     packages="mp_thompson_aers_in"/>

                <var name="pwif_gocart_clim" type="real"  dimensions="nMonths nGocartLevels nCells" units="Pa"
                     description="climatological Gocart pressure levels"
                     packages="mp_thompson_aers_in"/>

                <!-- description of the vertical grid structure -->
                <var name="hx" type="real" dimensions="nVertLevelsP1 nCells" units="m"
                     description="terrain influence in vertical coordinate, $h_s(x,y,\zeta)$ in Klemp (MWR 2011)"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="zgrid" type="real" dimensions="nVertLevelsP1 nCells" units="m MSL"
                     description="Geometric height of layer interfaces"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="rdzw" type="real" dimensions="nVertLevels" units="unitless"
                     description="Reciprocal dzw"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="dzu" type="real" dimensions="nVertLevels" units="unitless"
                     description="d(zeta) at w levels"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="rdzu" type="real" dimensions="nVertLevels" units="unitless"
                     description="Reciprocal dzu"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="fzm" type="real" dimensions="nVertLevels" units="unitless"
                     description="Weight for linear interpolation to w(k) point for u(k) level variable"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="fzp" type="real" dimensions="nVertLevels" units="unitless"
                     description="Weight for linear interpolation to w(k) point for u(k-1) level variable"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="zxu" type="real" dimensions="nVertLevels nEdges" units="unitless"
                     description="dz/dx on horizontal coordinate surfaces at u levels"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="zz" type="real" dimensions="nVertLevels nCells" units="unitless"
                     description="d(zeta)/dz, vertical metric term"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="zb" type="real" dimensions="nVertLevelsP1 TWO nEdges" units="unitless"
                     description="Coefficients for contribution from u to omega diagnosis, edge-oriented"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="zb3" type="real" dimensions="nVertLevelsP1 TWO nEdges" units="unitless"
                     description="Coefficients for 3rd-order correction to contribution from u to omega diagnosis, edge-oriented"
                     packages="vertical_stage_out;met_stage_out"/>

                <!-- W-Rayleigh damping coefficient -->
                <var name="dss" type="real" dimensions="nVertLevels nCells" units="unitless"
                     description="w-damping coefficient"
                     packages="vertical_stage_out;met_stage_out"/>

                <var name="u_init" type="real" dimensions="nVertLevels" units="m s^{-1}"
                     description="u reference profile"
                     packages="met_stage_out"/>

                <var name="v_init" type="real" dimensions="nVertLevels" units="m s^{-1}"
                     description="v reference profile"
                     packages="met_stage_out"/>

                <var name="t_init" type="real" dimensions="nVertLevels nCells" units="K"
                     description="theta reference profile"
                     packages="met_stage_out"/>

                <var name="qv_init" type="real" dimensions="nVertLevels" units="kg kg^{-1}"
                     description="qv reference profile"
                     packages="met_stage_out"/>

                <!-- variables needed for advection -->
                <var name="deriv_two" type="real" dimensions="FIFTEEN TWO nEdges" units="unitless"
                     description="weights for cell-centered second derivative, normal to edge, for transport scheme"/>

                <var name="advCells" type="integer" dimensions="TWENTYONE nCells" units="-"
                     description="cells used in least squares fit of polynomial for advection"/>

                <var name="cell_gradient_coef_x" type="real" dimensions="maxEdges nCells" units="m^-1"
                     description="Coefficients for computing the x (zonal) derivative of a cell-centered variable"/>

                <var name="cell_gradient_coef_y" type="real" dimensions="maxEdges nCells" units="m^-1"
                     description="Coefficients for computing the y (meridional) derivative of a cell-centered variable"/>

                <!-- arrays required for reconstruction of velocity field -->
                <var name="coeffs_reconstruct" type="real" dimensions="R3 maxEdges nCells" units="unitless"
                     description="Coefficients to reconstruct velocity vectors at cell centers"/>

                <var name="bdyMaskCell" type="integer" dimensions="nCells" default_value="0" units="-"
                     description="Indicator of whether a cell is an interior cell, a relaxation-zone cell, or a specified-zone cell"/>

                <var name="bdyMaskEdge" type="integer" dimensions="nEdges" default_value="0" units="-"
                     description="Indicator of whether an edge is an interior edge, a relaxation-zone edge, or a specified-zone edge"/>

                <var name="bdyMaskVertex" type="integer" dimensions="nVertices" default_value="0" units="-"
                     description="Indicator of whether an vertex is an interior vertex, a relaxation-zone vertex, or a specified-zone vertex"/>
        </var_struct>

        <var_struct name="state" time_levs="1">
                <var name="initial_time" type="text" dimensions="" units="YYYY-MM-DD_hh:mm:ss"
                     description="Model initialization time"/>

                <var name="xtime" type="text" dimensions="Time" units="YYYY-MM-DD_hh:mm:ss"
                     description="Model valid time"/>

                <var name="Time" type="real" dimensions="Time" units="seconds since YYYY-MM-DD hh:mm:ss"
                     description="CF-compliant valid time" standard_name="time"/>

                <var name="u" type="real" dimensions="nVertLevels nEdges Time" units="m s^{-1}"
                     description="Horizontal normal velocity at edges"
                     packages="met_stage_out"/>

                <var name="w" type="real" dimensions="nVertLevelsP1 nCells Time" units="m s^{-1}"
                     description="Vertical velocity at vertical cell faces"
                     packages="met_stage_out"/>

                <var name="rho_zz" type="real" dimensions="nVertLevels nCells Time" units="kg m^{-3}"
                     description="Dry air density divided by d(zeta)/dz"
                     packages="met_stage_out"/>

                <var name="theta_m" type="real" dimensions="nVertLevels nCells Time" units="K"
                     description="Moist potential temperature: theta*(1+q_v*R_v/R_d)"
                     packages="met_stage_out"/>

                <var_array name="scalars" type="real" dimensions="nVertLevels nCells Time" packages="met_stage_out">
                        <var name="qv" array_group="moist" units="kg kg^{-1}"
                             description="Water vapor mixing ratio"/>

                        <var name="qc" array_group="moist" units="kg kg^{-1}"
                             description="Cloud water mixing ratio"/>

                        <var name="qr" array_group="moist" units="kg kg^{-1}"
                             description="Rain water mixing ratio"/>

                        <var name="qAB" array_group="passive" units="kg kg^{-1}"
                             description="Molecular AB mixing ratio"/>

                        <var name="qA" array_group="passive" units="kg kg^{-1}"
                             description="Atomic A mixing ratio"/>

                        <var name="qB" array_group="passive" units="kg kg^{-1}"
                             description="Atomic B mixing ratio"/>

                        <var name="nifa" array_group="number" units="nb kg^{-1}"
                             description="Gocart ice-friendly aerosol number concentration"
                             packages="mp_thompson_aers_in"/>

                        <var name="nwfa" array_group="number" units="nb kg^{-1}"
                             description="Gocart water-friendly aerosol number concentration"
                             packages="mp_thompson_aers_in"/>
                        <var name="tke" array_group="turbulence" units="m^2 s^{-2}"
                             description="Turbulent kinetic energy for the prognostic tke LES scheme"/>
                </var_array>

                <!-- Climatology for ocean mixed-layer depth -->
                <var name="h_oml_initial" type="real" dimensions="nCells Time" units="m"
                     description="Initial depth of ocean mix layer"
                     packages="met_stage_out"/>
        </var_struct>

        <var_struct name="lbc_state" time_levs="1">
                <var name="lbc_u" type="real" dimensions="nVertLevels nEdges Time" packages="lbcs"
                     units="m s^{-1}"
                     description="Horizontal normal velocity on domain lateral boundary edges"/>

                <var name="lbc_w" type="real" dimensions="nVertLevelsP1 nCells Time" packages="lbcs"
                     units="m s^{-1}"
                     description="Vertical velocity on domain lateral boundary vertical cell faces"/>

                <var name="lbc_rho" type="real" dimensions="nVertLevels nCells Time" packages="lbcs"
                     units="kg m^{-3}"
                     description="Dry air density on lateral boundary cells"/>

                <var name="lbc_theta" type="real" dimensions="nVertLevels nCells Time" packages="lbcs"
                     units="K"
                     description="Potential temperature on lateral boundary cells"/>

                <var_array name="lbc_scalars" type="real" dimensions="nVertLevels nCells Time" packages="lbcs">
                        <var name="lbc_qv" array_group="moist" units="kg kg^{-1}"
                             description="Water vapor mixing ratio on lateral boundary cells"/>

                        <var name="lbc_qc" array_group="moist" units="kg kg^{-1}"
                             description="Cloud water mixing ratio on lateral boundary cells"/>

                        <var name="lbc_qr" array_group="moist" units="kg kg^{-1}"
                             description="Rain water mixing ratio on lateral boundary cells"/>

                        <var name="lbc_nifa" array_group="number" units="nb kg^{-1}"
                             description="Gocart ice-friendly aerosol number concentration on lateral boundary cells"
                             packages="mp_thompson_aers_in"/>

                        <var name="lbc_nwfa" array_group="number" units="nb kg^{-1}"
                             description="Gocart water-friendly aerosol number concentration on lateral boundary cells"
                             packages="mp_thompson_aers_in"/>
                </var_array>

        </var_struct>

        <var_struct name="fg" time_levs="1">

                <!-- horizontally interpolated from first-guess data -->
                <var name="u_fg" name_in_code="u" type="real" dimensions="nFGLevels nEdges Time" units="m s^{-1}"
                     description="First guess zonal wind component"
                     packages="first_guess_field"/>

                <var name="v_fg" name_in_code="v" type="real" dimensions="nFGLevels nEdges Time" units="m s^{-1}"
                     description="First guess merdian wind component"
                     packages="first_guess_field"/>

                <var name="t_fg" name_in_code="t" type="real" dimensions="nFGLevels nCells Time" units="K"
                     description="First guess temperature"
                     packages="first_guess_field"/>

                <var name="p_fg" name_in_code="p" type="real" dimensions="nFGLevels nCells Time" units="Pa"
                     description="First guess pressure"
                     packages="first_guess_field"/>

                <var name="z_fg" name_in_code="z" type="real" dimensions="nFGLevels nCells Time" units="m"
                     description="First guess geopotential height"
                     packages="first_guess_field"/>

                <var name="rh_fg" name_in_code="rh" type="real" dimensions="nFGLevels nCells Time" units="percent"
                     description="First guess relative humidity with respect to liquid water"
                     packages="first_guess_field"/>

                <var name="sh_fg" name_in_code="sh" type="real" dimensions="nFGLevels nCells Time" units="kg kg^{-1}"
                     description="First guess specific humidity"
                     packages="first_guess_field"/>

                <var name="soilz_fg" name_in_code="soilz" type="real" dimensions="nCells Time" units="m"
                     description="First guess soil height"
                     packages="first_guess_field"/>

                <var name="psfc_fg" name_in_code="psfc" type="real" dimensions="nCells Time" units="Pa"
                     description="First guess surface pressure"
                     packages="first_guess_field"/>

                <var name="pmsl_fg" name_in_code="pmsl" type="real" dimensions="nCells Time" units="Pa"
                     description="First guess mean sea level pressure"
                     packages="first_guess_field"/>

                <var name="dz_fg" type="real" dimensions="nFGSoilLevels nCells Time" units="cm"
                     description="First guess depth of soil layer bottom"
                     packages="first_guess_field"/>

                <var name="dzs_fg" type="real" dimensions="nFGSoilLevels nCells Time" units="cm"
                     description="First guess soil layer thickness"
                     packages="first_guess_field"/>

                <var name="zs_fg" type="real" dimensions="nFGSoilLevels nCells Time" units="cm"
                     description="First guess depth of centers of soil levels"
                     packages="first_guess_field"/>

                <var name="st_fg" type="real" dimensions="nFGSoilLevels nCells Time" units="k"
                     description="First guess soil temperature"
                     packages="first_guess_field"/>

                <var name="sm_fg" type="real" dimensions="nFGSoilLevels nCells Time" units="m^{3} m^{-3}"
                     description="First guess soil moisture"
                     packages="first_guess_field"/>

                <!-- horizontally interpolated from first-guess data, and should be read in by model -->
                <var name="dz" type="real" dimensions="nSoilLevels nCells Time" units="m"
                     description="depth of soil layer bottom"
                     packages="met_stage_out"/>

                <var name="dzs" type="real" dimensions="nSoilLevels nCells Time" units="m"
                     description="soil layer thickness"
                     packages="met_stage_out"/>

                <var name="zs" type="real" dimensions="nSoilLevels nCells Time" units="m"
                     description="depth of centers of soil layers"
                     packages="met_stage_out"/>

                <var name="sh2o" type="real" dimensions="nSoilLevels nCells Time" units="m3 m^{-3}"
                     description="soil equivalent liquid water "
                     packages="met_stage_out"/>

                <var name="smois" type="real" dimensions="nSoilLevels nCells Time" units="m3 m^{-3}"
                     description="soil moisture"
                     packages="met_stage_out"/>

                <var name="tslb" type="real" dimensions="nSoilLevels nCells Time" units="K"
                     description="soil layer temperature"
                     packages="met_stage_out"/>

                <var name="smcrel" type="real" dimensions="nSoilLevels nCells Time" units="m3 m^{-3}"
                     description="soil moisture threshold below which transpiration begins to stress"
                     packages="met_stage_out"/>

                <var name="tmn" type="real" dimensions="nCells Time" units="K"
                     description="deep soil temperature"
                     packages="met_stage_out"/>

                <var name="skintemp" type="real" dimensions="nCells Time" units="K"
                     description="ground or water surface temperature"
                     packages="met_stage_out"/>

                <var name="sst" type="real" dimensions="nCells Time" units="K"
                     description="sea-surface temperature"
                     packages="met_stage_out;sfc_update"/>

                <var name="snow" type="real" dimensions="nCells Time" units="kg m^{-2}"
                     description="snow water equivalent"
                     packages="met_stage_out"/>

                <var name="snowc" type="real" dimensions="nCells Time" units="unitless"
                     description="flag for snow on ground (=0 no snow; =1,otherwise"
                     packages="met_stage_out"/>

                <var name="snowh" type="real" dimensions="nCells Time" units="m"
                     description="physical snow depth"
                     packages="met_stage_out"/>

                <var name="xice" type="real" dimensions="nCells Time" units="unitless"
                     description="fractional area coverage of sea-ice"
                     packages="met_stage_out;sfc_update"/>

                <var name="seaice" type="real" dimensions="nCells Time" units="unitless"
                     description="sea-ice flag (0=no seaice; =1 otherwise)"
                     packages="met_stage_out"/>

                <var name="gfs_z" type="real" dimensions="nVertLevels nCells Time" units="m"
                     description="geopotential height vertically interpolated from first guess"
                     packages="met_stage_out"/>

                <var name="vegfra" type="real" dimensions="nCells Time" units="percent"
                     description="vegetation fraction"
                     packages="met_stage_out"/>

                <var name="sfc_albbck" type="real" dimensions="nCells Time" units="unitless"
                     description="background surface albedo"
                     packages="met_stage_out"/>

                <var name="xland" type="real" dimensions="nCells Time" units="unitless"
                     description="land-ocean mask (1=land including sea-ice ; 2=ocean)"
                     packages="met_stage_out"/>

                <!-- diagnostic fields that we get from first-guess data, and which we want to have at model initial output time -->
                <var name="u10" type="real" dimensions="nCells Time" units="m s^{-1}"
                     description="10-meter zonal wind"
                     packages="met_stage_out"/>

                <var name="v10" type="real" dimensions="nCells Time" units="m s^{-1}"
                     description="10-meter meridional wind"
                     packages="met_stage_out"/>

                <var name="q2" type="real" dimensions="nCells Time" units="kg kg^{-1}"
                     description="2-meter specific humidity"
                     packages="met_stage_out"/>

                <var name="rh2" type="real" dimensions="nCells Time" units="percent"
                     description="2-meter relative humidity"
                     packages="met_stage_out"/>

                <var name="t2m" type="real" dimensions="nCells Time" units="K"
                     description="2-meter temperature"
                     packages="met_stage_out"/>

        </var_struct>

        <var_struct name="diag" time_levs="1">
                <var name="pressure_p" type="real" dimensions="nVertLevels nCells Time" units="Pa"
                     description="Perturbation pressure"
                     packages="met_stage_out"/>

                <var name="rho" type="real" dimensions="nVertLevels nCells Time" units="kg m^{-3}"
                     description="Dry air density"
                     packages="met_stage_out"/>

                <var name="theta" type="real" dimensions="nVertLevels nCells Time" units="K"
                     description="Potential temperature"
                     packages="met_stage_out"/>

                <var name="v" type="real" dimensions="nVertLevels nEdges Time" units="m s^{-1}"
                     description="Horizontal tangential velocity at edges"
                     packages="met_stage_out"/>

                <var name="relhum" type="real" dimensions="nVertLevels nCells Time" units="percent"
                     description="Relative humidity"
                     packages="met_stage_out"/>

                <var name="spechum" type="real" dimensions="nVertLevels nCells Time" units="kg kg^{-1}"
                     description="Specific humidity"
                     packages="met_stage_out"/>

                <!-- reconstructed horizontal velocity vectors at cell centers -->
                <var name="uReconstructX" type="real" dimensions="nVertLevels nCells Time" units="m s^{-1}"
                     description="Cartesian x-component of reconstructed horizontal velocity at cell centers"
                     packages="met_stage_out"/>

                <var name="uReconstructY" type="real" dimensions="nVertLevels nCells Time" units="m s^{-1}"
                     description="Cartesian y-component of reconstructed horizontal velocity at cell centers"
                     packages="met_stage_out"/>

                <var name="uReconstructZ" type="real" dimensions="nVertLevels nCells Time" units="m s^{-1}"
                     description="Cartesian z-component of reconstructed horizontal velocity at cell centers"
                     packages="met_stage_out"/>

                <var name="uReconstructZonal" type="real" dimensions="nVertLevels nCells Time" units="m s^{-1}"
                     description="Zonal component of reconstructed horizontal velocity at cell centers"
                     packages="met_stage_out"/>

                <var name="uReconstructMeridional" type="real" dimensions="nVertLevels nCells Time" units="m s^{-1}"
                     description="Meridional component of reconstructed horizontal velocity at cell centers"
                     packages="met_stage_out"/>

                <var name="exner" type="real" dimensions="nVertLevels nCells Time" units="unitless"
                     description="Exner function"
                     packages="met_stage_out"/>

                <var name="exner_base" type="real" dimensions="nVertLevels nCells Time" units="unitless"
                     description="Base-state Exner function"
                     packages="met_stage_out"/>

                <var name="rtheta_base" type="real" dimensions="nVertLevels nCells Time" units="kg K m^{-3}"
                     description="reference state rho*theta/zz"
                     packages="met_stage_out"/>

                <var name="pressure" type="real" dimensions="nVertLevels nCells Time" units="Pa"
                     description="Pressure"
                     packages="met_stage_out"/>

                <var name="pressure_base" type="real" dimensions="nVertLevels nCells Time" units="Pa"
                     description="Base state pressure"
                     packages="met_stage_out"/>

                <var name="rho_base" type="real" dimensions="nVertLevels nCells Time" units="kg m^{-3}"
                     description="Base state dry air density"
                     packages="met_stage_out"/>

                <var name="theta_base" type="real" dimensions="nVertLevels nCells Time" units="K"
                     description="Base state potential temperature"
                     packages="met_stage_out"/>

                <var name="cqw" type="real" dimensions="nVertLevels nCells Time" units="unitless"
                     description="rho_d/rho_m at w points"
                     packages="met_stage_out"/>

                <var name="surface_pressure" type="real" dimensions="nCells Time" units="Pa"
                     description="Diagnosed surface pressure"
                     packages="met_stage_out"/>

                <!-- coupled variables needed by the solver, but not output -->
                <var name="ru" type="real" dimensions="nVertLevels nEdges Time" units="kg m^{-2} s^{-1}"
                     description="horizontal momentum at cell edge (rho*u/zz)"
                     packages="met_stage_out"/>

                <var name="rw" type="real" dimensions="nVertLevelsP1 nCells Time" units="kg m^{-2} s^{-1}"
                     description="rho*omega/zz carried at w points"
                     packages="met_stage_out"/>

                <var name="rtheta_p" type="real" dimensions="nVertLevels nCells Time" units="kg K m^{-3}"
                     description="rho*theta_m/zz perturbation from the reference state value"
                     packages="met_stage_out"/>

                <var name="rho_p" type="real" dimensions="nVertLevels nCells Time" units="kg m^{-3}"
                     description="rho/zz perturbation from the reference state value, advanced over acoustic steps"
                     packages="met_stage_out"/>
        </var_struct>

        <var_struct name="diag_physics" time_levs="1">
                <var name="precipw" type="real" dimensions="nCells Time" units="kg m^{-2}"
                     description="precipitable water"
                     packages="met_stage_out"/>
        </var_struct>
</registry>