<!-- ================================================================================== -->
<!-- field_def_physics.xml                                                              -->
<!-- Definition of all existing variables that can be outputed from Martian physics     -->
<!-- ================================================================================== -->

    <field_definition prec="4"
                      freq_op="1ts"
                      enabled="true">

        <!-- 0D variables -->
        <field_group id="fields_0D" grid_ref="scalarpoint">
            <field id="Ls"
                   long_name="Solar longitude"
                   unit="degrees" />
            <field id="Sols"
                   long_name="Time"
                   unit="sols" />

            <!-- planetary escape fluxes -->
            <field id="PhiH"
                   long_name="H escape flux"
                   unit="s-1" />
            <field id="PhiH2"
                   long_name="H2 escape flux"
                   unit="s-1" />
            <field id="PhiD"
                   long_name="D escape flux"
                   unit="s-1" />

            <field id="co2conservation"
                   long_name="Total CO2 mass conservation in physic"
                   unit="kg" />

            <!-- tracer total mass (chemistry) -->
            <field id="mass_co2"
                   long_name="CO2 global mass"
                   unit="g" />
            <field id="mass_co"
                   long_name="CO global mass"
                   unit="g" />
            <field id="mass_o"
                   long_name="O global mass"
                   unit="g" />
            <field id="mass_o1d"
                   long_name="O1d global mass"
                   unit="g" />
            <field id="mass_o2"
                   long_name="O2 global mass"
                   unit="g" />
            <field id="mass_o3"
                   long_name="O3 global mass"
                   unit="g" />
            <field id="mass_h"
                   long_name="H global mass"
                   unit="g" />
            <field id="mass_h2"
                   long_name="H2 global mass"
                   unit="g" />
            <field id="mass_oh"
                   long_name="OH global mass"
                   unit="g" />
            <field id="mass_ho2"
                   long_name="HO2 global mass"
                   unit="g" />
            <field id="mass_h2o2"
                   long_name="H2O2 global mass"
                   unit="g" />
            <field id="mass_n2"
                   long_name="N2 global mass"
                   unit="g" />
            <field id="mass_ar"
                   long_name="Ar global mass"
                   unit="g" />
            <field id="mass_h2o_ice"
                   long_name="H2O_ice global mass"
                   unit="g" />
            <field id="mass_h2o_vap"
                   long_name="H2O_vap global mass"
                   unit="g" />
            <field id="mass_n"
                   long_name="N global mass"
                   unit="g" />
            <field id="mass_no"
                   long_name="NO global mass"
                   unit="g" />
            <field id="mass_no2"
                   long_name="NO2 global mass"
                   unit="g" />
            <field id="mass_n2d"
                   long_name="N2d global mass"
                   unit="g" />
            <field id="mass_ch4"
                   long_name="CH4 global mass"
                   unit="g" />
            <field id="mass_he"
                   long_name="He global mass"
                   unit="g" />

            <!-- tracer total mass (ions) -->
            <field id="mass_co2plus"
                   long_name="CO2+ global mass"
                   unit="g" />
            <field id="mass_oplus"
                   long_name="O+ global mass"
                   unit="g" />
            <field id="mass_o2plus"
                   long_name="O2+ global mass"
                   unit="g" />
            <field id="mass_coplus"
                   long_name="CO+ global mass"
                   unit="g" />
            <field id="mass_o2plus"
                   long_name="O2+ global mass"
                   unit="g" />
            <field id="mass_coplus"
                   long_name="CO+ global mass"
                   unit="g" />
            <field id="mass_cplus"
                   long_name="C+ global mass"
                   unit="g" />
            <field id="mass_nplus"
                   long_name="N+ global mass"
                   unit="g" />
            <field id="mass_noplus"
                   long_name="NO+ global mass"
                   unit="g" />
            <field id="mass_n2plus"
                   long_name="N2+ global mass"
                   unit="g" />
            <field id="mass_hplus"
                   long_name="H+ global mass"
                   unit="g" />
            <field id="mass_hco2plus"
                   long_name="HCO2+ global mass"
                   unit="g" />
            <field id="mass_hcoplus"
                   long_name="HCO+ global mass"
                   unit="g" />
            <field id="mass_h2oplus"
                   long_name="H2O+ global mass"
                   unit="g" />
            <field id="mass_h3oplus"
                   long_name="H3O+ global mass"
                   unit="g" />
            <field id="mass_ohplus"
                   long_name="OH+ global mass"
                   unit="g" />
            <field id="mass_elec"
                   long_name="electron global mass"
                   unit="g" />

        </field_group>

        <field_group id="fields_controle" grid_ref="controle_grid">
            <field id="controle"
                   long_name="controle" />
        </field_group>

        <field_group id="fields_altitude" grid_ref="altitude_grid">
            <field id="aps"
                   long_name="aps" />
            <field id="bps"
                   long_name="bps" />
        </field_group>

        <field_group id="fields_interlayer" grid_ref="interlayer_grid">
            <field id="ap"
                   long_name="ap" />
            <field id="bp"
                   long_name="bp" />
        </field_group>


        <!-- 2D variables -->
        <field_group id="fields_2D" domain_ref="dom_glo">
            <field id="area"
                   long_name="Mesh area"
                   unit="m2" />
            <field id="phisinit"
                   long_name="Surface geopotential"
                   unit="m2.s-2" />
            <field id="emis"
                   long_name="Surface emissivity"
                   unit="" />
            <field id="emis_slope01"
                   long_name="Surface emissivity of slope 01"
                   unit="" />
            <field id="emis_slope02"
                   long_name="Surface emissivity of slope 02"
                   unit="" />
            <field id="emis_slope03"
                   long_name="Surface emissivity of slope 03"
                   unit="" />
            <field id="emis_slope04"
                   long_name="Surface emissivity of slope 04"
                   unit="" />
            <field id="emis_slope05"
                   long_name="Surface emissivity of slope 05"
                   unit="" />
            <field id="emis_slope06"
                   long_name="Surface emissivity of slope 06"
                   unit="" />
            <field id="emis_slope07"
                   long_name="Surface emissivity of slope 07"
                   unit="" />
            <field id="albedo"
                   long_name="Albedo of the surface"
                   unit="" />
            <field id="albedo_slope01"
                   long_name="Albedo of the surface for slope 01"
                   unit="" />
            <field id="albedo_slope02"
                   long_name="Albedo of the surface for slope 02"
                   unit="" />
            <field id="albedo_slope03"
                   long_name="Albedo of the surface for slope 03"
                   unit="" />
            <field id="albedo_slope04"
                   long_name="Albedo of the surface for slope 04"
                   unit="" />
            <field id="albedo_slope05"
                   long_name="Albedo of the surface for slope 05"
                   unit="" />
            <field id="albedo_slope06"
                   long_name="Albedo of the surface for slope 06"
                   unit="" />
            <field id="albedo_slope07"
                   long_name="Albedo of the surface for slope 07"
                   unit="" />
            <field id="local_time"
                   long_name="Local time"
                   unit="sol" />

            <field id="ps"
                   long_name="Surface Pressure"
                   unit="Pa" />
            <field id="tsurf"
                   long_name="Surface Temperature"
                   unit="K" />
            <field id="tsurf_slope01"
                   long_name="Surface Temperature of slope 01"
                   unit="K" />
            <field id="tsurf_slope02"
                   long_name="Surface Temperature of slope 02"
                   unit="K" />
            <field id="tsurf_slope03"
                   long_name="Surface Temperature of slope 03"
                   unit="K" />
            <field id="tsurf_slope04"
                   long_name="Surface Temperature of slope 04"
                   unit="K" />
            <field id="tsurf_slope05"
                   long_name="Surface Temperature of slope 05"
                   unit="K" />
            <field id="tsurf_slope06"
                   long_name="Surface Temperature of slope 06"
                   unit="K" />
            <field id="tsurf_slope07"
                   long_name="Surface Temperature of slope 07"
                   unit="K" />
            <field id="temp_layer1"
                   long_name="Temperature in layer 1"
                   unit="K" />
            <field id="temp7"
                   long_name="Temperature in layer 7"
                   unit="K" />

            <!-- Radiation -->
            <field id="fluxsurf_lw"
                   long_name="Longwave radiation at the surface"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope01"
                   long_name="Longwave radiation at the surface on slope 01"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope02"
                   long_name="Longwave radiation at the surface on slope 02"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope03"
                   long_name="Longwave radiation at the surface on slope 03"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope04"
                   long_name="Longwave radiation at the surface on slope 04"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope05"
                   long_name="Longwave radiation at the surface on slope 05"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope06"
                   long_name="Longwave radiation at the surface on slope 06"
                   unit="W.m-2" />
            <field id="fluxsurf_lw_slope07"
                   long_name="Longwave radiation at the surface on slope 07"
                   unit="W.m-2" />

            <field id="fluxtop_lw"
                   long_name="Longwave radiation at the top of the atmosphere"
                   unit="W.m-2" />
            <field id="fluxtop_dn_sw"
                   long_name="Incoming shortwave radiation at the top of the atmosphere"
                   unit="W.m-2" />
            <field id="fluxtop_up_sw"
                   long_name="Upward shortwave radiation at the top of the atmosphere"
                   unit="W.m-2" />

            <field id="fluxsurf_dn_sw"
                   long_name="Incoming shortwave radiation at the surface"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope01"
                   long_name="Incoming shortwave radiation at the surface on slope 01"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope02"
                   long_name="Incoming shortwave radiation at the surface on slope 02"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope03"
                   long_name="Incoming shortwave radiation at the surface on slope 03"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope04"
                   long_name="Incoming shortwave radiation at the surface on slope 04"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope05"
                   long_name="Incoming shortwave radiation at the surface on slope 05"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope06"
                   long_name="Incoming shortwave radiation at the surface on slope 06"
                   unit="W.m-2" />
            <field id="fluxsurf_dn_sw_slope07"
                   long_name="Incoming shortwave radiation at the surface on slope 07"
                   unit="W.m-2" />
            <field id="fluxsurf_up_sw"
                   long_name="Reflected solar radiative flux from surface"
                   unit="W.m-2" />

            <!-- dust cycle -->
            <field id="tau_pref_scenario"
                   long_name="Prescribed visible dust optical depth at 610Pa"
                   unit="" />
            <field id="tau_pref_gcm"
                   long_name="Visible dust optical depth at 610Pa in the GCM"
                   unit="" />
            <field id="IRtoVIScoef"
                   long_name="Conversion coeff for dust tau from abs9.3um to ext0.67um"
                   unit="" />
            <field id="IRtoVIScoef"
                   long_name="Conversion coeff for dust tau from abs9.3um to ext0.67um"
                   unit="" />
            <field id="dust_rad_adjust"
                   long_name="radiative adjustment coefficient for dust"
                   unit="" />

            <field id="dqsdust"
                   long_name="deposited surface dust mass"
                   unit="kg.m-2.s-1" />
            <field id="dqndust"
                   long_name="deposited surface dust number"
                   unit="number.m-2.s-1" />
            <field id="dqsdifdustq"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope01"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope02"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope03"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope04"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope05"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope06"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifdustq_slope07"
                   long_name="tendency due to vertical diffusion of background dust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope01"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope02"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope03"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope04"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope05"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope06"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />
            <field id="dqsdifrdsq_slope07"
                   long_name="tendency due to vertical diffusion of stormdust on surface"
                   unit="kg.m-2.s-1" />


            <field id="tau_pref_target"
                   long_name="target visible dust opacity column at 610Pa"
                   unit="" />
            <field id="dtau"
                   long_name="opacity difference wrt scenario"
                   unit="" />
            <field id="dustliftday"
                   long_name="dust injection rate"
                   unit="s-1" />

            <field id="stormfract"
                   long_name="fraction of the mesh with stormdust"
                   unit="none" />
            <field id="mstormdtot"
                   long_name="total mass of stormdust in the column"
                   unit="kg.m-2" />
            <field id="mdusttot"
                   long_name="total mass of background dust in the column"
                   unit="kg.m-2" />
            <field id="rdsdqsdust"
                   long_name="deposited surface stormdust mass rate"
                   unit="kg.m-2.s-1" />
            <field id="zdqssed_dust_mass"
                   long_name="sedimentation tendency of background dust on surface"
                   unit="kg.m-2.s-1" />

            <!-- water cycle -->
            <field id="watercaptag"
                   long_name="Perennial water ice location"
                   unit="T/F" />
            <field id="watercap"
                   long_name="Perennial water ice thickness"
                   unit="kg.m-2" />
            <field id="watercap_slope01"
                   long_name="Perennial water ice thickness of slope 01"
                   unit="kg.m-2" />
            <field id="watercap_slope02"
                   long_name="Perennial water ice thickness of slope 02"
                   unit="kg/m2" />
            <field id="watercap_slope03"
                   long_name="Perennial water ice thickness of slope 03"
                   unit="kg/m2" />
            <field id="watercap_slope04"
                   long_name="Perennial water ice thickness of slope 04"
                   unit="kg/m2" />
            <field id="watercap_slope05"
                   long_name="Perennial water ice thickness of slope 05"
                   unit="kg/m2" />
            <field id="watercap_slope06"
                   long_name="Perennial water ice thickness of slope 06"
                   unit="kg/m2" />
            <field id="watercap_slope07"
                   long_name="Perennial water ice thickness of slope 07"
                   unit="kg/m2" />

            <field id="surf_h2o_lh"
                   long_name="Ground ice latent heat flux"
                   unit="W.m-2" />
            <field id="tauTES"
                   long_name="Water ice cloud opacity at TES abs 825 cm-1"
                   unit="" />
            <field id="tauTESap"
                   long_name="Water ice cloud opacity at TES abs 825 cm-1"
                   unit="" />
            <field id="mtot"
                   long_name="Total mass of water vapor in the column"
                   unit="kg.m-2" />
            <field id="icetot"
                   long_name="Total mass of water ice in the atmospheric column"
                   unit="kg.m-2" />
            <field id="reffice"
                   long_name="Mean effective radius of water ice in the atmospheric column"
                   unit="m" />
            <field id="Nccntot"
                   long_name="Number of water condensation nuclei in the atmospheric column"
                   unit="nbr.m-2" />
            <field id="Mccntot"
                   long_name="Mass of water condensation nuclei in the atmospheric column"
                   unit="kg.m-2" />
            <field id="h2o_ice_s"
                   long_name="Mass of water ice on the surface"
                   unit="kg.m-2" />
            <field id="h2o_ice_s"
                   long_name="Mass of water ice on the surface"
                   unit="kg.m-2" />
            <field id="h2o_ice_s_slope01"
                   long_name="Mass of water ice on the surface of slope 01"
                   unit="kg.m-2" />
            <field id="h2o_ice_s_slope02"
                   long_name="Mass of water ice on the surface of slope 02"
                   unit="kg/m2" />
            <field id="h2o_ice_s_slope03"
                   long_name="Mass of water ice on the surface of slope 03"
                   unit="kg/m2" />
            <field id="h2o_ice_s_slope04"
                   long_name="Mass of water ice on the surface of slope 04"
                   unit="kg/m2" />
            <field id="h2o_ice_s_slope05"
                   long_name="Mass of water ice on the surface of slope 05"
                   unit="kg/m2" />
            <field id="h2o_ice_s_slope06"
                   long_name="Mass of water ice on the surface of slope 06"
                   unit="kg/m2" />
            <field id="h2o_ice_s_slope07"
                   long_name="Mass of water ice on the surface of slope 07"
                   unit="kg/m2" />
            <field id="waterdensity_surface"
                   long_name="waterdensity_surface"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope01"
                   long_name="waterdensity_surface of slope 01"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope02"
                   long_name="waterdensity_surface of slope 02"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope03"
                   long_name="waterdensity_surface of slope 03"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope04"
                   long_name="waterdensity_surface of slope 04"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope05"
                   long_name="waterdensity_surface of slope 05"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope06"
                   long_name="waterdensity_surface of slope 06"
                   unit="kg.m-3" />
            <field id="waterdensity_surface_slope07"
                   long_name="waterdensity_surface of slope 07"
                   unit="kg.m-3" />
            <field id="h2o_layer1"
                   long_name="h2o in the first layer"
                   unit="kg/kg" />
            <field id="co2_layer1"
                   long_name="co2 in the first layer"
                   unit="kg/kg" />

            <!-- Subsurface ice interactions -->
            <field id="zdqsdif_ssi_frost_tot"
                   long_name="Flux between frost and subsurface ice"
                   unit="kg.m-2.s-1" />
            <field id="zdqsdif_ssi_atm_tot"
                   long_name="Flux between atmosphere and subsurface ice"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice"
                   long_name="Total flux exchanged with subsurface ice"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope01"
                   long_name="Total flux exchanged with subsurface ice of slope 01"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope02"
                   long_name="Total flux exchanged with subsurface ice of slope 02"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope03"
                   long_name="Total flux exchanged with subsurface ice of slope 03"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope04"
                   long_name="Total flux exchanged with subsurface ice of slope 04"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope05"
                   long_name="Total flux exchanged with subsurface ice of slope 05"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope06"
                   long_name="Total flux exchanged with subsurface ice of slope 06"
                   unit="kg.m-2.s-1" />
             <field id="flux_ssice_slope07"
                   long_name="Total flux exchanged with subsurface ice of slope 07"
                    unit="kg.m-2.s-1" />

            <!-- Water ice sublimation parameters -->
            <field id="rib_dry_vdif_cd"
                   long_name="Dry Richardson number in vdif_cd"
                   unit="1" />
            <field id="rib_vdif_cd"
                   long_name="Richardson number in vdif_cd"
                   unit="1" />
            <field id="fm_vdif_cd"
                   long_name="Momentum stability function  in vdif_cd"
                   unit="1" />
            <field id="fh_vdif_cd"
                   long_name="Heat stability function  in vdif_cd"
                   unit="1" />
            <field id="z0t_vdif_cd"
                   long_name="Thermal roughness length  in vdif_cd"
                   unit="1" />
            <field id="z0_vdif_cd"
                   long_name="Momentum roughness length  in vdif_cd"
                   unit="1" />
            <field id="Reynolds_vdif_cd"
                   long_name="Reynolds number in vdif_cd"
                   unit="1" />

            <!-- CO2 cycle -->
            <field id="co2ice"
                   long_name="CO2 ice thickness"
                   unit="kg.m-2" />
            <field id="co2ice_slope01"
                   long_name="CO2 ice thickness on slope 01"
                   unit="kg.m-2" />
            <field id="co2ice_slope02"
                   long_name="CO2 ice thickness on slope 02"
                   unit="kg/m2" />
            <field id="co2ice_slope03"
                   long_name="CO2 ice thickness on slope 03"
                   unit="kg/m2" />
            <field id="co2ice_slope04"
                   long_name="CO2 ice thickness on slope 04"
                   unit="kg/m2" />
            <field id="co2ice_slope05"
                   long_name="CO2 ice thickness on slope 05"
                   unit="kg/m2" />
            <field id="co2ice_slope06"
                   long_name="CO2 ice thickness on slope 06"
                   unit="kg/m2" />
            <field id="co2ice_slope07"
                   long_name="CO2 ice thickness on slope 07"
                   unit="kg/m2" />

            <field id="perennial_co2ice"
                   long_name="Perennial co2 ice thickness"
                   unit="kg.m-2" />
            <field id="perennial_co2ice_slope01"
                   long_name="Perennial co2 ice thickness of slope 01"
                   unit="kg.m-2" />
            <field id="perennial_co2ice_slope02"
                   long_name="Perennial co2 ice thickness of slope 02"
                   unit="kg/m2" />
            <field id="perennial_co2ice_slope03"
                   long_name="Perennial co2 ice thickness of slope 03"
                   unit="kg/m2" />
            <field id="perennial_co2ice_slope04"
                   long_name="Perennial co2 ice thickness of slope 04"
                   unit="kg/m2" />
            <field id="perennial_co2ice_slope05"
                   long_name="Perennial co2 ice thickness of slope 05"
                   unit="kg/m2" />
            <field id="perennial_co2ice_slope06"
                   long_name="Perennial co2 ice thickness of slope 06"
                   unit="kg/m2" />
            <field id="perennial_co2ice_slope07"
                   long_name="Perennial co2 ice thickness of slope 07"
                   unit="kg/m2" />

            <field id="co2condens_zfallice"
                   long_name="CO2 ice tendency on the surface"
                   unit="kg.m-2.s-1" />

            <!-- tracer columns (chemistry) -->
            <field id="c_co2"
                   long_name="CO2 column"
                   unit="mol.cm-2" />
            <field id="c_co"
                   long_name="CO column"
                   unit="mol.cm-2" />
            <field id="c_o"
                   long_name="O column"
                   unit="mol.cm-2" />
            <field id="c_o1d"
                   long_name="O1d column"
                   unit="mol.cm-2" />
            <field id="c_o2"
                   long_name="O2 column"
                   unit="mol.cm-2" />
            <field id="c_o3"
                   long_name="O3 column"
                   unit="mol.cm-2" />
            <field id="c_h"
                   long_name="H column"
                   unit="mol.cm-2" />
            <field id="c_h2"
                   long_name="H2 column"
                   unit="mol.cm-2" />
            <field id="c_oh"
                   long_name="OH column"
                   unit="mol.cm-2" />
            <field id="c_ho2"
                   long_name="HO2 column"
                   unit="mol.cm-2" />
            <field id="c_h2o2"
                   long_name="H2 column"
                   unit="mol.cm-2" />
            <field id="c_n2"
                   long_name="N2 column"
                   unit="mol.cm-2" />
            <field id="c_ar"
                   long_name="Ar column"
                   unit="mol.cm-2" />
            <field id="c_h2o_ice"
                   long_name="H2O_ice column"
                   unit="mol.cm-2" />
            <field id="c_h2o_vap"
                   long_name="H2O_vap column"
                   unit="mol.cm-2" />
            <field id="c_n"
                   long_name="N column"
                   unit="mol.cm-2" />
            <field id="c_no"
                   long_name="NO column"
                   unit="mol.cm-2" />
            <field id="c_no2"
                   long_name="NO2 column"
                   unit="mol.cm-2" />
            <field id="c_n2d"
                   long_name="N2d column"
                   unit="mol.cm-2" />
            <field id="c_ch4"
                   long_name="CH4 column"
                   unit="mol.cm-2" />
            <field id="c_he"
                   long_name="He column"
                   unit="mol.cm-2" />

            <!-- ions -->
            <field id="c_co2plus"
                   long_name="CO2+ column"
                   unit="mol.cm-2" />
            <field id="c_oplus"
                   long_name="O+ column"
                   unit="mol.cm-2" />
            <field id="c_o2plus"
                   long_name="O2+ column"
                   unit="mol.cm-2" />
            <field id="c_coplus"
                   long_name="CO+ column"
                   unit="mol.cm-2" />
            <field id="c_coplus"
                   long_name="CO+ column"
                   unit="mol.cm-2" />
            <field id="c_cplus"
                   long_name="C+ column"
                   unit="mol.cm-2" />
            <field id="c_nplus"
                   long_name="N+ column"
                   unit="mol.cm-2" />
            <field id="c_noplus"
                   long_name="NO+ column"
                   unit="mol.cm-2" />
            <field id="c_n2plus"
                   long_name="N2+ column"
                   unit="mol.cm-2" />
            <field id="c_hplus"
                   long_name="H+ column"
                   unit="mol.cm-2" />
            <field id="c_hco2plus"
                   long_name="HCO2+ column"
                   unit="mol.cm-2" />
            <field id="c_hcoplus"
                   long_name="HCO+ column"
                   unit="mol.cm-2" />
            <field id="c_h2oplus"
                   long_name="H2O+ column"
                   unit="mol.cm-2" />
            <field id="c_h3oplus"
                   long_name="H3O+ column"
                   unit="mol.cm-2" />
            <field id="c_ohplus"
                   long_name="OH+ column"
                   unit="mol.cm-2" />
            <field id="c_elec"
                   long_name="electron column"
                   unit="mol.cm-2" />

            <!-- Thermals parametrization -->
            <field id="lmax_th"
                   long_name="index of vertical extension of thermals"
                   unit="grid level" />
            <field id="zmax_th"
                   long_name="vertical extension of thermals"
                   unit="m" />
            <field id="hfmax_th"
                   long_name="maximum heat flux in thermals"
                   unit="K.m/s" />
            <field id="wstar"
                   long_name="maximum vertical velocity in thermals"
                   unit="m/s" />

            <!-- Subsurface tracers  (adsorption) -->

            <field id="flux_rego"
                   long_name="flux of water from the regolith"
                   unit="kg/m^2" />
            <field id="mass_h2o_soil"
                   long_name="Mass of subsurface water column at each point"
                   unit="kg m-2" />
            <field id="mass_ice_soil"
                   long_name="Mass of subsurface ice at each point"
                   unit="kg m-2" />
            <field id="nsurf"
                   long_name="surface water vapor density"
                   unit="kg/m^3" />


        </field_group>

        <!-- 3D variables -->
        <field_group id="fields_3D" domain_ref="dom_glo" axis_ref="altitude">
        <!--<field_group id="fields_3D" grid_ref="grid_3d">-->

            <field id="zzlay"
                   long_name="Midlayer altitude"
                   unit="m" />
            <field id="zzlev"
                   long_name="Interlayer altitude"
                   unit="m" />
            <field id="pphi"
                   long_name="Geopotential"
                   unit="m2.s-2" />

            <!-- main meteorological fields -->
            <field id="temp"
                   long_name="Atmospheric temperature"
                   unit="K" />
            <field id="u"
                   long_name="Eastward Zonal Wind"
                   unit="m.s-1" />
            <field id="v"
                   long_name="Northward Meridional Wind"
                   unit="m.s-1" />
            <field id="w"
                   long_name="Vertical Wind (positive when downward)"
                   unit="m.s-1" />
            <field id="rho"
                   long_name="Atmospheric density"
                   unit="kg.m-3" />
            <field id="pressure"
                   long_name="Atmospheric pressure"
                   unit="Pa" />
            <field id="zplev"
                   long_name="Interlayer pressure"
                   unit="Pa" />

            <!-- heating rates -->
            <field id="dtrad"
                   long_name="total radiative heating rate"
                   unit="K.s-1" />
            <field id="sw_htrt"
                   long_name="Shortwave heating rate"
                   unit="K.s-1" />
            <field id="lw_htrt"
                   long_name="Longwave heating rate"
                   unit="K.s-1" />
            <field id="q15um"
                   long_name="15 um cooling"
                   unit="K.s-1" />
            <field id="qnir"
                   long_name="NIR heating rate"
                   unit="K.s-1" />
            <field id="quv"
                   long_name="EUV heating rate"
                   unit="K.s-1" />
            <field id="cond"
                   long_name="Thermal conduction heating rate"
                   unit="K.s-1" />

            <!-- tracers (CO2 and water) -->
            <field id="co2"
                   long_name="CO2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_co2"
                   long_name="CO2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_co2"
                   long_name="CO2 number density"
                   unit="cm-3" />
            <field id="h2o_ice"
                   long_name="water ice mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_h2oice"
                   long_name="water ice volume mixing ratio"
                   unit="mol/mol" />
            <field id="vmr_h2o_ice"
                   long_name="water ice volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_h2o_ice"
                   long_name="water ice number density"
                   unit="cm-3" />
            <field id="h2o_vap"
                   long_name="water vapor mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_h2ovap"
                   long_name="water vapor volume mixing ratio"
                   unit="mol/mol" />
            <field id="vmr_h2o_vap"
                   long_name="water vapor volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_h2o_vap"
                   long_name="water vapor number density"
                   unit="cm-3" />

            <!-- tracers (chemistry) -->
            <field id="co"
                   long_name="CO mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_co"
                   long_name="CO volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_co"
                   long_name="CO number density"
                   unit="cm-3" />
            <field id="o"
                   long_name="O mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_o"
                   long_name="O volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_o"
                   long_name="O number density"
                   unit="cm-3" />
            <field id="o1d"
                   long_name="O1d mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_o1d"
                   long_name="O1d volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_o1d"
                   long_name="O1d number density"
                   unit="cm-3" />
            <field id="o2"
                   long_name="O2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_o2"
                   long_name="O2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_o2"
                   long_name="O2 number density"
                   unit="cm-3" />
            <field id="o3"
                   long_name="O3 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_o3"
                   long_name="O3 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_o3"
                   long_name="O3 number density"
                   unit="cm-3" />
            <field id="h"
                   long_name="H mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_h"
                   long_name="H volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_h"
                   long_name="H number density"
                   unit="cm-3" />
            <field id="h2"
                   long_name="H2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_h2"
                   long_name="H2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_h2"
                   long_name="H2 number density"
                   unit="cm-3" />
            <field id="oh"
                   long_name="OH mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_oh"
                   long_name="OH volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_oh"
                   long_name="OH number density"
                   unit="cm-3" />
            <field id="ho2"
                   long_name="HO2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_ho2"
                   long_name="HO2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_ho2"
                   long_name="HO2 number density"
                   unit="cm-3" />
            <field id="h2o2"
                   long_name="H2O2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_h2o2"
                   long_name="H2O2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_h2o2"
                   long_name="H2O2 number density"
                   unit="cm-3" />
            <field id="n2"
                   long_name="N2 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_n2"
                   long_name="N2 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_n2"
                   long_name="N2 number density"
                   unit="cm-3" />
            <field id="ar"
                   long_name="Ar mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_ar"
                   long_name="Ar volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_ar"
                   long_name="Ar number density"
                   unit="cm-3" />
            <field id="num_n"
                   long_name="N number density"
                   unit="cm-3" />
            <field id="num_no"
                   long_name="NO number density"
                   unit="cm-3" />
            <field id="num_no2"
                   long_name="NO2 number density"
                   unit="cm-3" />
            <field id="num_n2d"
                   long_name="N2d number density"
                   unit="cm-3" />
            <field id="ch4"
                   long_name="CH4 mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_ch4"
                   long_name="CH4 volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_ch4"
                   long_name="CH4 number density"
                   unit="cm-3" />
            <field id="he"
                   long_name="He mass mixing ratio"
                   unit="kg/kg" />
            <field id="vmr_he"
                   long_name="He volume mixing ratio"
                   unit="mol/mol" />
            <field id="num_he"
                   long_name="He number density"
                   unit="cm-3" />

            <!-- ions -->
            <field id="num_co2plus"
                   long_name="CO2+ number density"
                   unit="cm-3" />
            <field id="num_oplus"
                   long_name="O+ number density"
                   unit="cm-3" />
            <field id="num_o2plus"
                   long_name="O2+ number density"
                   unit="cm-3" />
            <field id="num_coplus"
                   long_name="CO+ number density"
                   unit="cm-3" />
            <field id="num_coplus"
                   long_name="CO+ number density"
                   unit="cm-3" />
            <field id="num_cplus"
                   long_name="C+ number density"
                   unit="cm-3" />
            <field id="num_nplus"
                   long_name="N+ number density"
                   unit="cm-3" />
            <field id="num_noplus"
                   long_name="NO+ number density"
                   unit="cm-3" />
            <field id="num_n2plus"
                   long_name="N2+ number density"
                   unit="cm-3" />
            <field id="num_hplus"
                   long_name="H+ number density"
                   unit="cm-3" />
            <field id="num_hco2plus"
                   long_name="HCO2+ number density"
                   unit="cm-3" />
            <field id="num_hcoplus"
                   long_name="HCO+ number density"
                   unit="cm-3" />
            <field id="num_h2oplus"
                   long_name="H2O+ number density"
                   unit="cm-3" />
            <field id="num_h3oplus"
                   long_name="H3O+ number density"
                   unit="cm-3" />
            <field id="num_ohplus"
                   long_name="OH+ number density"
                   unit="cm-3" />
            <field id="num_elec"
                   long_name="electron number density"
                   unit="cm-3" />

            <!-- Dust cycle parametrization -->
            <field id="reffdust"
                   long_name="Dust effective radius"
                   unit="m" />
            <field id="rdust"
                   long_name="Dust radius"
                   unit="m" />
            <field id="dustq"
                   long_name="Dust mass mixing ratio"
                   unit="kg/kg" />
            <field id="dustN"
                   long_name="Dust number mixing ratio"
                   unit="part/kg" />
            <field id="dsodust_TES"
                   long_name="Density scaled extinction opacity of std dust at 9.3um(TES)"
                   unit="m2.kg-1" />
            <field id="dso_TES"
                   long_name="Density scaled extinction opacity of all dust at 9.3um(TES)"
                   unit="m2.kg-1" />
            <field id="dsodust"
                   long_name="Density scaled extinction opacity of std dust at 21.6um(MCS)"
                   unit="m2.kg-1" />
            <field id="dso"
                   long_name="Density scaled extinction opacity of all dust at 21.6um(MCS)"
                   unit="m2.kg-1" />
            <field id="aerosol_dust"
                   long_name="Opacity of background dust"
                   unit="" />

            <field id="reffstormdust"
                   long_name="Stormdust dust effective radius"
                   unit="m" />
            <field id="rstormdust"
                   long_name="Stormdust dust radius"
                   unit="m" />
            <field id="dsords_TES"
                   long_name="Density scaled extinction opacity of stormdust at 9.3um(TES)"
                   unit="m2.kg-1" />
            <field id="dsords"
                   long_name="Density scaled extinction opacity of stormdust at 21.6um(MCS)"
                   unit="m2.kg-1" />
            <field id="rdsdustq"
                   long_name="Stormdust dust mass mixing ratio"
                   unit="kg/kg" />
            <field id="rdsdustqmodel"
                   long_name="Stormdust dust mass mixing ratio"
                   unit="kg/kg" />
            <field id="rdsdustN"
                   long_name="Stormdust dust number mixing ratio"
                   unit="part/kg" />
            <field id="rdsdustN"
                   long_name="Stormdust dust number mixing ratio"
                   unit="part/kg" />
            <field id="wspeed_stormdust"
                   long_name="vertical velocity of stormdust"
                   unit="m/s" />
            <field id="zdqsed_dust_mass"
                   long_name="sedimentation tendency of background dust mmr"
                   unit="kg/kg.s-1" />
            <field id="zdqsed_dust_number"
                   long_name="sedimentation tendency of background dust number"
                   unit="nbr/kg.s-1" />
            <field id="zdqsed_stormdust_mass"
                   long_name="sedimentation tendency of stormdust dust mmr"
                   unit="kg/kg.s-1" />
            <field id="aerosol_stormdust"
                   long_name="Opacity of stormdust dust"
                   unit="" />
            <field id="rds_lapserate"
                   long_name="lapse rate in the rocket dust storm"
                   unit="K/m" />
            <field id="rds_deltahr"
                   long_name="extra heating rate in the rocket dust storm"
                   unit="K/s" />

            <field id="topdustq"
                   long_name="Topdust dust mass mixing ratio"
                   unit="kg/kg" />
            <field id="topdustN"
                   long_name="Topdust dust number mixing ratio"
                   unit="part/kg" />
            <field id="refftopdust"
                   long_name="Topdust dust effective radius"
                   unit="m" />
            <field id="dsotop_TES"
                   long_name="Density scaled extinction opacity of topdust at 9.3um(TES)"
                   unit="m2.kg-1" />
            <field id="dsotop"
                   long_name="Density scaled extinction opacity of topdust at 21.6um(MCS)"
                   unit="m2.kg-1" />


            <!-- Water cycle parametrization -->
            <field id="rice"
                   long_name="Water ice particle size"
                   unit="m" />
            <field id="h2o_saturation"
                   long_name="h2o vapor saturation ratio"
                   unit="" />
            <field id="watercloud_pdqh2oice"
                   long_name="pdqcloud_h2o_ice after microphysics"
                   unit="kg/kg.s-1" />
            <field id="watercloud_pdqh2ovap"
                   long_name="pdqcloud_h2o_vap after microphysics"
                   unit="kg/kg.s-1" />
            <field id="watercloud_pdqhdoice"
                   long_name="pdqcloud_hdo_ice after microphysics"
                   unit="kg/kg.s-1" />
            <field id="watercloud_pdqhdovap"
                   long_name="pdqcloud_hdo_vap after microphysics"
                   unit="kg/kg.s-1" />

            <field id="zpotcond_inst"
                   long_name="zpotcond_inst microphysics"
                   unit="(kg/kg)" />
            <field id="zpotcond_full"
                   long_name="zpotcond_full microphysics"
                   unit="(kg/kg)" />
            <field id="zpotcond"
                   long_name="zpotcond microphysics"
                   unit="(kg/kg)" />
            <field id="count_micro"
                   long_name="count_micro after microphysics"
                   unit="integer" />

            <field id="zdqcloud_ice"
                   long_name="cloud ice"
                   unit="kg.m-2.s-1" />
            <field id="zdqcloud_vap"
                   long_name="cloud vap"
                   unit="kg.m-2.s-1" />
            <field id="zdqcloud"
                   long_name="cloud"
                   unit="kg.m-2.s-1" />

            <!-- CO2 condensation parametrization -->
            <field id="co2condens_pdtc"
                   long_name="Temperature tendency due to CO2 condensation"
                   unit="K.s-1" />

            <!-- Non-orographic GW parametrization -->
            <field id="nonoro_bv"
                   long_name="Brunt Vaisala frequency in nonoro"
                   unit="Hz" />
            <field id="nonoro_u_epflux_tot"
                   long_name="Total EP Flux along U in nonoro"
                   unit="" />
            <field id="nonoro_v_epflux_tot"
                   long_name="Total EP Flux along V in nonoro"
                   unit="" />
            <field id="du_nonoro_gwd"
                   long_name="Tendency on U due to nonoro GW"
                   unit="m.s-2" />
            <field id="dv_nonoro_gwd"
                   long_name="Tendency on V due to nonoro GW"
                   unit="m.s-2" />

        </field_group>

        <!-- soil variables -->
        <field_group id="fields_soil" domain_ref="dom_glo" axis_ref="soil_layers">
            <field id="soiltemp"
                   long_name="Soil temperature"
                   unit="K" />
            <field id="soiltemp_slope01"
                   long_name="Soil temperature for slope 01"
                   unit="K" />
            <field id="soiltemp_slope02"
                   long_name="Soil temperature for slope 02"
                   unit="K" />
            <field id="soiltemp_slope03"
                   long_name="Soil temperature for slope 03"
                   unit="K" />
            <field id="soiltemp_slope04"
                   long_name="Soil temperature for slope 04"
                   unit="K" />
            <field id="soiltemp_slope05"
                   long_name="Soil temperature for slope 05"
                   unit="K" />
            <field id="soiltemp_slope06"
                   long_name="Soil temperature for slope 06"
                   unit="K" />
            <field id="soiltemp_slope07"
                   long_name="Soil temperature for slope 07"
                   unit="K" />
            <field id="inertiedat"
                   long_name="Soil thermal inertia"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope01"
                   long_name="Soil thermal inertia for slope 01"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope02"
                   long_name="Soil thermal inertia for slope 02"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope03"
                   long_name="Soil thermal inertia for slope 03"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope04"
                   long_name="Soil thermal inertia for slope 04"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope05"
                   long_name="Soil thermal inertia for slope 05"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope06"
                   long_name="Soil thermal inertia for slope 06"
                   unit="J/kg/K" />
            <field id="inertiesoil_slope07"
                   long_name="Soil thermal inertia for slope 07"
                   unit="J/kg/K" />
            <field id="waterdensity_soil"
                   long_name="waterdensity_soil"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope01"
                   long_name="waterdensity_soil of slope 01"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope02"
                   long_name="waterdensity_soil of slope 02"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope03"
                   long_name="waterdensity_soil of slope 03"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope04"
                   long_name="waterdensity_soil of slope 04"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope05"
                   long_name="waterdensity_soil of slope 05"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope06"
                   long_name="waterdensity_soil of slope 06"
                   unit="kg.m-3" />
            <field id="waterdensity_soil_slope07"
                   long_name="waterdensity_soil of slope 07"
                   unit="kg.m-3" />
            <field id="nsubtimestep"
                   long_name="vdifc substimestep number"
                   unit="s" />
            <field id="zdq_subtimestep"
                   long_name="Actual flux zdqsdif*subtimestep"
                   unit="kg.m-2" />
            <field id="zdq_end"
                   long_name="Flux after all contributions"
                   unit="kg.m-2.s-1" />
            <field id="flux_soillayer"
                   long_name="flux of water between the soil layers"
                   unit="kg.m-2.s-1" />
            <field id="ice_saturation_soil"
                   long_name="Water ice saturation in the soil layers"
                   unit="Percent" />
            <field id="znsoil"
                   long_name="Water vapor soil concentration"
                   unit="kg m-3 of pore air" />
            <field id="nsatsoil"
                   long_name="subsurface water vapor saturation density"
                   unit="kg/m^3" />
            <field id="adswater"
                   long_name="subsurface adsorbed water"
                   unit="kg/m^3" />
            <field id="coeff_diffusion_soil"
                   long_name="interlayer diffusion coefficient"
                   unit="m^2/s" />

        </field_group>

    </field_definition>
