New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 9299 for branches/UKMO/r6232_tracer_advection/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2018-02-01T13:12:00+01:00 (6 years ago)
Author:
jcastill
Message:

First implementation of tracers - it has not been properly tested yet

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/r6232_tracer_advection/NEMOGCM/CONFIG/SHARED/field_def.xml

    r5517 r9299  
    2323      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    2424         <field id="e3t"          long_name="T-cell thickness"   standard_name="cell_thickness"   unit="m"   grid_ref="grid_T_3D"/> 
     25         <field id="e3t_surf"     long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"   unit="m"   grid_ref="grid_T_SFC"/> 
    2526 
    2627         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
     
    613614    --> 
    614615 
    615     <field_group id="trendT" grid_ref="grid_T_3D"> 
    616616      <!-- variables available with ln_tra_trd --> 
    617       <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                                                                                          unit="degC/s"                        /> 
     617 
     618      <!-- Asselin trends  calculated on odd time steps-->  
     619      <field_group id="trendT_odd"  grid_ref="grid_T_3D">  
     620       <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"       unit="degree_C/s" />  
     621       <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"       unit="0.001/s" />  
     622       <!-- Thickness weighted versions: -->  
     623       <field id="ttrd_atf_e3t"      unit="degC/s * m"  >  ttrd_atf * e3t </field>  
     624       <field id="strd_atf_e3t"      unit="1e-3/s * m"  >  strd_atf * e3t </field>  
     625       <!-- OMIP  layer-integrated trends -->  
     626       <field id="ttrd_atf_li"      long_name="layer integrated heat-trend: asselin time filter "       unit="W/m^2" > ttrd_atf_e3t * 1026.0 * 3991.86795711963  </field>  
     627       <field id="strd_atf_li"      long_name="layer integrated salt   -trend: asselin time filter "       unit="kg/(m^2 s)" > strd_atf_e3t * 1026.0 * 0.001 </field>  
     628       </field_group>  
     629 
     630    <!-- Other trends  calculated on even time steps-->  
     631  <field_group id="trendT_even" grid_ref="grid_T_3D">  
     632      <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                                                                                          unit="degC/s"                        />  
    618633      <field id="strd_xad"      long_name="salinity   -trend: i-advection"                                                                                          unit="1e-3/s"                        /> 
    619634      <field id="ttrd_yad"      long_name="temperature-trend: j-advection"                                                                                          unit="degC/s"                        /> 
     
    630645      <field id="strd_zdf"      long_name="salinity   -trend: vertical diffusion"      standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing"        unit="1e-3/s"                        /> 
    631646 
     647      <field id="strd_evd_e3t"      unit="1e-3/s * m" >  strd_evd * e3t </field> 
     648 
    632649      <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 
    633650      <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"   unit="degC/s" /> 
     
    641658      <field id="ttrd_npc"      long_name="temperature-trend: non-penetrative conv."     unit="degC/s" /> 
    642659      <field id="strd_npc"      long_name="salinity   -trend: non-penetrative conv."     unit="1e-3/s" /> 
    643       <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"   unit="degC/s" /> 
    644       <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"   unit="degC/s" /> 
     660      <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"   unit="degC/s" grid_ref="grid_T_2D" /> 
     661      <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"   unit="degC/s" grid_ref="grid_T_2D" /> 
    645662      <field id="ttrd_qsr"      long_name="temperature-trend: solar penetr. heating"     unit="degC/s" /> 
    646663      <field id="ttrd_bbc"      long_name="temperature-trend: geothermal heating"        unit="degC/s" /> 
    647       <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"       unit="degC/s" /> 
    648       <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"       unit="1e-3/s" /> 
    649  
     664 
     665      <!-- OMIP  layer-integrated trends -->  
     666      <field id="ttrd_totad_li"    long_name="layer integrated heat-trend : total advection"       unit="W/m^2"     > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field>  
     667      <field id="strd_totad_li"    long_name="layer integrated salt   -trend : total advection"      unit="kg/(m^2 s)"    > strd_totad_e3t * 1026.0 * 0.001  </field>  
     668      <field id="ttrd_evd_li"      long_name="layer integrated heat-trend : EVD convection"         unit="W/m^2"    > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>  
     669      <field id="strd_evd_li"      long_name="layer integrated salt   -trend : EVD convection"      unit="kg/(m^2 s)"  > strd_evd_e3t * 1026.0 * 0.001  </field>  
     670      <field id="ttrd_iso_li"      long_name="layer integrated heat-trend : isopycnal diffusion"    unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field>  
     671      <field id="strd_iso_li"      long_name="layer integrated salt   -trend : isopycnal diffusion"   unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001  </field>  
     672      <field id="ttrd_zdfp_li"     long_name="layer integrated heat-trend : pure vert. diffusion"   unit="W/m^2" > ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field>  
     673      <field id="strd_zdfp_li"     long_name="layer integrated salt   -trend : pure vert. diffusion"   unit="kg/(m^2 s)" > strd_zdfp_e3t * 1026.0 * 0.001  </field>  
     674      <field id="ttrd_qns_li"      long_name="layer integrated heat-trend : non-solar flux + runoff"   unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field>  
     675      <field id="ttrd_qsr_li"      long_name="layer integrated heat-trend : solar flux"   unit="W/m^2"  grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field>  
     676      <field id="ttrd_bbl_li"      long_name="layer integrated heat-trend: bottom boundary layer "     unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field>  
     677      <field id="strd_bbl_li"      long_name="layer integrated salt   -trend: bottom boundary layer "     unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001  </field>  
     678      <field id="ttrd_evd_li"      long_name="layer integrated heat -trend: evd convection "       unit="W/m^2" >ttrd_evd_e3t * 1026.0 * 3991.86795711963  </field>  
     679      <field id="strd_evd_li"      long_name="layer integrated salt -trend: evd convection "       unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001  </field>  
     680     </field_group>  
     681            
     682     <!--  Total trends calculated every time step-->  
     683     <field_group id="trendT" grid_ref="grid_T_3D">  
     684      <field id="ttrd_tot"      long_name="temperature-trend: total model trend"         unit="degC/s" />  
     685      <field id="strd_tot"      long_name="salinity   -trend: total model trend"         unit="1e-3/s" />  
     686      <!-- Thickness weighted versions: -->  
     687      <!-- OMIP  layer-integrated total trends -->  
     688      <field id="ttrd_tot_li"      long_name="layer integrated heat-trend: total model trend :"         unit="W/m^2" > ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field>  
     689      <field id="strd_tot_li"      long_name="layer integrated salt   -trend: total model trend :"         unit="kg/(m^2 s)" > strd_tot_e3t * 1026.0 * 0.001  </field>  
     690                     
     691      <!-- **** these trends have not been apportioned to all/even/odd ts yet **** -->  
    650692      <!-- variables available with ln_KE_trd --> 
    651693      <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"          unit="W/s^3"                        /> 
     
    946988    --> 
    947989 
     990     <field_group id="TRD" >  
     991          <field field_ref="ttrd_totad_li"   name="opottempadvect"  />  
     992          <field field_ref="ttrd_iso_li"     name="opottemppmdiff"  />  
     993          <field field_ref="ttrd_zdfp_li"    name="opottempdiff"  />  
     994          <field field_ref="ttrd_evd_li"     name="opottempevd" />  
     995          <field field_ref="strd_evd_li"     name="osaltevd" />  
     996          <field field_ref="ttrd_qns_li"     name="opottempqns"  />  
     997          <field field_ref="ttrd_qsr_li"     name="rsdoabsorb" operation="accumulate" />  
     998          <field field_ref="strd_totad_li"   name="osaltadvect" />  
     999          <field field_ref="strd_iso_li"     name="osaltpmdiff"  />  
     1000          <field field_ref="strd_zdfp_li"    name="osaltdiff" />  
     1001    </field_group>  
     1002      
    9481003    <field_group id="mooring" > 
    9491004      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"      /> 
Note: See TracChangeset for help on using the changeset viewer.