<!-- Context LMDZ -->
<context id="LMDZ">
<!--  <calendar type="D360" start_date="1980-01-01 00:00:00" /> -->

  <!-- Define available variables -->
  <field_definition src="./field_def_lmdz.xml"/>

  <!-- Define output files
              Each file contains the list of variables and their output levels -->

  <file_definition src="./file_def_histday_lmdz.xml"/>
  <file_definition src="./file_def_histhf_lmdz.xml"/>
  <file_definition src="./file_def_histins_lmdz.xml"/>
<!--  <file_definition src="./file_def_histLES_lmdz.xml"/> -->
  <file_definition src="./file_def_histmth_lmdz.xml"/>
<!--  <file_definition src="./file_def_histstn_lmdz.xml"/> -->
  <file_definition src="./file_def_histmthNMC_lmdz.xml"/>
  <file_definition src="./file_def_histdayNMC_lmdz.xml"/>
  <file_definition src="./file_def_histhfNMC_lmdz.xml"/>
<!--  <file_definition src="./file_def_histmthCOSP_lmdz.xml"/> -->
<!--  <file_definition src="./file_def_histdayCOSP_lmdz.xml"/>  -->
<!--  <file_definition src="./file_def_histhfCOSP_lmdz.xml"/> --> 
  <file_definition src="./file_def_histstrataer_lmdz.xml"/>
 




  <field_definition>

    <field_group id="etat0_limit_read" read_access="false" enabled="false" operation="instant">
       <field id="ts"           field_ref="ts_reg"      domain_ref="from_ecphy"     />
       <field id="qs"           field_ref="qs_reg"      domain_ref="from_ecphy"     />
       <field id="landice"      field_ref="landice_reg" domain_ref="from_landiceref"  />
       <field id="mask"         field_ref="mask_reg"    domain_ref="from_relief" />
       <field id="sic_limit"    field_ref="sic_reg"     domain_ref="from_limit_amip"  />
       <field id="sst_limit"    field_ref="sst_reg"     domain_ref="from_limit_amip" />
       <field id="rugos_limit"  field_ref="rugos_reg"   domain_ref="from_limit" > exp(this) </field>
       <field id="albedo_limit" field_ref="albedo_reg"  domain_ref="from_limit"  />

       <field id="zmea" field_ref="zmea_reg"  domain_ref="from_orography"  />
       <field id="zstd" field_ref="zstd_reg"  domain_ref="from_orography"  />
       <field id="zsig" field_ref="zsig_reg"  domain_ref="from_orography"  />
       <field id="zgam" field_ref="zgam_reg"  domain_ref="from_orography"  />
       <field id="zthe" field_ref="zthe_reg"  domain_ref="from_orography"  />

    </field_group>
   
   
    <field_group id="ozone_read" operation="instant " >
      <field id="tro3_read"             field_ref="tro3_reg"             domain_ref="from_ozone" axis_ref="axis_ozone"  enabled="false"/>
      <field id="tro3_daylight_read"    field_ref="tro3_daylight_reg"    domain_ref="from_ozone" axis_ref="axis_ozone"  enabled="false" />
    </field_group>


    <field id="tro3_in"   grid_ref="grid_climoz" enabled="false"/>
    <field id="tro3_out"  read_access="true"   field_ref="tro3_in"  grid_ref="grid_from_climoz" enabled="false"/>
 
  </field_definition>





  <file_definition>
    
    <file id="test_o3" output_freq="1ts" type="one_file" enabled="false">
      <field id="o3_out" grid_ref="grid_from_climoz" operation="instant" />
     </file>
      
      
    <file_group id="etat0_limit_read" type="one_file" mode="read" enabled=".FALSE." output_level="10" output_freq="1y" >
      
      <file id="ECPHY"  name="ECPHY" >
        <field id="ts_reg"   name="ST"      domain_ref="domain_ecphy" operation="instant"  freq_offset="1ts"/>
        <field id="qs_reg"   name="CDSW"    domain_ref="domain_ecphy" operation="instant"  freq_offset="1ts"/>
      </file>

      <file id="landiceref"  name="landiceref" >
        <field id="landice_reg"   name="landice"   domain_ref="domain_landiceref" operation="instant" freq_offset="1ts"/>
      </file>

      <file id="relief" name="Relief" >
        <field id="mask_reg" name="RELIEF" domain_ref="domain_relief" operation="instant" freq_offset="1ts"> this &gt; 1e-5 </field>
      </file>
 
      <file id="amipbc_sic_1x1" name="amipbc_sic_1x1" >
        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_month" operation="instant" freq_offset="1ts"/>
      </file>
 
      <file id="amipbc_sst_1x1" name="amipbc_sst_1x1" >
         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_month" operation="instant" freq_offset="1ts"/> 
      </file>

      <file id="rugos" name="Rugos" >
        <field id="rugos_reg" name="RUGOS" domain_ref="domain_limit" axis_ref="time_month" operation="instant" freq_offset="1ts"> log(this) </field> 
      </file>

      <file id="albedo" name="Albedo" >
        <field id="albedo_reg" name="ALBEDO" domain_ref="domain_limit" axis_ref="time_month" operation="instant" freq_offset="1ts"/> 
      </file>
     
      <file id="orography" name="orography_2_00" >

         <field_group domain_ref="domain_orography" operation="instant" freq_offset="1ts"> 
            <field id="zmea_reg" name="phi" />
            <field id="zstd_reg" name="mu" />
            <field id="zsig_reg" name="sigma" />
            <field id="zgam_reg" name="gama" />
            <field id="zthe_reg" name="theta" />
         </field_group >

      </file>

     </file_group>
     

     <file id="limit_write" name="limit" type="one_file" output_freq="1ts" sync_freq="1ts" enabled="false">         
       <field_group domain_ref="dom_glo" axis_ref="time_year" operation="once">
          <field id="foce_limout"  name="FOCE" long_name="Fraction ocean"/>
          <field id="fsic_limout"  name="FSIC" long_name="Fraction glace de mer"/>
          <field id="fter_limout"  name="FTER" long_name="Fraction terre"/>
          <field id="flic_limout"  name="FLIC" long_name="Fraction land ice"/>
          <field id="sst_limout"   name="SST"  long_name="Temperature superficielle de la mer"/>
          <field id="bils_limout"  name="BILS" long_name="Reference flux de chaleur au sol"/> 
          <field id="alb_limout"   name="ALB"  long_name="Albedo a la surface"/>
          <field id="rug_limout"   name="RUG"  long_name="Rugosite" />
       </field_group>
     </file>

     <file id="limit_read" name="limit" time_counter_name="time_year" type="one_file" output_freq="1d" mode="read" cyclic="true" enabled="false">         
       <field_group domain_ref="dom_glo" freq_offset="1ts" operation="instant">
          <field id="foce_limin"  name="FOCE" long_name="Fraction ocean"/>
          <field id="fsic_limin"  name="FSIC" long_name="Fraction glace de mer"/>
          <field id="fter_limin"  name="FTER" long_name="Fraction terre"/>
          <field id="flic_limin"  name="FLIC" long_name="Fraction land ice"/>
          <field id="sst_limin"   name="SST"  long_name="Temperature superficielle de la mer"/>
          <field id="bils_limin"  name="BILS" long_name="Reference flux de chaleur au sol" enabled="false"/> 
          <field id="alb_limin"   name="ALB"  long_name="Albedo a la surface"/>
          <field id="rug_limin"   name="RUG"  long_name="Rugosite" />
       </field_group>
     </file>
     

     <file id="climoz_LMDZ" type="one_file" output_freq="1d" mode="read" >
        <field id="tro3_reg" name="tro3" domain_ref="domain_ozone" axis_ref="axis_ozone" operation="instant" freq_offset="1ts" enabled="false"/> 
        <field id="tro3_daylight_reg" name="tro3_daylight" domain_ref="domain_ozone" axis_ref="axis_ozone" operation="instant" freq_offset="1ts" enabled="false"/> 
     </file>         
    
  </file_definition> 






  <!-- Define domains and groups of domains -->
  <domain_definition>

    <domain_group id="dom_glo" data_dim="1" >
  
      <domain id="dom_glo"  />
      <domain id="greordered"  domain_ref="dom_glo">
        <reorder_domain invert_lat="true" shift_lon_fraction="0.5" min_lon="0" max_lon="360" />
      </domain> 
    
      <domain id="from_ecphy" domain_ref="dom_glo">
        <interpolate_domain write_weight="true" mode="read_or_compute"/>
      </domain>  

      <domain id="from_landiceref" domain_ref="dom_glo">
        <interpolate_domain write_weight="true" mode="read_or_compute"/>
      </domain>

      <domain id="from_relief" domain_ref="dom_glo">
       <interpolate_domain write_weight="true" mode="read_or_compute" />
      </domain>

      <domain id="from_orography" domain_ref="dom_glo">
        <interpolate_domain order="1" write_weight="true" mode="read_or_compute"/>
      </domain>

      <domain id="from_limit" domain_ref="dom_glo">
        <interpolate_domain write_weight="true" mode="read_or_compute"/>
      </domain>

      <domain id="from_limit_amip"domain_ref="dom_glo" >
        <interpolate_domain write_weight="true" mode="read_or_compute"/>
      </domain>

      <domain id="from_ozone" domain_ref="dom_glo">
        <interpolate_domain order="1" write_weight="true" mode="read_or_compute" />
      </domain>

      <domain id="from_climoz" domain_ref="dom_glo">
        <interpolate_domain order="1" write_weight="true" mode="read_or_compute"/>
      </domain>
   
    </domain_group>
    
    
    <domain id="dom_regular" ni_glo="144" nj_glo="142" type="rectilinear"  >
      <generate_rectilinear_domain/>
      <interpolate_domain order="1"/>
    </domain>

    <domain id="domain_ecphy"  type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>

    <domain id="domain_landiceref"  type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>

    <domain id="domain_relief"  type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>
 
    <domain id="domain_limit"  type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>

    <domain id="domain_limit_amip"  type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>

    <domain id="domain_ozone" type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>

    <domain id="domain_orography" type="rectilinear" >
       <generate_rectilinear_domain /> 
    </domain>
    
    <domain id="dom_out" domain_ref="dom_regular"/>

    <domain id="domain_climoz" type="rectilinear" />
    
  </domain_definition>
  
  <!-- Need to define a grid? Do it here -->
<!--  <grid_definition>
    <grid_group id="vertical" axis_ref="presnivs" />
  </grid_> -->

<!-- Define Scalar grid for GHG, orbital parameters and solar constants -->
  <grid_definition>
    <grid id="grid_scalar" >
    <scalar/>
    </grid>
  </grid_definition>
 
  <!-- Define groups of vertical axes -->
  <axis_definition>
    <axis id="time_month" n_glo="12" value="(0,11) [1 2 3 4 5 6 7 8 9 10 11 12]"/>
    <axis id="time_year" unit="day" />
    <axis id="axis_ozone" unit="millibar" axis_type="Z"/>
    <axis id="axis_climoz"  unit="millibar" axis_type="Z"/>
    <axis id="tr_climoz"  />
    <axis id="time_axis_climoz" />


    <axis id="presnivs" standard_name="Vertical levels" unit="Pa">
    </axis>
    <axis id="Ahyb" standard_name="Ahyb comp of Hyb Cord" unit="Pa">
    </axis> 
    <axis id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit="">
    </axis>
    <axis id="Ahyb_inter" standard_name="A comp of Hyb Cord at interface" unit="Pa">
    </axis> 
    <axis id="Bhyb_inter" standard_name="B comp of Hyb Cord at interface" unit="">
    </axis>
    <axis id="Alt" standard_name="Height approx for scale heigh of 8km at levels" unit="km">
    </axis>
    <axis id="plev" standard_name="model_level_number" unit="Pa">
    </axis>
    <axis id="klev"  prec="8" long_name="number of layers"
          standard_name="number of layers" unit="1">
    </axis>
    <axis id="klevp1"  prec="8" long_name="number of layer interfaces"
          standard_name="number of layer interfaces" unit="1">
    </axis>
    <axis id="bnds" standard_name="bounds" unit="1" >
    </axis>
    <axis id="spectband" standard_name="Sensor Band Central Radiation Wavenumber" unit="m-1">
    </axis>
    <axis id="axis_lat" standard_name="Latitude axis">
           <reduce_domain operation="average" direction="iDir" />
    </axis>

<!-- Cosp axis definitions-->
    <axis id="height" standard_name="Cosp levels" unit="m"> 
    </axis>
    <axis id="height_mlev" standard_name="height_mlev" unit="m"> 
    </axis>
    <axis id="sza" standard_name="solar_zenith_angle" unit="degrees">
    </axis>
    <axis id="pressure2" standard_name="pressure" unit="mb">
    </axis>
    <axis id="column" standard_name="column" unit="count"> 
    </axis> 
    <axis id="temp" standard_name="temperature" unit="K"> 
    </axis> 
    <axis id="cth" standard_name="altitude" unit="m"> 
    </axis> 
    <axis id="ReffIce" standard_name="ReffIce" unit="microne" >
    </axis>
    <axis id="ReffLiq" standard_name="ReffLiq" unit="microne" >
    </axis>
    <axis id="scatratio" standard_name="scatratio" unit="1" >
    </axis>
    <axis id="dbze" standard_name="dbze" unit="dBZ" >
    </axis>
    <axis id="tau" standard_name="tau" unit="1" >
    </axis>
  </axis_definition>

  <grid_definition>

    <grid id="klev_bnds"> <axis axis_ref="klev" /> <axis axis_ref="bnds" /> </grid>
    <grid id="klevp1_bnds"> <axis axis_ref="klevp1" /> <axis axis_ref="bnds" /> </grid>

     <grid id="grid_glo">
        <domain domain_ref="dom_glo" />
     </grid>

     <grid id="grid_glo_presnivs">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="presnivs" />
     </grid>


<!------------------------------>
      <grid id="grid_out">
       <domain id="dom_out" />
      </grid>

     <grid id="grid_out_presnivs">
        <domain id="dom_out" />
        <axis id="presnivs" />
     </grid>
<!--------------------------------------------->



     <grid id="grid_glo_plev">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="plev" />
     </grid>

     <grid id="grid_glo_spectband">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="spectband" />
     </grid>

     <grid id="grid_glo_height">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="height" />
     </grid>

     <grid id="grid_glo_heightmlev">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="height_mlev" />
     </grid>

     <grid id="grid_glo_temp">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="temp" />
     </grid>

     <grid id="grid_glo_sza">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="sza" />
     </grid>

     <grid id="grid_glo_column">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="column" />
     </grid>


<!-- Define 4D grids for Cosp simulator -->
     <grid id="grid4Dcol">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="column" />
        <axis axis_ref="height_mlev" />
     </grid>

     <grid id="grid4Dsrbin">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="height" />
        <axis axis_ref="scatratio" />
     </grid>

     <grid id="grid4Ddbze">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="height" />
        <axis axis_ref="dbze" />
     </grid>

     <grid id="grid4Dtau">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="tau" />
        <axis axis_ref="pressure2" />
     </grid>

     <grid id="grid4Dmisr">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="cth" />
        <axis axis_ref="tau" />
     </grid>

     <grid id="grid4Dreffi">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="tau" />
        <axis axis_ref="ReffIce" />
     </grid>

     <grid id="grid4Dreffl">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="tau" />
        <axis axis_ref="ReffLiq" />
     </grid>
     <grid id="grid4Dcol2">
        <domain domain_ref="dom_glo" />
        <axis axis_ref="height" />
        <axis axis_ref="column" />
     </grid>

 <!-- Grid definitions to allow summing of a 3D varaible -->    
      <grid id="grid_3D" >
         <domain domain_ref="dom_glo" />
         <axis axis_ref="lev" />
       </grid>

       <grid id="grid_sum_axis">
         <domain domain_ref="dom_glo" />
         <scalar>
            <reduce_axis operation="sum" />
         </scalar>
       </grid>

       <grid id="grid_sum">
         <scalar >
        <reduce_domain operation="sum" local="true" />
        <reduce_scalar operation="sum" />
         </scalar>
         <scalar>
         </scalar>
       </grid>
       
       <grid id="grid_climoz">
         <domain domain_ref="domain_climoz" />
         <axis axis_ref="axis_climoz" />
         <axis axis_ref="time_axis_climoz" />
         <axis axis_ref="tr_climoz" />
       </grid >
       
       <grid id="grid_from_climoz">
         <domain domain_ref="from_climoz" />
         <axis axis_ref="axis_climoz" />
         <axis axis_ref="time_axis_climoz" />
         <axis axis_ref="tr_climoz" />
       </grid>
       

  </grid_definition>

</context>
