source: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/context_pcm_physics.xml

<context id="LMDZ">

    <!-- Definition of a calendar is done in code via wxios_set_cal -->

    <!-- Definition of domains and groups of domains -->
    <domain_definition>
        <domain id="dom_glo" data_dim="2" />

        <domain id="dom_32_24" type="rectilinear" ni_glo="32" nj_glo="24"   >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>
        
        <domain id="dom_64_48" type="rectilinear" ni_glo="64" nj_glo="48"   >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_96_95" type="rectilinear"  ni_glo="96" nj_glo="95"  >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_144_142" type="rectilinear" ni_glo="144" nj_glo="142"   >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_512_360" type="rectilinear" ni_glo="512" nj_glo="360"   >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_720_360" type="rectilinear" ni_glo="720" nj_glo="360" >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_128_96" type="rectilinear" ni_glo="128" nj_glo="96" >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_256_192" type="rectilinear" ni_glo="256" nj_glo="192" >
          <generate_rectilinear_domain/>
          <interpolate_domain order="1"/>
        </domain>

        <domain id="dom_out" domain_ref="dom_32_24"/>
    </domain_definition>
    <!-- <domain_definition src="./domain_def.xml"/> -->

    <!-- Definition of vertical axes
                           NB: these must be defined and set in code via
         xios_output_mod -->
    <axis_definition>
        <axis id="presnivs"
              standard_name="Pseudo-pressure of model vertical levels"
              unit="Pa">
        </axis>
        <axis id="altitude"
              standard_name="Pseudo-altitude of model vertical levels"
              unit="km">
        </axis>
        <axis id="soil_layers"
              standard_name="soil layer depth at intermediate level"
              unit="m">
        </axis>
        <axis id="interlayer"
              standard_name="Levels between layers (altitude+1)"
              unit="km">
        </axis>
        <axis id="controle_axe"
              standard_name="Controle"
              n_glo="100">
        </axis>
        <!--Added for diagspecIR-->
        <axis id="IR_Wavenumber"
               standard_name="Band mid frequency in the infrared"
               unit="cm-1">
        </axis>
        <!--Added for diagspecVI-->
        <axis id="VI_Wavenumber"
              standard_name="Band mid frequency in the visible"
              unit=
            "cm-1">
        </axis>

    </axis_definition>

    <scalar_definition>
           <scalar id="scalar" />
    </scalar_definition>

    <grid_definition>
        <!-- toggle axis id below to change output vertical axis -->
        <grid id="grid_2D">
            <domain domain_ref="dom_glo" />
        </grid>
        <grid id="grid_3D">
            <domain domain_ref="dom_glo" />
            <axis axis_ref="altitude" />
        </grid>
        <grid id="grid_soil">
            <domain domain_ref="dom_glo" />
            <axis axis_ref="soil_layers" />
        </grid>       
        <!--Added for diagspecIR-->
        <grid id="grid_IR_3D">
            <domain domain_ref="dom_glo" />
            <axis axis_ref='IR_Wavenumber'/>
        </grid>
        <!-- Added for diagspecVI-->
        <grid id="grid_VI_3D">
          <domain domain_ref="dom_glo" />
          <axis axis_ref="VI_Wavenumber"/>
        </grid>
        <!-- grid for interpolated lon-lat outputs -->
        <!-- to output scalars -->
        <grid id="scalarpoint" >
            <scalar scalar_ref="scalar" />
        </grid>
        <grid id="altitude_grid" >
            <axis axis_ref="altitude" />
        </grid>
        <!-- output grids -->
        <grid id="grid_3D_out">
            <domain domain_ref="dom_out" />
            <axis axis_ref="altitude" />
        </grid>
        <grid id="grid_2D_out">
            <domain domain_ref="dom_out" />
        </grid>
        <!--added for diagspecIR-->
        <grid id="grid_IR_3D_out">
            <domain domain_ref="dom_out"/>
            <axis axis_ref="IR_Wavenumber"/>
        </grid>
        <grid id="grid_IR">
            <axis axis_ref="IR_Wavenumber"/>
        </grid>
        <!--Added for diagspecVI-->
        <grid id="grid_VI_3D_out">
            <domain domain_ref="dom_out"/>
            <axis axis_ref="VI_Wavenumber"/>
        </axis>
        <grid id="grid_VI">
          <axis axis_ref="VI_Wavenumber"/>
        </grid>
    </grid_definition>

    <!-- Definition of model variables -->
    <field_definition src="./field_def_physics.xml" />

    <!-- Definition of output files -->
    <file_definition src="./file_def_physics.xml" />
</context>