Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/domain_def.xml
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/domain_def.xml	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/domain_def.xml	(revision 6534)
@@ -0,0 +1,1 @@
+link ../../SHARED/domain_def.xml
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/field_def.xml
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/field_def.xml	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/field_def.xml	(revision 6534)
@@ -0,0 +1,1072 @@
+<!--
+<?xml version="1.0"?>
+       -->
+ 
+    <!-- $id$ -->
+    
+    <!-- 
+============================================================================================================
+=                                  definition of all existing variables                                    =
+=                                            DO NOT CHANGE                                                 =
+============================================================================================================
+    -->
+   <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->
+
+    <!-- 
+============================================================================================================
+                                  Physical ocean and ice model variables
+============================================================================================================
+    -->
+
+      <!-- T grid -->
+      
+      <field_group id="grid_T" grid_ref="grid_T_2D" >
+         <field id="e3t"          long_name="T-cell thickness"   standard_name="cell_thickness"   unit="m"   grid_ref="grid_T_3D"/>
+         <field id="e3t_0"        long_name="Initial T-cell thickness"   standard_name="ref_cell_thickness"   unit="m"   grid_ref="grid_T_3D"/>
+
+         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/>
+         <field id="toce_e3t"     long_name="temperature * e3t"                                                     unit="degC*m"   grid_ref="grid_T_3D" > toce * e3t </field >
+         <field id="soce"         long_name="salinity"            standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/>
+         <field id="soce_e3t"     long_name="salinity * e3t"                                                        unit="1e-3*m"   grid_ref="grid_T_3D" > soce * e3t </field >
+
+         <field id="sst"          long_name="sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     />
+         <field id="sst2"         long_name="square of sea surface temperature"   standard_name="square_of_sea_surface_temperature"   unit="degC2"     > sst * sst </field >
+         <field id="sstmax"       long_name="max of sea surface temperature"   field_ref="sst"   operation="maximum"                                  />
+         <field id="sstmin"       long_name="min of sea surface temperature"   field_ref="sst"   operation="minimum"                                  />
+         <field id="sstgrad"      long_name="module of sst gradient"                                                                  unit="degC/m"   />
+         <field id="sstgrad2"     long_name="square of module of sst gradient"                                                        unit="degC2/m2" />
+         <field id="sbt"          long_name="sea bottom temperature"                                                                  unit="degC"     />
+	 
+         <field id="sss"          long_name="sea surface salinity"             standard_name="sea_surface_salinity"   unit="1e-3" />
+         <field id="sss2"         long_name="square of sea surface salinity"                                          unit="1e-6"  > sss * sss </field >
+         <field id="sssmax"       long_name="max of sea surface salinity"   field_ref="sss"   operation="maximum"                 />
+         <field id="sssmin"       long_name="min of sea surface salinity"   field_ref="sss"   operation="minimum"                 />
+         <field id="sbs"          long_name="sea bottom salinity"                                                     unit="1e-3" />
+
+         <field id="taubot"       long_name="bottom stress module"                                                    unit="N/m2" /> 
+
+         <field id="ssh"          long_name="sea surface height"             standard_name="sea_surface_height_above_geoid"             unit="m" />
+         <field id="ssh2"         long_name="square of sea surface height"   standard_name="square_of_sea_surface_height_above_geoid"   unit="m2" > ssh * ssh </field >
+         <field id="sshmax"       long_name="max of sea surface height"   field_ref="ssh"   operation="maximum"                                  />
+
+         <field id="mldkz5"       long_name="Turbocline depth (Kz = 5e-4)"                       standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity"                unit="m"          />
+         <field id="mldr10_1"     long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)"          standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"          />
+         <field id="mldr10_1max"  long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="maximum"                                                                          />
+         <field id="mldr10_1min"  long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="minimum"                                                                          />
+         <field id="heatc"        long_name="Heat content vertically integrated"                 standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"       />
+         <field id="saltc"        long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3*kg/m2" />
+
+         <!-- EOS -->
+         <field id="alpha"        long_name="thermal expansion"                                                         unit="degC-1" grid_ref="grid_T_3D" />
+         <field id="beta"         long_name="haline contraction"                                                        unit="1e3"    grid_ref="grid_T_3D" />
+         <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" />
+
+         <!-- Energy - horizontal divergence -->
+         <field id="eken"         long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" />
+         <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" />
+
+         <!-- variables available with MLE -->
+         <field id="Lf_NHpf"      long_name="MLE: Lf = N H / f"   unit="m" />
+
+         <!-- next variables available with key_diahth -->
+         <field id="mlddzt"       long_name="Thermocline Depth (depth of max dT/dz)"         standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature"             unit="m"                         />
+         <field id="mldr10_3"     long_name="Mixed Layer Depth (dsigma = 0.03 wrt 10m)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
+         <field id="mldr0_1"      long_name="Mixed Layer Depth (dsigma = 0.01 wrt sfc)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
+         <field id="mldr0_3"      long_name="Mixed Layer Depth (dsigma = 0.03 wrt sfc)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
+         <field id="mld_dt02"     long_name="Mixed Layer Depth (|dT| = 0.2 wrt 10m)"         standard_name="ocean_mixed_layer_thickness_defined_by_temperature"                                unit="m"                         />
+         <field id="topthdep"     long_name="Top of Thermocline Depth (dT = -0.2 wrt 10m)"   standard_name="ocean_mixed_layer_thickness_defined_by_temperature"                                unit="m"                         />
+         <field id="pycndep"      long_name="Pycnocline Depth (dsigma[dT=-0.2] wrt 10m)"     standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
+         <field id="BLT"          long_name="Barrier Layer Thickness"                                                                                                                          unit="m"                          > topthdep - pycndep </field>
+         <field id="tinv"         long_name="Max of vertical invertion of temperature"                                                                                                         unit="degC"                      />
+         <field id="depti"        long_name="Depth of max. vert. inv. of temperature"                                                                                                          unit="m"                         />
+         <field id="20d"          long_name="Depth of 20C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"      axis_ref="iax_20C" />
+         <field id="28d"          long_name="Depth of 28C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"      axis_ref="iax_28C" />
+         <field id="hc300"        long_name="Heat content 0-300m"                            standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"                      />
+
+         <!-- variables available with key_diaar5 -->
+         <field id="botpres"      long_name="Pressure at sea floor"   standard_name="sea_water_pressure_at_sea_floor"   unit="dbar" />
+
+         <!-- variables available with key_vvl -->
+         <field id="tpt_dep"      long_name="T-point depth"                  standard_name="depth_below_geoid"   unit="m"   grid_ref="grid_T_3D" />
+         <field id="e3tdef"       long_name="T-cell thickness deformation"                                       unit="%"   grid_ref="grid_T_3D" />
+      </field_group>
+
+      <!-- Tides -->
+
+      <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" >
+         <!-- tidal composante -->
+         <field id="M2x"          long_name="M2 Elevation harmonic real part "                             unit="m"        />
+         <field id="M2y"          long_name="M2 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="S2x"          long_name="S2 Elevation harmonic real part "                             unit="m"        />
+         <field id="S2y"          long_name="S2 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="N2x"          long_name="N2 Elevation harmonic real part "                             unit="m"        />
+         <field id="N2y"          long_name="N2 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="K1x"          long_name="K1 Elevation harmonic real part "                             unit="m"        />
+         <field id="K1y"          long_name="K1 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="O1x"          long_name="O1 Elevation harmonic real part "                             unit="m"        />
+         <field id="O1y"          long_name="O1 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="Q1x"          long_name="Q1 Elevation harmonic real part "                             unit="m"        />
+         <field id="Q1y"          long_name="Q1 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="M4x"          long_name="M4 Elevation harmonic real part "                             unit="m"        />
+         <field id="M4y"          long_name="M4 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="K2x"          long_name="K2 Elevation harmonic real part "                             unit="m"        />
+         <field id="K2y"          long_name="K2 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="P1x"          long_name="P1 Elevation harmonic real part "                             unit="m"        />
+         <field id="P1y"          long_name="P1 Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="Mfx"          long_name="Mf Elevation harmonic real part "                             unit="m"        />
+         <field id="Mfy"          long_name="Mf Elevation harmonic imaginary part"                         unit="m"        />
+         <field id="Mmx"          long_name="Mm Elevation harmonic real part "                             unit="m"        />
+         <field id="Mmy"          long_name="Mm Elevation harmonic imaginary part"                         unit="m"        />
+      </field_group>
+	 
+      <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once" >
+         <field id="M2x_u"        long_name="M2 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="M2y_u"        long_name="M2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="S2x_u"        long_name="S2 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="S2y_u"        long_name="S2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="N2x_u"        long_name="N2 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="N2y_u"        long_name="N2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="K1x_u"        long_name="K1 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="K1y_u"        long_name="K1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="O1x_u"        long_name="O1 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="O1y_u"        long_name="O1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Q1x_u"        long_name="Q1 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="Q1y_u"        long_name="Q1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="M4x_u"        long_name="M4 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="M4y_u"        long_name="M4 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="K2x_u"        long_name="K2 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="K2y_u"        long_name="K2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="P1x_u"        long_name="P1 current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="P1y_u"        long_name="P1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Mfx_u"        long_name="Mf current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="Mfy_u"        long_name="Mf current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Mmx_u"        long_name="Mm current barotrope along i-axis harmonic real part "        unit="m/s"      />
+         <field id="Mmy_u"        long_name="Mm current barotrope along i-axis harmonic imaginary part "   unit="m/s"      />
+      </field_group>
+	 
+      <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once" >
+         <field id="M2x_v"        long_name="M2 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="M2y_v"        long_name="M2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="S2x_v"        long_name="S2 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="S2y_v"        long_name="S2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="N2x_v"        long_name="N2 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="N2y_v"        long_name="N2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="K1x_v"        long_name="K1 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="K1y_v"        long_name="K1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="O1x_v"        long_name="O1 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="O1y_v"        long_name="O1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Q1x_v"        long_name="Q1 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="Q1y_v"        long_name="Q1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="M4x_v"        long_name="M4 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="M4y_v"        long_name="M4 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="K2x_v"        long_name="K2 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="K2y_v"        long_name="K2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="P1x_v"        long_name="P1 current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="P1y_v"        long_name="P1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Mfx_v"        long_name="Mf current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="Mfy_v"        long_name="Mf current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />
+         <field id="Mmx_v"        long_name="Mm current barotrope along j-axis harmonic real part "        unit="m/s"      />
+         <field id="Mmy_v"        long_name="Mm current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />	 
+      </field_group>
+
+      <!-- SBC -->
+      
+      <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc -->
+         <field id="empmr"        long_name="Net Upward Water Flux"                standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   />
+         <field id="empbmr"       long_name="Net Upward Water Flux at pre. tstep"  standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   />
+         <field id="emp_oce"      long_name="Evap minus Precip over ocean"         standard_name="evap_minus_precip_over_sea_water"                                     unit="kg/m2/s"   />
+         <field id="emp_ice"      long_name="Evap minus Precip over ice"           standard_name="evap_minus_precip_over_sea_ice"                                       unit="kg/m2/s"   />
+         <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="1e-3/m2/s" />
+         <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   />
+         <field id="snowpre"      long_name="Snow precipitation"                   standard_name="snowfall_flux"                                                        unit="kg/m2/s"   />
+         <field id="runoffs"      long_name="River Runoffs"                        standard_name="water_flux_into_sea_water_from_rivers"                                unit="kg/m2/s"   />
+         <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   />
+ 
+         <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           />
+         <field id="qns"          long_name="non solar Downward Heat Flux"                                                                                               unit="W/m2"                           />
+         <field id="qsr"          long_name="Shortwave Radiation"                   standard_name="net_downward_shortwave_flux_at_sea_water_surface"                     unit="W/m2"                           />
+         <field id="qsr3d"        long_name="Shortwave Radiation 3D distribution"   standard_name="downwelling_shortwave_flux_in_sea_water"                              unit="W/m2"      grid_ref="grid_T_3D" />
+         <field id="qrp"          long_name="Surface Heat Flux: Damping"            standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation"                 unit="W/m2"                           />
+         <field id="erp"          long_name="Surface Water Flux: Damping"           standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation"              unit="kg/m2/s"                        />
+         <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           />
+         <field id="wspd"         long_name="wind speed module"                     standard_name="wind_speed"                                                           unit="m/s"                            />
+         
+         <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn -->
+         <field id="ssh_ib"       long_name="Inverse barometer sea surface height"  standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency"   unit="m"        />
+
+         <!-- * variable related to ice shelf forcing * -->
+         <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  />
+         <field id="qisf"         long_name="Ice Shelf Heat Flux"                           unit="W/m2"     />
+         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature)"   unit="m/s"      />
+         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)"      unit="m/s"      />
+         <field id="stbl"         long_name="salinity in the Losh tbl"                      unit="1e-3"     />
+         <field id="ttbl"         long_name="temperature in the Losh tbl"                   unit="degC"     />
+
+         <!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
+         <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"  standard_name="surface_net_downward_longwave_flux"   unit="W/m2"  />
+         <field id="qsb_oce"      long_name="Sensible Downward Heat Flux over open ocean"  standard_name="surface_downward_sensible_heat_flux"  unit="W/m2"  />
+         <field id="qla_oce"      long_name="Latent Downward Heat Flux over open ocean"    standard_name="surface_downward_latent_heat_flux"    unit="W/m2"  />
+         <field id="qemp_oce"     long_name="Downward Heat Flux from E-P over open ocean"                                                       unit="W/m2"  />
+         <field id="taum_oce"     long_name="wind stress module over open ocean"           standard_name="magnitude_of_surface_downward_stress" unit="N/m2"  />
+
+         <!-- available key_oasis3 -->
+         <field id="snow_ao_cea"  long_name="Snow over ice-free ocean (cell average)"   standard_name="snowfall_flux"                             unit="kg/m2/s"  />
+         <field id="snow_ai_cea"  long_name="Snow over sea-ice (cell average)"          standard_name="snowfall_flux"                             unit="kg/m2/s"  />
+         <field id="subl_ai_cea"  long_name="Sublimation over sea-ice (cell average)"   standard_name="surface_snow_and_ice_sublimation_flux"     unit="kg/m2/s"  />
+         <field id="icealb_cea"   long_name="Ice albedo (cell average)"                 standard_name="sea_ice_albedo"                            unit="1"        />
+         <field id="calving_cea"  long_name="Calving"                                   standard_name="water_flux_into_sea_water_from_icebergs"   unit="kg/m2/s"  />
+
+         <!-- available if key_oasis3 + conservative method -->
+         <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  />
+         <field id="evap_ao_cea"   long_name="Evaporation over ice-free ocean (cell average)"           standard_name="water_evaporation_flux"                                                        unit="kg/m2/s"  />
+         <field id="isnwmlt_cea"   long_name="Snow over Ice melting (cell average)"                     standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  />
+         <field id="fsal_virt_cea" long_name="Virtual salt flux due to ice formation (cell average)"    standard_name="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics"                unit="kg/m2/s"  />
+         <field id="fsal_real_cea" long_name="Real salt flux due to ice formation (cell average)"       standard_name="downward_sea_ice_basal_salt_flux"                                              unit="kg/m2/s"  />
+         <field id="hflx_rain_cea" long_name="heat flux due to rainfall"                                standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"        unit="W/m2"     />
+         <field id="hflx_evap_cea" long_name="heat flux due to evaporation"                             standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water"   unit="W/m2"     />
+         <field id="hflx_snow_cea" long_name="heat flux due to snow falling over ice-free ocean"        standard_name="heat_flux_into_sea_water_due_to_snow_thermodynamics"                           unit="W/m2"     />
+         <field id="hflx_ice_cea"  long_name="heat flux due to ice thermodynamics"                      standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics"                        unit="W/m2"     />
+         <field id="hflx_rnf_cea"  long_name="heat flux due to runoffs"                                 standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water"          unit="W/m2"     />
+         <field id="hflx_cal_cea"  long_name="heat flux due to calving"                                 standard_name="heat_flux_into_sea_water_due_to_iceberg_thermodynamics"                        unit="W/m2"     />
+         <field id="bicemel_cea"   long_name="Rate of Melt at Sea Ice Base (cell average)"              standard_name="tendency_of_sea_ice_amount_due_to_basal_melting"                               unit="kg/m2/s"  />
+         <field id="licepro_cea"   long_name="Lateral Sea Ice Growth Rate (cell average)"               standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes"                 unit="kg/m2/s"  />
+         <field id="snowmel_cea"   long_name="Snow Melt Rate (cell average)"                            standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  />
+         <field id="sntoice_cea"   long_name="Snow-Ice Formation Rate (cell average)"                   standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion"                             unit="kg/m2/s"  />
+         <field id="ticemel_cea"   long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"  standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"                             unit="kg/m2/s"  />
+
+         <!-- ice fields -->
+         <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            />
+
+         <field id="qsr_ai_cea"   long_name="Air-Ice downward solar heat flux (cell average)"              standard_name="surface_downwelling_shortwave_flux_in_air"          unit="W/m2"         />
+         <field id="qns_ai_cea"   long_name="Air-Ice downward non-solar heat flux (cell average)"                                                                             unit="W/m2"         />
+         <field id="qla_ai_cea"   long_name="Air-Ice downward Latent heat flux (cell average)"             standard_name="surface_downward_latent_heat_flux"                  unit="W/m2"         />
+         
+         <field id="qsr_io_cea"   long_name="Ice-Oce downward solar heat flux (cell average)"              standard_name="net_downward_shortwave_flux_at_sea_water_surface"   unit="W/m2"         />
+         <field id="qns_io_cea"   long_name="Ice-Oce downward non-solar heat flux (cell average)"                                                                             unit="W/m2"         />
+         
+         <field id="snowthic_cea" long_name="Snow thickness (cell average)"                                standard_name="surface_snow_thickness"                             unit="m"            />
+         <field id="icethic_cea"  long_name="Ice thickness (cell average)"                                 standard_name="sea_ice_thickness"                                  unit="m"            />
+         <field id="iceprod_cea"  long_name="Ice production (cell average)"                                                                                                   unit="m/s"          />
+         <field id="iiceconc"     long_name="Ice concentration"                                            standard_name="sea_ice_area_fraction"                              unit="1"            />
+         
+         <field id="ice_pres"     long_name="Ice presence"                                                                                                                    unit=""             />
+         <field id="ist_cea"      long_name="Ice surface temperature (cell average)"                       standard_name="surface_temperature"                                unit="degC"         />
+         <field id="ist_ipa"      long_name="Ice surface temperature (ice presence average)"               standard_name="surface_temperature"                                unit="degC"         />      
+         <field id="uice_ipa"     long_name="Ice velocity along i-axis at I-point (ice presence average)"  standard_name="sea_ice_x_velocity"                                 unit="m/s"          />      
+         <field id="vice_ipa"     long_name="Ice velocity along j-axis at I-point (ice presence average)"  standard_name="sea_ice_y_velocity"                                 unit="m/s"          />      
+         
+         <field id="utau_ice"     long_name="Wind stress along i-axis over the ice at i-point"             standard_name="surface_downward_x_stress"                          unit="N/m2"         />
+         <field id="vtau_ice"     long_name="Wind stress along j-axis over the ice at i-point"             standard_name="surface_downward_y_stress"                          unit="N/m2"         />
+         
+         <field id="u_imasstr"    long_name="Sea-ice mass transport along i-axis"                          standard_name="sea_ice_x_transport"                                unit="kg/s"         />
+         <field id="v_imasstr"    long_name="Sea-ice mass transport along j-axis"                          standard_name="sea_ice_y_transport"                                unit="kg/s"         />
+         <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                                                                                              unit="kg*degC/m2/s" />
+         <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" />        
+         <field id="rnf_x_sst"    long_name="Runoff term on SST"                                                                                                              unit="kg*degC/m2/s" />
+         <field id="rnf_x_sss"    long_name="Runoff term on SSS"                                                                                                              unit="kg*1e-3/m2/s" />
+       
+         <field id="iceconc"      long_name="ice concentration"                                            standard_name="sea_ice_area_fraction"                              unit="%"            />
+       	 <field id="isst"         long_name="sea surface temperature"                                      standard_name="sea_surface_temperature"                            unit="degC"         />
+         <field id="isss"         long_name="sea surface salinity"                                         standard_name="sea_surface_salinity"                               unit="1e-3"         /> 
+         <field id="qt_oce"       long_name="total flux at ocean surface"                                  standard_name="surface_downward_heat_flux_in_sea_water"            unit="W/m2"         />
+         <field id="qsr_oce"      long_name="solar heat flux at ocean surface"                             standard_name="net_downward_shortwave_flux_at_sea_water_surface"   unit="W/m2"         />
+         <field id="qns_oce"      long_name="non-solar heat flux at ocean surface"                                                                                            unit="W/m2"         />
+         <field id="qt_ice"       long_name="total heat flux at ice surface: sum over categories"          standard_name="surface_downward_heat_flux_in_air"                  unit="W/m2"         />
+         <field id="qsr_ice"      long_name="solar heat flux at ice surface: sum over categories"          standard_name="surface_downwelling_shortwave_flux_in_air"          unit="W/m2"         />
+         <field id="qns_ice"      long_name="non-solar heat flux at ice surface: sum over categories"                                                                         unit="W/m2"         />
+         <field id="qtr_ice"      long_name="solar heat flux transmitted through ice: sum over categories"                                                                    unit="W/m2"         />
+         <field id="qemp_ice"     long_name="Downward Heat Flux from E-P over ice"                                                                                            unit="W/m2"         />
+         <field id="micesalt"     long_name="Mean ice salinity"                                                                                                               unit="1e-3"         />
+         <field id="miceage"      long_name="Mean ice age"                                                                                                                    unit="years"        />
+         <field id="alb_ice"      long_name="Mean albedo over sea ice"                                                                                                        unit=""             />
+         <field id="albedo"       long_name="Mean albedo over sea ice and ocean"                                                                                              unit=""             />
+
+         <field id="iceage_cat"   long_name="Ice age for categories"                                       unit="days"   axis_ref="ncatice" />
+         <field id="iceconc_cat"  long_name="Ice concentration for categories"                             unit="%"      axis_ref="ncatice" />
+         <field id="icethic_cat"  long_name="Ice thickness for categories"                                 unit="m"      axis_ref="ncatice" />
+         <field id="snowthic_cat" long_name="Snow thicknessi for categories"                               unit="m"      axis_ref="ncatice" />
+         <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories"                         unit="g/kg"   axis_ref="ncatice" />
+         <field id="brinevol_cat" long_name="Brine volume for categories"                                  unit="%"      axis_ref="ncatice" />
+         <field id="icetemp_cat"  long_name="Ice temperature for categories"                               unit="degC"   axis_ref="ncatice" />
+         <field id="snwtemp_cat"  long_name="Snow temperature for categories"                              unit="degC"   axis_ref="ncatice" />
+
+         <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     />
+         <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"    /> 
+         <field id="isnowhc"      long_name="snow total heat content"                                      unit="10^9J"    />
+         <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     />
+         <field id="ibrinv"       long_name="brine volume"                                                 unit="%"        />
+         <field id="icecolf"      long_name="frazil ice collection thickness"                              unit="m"        />
+         <field id="icestr"       long_name="ice strength"                                                 unit="N/m"      />
+         <field id="icevel"       long_name="ice velocity"                                                 unit="m/s"      />
+         <field id="idive"        long_name="divergence"                                                   unit="1e-8s-1"  />
+         <field id="ishear"       long_name="shear"                                                        unit="1e-8s-1"  />
+         <field id="icevolu"      long_name="ice volume"                                                   unit="m"        />
+         <field id="snowvol"      long_name="snow volume"                                                  unit="m"        />
+
+         <field id="icetrp"       long_name="ice volume transport"                                         unit="m/day"          />
+         <field id="snwtrp"       long_name="snw volume transport"                                         unit="m/day"          />
+         <field id="saltrp"       long_name="salt content transport"                                       unit="1e-3*kg/m2/day" />
+         <field id="deitrp"       long_name="advected ice enthalpy"                                        unit="W/m2"           />
+         <field id="destrp"       long_name="advected snw enthalpy"                                        unit="W/m2"           />
+
+         <field id="sfxbri"       long_name="brine salt flux"                                              unit="1e-3*kg/m2/day" />
+         <field id="sfxdyn"       long_name="salt flux from ridging rafting"                               unit="1e-3*kg/m2/day" />
+         <field id="sfxres"       long_name="salt flux from lipupdate (resultant)"                         unit="1e-3*kg/m2/day" />
+         <field id="sfxbog"       long_name="salt flux from bot growth"                                    unit="1e-3*kg/m2/day" />
+         <field id="sfxbom"       long_name="salt flux from bot melt"                                      unit="1e-3*kg/m2/day" />
+         <field id="sfxsum"       long_name="salt flux from surf melt"                                     unit="1e-3*kg/m2/day" />
+         <field id="sfxsni"       long_name="salt flux from snow-ice formation"                            unit="1e-3*kg/m2/day" />
+         <field id="sfxopw"       long_name="salt flux from open water ice formation"                      unit="1e-3*kg/m2/day" />
+         <field id="sfxsub"       long_name="salt flux from sublimation"                                   unit="1e-3*kg/m2/day" />
+         <field id="sfx"          long_name="salt flux total"                                              unit="1e-3*kg/m2/day" />
+
+         <field id="vfxbog"       long_name="daily bottom thermo ice prod."                                unit="m/day"   />
+         <field id="vfxdyn"       long_name="daily  dynamic ice prod."                                     unit="m/day"   />
+         <field id="vfxopw"       long_name="daily lateral thermo ice prod."                               unit="m/day"   />
+         <field id="vfxsni"       long_name="daily snowice ice prod."                                      unit="m/day"   />
+         <field id="vfxsum"       long_name="surface melt"                                                 unit="m/day"   />
+         <field id="vfxbom"       long_name="bottom melt"                                                  unit="m/day"   />
+         <field id="vfxres"       long_name="daily resultant ice prod./melting from limupdate"             unit="m/day"   />
+         <field id="vfxice"       long_name="ice melt/growth"                                              unit="m/day"   />
+         <field id="vfxsnw"       long_name="snw melt/growth"                                              unit="m/day"   />
+         <field id="vfxsub"       long_name="snw sublimation"                                              unit="m/day"   />
+         <field id="vfxspr"       long_name="snw precipitation on ice"                                     unit="m/day"   />
+         <field id="vfxthin"      long_name="daily thermo ice prod. for thin ice(<20cm) + open water"      unit="m/day"   />
+
+         <field id="afxtot"       long_name="area tendency (total)"                                        unit="day-1"   />
+         <field id="afxdyn"       long_name="area tendency (dynamics)"                                     unit="day-1"   />
+         <field id="afxthd"       long_name="area tendency (thermo)"                                       unit="day-1"   />
+
+         <field id="hfxsum"       long_name="heat fluxes causing surface ice melt"                         unit="W/m2"  />
+         <field id="hfxbom"       long_name="heat fluxes causing bottom ice melt"                          unit="W/m2"  />
+         <field id="hfxbog"       long_name="heat fluxes causing bottom ice growth"                        unit="W/m2"  />
+         <field id="hfxdif"       long_name="heat fluxes causing ice temperature change"                   unit="W/m2"  />
+         <field id="hfxopw"       long_name="heat fluxes causing open water ice formation"                 unit="W/m2"  />
+         <field id="hfxsnw"       long_name="heat fluxes causing snow melt"                                unit="W/m2"  />
+         <field id="hfxerr"       long_name="heat fluxes error after heat diffusion"                       unit="W/m2"  />
+         <field id="hfxerr_rem"   long_name="heat fluxes error after remapping"                            unit="W/m2"  />
+         <field id="hfxout"       long_name="total heat fluxes received by the ocean"                      unit="W/m2"  />
+         <field id="hfxin"        long_name="total heat fluxes at the ice/ocean surface"                   unit="W/m2"  />
+
+         <!-- heat flux associated with mass exchange -->
+         <field id="hfxthd"       long_name="heat fluxes from ice-ocean mass exchange during thermo"       unit="W/m2"  />
+         <field id="hfxdyn"       long_name="heat fluxes from ice-ocean mass exchange during dynamic"      unit="W/m2"  />
+         <field id="hfxres"       long_name="heat fluxes from ice-ocean mass exchange during resultant"    unit="W/m2"  />
+         <field id="hfxsub"       long_name="heat fluxes from ice-atm. mass exchange during sublimation"   unit="W/m2"  />
+         <field id="hfxspr"       long_name="heat fluxes from ice-atm. mass exchange during snow precip"   unit="W/m2"  />
+
+         <!-- diags -->
+         <field id="hfxdhc"       long_name="Heat content variation in snow and ice"                       unit="W/m2" />
+         <field id="hfxtur"       long_name="turbulent heat flux at the ice base"                          unit="W/m2" />
+	 <!-- sbcssm variables -->
+         <field id="sst_m"    unit="degC" />
+         <field id="sss_m"    unit="psu"  />
+         <field id="ssu_m"    unit="m/s"  />
+         <field id="ssv_m"    unit="m/s"  />
+         <field id="ssh_m"    unit="m"    />
+         <field id="e3t_m"    unit="m"    />
+         <field id="frq_m"    unit="-"    />
+
+      </field_group>
+
+      <!-- U grid -->
+      
+      <field_group id="grid_U"   grid_ref="grid_U_2D">
+         <field id="e3u"          long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" />
+         <field id="e3u_0"        long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/>
+         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            />
+         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" />
+         <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"                                                                   unit="m2/s"       grid_ref="grid_U_3D"  > uoce * e3u </field>
+         <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             />
+         <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             />
+         <field id="ubar"         long_name="ocean barotropic current along i-axis"                                                              unit="m/s"                             />
+         <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"                 standard_name="ocean_volume_x_transport"    unit="m3/s"       grid_ref="grid_U_3D" />
+         <field id="uocet"        long_name="ocean transport along i-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_U_3D" />
+         <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" />
+
+         <!-- variables available with MLE -->
+         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" />
+
+         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv -->
+         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" />
+
+         <!-- uoce_eiv: available with key_trabbl -->
+         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"    unit="m/s"  />
+         <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"   unit="m3/s" />
+
+         <!-- variable for ice shelves -->
+         <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" />
+
+         <!-- variables available with key_diaar5 -->
+         <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"    standard_name="ocean_mass_x_transport"                          unit="kg/s"        grid_ref="grid_U_3D" />
+         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"    standard_name="ocean_heat_x_transport"                          unit="W"                                />
+         <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_x_transport"                          unit="1e-3*kg/s"                        />
+         <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"       standard_name="ocean_heat_x_transport_due_to_bolus_advection"   unit="W"                                />
+         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis"   standard_name="ocean_heat_x_transport_due_to_diffusion"         unit="W"                                />
+      </field_group>
+      
+      <!-- V grid -->
+      
+      <field_group id="grid_V"   grid_ref="grid_V_2D">
+         <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" />
+         <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D"/>
+         <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            />
+         <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" />
+         <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"                                                                   unit="m2/s"       grid_ref="grid_V_3D"  > voce * e3v </field>
+         <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             />
+         <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"                             />
+         <field id="vbar"         long_name="ocean barotropic current along j-axis"                                                              unit="m/s"                             />
+         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"                 standard_name="ocean_volume_y_transport"    unit="m3/s"       grid_ref="grid_V_3D" />
+         <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" />
+         <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" />
+
+         <!-- variables available with MLE -->
+         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" />
+
+         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv -->
+         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" />
+
+         <!-- voce_eiv: available with key_trabbl -->
+         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"    unit="m/s"  />
+         <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"   unit="m3/s" />
+
+         <!-- variable for ice shelves -->
+         <field id="vtbl"         long_name="meridional current in the Losh tbl"   unit="m/s" />
+
+         <!-- variables available with key_diaar5 -->
+         <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"    standard_name="ocean_mass_y_transport"                          unit="kg/s"        grid_ref="grid_V_3D" />
+         <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"    standard_name="ocean_heat_y_transport"                          unit="W"                                />
+         <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"                        />
+         <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"                                />
+         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis"   standard_name="ocean_heat_y_transport_due_to_diffusion"         unit="W"                                />
+      </field_group>
+      
+      <!-- W grid -->
+      
+      <field_group id="grid_W" grid_ref="grid_W_3D">
+        <field id="e3w"          long_name="W-cell thickness"                     standard_name="cell_thickness"              unit="m"    />
+        <field id="woce"         long_name="ocean vertical velocity"              standard_name="upward_sea_water_velocity"   unit="m/s"  />
+        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                               unit="m3/s" />
+
+        <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv -->
+        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" />
+
+        <field id="avt"          long_name="vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />
+        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />
+        <field id="avm"          long_name="vertical eddy viscosity"     standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" />
+
+        <!-- avs: available with key_zdfddm -->
+        <field id="avs"          long_name="salt vertical eddy diffusivity"   standard_name="ocean_vertical_salt_diffusivity"   unit="m2/s" />
+        <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />
+
+        <!-- avt_evd and avm_evd: available with ln_zdfevd -->
+        <field id="avt_evd"      long_name="convective enhancement to vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" />
+        <field id="avm_evd"      long_name="convective enhancement to vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" />
+
+        <!-- avt_tide: available with key_zdftmx -->
+        <field id="av_tide"      long_name="tidal vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_tides"   unit="m2/s" />
+
+       <!-- variables available with key_zdftmx_new -->
+        <field id="av_ratio"     long_name="S over T diffusivity ratio"            standard_name="salinity_over_temperature_diffusivity_ratio"                     unit="1"    />
+        <field id="av_wave"      long_name="wave-induced vertical diffusivity"     standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves"         unit="m2/s" />
+        <field id="bn2"          long_name="squared Brunt-Vaisala frequency"       standard_name="squared_brunt_vaisala_frequency"                                 unit="s-1"  />
+        <field id="bflx_tmx"     long_name="wave-induced buoyancy flux"            standard_name="buoyancy_flux_due_to_internal_waves"                             unit="W/kg" />
+        <field id="pcmap_tmx"    long_name="power consumed by wave-driven mixing"  standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing"   unit="W/m2"      grid_ref="grid_W_2D" />
+        <field id="emix_tmx"     long_name="power density available for mixing"    standard_name="power_available_for_mixing_from_breaking_internal_waves"         unit="W/kg" />
+
+        <!-- variables available with key_diaar5 -->   
+        <field id="w_masstr"     long_name="vertical mass transport"             standard_name="upward_ocean_mass_transport"             unit="kg/s"   />
+        <field id="w_masstr2"    long_name="square of vertical mass transport"   standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" />
+
+        <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d -->
+        <field id="aht2d"        long_name="lateral eddy diffusivity"       standard_name="ocean_tracer_xy_laplacian_diffusivity"      unit="m2/s"   grid_ref="grid_W_2D" />
+        <field id="aht2d_eiv"    long_name="EIV lateral eddy diffusivity"   standard_name="ocean_tracer_bolus_laplacian_diffusivity"   unit="m2/s"   grid_ref="grid_W_2D" />
+      </field_group>
+          
+      <!-- scalar variables available with key_diaar5 -->
+
+      <field_group id="scalar"  domain_ref="1point" >
+         <field id="voltot"     long_name="global total volume"                          standard_name="sea_water_volume"                               unit="m3"   />
+         <field id="sshtot"     long_name="global mean ssh"                              standard_name="global_average_sea_level_change"                unit="m"    />
+         <field id="sshsteric"  long_name="global mean ssh steric"                       standard_name="global_average_steric_sea_level_change"         unit="m"    />
+         <field id="sshthster"  long_name="global mean ssh thermosteric"                 standard_name="global_average_thermosteric_sea_level_change"   unit="m"    />
+         <field id="masstot"    long_name="global total mass"                            standard_name="sea_water_mass"                                 unit="kg"   />
+         <field id="temptot"    long_name="global mean temperature"                      standard_name="sea_water_potential_temperature"                unit="degC" />
+         <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" />
+         <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait"   standard_name="sea_ice_transport_across_line"                  unit="kg/s" />
+
+      	 <!-- available with ln_diahsb -->
+    	 <field id="bgtemper"     long_name="drift in global mean temperature wrt timestep 1"                 standard_name="change_over_time_in_sea_water_potential_temperature"   unit="degC"     />
+    	 <field id="bgsaline"     long_name="drift in global mean salinity wrt timestep 1"                    standard_name="change_over_time_in_sea_water_practical_salinity"      unit="1e-3"     />
+    	 <field id="bgheatco"     long_name="drift in global mean heat content wrt timestep 1"                                                                                      unit="10^9J"    />
+    	 <field id="bgsaltco"     long_name="drift in global mean salt content wrt timestep 1"                                                                                      unit="1e-3*m3"  />
+    	 <field id="bgvolssh"     long_name="drift in global mean ssh volume wrt timestep 1"                                                                                        unit="km3"      />
+         <field id="bgvole3t"     long_name="drift in global mean volume variation (e3t) wrt timestep 1"                                                                            unit="km3"      />
+    	 <field id="bgvoltot"     long_name="drift in global mean volume wrt timestep 1"                                                                                            unit="km3"      />
+         <!-- NOTE: No matching iom_put call -->
+    	 <field id="bgsshtot"     long_name="drift in global mean ssh wrt timestep 1"                         standard_name="global_average_sea_level_change"                       unit="m"        />
+    	 <field id="bgfrcvol"     long_name="drift in global mean volume from forcing wrt timestep 1"                                                                               unit="km3"      />
+    	 <field id="bgfrctem"     long_name="drift in global mean heat content from forcing wrt timestep 1"                                                                         unit="10^9J"    />
+    	 <field id="bgfrcsal"     long_name="drift in global mean salt content from forcing wrt timestep 1"                                                                         unit="1e-3*km3" />
+    	 <field id="bgmistem"     long_name="global mean temperature error due to free surface"                                                                                     unit="degC"     />
+    	 <field id="bgmissal"     long_name="global mean salinity error due to free surface"                                                                                        unit="1e-3"     />
+      </field_group>
+
+      <!-- LIM3 scalar variables -->
+
+      <field_group id="SBC_scalar"  domain_ref="1point" >
+         <!-- available with ln_limdiaout -->
+         <field id="ibgvoltot"    long_name="global mean ice volume"                                 unit="km3"        />
+         <field id="sbgvoltot"    long_name="global mean snow volume"                                unit="km3"        />
+         <field id="ibgarea"      long_name="global mean ice area"                                   unit="km2"        />
+         <field id="ibgsaline"    long_name="global mean ice salinity"                               unit="1e-3"       />
+         <field id="ibgtemper"    long_name="global mean ice temperature"                            unit="degC"       />
+         <field id="ibgheatco"    long_name="global mean ice heat content"                           unit="10^20J"     />
+         <field id="sbgheatco"    long_name="global mean snow heat content"                          unit="10^20J"     />
+         <field id="ibgsaltco"    long_name="global mean ice salt content"                           unit="1e-3*km3"   />
+
+         <field id="ibgvfx"       long_name="global mean volume flux (emp)"                          unit="m/day"      />
+         <field id="ibgvfxbog"    long_name="global mean volume flux (bottom growth)"                unit="m/day"      />
+         <field id="ibgvfxopw"    long_name="global mean volume flux (open water growth)"            unit="m/day"      />
+         <field id="ibgvfxsni"    long_name="global mean volume flux (snow-ice growth)"              unit="m/day"      />
+         <field id="ibgvfxdyn"    long_name="global mean volume flux (dynamic growth)"               unit="m/day"      />
+         <field id="ibgvfxbom"    long_name="global mean volume flux (bottom melt)"                  unit="m/day"      />
+         <field id="ibgvfxsum"    long_name="global mean volume flux (surface melt)"                 unit="m/day"      />
+         <field id="ibgvfxres"    long_name="global mean volume flux (resultant)"                    unit="m/day"      />
+         <field id="ibgvfxspr"    long_name="global mean volume flux (snow precip)"                  unit="m/day"      />
+         <field id="ibgvfxsnw"    long_name="global mean volume flux (snow melt)"                    unit="m/day"      />
+         <field id="ibgvfxsub"    long_name="global mean volume flux (snow sublimation)"             unit="m/day"      />
+
+         <field id="ibgsfx"       long_name="global mean salt flux (total)"                          unit="1e-3*m/day" />
+         <field id="ibgsfxbri"    long_name="global mean salt flux (brines)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxdyn"    long_name="global mean salt flux (dynamic)"                        unit="1e-3*m/day" />
+         <field id="ibgsfxres"    long_name="global mean salt flux (resultant)"                      unit="1e-3*m/day" />
+         <field id="ibgsfxbog"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxopw"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxsni"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxbom"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxsum"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+         <field id="ibgsfxsub"    long_name="global mean salt flux (thermo)"                         unit="1e-3*m/day" />
+
+         <field id="ibghfxdhc"    long_name="Heat content variation in snow and ice"                 unit="W"          />
+         <field id="ibghfxspr"    long_name="Heat content of snow precip"                            unit="W"          />
+
+         <field id="ibghfxthd"    long_name="heat fluxes from ice-ocean exchange during thermo"      unit="W"          />
+         <field id="ibghfxsum"    long_name="heat fluxes causing surface ice melt"                   unit="W"          />
+         <field id="ibghfxbom"    long_name="heat fluxes causing bottom ice melt"                    unit="W"          />
+         <field id="ibghfxbog"    long_name="heat fluxes causing bottom ice growth"                  unit="W"          />
+         <field id="ibghfxdif"    long_name="heat fluxes causing ice temperature change"             unit="W"          />
+         <field id="ibghfxopw"    long_name="heat fluxes causing open water ice formation"           unit="W"          />
+         <field id="ibghfxdyn"    long_name="heat fluxes from ice-ocean exchange during dynamic"     unit="W"          />
+         <field id="ibghfxres"    long_name="heat fluxes from ice-ocean exchange during resultant"   unit="W"          />
+         <field id="ibghfxsub"    long_name="heat fluxes from sublimation"                           unit="W"          />
+         <field id="ibghfxsnw"    long_name="heat fluxes from snow-ocean exchange"                   unit="W"          />
+         <field id="ibghfxout"    long_name="non solar heat fluxes received by the ocean"            unit="W"          />
+         <field id="ibghfxin"     long_name="total heat fluxes at the ice surface"                   unit="W"          />
+
+         <field id="ibgfrcvol"    long_name="global mean forcing volume (emp)"                       unit="km3"        />
+         <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"                       unit="1e-3*km3"   />
+         <field id="ibgvolgrm"    long_name="global mean ice growth+melt volume"                     unit="km3"        />
+      </field_group>
+  
+      <!-- variables available with key_float -->
+
+      <field_group id="floatvar" domain_ref="1point" axis_ref="nfloat" operation="instant" >
+        <field id="traj_lon"      long_name="floats longitude"                                                           unit="degrees_east"  />
+        <field id="traj_lat"      long_name="floats latitude"                                                            unit="degrees_north" />
+        <field id="traj_dep"      long_name="floats depth"                                                               unit="m"             />
+        <field id="traj_temp"     long_name="floats temperature"       standard_name="sea_water_potential_temperature"   unit="degC"          />
+        <field id="traj_salt"     long_name="floats salinity"          standard_name="sea_water_practical_salinity"      unit="1e-3"          />
+        <field id="traj_dens"     long_name="floats in-situ density"   standard_name="sea_water_density"                 unit="kg/m3"         />
+        <field id="traj_group"    long_name="floats group"                                                               unit="1"             />
+      </field_group>
+
+      <!-- variables available with iceberg trajectories -->
+
+      <field_group id="icbvar" domain_ref="grid_T"  > 
+        <field id="berg_melt"          long_name="icb melt rate of icebergs"                       unit="kg/m2/s"                    />
+        <field id="berg_buoy_melt"     long_name="icb buoyancy component of iceberg melt rate"     unit="kg/m2/s"                    />
+        <field id="berg_eros_melt"     long_name="icb erosion component of iceberg melt rate"      unit="kg/m2/s"                    />
+        <field id="berg_conv_melt"     long_name="icb convective component of iceberg melt rate"   unit="kg/m2/s"                    />
+        <field id="berg_virtual_area"  long_name="icb virtual coverage by icebergs"                unit="m2"                         />
+        <field id="bits_src"           long_name="icb mass source of bergy bits"                   unit="kg/m2/s"                    />
+        <field id="bits_melt"          long_name="icb melt rate of bergy bits"                     unit="kg/m2/s"                    />
+        <field id="bits_mass"          long_name="icb bergy bit density field"                     unit="kg/m2"                      />
+        <field id="berg_mass"          long_name="icb iceberg density field"                       unit="kg/m2"                      />
+        <field id="calving"            long_name="icb calving mass input"                          unit="kg/s"                       />
+        <field id="berg_floating_melt" long_name="icb melt rate of icebergs + bits"                unit="kg/m2/s"                    />
+        <field id="berg_real_calving"  long_name="icb calving into iceberg class"                  unit="kg/s"     axis_ref="icbcla" />
+        <field id="berg_stored_ice"    long_name="icb accumulated ice mass by class"               unit="kg"       axis_ref="icbcla" />
+      </field_group>
+
+      <!-- Poleward transport : ptr -->     
+      <field_group id="diaptr" domain_ref="ptr" > 
+        <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"       grid_ref="gznl_W_3D" />
+        <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"       grid_ref="gznl_W_3D" />
+        <field id="zomsfpac"          long_name="Meridional Stream-Function: Pacific"          unit="Sv"       grid_ref="gznl_W_3D" />
+        <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"       grid_ref="gznl_W_3D" />
+        <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"       grid_ref="gznl_W_3D" />
+        <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="degC"     grid_ref="gznl_T_3D" />
+        <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="degC"     grid_ref="gznl_T_3D" />
+        <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="degC"     grid_ref="gznl_T_3D" />
+        <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="degC"     grid_ref="gznl_T_3D" />
+        <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="degC"     grid_ref="gznl_T_3D" />
+        <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="1e-3"     grid_ref="gznl_T_3D" />
+        <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="1e-3"     grid_ref="gznl_T_3D" />
+        <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="1e-3"     grid_ref="gznl_T_3D" />
+        <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="1e-3"     grid_ref="gznl_T_3D" />
+        <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="1e-3"     grid_ref="gznl_T_3D" />
+        <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"       grid_ref="gznl_T_3D" />
+        <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"       grid_ref="gznl_T_3D" />
+        <field id="zosrfpac"          long_name="Zonal Mean Surface : Pacific"                 unit="m2"       grid_ref="gznl_T_3D" />
+        <field id="zosrfind"          long_name="Zonal Mean Surface : Indian"                  unit="m2"       grid_ref="gznl_T_3D" />
+        <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"       grid_ref="gznl_T_3D" />
+        <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" />
+        <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" />
+        <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" />
+        <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" />
+      </field_group>
+
+    <!-- 
+============================================================================================================
+                  Physical ocean model trend diagnostics : temperature, KE, PE, momentum
+============================================================================================================
+    -->
+
+    <field_group id="trendT" grid_ref="grid_T_3D">
+      <!-- variables available with ln_tra_trd -->
+      <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                                                                                          unit="degC/s"                        />
+      <field id="strd_xad"      long_name="salinity   -trend: i-advection"                                                                                          unit="1e-3/s"                        />
+      <field id="ttrd_yad"      long_name="temperature-trend: j-advection"                                                                                          unit="degC/s"                        />
+      <field id="strd_yad"      long_name="salinity   -trend: j-advection"                                                                                          unit="1e-3/s"                        />
+      <field id="ttrd_zad"      long_name="temperature-trend: k-advection"                                                                                          unit="degC/s"                        />
+      <field id="strd_zad"      long_name="salinity   -trend: k-advection"                                                                                          unit="1e-3/s"                        />
+      <field id="ttrd_ad"       long_name="temperature-trend: advection"               standard_name="tendency_of_sea_water_temperature_due_to_advection"           unit="degC/s"                         > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field>
+      <field id="strd_ad"       long_name="salinity   -trend: advection"               standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                         > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field>
+      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (no-vvl)"                                                                                unit="degC/s"   grid_ref="grid_T_2D" />
+      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (no-vvl)"                                                                                unit="1e-3/s"   grid_ref="grid_T_2D" />
+      <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing"   unit="degC/s"                        />
+      <field id="strd_ldf"      long_name="salinity   -trend: lateral  diffusion"      standard_name="tendency_of_sea_water_salinity_due_to_horizontal_mixing"      unit="1e-3/s"                        />
+      <field id="ttrd_zdf"      long_name="temperature-trend: vertical diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing"     unit="degC/s"                        />
+      <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"                        />
+
+      <!-- ln_traldf_iso=T only (iso-neutral diffusion) -->
+      <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"   unit="degC/s" />
+      <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"   unit="1e-3/s" />
+
+      <!-- -->
+      <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"        unit="degC/s" />
+      <field id="strd_dmp"      long_name="salinity   -trend: interior restoring"        unit="1e-3/s" />
+      <field id="ttrd_bbl"      long_name="temperature-trend: bottom boundary layer"     unit="degC/s" />
+      <field id="strd_bbl"      long_name="salinity   -trend: bottom boundary layer"     unit="1e-3/s" />
+      <field id="ttrd_npc"      long_name="temperature-trend: non-penetrative conv."     unit="degC/s" />
+      <field id="strd_npc"      long_name="salinity   -trend: non-penetrative conv."     unit="1e-3/s" />
+      <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"   unit="degC/s" />
+      <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"   unit="degC/s" />
+      <field id="ttrd_qsr"      long_name="temperature-trend: solar penetr. heating"     unit="degC/s" />
+      <field id="ttrd_bbc"      long_name="temperature-trend: geothermal heating"        unit="degC/s" />
+      <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"       unit="degC/s" />
+      <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"       unit="1e-3/s" />
+
+      <!-- variables available with ln_KE_trd -->
+      <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"          unit="W/s^3"                        />
+      <field id="ketrd_spg"     long_name="ke-trend: surface     pressure gradient"          unit="W/s^3"                        />
+      <field id="ketrd_spgexp"  long_name="ke-trend: surface pressure gradient (explicit)"   unit="W/s^3"                        />
+      <field id="ketrd_spgflt"  long_name="ke-trend: surface pressure gradient (filter)"     unit="W/s^3"                        />
+      <field id="ssh_flt"       long_name="filtered contribution to ssh (dynspg_flt)"        unit="m"       grid_ref="grid_T_2D" />
+      <field id="w0"            long_name="surface vertical velocity"                        unit="m/s"     grid_ref="grid_T_2D" />
+      <field id="pw0_exp"       long_name="surface pressure flux due to ssh"                 unit="W/s^2"   grid_ref="grid_T_2D" />
+      <field id="pw0_flt"       long_name="surface pressure flux due to filtered ssh"        unit="W/s^2"   grid_ref="grid_T_2D" />
+      <field id="ketrd_keg"     long_name="ke-trend: KE gradient         or hor. adv."       unit="W/s^3"                        />
+      <field id="ketrd_rvo"     long_name="ke-trend: relative  vorticity or metric term"     unit="W/s^3"                        />
+      <field id="ketrd_pvo"     long_name="ke-trend: planetary vorticity"                    unit="W/s^3"                        />
+      <field id="ketrd_zad"     long_name="ke-trend: vertical  advection"                    unit="W/s^3"                        />
+      <field id="ketrd_udx"     long_name="ke-trend: U.dx[U]"                                unit="W/s^3"                        />
+      <field id="ketrd_ldf"     long_name="ke-trend: lateral   diffusion"                    unit="W/s^3"                        />
+      <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                    unit="W/s^3"                        />
+      <field id="ketrd_tau"     long_name="ke-trend: wind stress "                           unit="W/s^3"   grid_ref="grid_T_2D" />
+      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        />   
+      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        />   
+      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        />  
+      <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        />
+      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />   
+
+      <!-- variables available with ln_PE_trd -->
+      <field id="petrd_xad"     long_name="pe-trend: i-advection"                unit="W/m^3"                        />
+      <field id="petrd_yad"     long_name="pe-trend: j-advection"                unit="W/m^3"                        />
+      <field id="petrd_zad"     long_name="pe-trend: k-advection"                unit="W/m^3"                        />
+      <field id="petrd_sad"     long_name="pe-trend: surface adv. (no-vvl)"      unit="W/m^3"   grid_ref="grid_T_2D" />
+      <field id="petrd_ldf"     long_name="pe-trend: lateral  diffusion"         unit="W/m^3"                        />
+      <field id="petrd_zdf"     long_name="pe-trend: vertical diffusion"         unit="W/m^3"                        />
+      <field id="petrd_zdfp"    long_name="pe-trend: pure vert. diffusion"       unit="W/m^3"                        />
+      <field id="petrd_dmp"     long_name="pe-trend: interior restoring"         unit="W/m^3"                        />
+      <field id="petrd_bbl"     long_name="pe-trend: bottom boundary layer"      unit="W/m^3"                        />
+      <field id="petrd_npc"     long_name="pe-trend: non-penetrative conv."      unit="W/m^3"                        />
+      <field id="petrd_nsr"     long_name="pe-trend: surface forcing + runoff"   unit="W/m^3"                        />
+      <field id="petrd_qsr"     long_name="pe-trend: solar penetr. heating"      unit="W/m^3"                        />
+      <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"         unit="W/m^3"                        />
+      <field id="petrd_atf"     long_name="pe-trend: asselin time filter"        unit="W/m^3"                        />
+      <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            />   
+      <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       />   
+      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          />   
+    </field_group>
+
+    <field_group id="trendU" grid_ref="grid_U_3D">
+     <!-- variables available with ln_dyn_trd -->
+     <field id="utrd_hpg"       long_name="i-trend: hydrostatic pressure gradient"          unit="m/s^2"                        />
+     <field id="utrd_spg"       long_name="i-trend: surface     pressure gradient"          unit="m/s^2"                        />
+     <field id="utrd_spgexp"    long_name="i-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        />
+     <field id="utrd_spgflt"    long_name="i-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        />
+     <field id="utrd_keg"       long_name="i-trend: KE gradient         or hor. adv."       unit="m/s^2"                        />
+     <field id="utrd_rvo"       long_name="i-trend: relative  vorticity or metric term"     unit="m/s^2"                        />
+     <field id="utrd_pvo"       long_name="i-trend: planetary vorticity"                    unit="m/s^2"                        />
+     <field id="utrd_zad"       long_name="i-trend: vertical  advection"                    unit="m/s^2"                        />
+     <field id="utrd_udx"       long_name="i-trend: U.dx[U]"                                unit="m/s^2"                        />
+     <field id="utrd_ldf"       long_name="i-trend: lateral   diffusion"                    unit="m/s^2"                        />
+     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                    unit="m/s^2"                        />
+     <field id="utrd_tau"       long_name="i-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_U_2D" />
+     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        />   
+     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        />   
+     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        />   
+     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        />   
+    </field_group>
+
+    <field_group id="trendV" grid_ref="grid_V_3D">
+     <!-- variables available with ln_dyn_trd -->
+     <field id="vtrd_hpg"       long_name="j-trend: hydrostatic pressure gradient"          unit="m/s^2"                        />
+     <field id="vtrd_spg"       long_name="j-trend: surface     pressure gradient"          unit="m/s^2"                        />
+     <field id="vtrd_spgexp"    long_name="j-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        />
+     <field id="vtrd_spgflt"    long_name="j-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        />
+     <field id="vtrd_keg"       long_name="j-trend: KE gradient         or hor. adv."       unit="m/s^2"                        />
+     <field id="vtrd_rvo"       long_name="j-trend: relative  vorticity or metric term"     unit="m/s^2"                        />
+     <field id="vtrd_pvo"       long_name="j-trend: planetary vorticity"                    unit="m/s^2"                        />
+     <field id="vtrd_zad"       long_name="j-trend: vertical  advection"                    unit="m/s^2"                        />
+     <field id="vtrd_vdy"       long_name="i-trend: V.dx[V]"                                unit="m/s^2"                        />
+     <field id="vtrd_ldf"       long_name="j-trend: lateral   diffusion"                    unit="m/s^2"                        />
+     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                    unit="m/s^2"                        />
+     <field id="vtrd_tau"       long_name="j-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_V_2D" />
+     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        />   
+     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        />   
+     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        />   
+     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        />   
+    </field_group>
+
+    <!-- 
+============================================================================================================
+                                     Biogeochemistry model variables
+============================================================================================================
+    -->
+
+      <!-- ptrc on T grid -->
+
+     <field_group id="ptrc_T" grid_ref="grid_T_3D">
+       <field id="DIC"          long_name="Dissolved inorganic Concentration"        unit="mmol/m3" />
+       <field id="DIC_E3T"      long_name="DIC * E3T"                                unit="mmol/m2" > DIC * e3t </field >
+       <field id="Alkalini"     long_name="Total Alkalinity Concentration"           unit="mmol/m3" />
+       <field id="Alkalini_E3T" long_name="Alkalini * E3T"                           unit="mmol/m2"  > Alkalini * e3t </field >
+       <field id="O2"           long_name="Oxygen Concentration"                     unit="mmol/m3" />
+       <field id="O2_E3T"       long_name="O2 * E3T"                                 unit="mmol/m2"  > O2 * e3t </field >
+       <field id="CaCO3"        long_name="Calcite Concentration"                    unit="mmol/m3" />
+       <field id="CaCO3_E3T"    long_name="CaCO3 * E3T"                              unit="mmol/m2"  > CaCO3 * e3t </field >
+       <field id="PO4"          long_name="Phosphate Concentration"                  unit="mmol/m3" />
+       <field id="PO4_E3T"      long_name="PO4 * E3T"                                unit="mmol/m2"  > PO4 * e3t </field >
+       <field id="POC"          long_name="Small organic carbon Concentration"       unit="mmol/m3" />
+       <field id="POC_E3T"      long_name="POC * E3T"                                unit="mmol/m2"  > POC * e3t </field >
+       <field id="Si"           long_name="Silicate Concentration"                   unit="mmol/m3" />
+       <field id="Si_E3T"       long_name="Si * E3T"                                 unit="mmol/m2"  > Si * e3t </field >
+       <field id="PHY"          long_name="(Nano)Phytoplankton Concentration"        unit="mmol/m3" />
+       <field id="PHY_E3T"      long_name="PHY * E3T"                                unit="mmol/m2"  > PHY * e3t </field >
+       <field id="ZOO"          long_name="(Micro)Zooplankton Concentration"         unit="mmol/m3" />
+       <field id="ZOO_E3T"      long_name="ZOO2 * E3T"                               unit="mmol/m2"  > ZOO * e3t </field >
+       <field id="DOC"          long_name="Dissolved organic Concentration"          unit="mmol/m3" />
+       <field id="DOC_E3T"      long_name="DOC * E3T"                                unit="mmol/m2"  > DOC * e3t </field >
+       <field id="PHY2"         long_name="Diatoms Concentration"                    unit="mmol/m3" />
+       <field id="PHY2_E3T"     long_name="PHY2 * E3T"                               unit="mmol/m2"  > PHY2 * e3t </field >
+       <field id="ZOO2"         long_name="Mesozooplankton Concentration"            unit="mmol/m3" />
+       <field id="ZOO2_E3T"     long_name="ZOO2 * E3T"                               unit="mmol/m2"  > ZOO2 * e3t </field >
+       <field id="DSi"          long_name="Diatoms Silicate Concentration"           unit="mmol/m3" />
+       <field id="DSi_E3T"      long_name="Dsi * E3T"                                unit="mmol/m2"  > DSi * e3t </field >
+       <field id="Fer"          long_name="Dissolved Iron Concentration"             unit="mmol/m3" />
+       <field id="Fer_E3T"      long_name="Fer * E3T"                                unit="mmol/m2"  > Fer * e3t </field >
+       <field id="BFe"          long_name="Big iron particles Concentration"         unit="mmol/m3" />
+       <field id="BFe_E3T"      long_name="BFe * E3T"                                unit="mmol/m2"  > BFe * e3t </field >
+       <field id="GOC"          long_name="Big organic carbon Concentration"         unit="mmol/m3" />
+       <field id="GOC_E3T"      long_name="GOC * E3T"                                unit="mmol/m2"  > GOC * e3t </field >
+       <field id="SFe"          long_name="Small iron particles Concentration"       unit="mmol/m3" />
+       <field id="SFe_E3T"      long_name="SFe * E3T"                                unit="mmol/m2"  > SFe * e3t </field >
+       <field id="DFe"          long_name="Diatoms iron  Concentration"              unit="mmol/m3" />
+       <field id="DFe_E3T"      long_name="DFe * E3T"                                unit="mmol/m2"  > DFe * e3t </field >
+       <field id="GSi"          long_name="Sinking biogenic Silicate Concentration"  unit="mmol/m3" />
+       <field id="GSi_E3T"      long_name="GSi * E3T"                                unit="mmol/m2"  > GSi * e3t </field >
+       <field id="NFe"          long_name="Nano iron Concentration"                  unit="mmol/m3" />
+       <field id="NFe_E3T"      long_name="NFe * E3T"                                unit="mmol/m2"  > NFe * e3t </field >
+       <field id="NCHL"         long_name="Nano chlorophyl Concentration"            unit="mg/m3"   />
+       <field id="NCHL_E3T"     long_name="NCHL * E3T"                               unit="mmol/m2"  > NCHL * e3t </field >
+       <field id="DCHL"         long_name="Diatoms chlorophyl Concentration"         unit="mg/m3"   />
+       <field id="DCHL_E3T"     long_name="DCHL * E3T"                               unit="mmol/m2"  > DCHL * e3t </field >
+       <field id="NO3"          long_name="Nitrate Concentration"                    unit="mmol/m3" />
+       <field id="NO3_E3T"      long_name="NO3 * E3T"                                unit="mmol/m2"  > NO3 * e3t </field >
+       <field id="NH4"          long_name="Ammonium Concentration"                   unit="mmol/m3" />
+       <field id="NH4_E3T"      long_name="NH4 * E3T"                                unit="mmol/m2"  > NH4 * e3t </field >
+
+
+       <!-- PISCES with N15 parametisation : variables available with key_n15 -->
+       <field id="POC_15"       long_name="N15 Small organic carbon Concentration"   unit="mmol/m3" />
+       <field id="PHY_15"       long_name="N15 (Nano)Phytoplankton Concentration"    unit="mmol/m3" />
+       <field id="ZOO_15"       long_name="N15 (Micro)Zooplankton Concentration"     unit="mmol/m3" />
+       <field id="DOC_15"       long_name="N15 Dissolved organic Concentration"      unit="mmol/m3" />
+       <field id="PHY2_15"      long_name="N15 Diatoms Concentration"                unit="mmol/m3" />
+       <field id="ZOO2_15"      long_name="N15 Mesozooplankton Concentration"        unit="mmol/m3" />
+       <field id="GOC_15"       long_name="N15 Big organic carbon Concentration"     unit="mmol/m3" />
+       <field id="NO3_15"       long_name="N15 Nitrate Concentration"                unit="mmol/m3" />
+       <field id="NO3_15_E3T"   long_name="NO3_15 * E3T"                             unit="mmol/m2"  > NO3_N15 * e3t </field >
+       <field id="NH4_15"       long_name="N15 Ammonium Concentration"               unit="mmol/m3" />
+
+       <!-- PISCES with C13 parametisation : variables available with key_c13 -->
+       <field id="POC_13"       long_name="C13 Small organic carbon Concentration"   unit="mmol/m3" />
+       <field id="PHY_13"       long_name="C13 (Nano)Phytoplankton Concentration"    unit="mmol/m3" />
+       <field id="ZOO_13"       long_name="C13 (Micro)Zooplankton Concentration"     unit="mmol/m3" />
+       <field id="DOC_13"       long_name="C13 Dissolved organic Concentration"      unit="mmol/m3" />
+       <field id="PHY2_13"      long_name="C13 Diatoms Concentration"                unit="mmol/m3" />
+       <field id="ZOO2_13"      long_name="C13 Mesozooplankton Concentration"        unit="mmol/m3" />
+       <field id="GOC_13"       long_name="C13 Big organic carbon Concentration"     unit="mmol/m3" />
+       <field id="DIC_13"       long_name="C13 Dissolved Inorganic Carbon Conc"      unit="mmol/m3" />
+       <field id="DIC_13_E3T"   long_name="DIC_13 * E3T"                             unit="mmol/m2"  > DIC_13 * e3t </field >
+       <field id="CaCO3_13"     long_name="C13 Calcite Concentration"                unit="mmol/m3" />
+
+       <!-- PISCES with Kriest parametisation : variables available with key_kriest -->
+       <field id="Num"         long_name="Number of organic particles"              unit="1" />
+       <field id="Num_E3T"     long_name="Num * E3T"                                unit="m"  > Num * e3t </field >
+
+       <!-- PISCES light : variables available with key_pisces_reduced -->
+       <field id="DET"         long_name="Detritus"                                 unit="mmol-N/m3" />
+       <field id="DET_E3T"     long_name="DET * E3T"                                unit="mmol-N/m2"  > DET * e3t </field >
+       <field id="DOM"         long_name="Dissolved Organic Matter"                 unit="mmol-N/m3" />
+       <field id="DOM_E3T"     long_name="DOM * E3T"                                unit="mmol-N/m2"  > DOM * e3t </field >
+
+       <!-- CFC11 : variables available with key_cfc -->
+       <field id="CFC11"       long_name="CFC-11 Concentration"                     unit="umol/m3" />
+       <field id="CFC11_E3T"   long_name="CFC11 * E3T"                              unit="umol/m2"  > CFC11 * e3t </field >
+       <!-- Bomb C14 : variables available with key_c14b -->
+       <field id="C14B"     long_name="Bomb C14 Concentration"                      unit="1"         />
+       <field id="C14B_E3T"    long_name="C14B * E3T"                               unit="m"  > C14B * e3t </field >
+     </field_group>
+
+     <!-- PISCES additional diagnostics on T grid  -->
+
+     <field_group id="diad_T" grid_ref="grid_T_2D">
+       <field id="PH"          long_name="PH"                                      unit="1"          grid_ref="grid_T_3D" />
+       <field id="CO3"         long_name="Bicarbonates"                            unit="mol/m3"     grid_ref="grid_T_3D" />
+       <field id="CO3sat"      long_name="CO3 saturation"                          unit="mol/m3"     grid_ref="grid_T_3D" />
+       <field id="PAR"         long_name="Photosynthetically Available Radiation"  unit="W/m2"       grid_ref="grid_T_3D" />
+       <field id="PARDM"       long_name="Daily mean PAR"                          unit="W/m2"       grid_ref="grid_T_3D" />
+       <field id="PPPHY"       long_name="Primary production of nanophyto"         unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PPPHY2"      long_name="Primary production of diatoms"           unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PPNEWN"      long_name="New Primary production of nanophyto"     unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PPNEWD"      long_name="New Primary production of diatoms"       unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PBSi"        long_name="Primary production of Si diatoms"        unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PFeN"        long_name="Primary production of nano iron"         unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="PFeD"        long_name="Primary production of diatoms iron"      unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="xfracal"     long_name="Calcifying fraction"                     unit="1"          grid_ref="grid_T_3D" />
+       <field id="PCAL"        long_name="Calcite production"                      unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="DCAL"        long_name="Calcite dissolution"                     unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="GRAZ1"       long_name="Grazing by microzooplankton"             unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="GRAZ2"       long_name="Grazing by mesozooplankton"              unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="REMIN"       long_name="Oxic remineralization of OM"             unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="DENIT"       long_name="Anoxic remineralization of OM"           unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="Nfix"        long_name="Nitrogen fixation"                       unit="mol/m3/s"   grid_ref="grid_T_3D" />
+       <field id="Mumax"       long_name="Maximum growth rate"                     unit="s-1"        grid_ref="grid_T_3D" />
+       <field id="MuN"         long_name="Realized growth rate for nanophyto"      unit="s-1"        grid_ref="grid_T_3D" />
+       <field id="MuD"         long_name="Realized growth rate for diatomes"       unit="s-1"        grid_ref="grid_T_3D" />
+       <field id="LNnut"       long_name="Nutrient limitation term in Nanophyto"   unit=""           grid_ref="grid_T_3D" />
+       <field id="LDnut"       long_name="Nutrient limitation term in Diatoms"     unit=""           grid_ref="grid_T_3D" />
+       <field id="LNFe"        long_name="Iron limitation term in Nanophyto"       unit=""           grid_ref="grid_T_3D" />
+       <field id="LDFe"        long_name="Iron limitation term in Diatoms"         unit=""           grid_ref="grid_T_3D" />
+       <field id="LNlight"     long_name="Light limitation term in Nanophyto"      unit=""           grid_ref="grid_T_3D" />
+       <field id="LDlight"     long_name="Light limitation term in Diatoms"        unit=""           grid_ref="grid_T_3D" />
+       <field id="Fe2"         long_name="Iron II concentration"                   unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="Fe3"         long_name="Iron III concentration"                  unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="FeL1"        long_name="Complexed Iron concentration with L1"    unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="FeL2"        long_name="Complexed Iron concentration with L2"    unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="FeP"         long_name="Precipitated Iron III"                   unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="TL1"         long_name="Total L1 concentration"                  unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="TL2"         long_name="Total L2 concentration"                  unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="pdust"       long_name="dust concentration"                      unit="g/m3"                            />
+       <field id="Totlig"      long_name="Total ligand concentation"               unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="Biron"       long_name="Bioavailable iron"                       unit="nmol/m3"    grid_ref="grid_T_3D" />
+       <field id="Sdenit"      long_name="Nitrate reduction in the sediments"      unit="molN/m2/s"                       />
+       <field id="SedCal"      long_name="Calcite burial in the sediments"         unit="molC/m2/s"                       />
+       <field id="SedSi"       long_name="Silicon burial in the sediments"         unit="molSi/m2/s"                      />
+       <field id="SedC"        long_name="Organic C burial in the sediments"       unit="molC/m2/s"                       />
+       <field id="Ironice"     long_name="Iron input/uptake due to sea ice"        unit="mol/m2/s"                        />
+       <field id="HYDR"        long_name="Iron input from hydrothemal vents"       unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="EPC100"      long_name="Export of carbon particles at 100 m"     unit="mol/m2/s"                        />
+       <field id="EPFE100"     long_name="Export of biogenic iron at 100 m"        unit="mol/m2/s"                        />
+       <field id="EPSI100"     long_name="Export of Silicate at 100 m"             unit="mol/m2/s"                        />
+       <field id="EPCAL100"    long_name="Export of Calcite at 100 m"              unit="mol/m2/s"                        />
+       <field id="EXPC"        long_name="Export of carbon"                        unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="EXPFE"       long_name="Export of biogenic iron"                 unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="EXPSI"       long_name="Export of Silicate"                      unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="EXPCAL"      long_name="Export of Calcite"                       unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="Cflx"        long_name="DIC flux"                                unit="mol/m2/s"                        />
+       <field id="C13flx"      long_name="DIC13 flux"                              unit="mol/m2/s"                        />
+       <field id="Oflx"        long_name="Oxygen flux"                             unit="mol/m2/s"                        />
+       <field id="Kg"          long_name="Gas transfer"                            unit="mol/m2/s/uatm"                   />
+       <field id="Dpco2"       long_name="Delta CO2"                               unit="uatm"                            />
+       <field id="Dpo2"        long_name="Delta O2"                                unit="uatm"                            />
+       <field id="Heup"        long_name="Euphotic layer depth"                    unit="m"                               />
+       <field id="Irondep"     long_name="Iron deposition from dust"               unit="mol/m2/s"                        />
+       <field id="Ironsed"     long_name="Iron deposition from sediment"           unit="mol/m2/s"   grid_ref="grid_T_3D" />
+
+
+       <!-- PISCES with Kriest parametisation : variables available with key_kriest -->
+       <field id="EPN100"      long_name="Particulate number flux at 100 m"        unit="mol/m2/s"                        />
+       <field id="EXPN"        long_name="Particulate number flux"                 unit="mol/m2/s"   grid_ref="grid_T_3D" />
+       <field id="XNUM"        long_name="Number of particles in aggregats"        unit="1"          grid_ref="grid_T_3D" />
+       <field id="WSC"         long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" />
+       <field id="WSN"         long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" />
+
+       <!-- dbio_T on T grid : variables available with key_diaar5 -->
+       <field id="TPP"         long_name="Total Primary production of phyto"                   unit="mol/m3/s"  grid_ref="grid_T_3D" />
+       <field id="TPNEW"       long_name="New Primary production of phyto"                     unit="mol/m3/s"  grid_ref="grid_T_3D" />
+       <field id="TPBFE"       long_name="Total biogenic iron production"                      unit="mol/m3/s"  grid_ref="grid_T_3D" />
+       <field id="INTDIC"      long_name="DIC content"                                         unit="kg/m2"                          />
+       <field id="O2MIN"       long_name="Oxygen minimum concentration"                        unit="mol/m3"                          />
+       <field id="ZO2MIN"      long_name="Depth of oxygen minimum concentration"               unit="m"                              />
+       <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="mol/m2/s"                       />
+       <field id="INTPPPHY"    long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"                       />
+       <field id="INTPPPHY2"   long_name="Vertically integrated primary production by diatom"  unit="mol/m2/s"                       />
+       <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="mol/m2/s"                       />
+       <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="mol/m2/s"                       />
+       <field id="INTPBFE"     long_name="Vertically integrated of biogenic iron production"   unit="mol/m2/s"                       />
+       <field id="INTPBSI"     long_name="Vertically integrated of biogenic Si production"     unit="mol/m2/s"                       />
+       <field id="INTPCAL"     long_name="Vertically integrated of calcite production"         unit="mol/m2/s"                       />
+
+       <!-- PISCES light : variables available with key_pisces_reduced -->
+       <field id="FNO3PHY"     long_name="FNO3PHY"                                 unit=""          grid_ref="grid_T_3D" /> 
+       <field id="FNH4PHY"     long_name="FNH4PHY"                                 unit=""          grid_ref="grid_T_3D" /> 
+       <field id="FNH4NO3"     long_name="FNH4NO3"                                 unit=""          grid_ref="grid_T_3D" /> 
+       <field id="TNO3PHY"     long_name="TNO3PHY"                                 unit=""  /> 
+       <field id="TNH4PHY"     long_name="TNH4PHY"                                 unit=""  /> 
+       <field id="TPHYDOM"     long_name="TPHYDOM"                                 unit=""  /> 
+       <field id="TPHYNH4"     long_name="TPHYNH4"                                 unit=""  /> 
+       <field id="TPHYZOO"     long_name="TPHYZOO"                                 unit=""  /> 
+       <field id="TPHYDET"     long_name="TPHYDET"                                 unit=""  /> 
+       <field id="TDETZOO"     long_name="TDETZOO"                                 unit=""  /> 
+       <field id="TZOODET"     long_name="TZOODET"                                 unit=""  /> 
+       <field id="TZOOBOD"     long_name="TZOOBOD"                                 unit=""  /> 
+       <field id="TZOONH4"     long_name="TZOONH4"                                 unit=""  /> 
+       <field id="TZOODOM"     long_name="TZOODOM"                                 unit=""  /> 
+       <field id="TNH4NO3"     long_name="TNH4NO3"                                 unit=""  /> 
+       <field id="TDOMNH4"     long_name="TDOMNH4"                                 unit=""  /> 
+       <field id="TDETNH4"     long_name="TDETNH4"                                 unit=""  /> 
+       <field id="TPHYTOT"     long_name="TPHYTOT"                                 unit=""  /> 
+       <field id="TZOOTOT"     long_name="TZOOTOT"                                 unit=""  /> 
+       <field id="SEDPOC"      long_name="SEDPOC"                                  unit=""  /> 
+       <field id="TDETSED"     long_name="TDETSED"                                 unit=""  /> 
+
+       <!-- CFC11 : variables available with key_cfc -->
+       <field id="qtrCFC11"    long_name="Air-sea flux of CFC-11"                  unit="mol/m2/s"   />
+       <field id="qintCFC11"   long_name="Cumulative air-sea flux of CFC-11"       unit="mol/m2"     />
+
+       <!-- Bomb C14 : variables available with key_c14b -->
+       <field id="qtrC14b"     long_name="Air-sea flux of Bomb C14"                unit="mol/m2/s"   />
+       <field id="qintC14b"    long_name="Cumulative air-sea flux of Bomb C14"     unit="mol/m2"     />
+       <field id="fdecay"      long_name="Radiactive decay of Bomb C14"            unit="mol/m3"  grid_ref="grid_T_3D"  />
+     </field_group>
+
+     <field_group id="PISCES_scalar"  domain_ref="1point" >
+       <field id="pno3tot"         long_name="global mean nitrate concentration"                  unit="mol/m3"   />
+       <field id="ppo4tot"         long_name="global mean phosphorus concentration"               unit="mol/m3"   />
+       <field id="psiltot"         long_name="global mean silicate concentration"                 unit="mol/m3"   />
+       <field id="palktot"         long_name="global mean alkalinity concentration"               unit="mol/m3"   />
+       <field id="pfertot"         long_name="global mean iron concentration"                     unit="mol/m3"   />
+       <field id="tcflx"           long_name="total Flux of Carbon out of the ocean"              unit="mol/s"   />
+       <field id="tcflxcum"        long_name="cumulative total Flux of Carbon out of the ocean"   unit="mol/s"   />
+       <field id="tcexp"           long_name="total Carbon export at 100m"                        unit="mol/s"   />
+       <field id="tintpp"          long_name="global total integrated primary production"         unit="mol/s"   />
+       <field id="tnfix"           long_name="global total nitrogen fixation"                     unit="mol/s"   />
+       <field id="tdenit"          long_name="Total denitrification"                              unit="mol/s"   />
+     </field_group>
+
+    <!-- 
+============================================================================================================
+                                        Definitions for iodef_demo.xml
+============================================================================================================
+    -->
+
+    <field_group id="mooring" >
+      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"      />
+      <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                   />
+      <field field_ref="uoce"         name="uo"       long_name="sea_water_x_velocity"                 />
+      <field field_ref="voce"         name="vo"       long_name="sea_water_y_velocity"                 />
+      <field field_ref="woce"         name="wo"       long_name="sea_water_z_velocity"                 />
+      <field field_ref="avt"          name="difvho"   long_name="ocean_vertical_heat_diffusivity"      />
+      <field field_ref="avm"          name="difvmo"   long_name="ocean_vertical_momentum_diffusivity"  />
+   
+      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       />
+      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             />
+      <field field_ref="sstgrad"      name="tosgrad"  long_name="module_of_sea_surface_temperature_gradient"    />
+      <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"                          />
+      <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                />
+      <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     />
+      <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           />
+      <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          />
+      <field field_ref="taum"                                                                                   />
+      <field field_ref="20d"                                                                                    />
+      <field field_ref="mldkz5"                                                                                 />
+      <field field_ref="mldr10_1"                                                                               />
+      <field field_ref="mldr10_3"                                                                               />
+      <field field_ref="mldr0_1"                                                                                />
+      <field field_ref="mldr0_3"                                                                                />
+      <field field_ref="mld_dt02"                                                                               />
+      <field field_ref="topthdep"                                                                               />
+      <field field_ref="pycndep"                                                                                />
+      <field field_ref="tinv"                                                                                   />
+      <field field_ref="depti"                                                                                  />
+      <field field_ref="BLT"          name="blt"      long_name="barrier_layer_thickness"                       />
+      <field field_ref="utau"         name="tauuo"    long_name="surface_downward_x_stress"                     />
+      <field field_ref="vtau"         name="tauvo"    long_name="surface_downward_y_stress"                     />
+    </field_group>
+
+    <field_group id="groupT" >
+      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"               />
+      <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                            />
+      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       />
+      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             />
+      <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"                          />
+      <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                />
+      <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     />
+      <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           />
+      <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          />
+      <field field_ref="taum"                                                                                   />
+      <field field_ref="20d"                                                                                    />
+      <field field_ref="mldkz5"                                                                                 />
+      <field field_ref="mldr10_1"                                                                               />
+      <field field_ref="mldr10_3"                                                                               />
+      <field field_ref="mld_dt02"                                                                               />
+      <field field_ref="topthdep"                                                                               />
+      <field field_ref="pycndep"                                                                                />
+      <field field_ref="tinv"                                                                                   />
+      <field field_ref="depti"                                                                                  />
+      <field field_ref="BLT"          name="blt"      long_name="Barrier Layer Thickness"                       />
+    </field_group>
+    
+    <field_group id="groupU" >
+      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      />
+      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    />
+      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" />
+    </field_group>
+    
+    <field_group id="groupV" >
+      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      />
+      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    />
+      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" />
+    </field_group>
+    
+    <field_group id="groupW" >
+      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  />
+    </field_group>
+    
+
+   </field_definition>
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/iodef.xml
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/iodef.xml	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/iodef.xml	(revision 6534)
@@ -0,0 +1,237 @@
+<?xml version="1.0"?>
+<simulation> 
+
+ <context id="nemo" time_origin="1900-01-01 00:00:00" >
+    
+    <!-- $id$ -->
+    
+    <!-- 
+============================================================================================================
+=                                  definition of all existing variables                                    =
+=                                            DO NOT CHANGE                                                 =
+============================================================================================================
+    -->
+    <field_definition src="./field_def.xml"/>
+    <!-- 
+============================================================================================================
+=                                           output files definition                                        =
+=                                            Define your own files                                         =
+=                                         put the variables you want...                                    =
+============================================================================================================
+    -->
+
+   <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="10d" min_digits="4">
+    
+      <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
+      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> <!-- 1h files -->
+      <file_group id="2h" output_freq="2h"  output_level="10" enabled=".TRUE."/> <!-- 2h files -->
+      <file_group id="3h" output_freq="3h"  output_level="10" enabled=".TRUE."/> <!-- 3h files -->     
+      <file_group id="4h" output_freq="4h"  output_level="10" enabled=".TRUE."/> <!-- 4h files -->
+      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files -->
+     
+      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."> <!-- 1d files -->
+        <file id="file1" name_suffix="_bioscalar" description="pisces sms variables" >
+           <field field_ref="tdenit"   name="tdenit"   unit="TgN/yr" operation="instant" > tdenit * 14. * 86400. * 365. / 1e12 </field>
+           <field field_ref="tnfix"    name="tnfix"    unit="TgN/yr" operation="instant" > tnfix * 14. * 86400. * 365. / 1e12 </field>
+           <field field_ref="tcflx"    name="tcflx"    unit="PgC/yr" operation="instant" > tcflx * -1. * 12. * 86400. * 365. / 1e15 </field>
+           <field field_ref="tcflxcum" name="tcflxcum" unit="PgC"    operation="instant" > tcflxcum * -1. * 12. / 1e15 </field>
+           <field field_ref="tcexp"    name="tcexp"    unit="PgC/yr" operation="instant" > tcexp * 12. * 86400. * 365. / 1e15 </field>
+           <field field_ref="tintpp"   name="tintpp"   unit="PgC/yr" operation="instant" > tintpp * 12. * 86400. * 365. / 1e15 </field>
+           <field field_ref="pno3tot"  name="pno3tot"  unit="umolN"  > pno3tot * 16. / 122. * 1e6 </field>
+           <field field_ref="ppo4tot"  name="ppo4tot"  unit="umolP"  > ppo4tot * 1. / 122. * 1e6 </field>
+           <field field_ref="psiltot"  name="psiltot"  unit="umolC"  > psiltot * 1e6  </field>
+           <field field_ref="palktot"  name="palktot"  unit="umolC"  > palktot * 1e6  </field>
+           <field field_ref="pfertot"  name="pfertot"  unit="nmolFe" > pfertot * 1e9  </field>
+        </file>
+      </file_group>
+
+      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->    
+      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE."/>  <!-- 5d files -->   
+
+      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files -->
+
+	<file id="file2" name_suffix="_ptrc_T" description="pisces sms variables" >
+          <field field_ref="DIC"      />
+          <field field_ref="Alkalini" />
+          <field field_ref="O2"       />
+          <field field_ref="PO4"      />
+          <field field_ref="Si"       />
+          <field field_ref="Fer"      />
+          <field field_ref="NCHL"     />
+          <field field_ref="DCHL"     />
+          <field field_ref="NO3"      />
+          <field field_ref="DIC_13"   />
+          <field field_ref="NO3_15"   />
+	</file>
+	
+	<file id="file3" name_suffix="_diad_T" description="additional pisces diagnostics" >
+          <field field_ref="Cflx"     />
+          <field field_ref="Dpco2"    />
+          <field field_ref="C13flx"     />
+	</file>
+
+      </file_group>
+      <file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
+      <file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
+      <file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
+      <file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
+
+      <file_group id="1y"  output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files -->
+
+	<file id="file4" name_suffix="_ptrc_T" description="pisces sms variables" >
+          <field field_ref="DIC"      />
+          <field field_ref="Alkalini" />
+          <field field_ref="O2"       />
+          <field field_ref="CaCO3"    />
+          <field field_ref="PO4"      />
+          <field field_ref="POC"      />
+          <field field_ref="Si"       />
+          <field field_ref="PHY"      />
+          <field field_ref="ZOO"      />
+          <field field_ref="DOC"      />
+          <field field_ref="PHY2"     />
+          <field field_ref="ZOO2"     />
+          <field field_ref="DSi"      />
+          <field field_ref="Fer"      />
+          <field field_ref="BFe"      />
+          <field field_ref="GOC"      />
+          <field field_ref="SFe"      />
+          <field field_ref="DFe"      />
+          <field field_ref="GSi"      />
+          <field field_ref="NFe"      />
+          <field field_ref="NCHL"     />
+          <field field_ref="DCHL"     />
+          <field field_ref="NO3"      />
+          <field field_ref="NH4"      />
+          <field field_ref="DIC_13"   />
+          <field field_ref="CaCO3_13"      />
+          <field field_ref="NO3_15"   />
+          <field field_ref="PHY_15"   />
+          <field field_ref="PHY_13"   />
+          <field field_ref="POC_15"      />
+          <field field_ref="POC_13"      />
+          <field field_ref="ZOO_15"      />
+          <field field_ref="ZOO_13"      />
+          <field field_ref="DOC_15"      />
+          <field field_ref="DOC_13"      />
+          <field field_ref="PHY2_15"     />
+          <field field_ref="PHY2_13"     />
+          <field field_ref="ZOO2_15"     />
+          <field field_ref="ZOO2_13"     />
+          <field field_ref="GOC_15"      />
+          <field field_ref="GOC_13"      />
+          <field field_ref="NH4_15"      />
+	</file>
+
+	<file id="file5" name_suffix="_diad_T" description="additional pisces diagnostics" >
+          <field field_ref="PH"       />
+          <field field_ref="CO3"      />
+          <field field_ref="CO3sat"   />
+          <field field_ref="PAR"      />
+          <field field_ref="PPPHY"    />
+          <field field_ref="PPPHY2"   />
+          <field field_ref="PPNEWN"   />
+          <field field_ref="PPNEWD"   />
+          <field field_ref="PBSi"     />
+          <field field_ref="PFeN"     />
+          <field field_ref="PFeD"     />
+          <field field_ref="xfracal"  />
+          <field field_ref="PCAL"     />
+          <field field_ref="DCAL"     />
+          <field field_ref="GRAZ1"    />
+          <field field_ref="GRAZ2"    />
+          <field field_ref="EPC100"   />
+          <field field_ref="EPFE100"  />
+          <field field_ref="EPSI100"  />
+          <field field_ref="EPCAL100" />
+          <field field_ref="Cflx"     />
+          <field field_ref="C13flx"     />
+          <field field_ref="Oflx"     />
+          <field field_ref="Kg"       />
+          <field field_ref="Dpco2"    />
+          <field field_ref="Dpo2"     />
+          <field field_ref="Heup"     />
+          <field field_ref="Irondep"  />
+          <field field_ref="Ironsed"  />
+          <field field_ref="Ironice"  />
+          <field field_ref="Nfix"     />
+          <field field_ref="MuN"      />
+          <field field_ref="MuD"      />
+          <field field_ref="LNnut"    />
+          <field field_ref="LDnut"    />
+          <field field_ref="LNFe"     />
+          <field field_ref="LDFe"     />
+          <field field_ref="LNlight"  />
+          <field field_ref="LDlight"  />
+          <field field_ref="pdust"    />
+          <field field_ref="Fe2"      />
+          <field field_ref="Fe3"      />
+          <field field_ref="FeL1"     />
+          <field field_ref="FeL2"     />
+          <field field_ref="FeP"      />
+          <field field_ref="TL1"      />
+          <field field_ref="TL2"      />
+          <field field_ref="Sdenit"   />
+          <field field_ref="Totlig"   />
+	</file>
+
+      </file_group>
+      <file_group id="2y"  output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
+      <file_group id="5y"  output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
+      <file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
+
+   </file_definition>
+    
+    <!-- 
+============================================================================================================
+= grid definition = = DO NOT CHANGE =
+============================================================================================================
+    -->
+    
+   <axis_definition>  
+      <axis id="deptht"  long_name="Vertical T levels"  unit="m" positive="down" />
+      <axis id="depthu"  long_name="Vertical U levels"  unit="m" positive="down" />
+      <axis id="depthv"  long_name="Vertical V levels"  unit="m" positive="down" />
+      <axis id="depthw"  long_name="Vertical W levels"  unit="m" positive="down" />
+      <axis id="nfloat"  long_name="Float number"       unit="1"                 />
+      <axis id="icbcla"  long_name="Iceberg class"      unit="1"                 />
+      <axis id="ncatice" long_name="Ice category"       unit="1"                 />
+      <axis id="iax_20C" long_name="20 degC isotherm"   unit="degC"              />
+      <axis id="iax_28C" long_name="28 degC isotherm"   unit="degC"              />
+   </axis_definition> 
+    
+   <domain_definition src="./domain_def.xml"/>
+   
+   <grid_definition>    
+     <grid id="grid_T_2D" domain_ref="grid_T"/>
+     <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/>
+     <grid id="grid_U_2D" domain_ref="grid_U"/>
+     <grid id="grid_U_3D" domain_ref="grid_U" axis_ref="depthu"/>
+     <grid id="grid_V_2D" domain_ref="grid_V"/>
+     <grid id="grid_V_3D" domain_ref="grid_V" axis_ref="depthv"/>
+     <grid id="grid_W_2D" domain_ref="grid_W"/>
+     <grid id="grid_W_3D" domain_ref="grid_W" axis_ref="depthw"/>
+    </grid_definition>    
+  
+  </context>
+  
+
+  <context id="xios">
+
+      <variable_definition>
+	
+     <!-- 
+        We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size)
+-->
+	  <variable id="buffer_size"               type="integer">50000000</variable>
+	  <variable id="buffer_server_factor_size" type="integer">2</variable>
+	  <variable id="info_level"                type="integer">0</variable>
+	  <variable id="using_server"              type="boolean">false</variable>
+	  <variable id="using_oasis"               type="boolean">false</variable>
+	  <variable id="oasis_codes_id"            type="string" >oceanx</variable>
+	
+      </variable_definition>
+               
+  </context>
+  
+</simulation>
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_cfg
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_cfg	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_cfg	(revision 6534)
@@ -0,0 +1,136 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA  :  ORCA2_OFF_PISCES configuration namelist used ot overwrite SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!-----------------------------------------------------------------------
+&namrun        !   parameters of the run
+!-----------------------------------------------------------------------
+   cn_exp      =  "PISCES"  !  experience name 
+   nn_it000    =       1    !  first time step
+   nn_itend    =    1460    !  last  time step (std 5475)
+   nn_stock    =    1460    !  frequency of creation of a restart file (modulo referenced to 1)
+/
+!-----------------------------------------------------------------------
+&namcfg        !   parameters of the configuration
+!-----------------------------------------------------------------------
+   cp_cfg      =  "orca"               !  name of the configuration
+   jp_cfg      =       2               !  resolution of the configuration
+   jpidta      =     182               !  1st lateral dimension ( >= jpi )
+   jpjdta      =     149               !  2nd    "         "    ( >= jpj )
+   jpkdta      =      31               !  number of levels      ( >= jpk )
+   jpiglo      =     182               !  1st dimension of global domain --> i =jpidta
+   jpjglo      =     149               !  2nd    -                  -    --> j  =jpjdta
+   jpizoom     =       1               !  left bottom (i,j) indices of the zoom
+   jpjzoom     =       1               !  in data domain indices
+   jperio      =       4               !  lateral cond. type (between 0 and 6)
+/
+!-----------------------------------------------------------------------
+&namzgr        !   vertical coordinate
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namdom        !   space and time domain (bathymetry, mesh, timestep)
+!-----------------------------------------------------------------------
+   nn_msh      =    1      !  create (=1) a mesh file or not (=0)
+   rn_rdt      = 21600.    !  time step for the dynamics (and tracer if nn_acc=0)
+   rn_rdtmin   = 21600.          !  minimum time step on tracers (used if nn_acc=1)
+   rn_rdtmax   = 21600.          !  maximum time step on tracers (used if nn_acc=1)
+   jphgr_msh   =       0               !  type of horizontal mesh
+   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1)
+   ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1)
+   ppe1_deg    =  999999.0             !  zonal      grid-spacing (degrees)
+   ppe2_deg    =  999999.0             !  meridional grid-spacing (degrees)
+   ppe1_m      =  999999.0             !  zonal      grid-spacing (degrees)
+   ppe2_m      =  999999.0             !  meridional grid-spacing (degrees)
+   ppsur       =   -4762.96143546300   !  ORCA r4, r2 and r05 coefficients
+   ppa0        =     255.58049070440   ! (default coefficients)
+   ppa1        =     245.58132232490   !
+   ppkth       =      21.43336197938   !
+   ppacr       =       3.0             !
+   ppdzmin     =  999999.              !  Minimum vertical spacing
+   pphmax      =  999999.              !  Maximum depth
+   ldbletanh   =  .FALSE.              !  Use/do not use double tanf function for vertical coordinates
+   ppa2        =  999999.              !  Double tanh function parameters
+   ppkth2      =  999999.              !
+   ppacr2      =  999999.              !
+/
+!-----------------------------------------------------------------------
+&namsplit      !   time splitting parameters                            ("key_dynspg_ts")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namcrs        !   Grid coarsening for dynamics output and/or
+               !   passive tracer coarsened online simulations
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc        !   Surface Boundary Condition (surface module)
+!-----------------------------------------------------------------------
+   nn_fsbc     = 1         !  frequency of surface boundary condition computation 
+   ln_rnf      = .false.   !  runoffs
+/  
+!-----------------------------------------------------------------------
+&namtra_qsr    !   penetrative solar radiation
+!-----------------------------------------------------------------------
+   ln_traqsr   = .false.   !  Light penetration (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namsbc_ssr    !   surface boundary condition : sea surface restoring
+!-----------------------------------------------------------------------
+/      
+
+!-----------------------------------------------------------------------
+&namcla        !   cross land advection
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nambbl        !   bottom boundary layer scheme
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nameos        !   ocean physical parameters
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtra_ldf    !   lateral diffusion scheme for tracer 
+!-----------------------------------------------------------------------
+   ln_triad_iso     =  .true.  !  griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp")
+/
+!-----------------------------------------------------------------------
+&namdta_dyn        !   offline dynamics read in files                ("key_offline")
+!-----------------------------------------------------------------------
+!          !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
+!          !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
+   sn_tem  = 'dyna_grid_T' ,    120            , 'votemper' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_sal  = 'dyna_grid_T' ,    120            , 'vosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_mld  = 'dyna_grid_T' ,    120            , 'somixhgt' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_emp  = 'dyna_grid_T' ,    120            , 'sowaflup' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_fmf  = 'dyna_grid_T' ,    120            , 'iowaflup' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_ice  = 'dyna_grid_T' ,    120            , 'soicecov' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_qsr  = 'dyna_grid_T' ,    120            , 'soshfldo' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_wnd  = 'dyna_grid_T' ,    120            , 'sowindsp' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_uwd  = 'dyna_grid_U' ,    120            , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_vwd  = 'dyna_grid_V' ,    120            , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_wwd  = 'dyna_grid_W' ,    120            , 'vovecrtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_avt  = 'dyna_grid_W' ,    120            , 'voddmavs' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_ahu  = 'dyna_grid_U' ,    120            , 'vozoahtu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_ahv  = 'dyna_grid_V' ,    120            , 'vomeahtv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_ahw  = 'dyna_grid_W' ,    120            , 'voveahtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_eiu  = 'dyna_grid_U' ,    120            , 'vozoaeiu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_eiv  = 'dyna_grid_V' ,    120            , 'vomeaeiv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+   sn_eiw  = 'dyna_grid_W' ,    120            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , ''
+!
+   cn_dir      = './'      !  root directory for the location of the dynamical files
+   ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad")
+   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file
+   ln_dynbbl   =  .true.   !  bbl coef are in files, so read them - requires ("key_trabbl")
+/
+!-----------------------------------------------------------------------
+&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl        !   Control prints & Benchmark
+!-----------------------------------------------------------------------
+/
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_cfg
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_cfg	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_cfg	(revision 6534)
@@ -0,0 +1,87 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! PISCES  :   ORCA2_OFF_PISCES configuration namelsit used to overwrite SHARED/namelist_pisces_ref
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisext     !   air-sea exchange
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+   d13c_co2  =  -6.5      ! d13C of atm. CO2
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisatm     !  Atmospheric prrssure 
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisbio     !   biological parameters
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+   nrdttrc    =  4        ! time step frequency for biology
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampislim     !   parameters for nutrient limitations
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisopt     !   parameters for optics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/ 
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisprod     !   parameters for phytoplankton growth
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+    e15n_prod  = 5      ! N15 assimilation fractionation
+    e13c_min  = 19      ! C13 assimilation fractionation min
+    e13c_max =  25      ! C13 assimilation fractionation max
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismort     !   parameters for phytoplankton sinks
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+    e13c_cal = 2.    ! C13 phyto calcification fractionation
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismes     !   parameters for mesozooplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+    e15n_ex2 =  6.   ! N15 mesozoo excretion fractionation
+    e15n_ing2 = -1.   ! N15 mesozoo ingestion fractionation
+    e13c_cal2 = 2.    ! C13 mesozoo calcification fractionation
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiszoo     !   parameters for microzooplankton
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+    e15n_ex = 6.     ! N15 microzoo excretion fractionation
+    e15n_ing = -1.   ! N15 microzoo ingestion fractionation
+    e13c_calz = 2.    ! C13 microzoo calcification fractionation
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisfer     !   parameters for iron chemistry
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/  
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisrem     !   parameters for remineralization
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+   e15n_denit  =  30.   ! N15 denitrification fractionation
+   e15n_nit    =  0.    ! N15 nitrification fractionation
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampiscal     !   parameters for Calcite chemistry
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampissbc     !   parameters for inputs deposition
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+    ln_ironice   =  .false.   ! boolean for Fe input from sea ice
+    ln_ndepo     =  .false.   ! boolean for atmospheric deposition of N
+    e15n_bdenit  =  3   ! N15 benthic denitrification fractionation
+    d15n_nfix    =  -1.0    ! dN15 Nfix
+    d15n_river   =  2    ! dN15 river input
+    d15n_depo    =  -4    ! dN15 atm. epo
+    d13c_rivdic  =  0.    ! dC13 river dic
+    d13c_rivdoc  =  0.    ! dC13 river doc
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampisdmp     !  Damping 
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+   ln_pisdmp    =  .false.     !  Relaxation fo some tracers to a mean value
+   nn_pisdmp    =  1460       !  Frequency of Relaxation 
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+&nampismass     !  Mass conservation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+/
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_ref
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_ref	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_pisces_ref	(revision 6534)
@@ -0,0 +1,1 @@
+link ../../SHARED/namelist_pisces_ref
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_ref
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_ref	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_ref	(revision 6534)
@@ -0,0 +1,1 @@
+link ../../SHARED/namelist_ref
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_cfg
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_cfg	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_cfg	(revision 6534)
@@ -0,0 +1,108 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/TOP1 : ORCA2_OFF_PISCES configuration namelist used to overwrite SHARED/namelist_top
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!-----------------------------------------------------------------------
+&namtrc_run     !   run information
+!-----------------------------------------------------------------------
+   nn_writetrc   =  1460     !  time step frequency for sn_tracer outputs
+/
+!-----------------------------------------------------------------------
+&namtrc     !   tracers definition
+!-----------------------------------------------------------------------
+!
+!              !    name   !           title of the field              !   units    ! initial data ! save   !
+!              !           !                                           !            ! from file    ! or not ! 
+!              !           !                                           !            ! or not       !        !
+   sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true.
+   sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true.     ,  .true.
+   sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(25)  = 'POC_15  ' , 'N15 Small organic carbon Concentration ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(26)  = 'PHY_15  ' , 'N15 Nanophytoplankton Concentration    ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(27)  = 'ZOO_15  ' , 'N15 Microzooplankton Concentration     ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(28)  = 'DOC_15  ' , 'N15 Dissolved organic Concentration    ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(29)  = 'PHY2_15 ' , 'N15 Diatoms Concentration              ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(30)  = 'ZOO2_15 ' , 'N15 Mesozooplankton Concentration      ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(31)  = 'GOC_15  ' , 'N15 Big organic carbon Concentration   ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(32)  = 'NO3_15  ' , 'N15 Nitrates Concentration             ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(33)  = 'NH4_15  ' , 'N15 Ammonium Concentration             ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(34)  = 'DIC_13  ' , 'C13 Dissolved inorganic carbon Conc    ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(35)  = 'CaCO3_13' , 'C13 Calcite Concentration              ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(36)  = 'POC_13  ' , 'C13 Small organic carbon Concentration ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(37)  = 'PHY_13  ' , 'C13 Nanophytoplankton Concentration    ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(38)  = 'ZOO_13  ' , 'C13 Microzooplankton Concentration     ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(39)  = 'DOC_13  ' , 'C13 Dissolved organic Concentration    ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(40)  = 'PHY2_13 ' , 'C13 Diatoms Concentration              ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(41)  = 'ZOO2_13 ' , 'C13 Mesozooplankton Concentration      ',  'mol-C/L' ,  .false.    ,  .true.
+   sn_tracer(42)  = 'GOC_13  ' , 'C13 Big organic carbon Concentration   ',  'mol-C/L' ,  .false.    ,  .true.
+/
+!-----------------------------------------------------------------------
+&namtrc_dta      !    Initialisation from data input file
+!-----------------------------------------------------------------------
+!
+!          !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
+!          !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
+   sn_trcdta(1)  = 'data_DIC_nomask'        ,        -12        ,  'DIC'     ,    .false.   , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(2)  = 'data_Alkalini_nomask'   ,        -12        ,  'Alkalini',    .false.   , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(3)  = 'data_O2_nomask'      ,        -1         ,  'O2'      ,    .true.    , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(5)  = 'data_PO4_nomask'     ,        -1         ,  'PO4'     ,    .true.    , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(7)  = 'data_Si_nomask'      ,        -1         ,  'Si'      ,    .true.    , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(14) = 'data_Fer_nomask'        ,        -12        ,  'Fer'     ,    .false.   , .true. , 'yearly'  , ''       , ''  , ''
+   sn_trcdta(23) = 'data_NO3_nomask'     ,        -1         ,  'NO3'     ,    .true.    , .true. , 'yearly'  , ''       , ''  , ''
+!
+   cn_dir        =  './'      !  root directory for the location of the data files
+   rn_trfac(1)   =   1.0e-06  !  multiplicative factor
+   rn_trfac(2)   =   1.0e-06  !  -      -      -     -
+   rn_trfac(3)   =  44.6e-06  !  -      -      -     -
+   rn_trfac(5)   = 122.0e-06  !  -      -      -     -
+   rn_trfac(7)   =   1.0e-06  !  -      -      -     -
+   rn_trfac(14)  =   1.0      !  -      -      -     -
+   rn_trfac(23)  =   7.6e-06  !  -      -      -     -
+/
+!-----------------------------------------------------------------------
+&namtrc_adv    !   advection scheme for passive tracer 
+!-----------------------------------------------------------------------
+   ln_trcadv_tvd     =  .false.  !  TVD scheme
+   ln_trcadv_muscl   =  .true.   !  MUSCL scheme
+/
+!-----------------------------------------------------------------------
+&namtrc_ldf    !   lateral diffusion scheme for passive tracer 
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtrc_zdf        !   vertical physics
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtrc_rad        !  treatment of negative concentrations 
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtrc_dia       !   parameters for passive tracer additional diagnostics
+!----------------------------------------------------------------------
+   nn_writedia   =  1460     !  time step frequency for diagnostics
+/
+!----------------------------------------------------------------------
+&namtrc_bc        !   data for boundary conditions
+!-----------------------------------------------------------------------
+/
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_ref
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_ref	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/EXP00/namelist_top_ref	(revision 6534)
@@ -0,0 +1,1 @@
+link ../../SHARED/namelist_top_ref
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/cpp_ORCA2_OFF_PISCES_ISO.fcm
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/cpp_ORCA2_OFF_PISCES_ISO.fcm	(revision 6534)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/ORCA2_OFF_PISCES_ISO/cpp_ORCA2_OFF_PISCES_ISO.fcm	(revision 6534)
@@ -0,0 +1,1 @@
+bld::tool::fppkeys key_trabbl key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_n15 key_c13 key_iomput key_mpp_mpi
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/cfg.txt
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/cfg.txt	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/CONFIG/cfg.txt	(revision 6534)
@@ -7,4 +7,5 @@
 AMM12 OPA_SRC
 ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC
+ORCA2_OFF_PISCES_ISO OPA_SRC OFF_SRC TOP_SRC
 ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC
 ORCA2_LIM3 OPA_SRC LIM_SRC_3 NST_SRC
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90	(revision 6534)
@@ -39,4 +39,7 @@
    !                               !!** Namelist  nampisext  **
    REAL(wp)          ::  atcco2     !: pre-industrial atmospheric [co2] (ppm) 	
+#ifdef key_c13
+   REAL(wp)          ::  d13c_co2   !: d13c of atmospheric [co2] (ppm)   
+#endif
    LOGICAL           ::  ln_co2int  !: flag to read in a file and interpolate atmospheric pco2 or not
    CHARACTER(len=34) ::  clname     !: filename of pco2 values
@@ -55,4 +58,7 @@
 
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: oce_co2   !: ocean carbon flux 
+#if defined key_c13
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: oce_c13  !: ocean C13 flux
+#endif
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: satmco2   !: atmospheric pco2 
 
@@ -84,15 +90,23 @@
       !
       INTEGER  ::   ji, jj, jm, iind, iindm1
-      REAL(wp) ::   ztc, ztc2, ztc3, ztc4, zws, zkgwan
+      REAL(wp) ::   ztc, ztc2, ztc3, zws, zkgwan
       REAL(wp) ::   zfld, zflu, zfld16, zflu16, zfact
       REAL(wp) ::   zph, zah2, zbot, zdic, zalk, zsch_o2, zalka, zsch_co2
       REAL(wp) ::   zyr_dec, zdco2dt
       CHARACTER (len=25) :: charout
-      REAL(wp), POINTER, DIMENSION(:,:) :: zkgco2, zkgo2, zh2co3, zoflx, zw2d 
+      REAL(wp), POINTER, DIMENSION(:,:) :: zkgco2, zkgo2, zh2co3, z_co3, zoflx, zw2d 
+#ifdef key_c13
+      REAL(wp) ::   zft, zfco3
+      REAL(wp), POINTER, DIMENSION(:,:) :: za_dic, za_g
+#endif
       !!---------------------------------------------------------------------
       !
       IF( nn_timing == 1 )  CALL timing_start('p4z_flx')
       !
-      CALL wrk_alloc( jpi, jpj, zkgco2, zkgo2, zh2co3, zoflx )
+#ifdef key_c13 
+      CALL wrk_alloc( jpi, jpj, zkgco2, zkgo2, zh2co3, z_co3, zoflx,  za_dic, za_g )
+#else
+      CALL wrk_alloc( jpi, jpj, zkgco2, zkgo2, zh2co3, z_co3, zoflx )
+#endif
       !
 
@@ -133,4 +147,7 @@
                zph   = MAX( hi(ji,jj,1), 1.e-10 ) / zfact
                zalka = trb(ji,jj,1,jptal) / zfact
+
+               ! CALCULATE [ALK]([CO3--], [HCO3-])
+               zalk  = zalka - (  akw3(ji,jj,1) / zph - zph + zbot / ( 1.+ zph / akb3(ji,jj,1) )  )
 
                ! CALCULATE [ALK]([CO3--], [HCO3-])
@@ -178,4 +195,13 @@
             zkgco2(ji,jj) = zkgwan * SQRT( 660./ zsch_co2 )
             zkgo2 (ji,jj) = zkgwan * SQRT( 660./ zsch_o2 )
+#if defined key_c13
+            ! Compute fractionation factors for C13 from Zhang et al. 1995
+            zfco3 = MAX(0.05,(z_co3(ji,jj)/trb(ji,jj,1,jpdic)+rtrn))
+            zfco3 = MIN(0.2 , zfco3)
+            zft = MIN( 25., ztc )
+            zft = MAX(  5., zft )
+            za_g  (ji,jj) = 1. + ( -0.0049 * zft - 1.31 ) / 1000.
+            za_dic(ji,jj) = 1. + ( 0.014 * zft * zfco3  - 0.105 * zft + 10.53 ) / 1000.
+#endif
          END DO
       END DO
@@ -183,5 +209,5 @@
       DO jj = 1, jpj
          DO ji = 1, jpi
-            ! Compute CO2 flux for the sea and air
+             ! Compute CO2 flux for the sea and air
             zfld = satmco2(ji,jj) * patm(ji,jj) * tmask(ji,jj,1) * chemc(ji,jj) * zkgco2(ji,jj)   ! (mol/L) * (m/s)
             zflu = zh2co3(ji,jj) * tmask(ji,jj,1) * zkgco2(ji,jj)                                   ! (mol/L) (m/s) ?
@@ -195,4 +221,15 @@
             zoflx(ji,jj) = zfld16 - zflu16
             tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + zoflx(ji,jj) * rfact2 / fse3t(ji,jj,1)
+#if defined key_c13
+            tra(ji,jj,1,jp13dic) = tra(ji,jj,1,jp13dic) + ( zfld * (1 + d13c_co2 / 1000.) &
+                                         - zflu * trb(ji,jj,1,jp13dic) / ( trb(ji,jj,1,jpdic) + rtrn ) &
+                                       / ( za_dic(ji,jj) + rtrn) ) &
+                                         * 0.99919 * za_g(ji,jj) * rfact2 / fse3t(ji,jj,1)
+
+             oce_c13(ji,jj) = ( zfld * (1 + d13c_co2 / 1000.) &
+                              - zflu * trb(ji,jj,1,jp13dic) / ( trb(ji,jj,1,jpdic) + rtrn ) &
+                                 / ( za_dic(ji,jj) + rtrn) ) &
+                                 * 0.99919 * za_g(ji,jj)  * rfact2 * e1e2t(ji,jj) * tmask(ji,jj,1) * 1000.
+#endif
          END DO
       END DO
@@ -215,4 +252,11 @@
             CALL iom_put( "Cflx"     , zw2d ) 
          ENDIF
+#if defined key_c13
+         CALL wrk_alloc( jpi, jpj, zw2d )
+         IF( iom_use( "C13flx"  ) )  THEN
+            zw2d(:,:) = oce_c13(:,:) / e1e2t(:,:) * rfact2r
+            CALL iom_put( "C13flx"     , zw2d )
+         ENDIF
+#endif
          IF( iom_use( "Oflx"  ) )  THEN
             zw2d(:,:) =  zoflx(:,:) * 1000 * tmask(:,:,1)
@@ -224,9 +268,9 @@
          ENDIF
          IF( iom_use( "Dpco2" ) ) THEN
-           zw2d(:,:) = ( satmco2(:,:) * patm(:,:) - zh2co3(:,:) / ( chemc(:,:) + rtrn ) ) * tmask(:,:,1)
+           zw2d(:,:) = ( satmco2(:,:) * patm(:,:) - zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1)
            CALL iom_put( "Dpco2" ,  zw2d )
          ENDIF
          IF( iom_use( "Dpo2" ) )  THEN
-           zw2d(:,:) = ( atcox * patm(:,:) - atcox * trn(:,:,1,jpoxy) / ( chemo2(:,:,1) + rtrn ) ) * tmask(:,:,1)
+           zw2d(:,:) = ( atcox * patm(:,:) - trb(:,:,1,jpoxy) / ( chemc(:,:,2) + rtrn ) ) * tmask(:,:,1)
            CALL iom_put( "Dpo2"  , zw2d )
          ENDIF
@@ -240,9 +284,13 @@
             trc2d(:,:,jp_pcs0_2d + 1) = zoflx(:,:) * 1000 * tmask(:,:,1) 
             trc2d(:,:,jp_pcs0_2d + 2) = zkgco2(:,:) * tmask(:,:,1) 
-            trc2d(:,:,jp_pcs0_2d + 3) = ( satmco2(:,:) * patm(:,:) - zh2co3(:,:) / ( chemc(:,:) + rtrn ) ) * tmask(:,:,1) 
-         ENDIF
-      ENDIF
-      !
-      CALL wrk_dealloc( jpi, jpj, zkgco2, zkgo2, zh2co3, zoflx )
+            trc2d(:,:,jp_pcs0_2d + 3) = ( satmco2(:,:) * patm(:,:) - zh2co3(:,:) / ( chemc(:,:,1) + rtrn ) ) * tmask(:,:,1) 
+         ENDIF
+      ENDIF
+      !
+#ifdef key_c13 
+      CALL wrk_dealloc( jpi, jpj, zkgco2, zkgo2, zh2co3, z_co3, zoflx,  za_dic, za_g )
+#else
+      CALL wrk_dealloc( jpi, jpj, zkgco2, zkgo2, zh2co3, z_co3, zoflx )
+#endif
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_flx')
@@ -261,5 +309,9 @@
       !! ** input   :   Namelist nampisext
       !!----------------------------------------------------------------------
-      NAMELIST/nampisext/ln_co2int, atcco2, clname, nn_offset
+      NAMELIST/nampisext/ln_co2int, atcco2, &
+#if defined key_c13
+                         &  d13c_co2,       &
+#endif
+                         &  clname, nn_offset
       INTEGER :: jm
       INTEGER :: ios                 ! Local integer output status for namelist read
@@ -286,4 +338,7 @@
          IF(lwp) THEN                         ! control print
             WRITE(numout,*) '    Constant Atmospheric pCO2 value  atcco2    =', atcco2
+#if defined key_c13
+            WRITE(numout,*) '    Constant Atm. d13C of pCO2 value  d13c_co2 =', d13c_co2
+#endif
             WRITE(numout,*) ' '
          ENDIF
@@ -316,4 +371,7 @@
       t_oce_co2_flx = 0._wp
       t_atm_co2_flx = 0._wp
+#ifdef key_c13
+       oce_c13(:,:)  = 0._wp
+#endif
       !
       CALL p4z_patm( nit000 )
@@ -384,5 +442,5 @@
       !!                     ***  ROUTINE p4z_flx_alloc  ***
       !!----------------------------------------------------------------------
-      ALLOCATE( oce_co2(jpi,jpj), satmco2(jpi,jpj), patm(jpi,jpj), STAT=p4z_flx_alloc )
+      ALLOCATE( oce_co2(jpi,jpj), satmco2(jpi,jpj), patm(jpi,jpj), oce_c13(jpi,jpj), STAT=p4z_flx_alloc )
       !
       IF( p4z_flx_alloc /= 0 )   CALL ctl_warn('p4z_flx_alloc : failed to allocate arrays')
@@ -402,3 +460,3 @@
 
    !!======================================================================
-END MODULE p4zflx
+END MODULE  p4zflx
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90	(revision 6534)
@@ -115,10 +115,10 @@
                ! Michaelis-Menten Limitation term for nutrients Small bacteria
                ! -------------------------------------------------------------
-               zdenom = 1. /  ( concbno3 * concbnh4 + concbnh4 * trb(ji,jj,jk,jpno3) + concbno3 * trb(ji,jj,jk,jpnh4) )
-               xnanono3(ji,jj,jk) = trb(ji,jj,jk,jpno3) * concbnh4 * zdenom
-               xnanonh4(ji,jj,jk) = trb(ji,jj,jk,jpnh4) * concbno3 * zdenom
+               xnanonh4(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-8 ) / (MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-8 ) + concbnh4)
+               xnanono3(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpno3) - 1.E-8 ) / (MAX(0., trb(ji,jj,jk,jpno3) - 1.E-8 )  &
+                          &                                                  + concbno3) * (1. - xnanonh4(ji,jj,jk) ) 
                !
                zlim1    = xnanono3(ji,jj,jk) + xnanonh4(ji,jj,jk)
-               zlim2    = trb(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + concbnh4 )
+               zlim2    = MAX(0., trb(ji,jj,jk,jppo4) - 1E-9 )/ ( trb(ji,jj,jk,jppo4) + concbnh4 )
                zlim3    = trb(ji,jj,jk,jpfer) / ( concbfe + trb(ji,jj,jk,jpfer) )
                zlim4    = trb(ji,jj,jk,jpdoc) / ( xkdoc   + trb(ji,jj,jk,jpdoc) )
@@ -128,10 +128,10 @@
                ! Michaelis-Menten Limitation term for nutrients Small flagellates
                ! -----------------------------------------------
-               zdenom = 1. /  ( zconc0n * zconc0nnh4 + zconc0nnh4 * trb(ji,jj,jk,jpno3) + zconc0n * trb(ji,jj,jk,jpnh4) )
-               xnanono3(ji,jj,jk) = trb(ji,jj,jk,jpno3) * zconc0nnh4 * zdenom
-               xnanonh4(ji,jj,jk) = trb(ji,jj,jk,jpnh4) * zconc0n    * zdenom
+               xnanonh4(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-7 ) / (MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-7 ) + zconc0nnh4)
+               xnanono3(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpno3) - 1.E-7 ) / (MAX(0., trb(ji,jj,jk,jpno3) - 1.E-7 )  &
+                          &                      + zconc0n) * (1. - xnanonh4(ji,jj,jk) )   
                !
                zlim1    = xnanono3(ji,jj,jk) + xnanonh4(ji,jj,jk)
-               zlim2    = trb(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + zconc0nnh4 )
+               zlim2    = MAX(0., trb(ji,jj,jk,jppo4) - 1.E-7 ) / ( trb(ji,jj,jk,jppo4) + zconc0nnh4 )
                zratio   = trb(ji,jj,jk,jpnfe) * z1_trbphy 
                zironmin = xcoef1 * trb(ji,jj,jk,jpnch) * z1_trbphy + xcoef2 * zlim1 + xcoef3 * xnanono3(ji,jj,jk)
@@ -143,10 +143,10 @@
                !   Michaelis-Menten Limitation term for nutrients Diatoms
                !   ----------------------------------------------
-               zdenom   = 1. / ( zconc1d * zconc1dnh4 + zconc1dnh4 * trb(ji,jj,jk,jpno3) + zconc1d * trb(ji,jj,jk,jpnh4) )
-               xdiatno3(ji,jj,jk) = trb(ji,jj,jk,jpno3) * zconc1dnh4 * zdenom
-               xdiatnh4(ji,jj,jk) = trb(ji,jj,jk,jpnh4) * zconc1d    * zdenom
+               xdiatnh4(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-7 ) / (MAX(0., trb(ji,jj,jk,jpnh4) - 1.E-7 ) + zconc1dnh4)
+               xdiatno3(ji,jj,jk) = MAX(0., trb(ji,jj,jk,jpno3) - 1.E-7 ) / (MAX(0., trb(ji,jj,jk,jpno3) - 1.E-7 ) &
+                          &                   + zconc1d ) * (1- xdiatnh4(ji,jj,jk)) 
                !
                zlim1    = xdiatno3(ji,jj,jk) + xdiatnh4(ji,jj,jk)
-               zlim2    = trb(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + zconc1dnh4  )
+               zlim2    = MAX(0., trb(ji,jj,jk,jppo4) - 1.E-7 ) / ( trb(ji,jj,jk,jppo4) + zconc1dnh4  )
                zlim3    = trb(ji,jj,jk,jpsil) / ( trb(ji,jj,jk,jpsil) + xksi(ji,jj) )
                zratio   = trb(ji,jj,jk,jpdfe) * z1_trbdia
@@ -265,3 +265,3 @@
 
    !!======================================================================
-END MODULE p4zlim
+END MODULE  p4zlim
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90	(revision 6534)
@@ -142,4 +142,11 @@
               tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) -      zcaldiss(ji,jj,jk)
               tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) +      zcaldiss(ji,jj,jk)
+
+#if defined key_c13
+              zaldi  = zcaldiss(ji,jj,jk) *  trb(ji,jj,jk,jp13cal) / ( trb(ji,jj,jk,jpcal) + rtrn )
+              tra(ji,jj,jk,jp13cal) = tra(ji,jj,jk,jp13cal) - zaldi
+              tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) + zaldi 
+#endif
+
             END DO
          END DO
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90	(revision 6534)
@@ -49,4 +49,11 @@
    REAL(wp), PUBLIC ::  epsher2      !: half sturation constant for grazing 2
    REAL(wp), PUBLIC ::  grazflux     !: mesozoo flux feeding rate
+#if defined key_n15
+   REAL(wp), PUBLIC ::  e15n_ex2     !: N15 mesozoo excretion fractionation
+   REAL(wp), PUBLIC ::  e15n_ing2    !: N15 mesozoo ingestion fractionation
+#endif
+#if defined key_c13
+   REAL(wp), PUBLIC ::  e13c_cal2    !: C13 mesozoo calcification fractionation
+#endif
 
    !!* Substitution
@@ -82,4 +89,12 @@
       REAL(wp) :: zgrazn, zgrazpoc, zgraznf, zgrazf
       REAL(wp) :: zgrazfffp, zgrazfffg, zgrazffep, zgrazffeg
+#if defined key_n15
+      REAL(wp) :: zgrazn15,zgrazd15,zgrazpo15,zgrazff15p,zgrazff15g,zgraztot15
+      REAL(wp) :: zgrarem2_15,zgrarsig_15,zgrarsigex_15,zmortzgoc_15,zgrapoc2_15,zfrac15
+#endif
+#if defined key_c13
+      REAL(wp) :: zgrazn13,zgrazd13,zgrazpo13,zgrazff13p,zgrazff13g,zgraztot13
+      REAL(wp) :: zgrarem2_13,zgrarsig_13,zmortzgoc_13,zgrapoc2_13,zfrac13
+#endif
       CHARACTER (len=25) :: charout
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d
@@ -139,4 +154,16 @@
                zgrazpof  = zgrazpoc * trb(ji,jj,jk,jpsfe) / ( trb(ji,jj,jk,jppoc) + rtrn)
 
+#ifdef key_n15
+               zgrazn15   = zgrazn   * trb(ji,jj,jk,jp15phy) / ( trb(ji,jj,jk,jpphy) + rtrn)
+               zgrazd15   = zgrazd   * trb(ji,jj,jk,jp15dia) / ( trb(ji,jj,jk,jpdia) + rtrn)
+               zgrazpo15  = zgrazpoc * trb(ji,jj,jk,jp15poc) / ( trb(ji,jj,jk,jppoc) + rtrn)
+#endif
+#ifdef key_c13
+               zgrazn13   = zgrazn   * trb(ji,jj,jk,jp13phy) / ( trb(ji,jj,jk,jpphy) + rtrn)
+               zgrazd13   = zgrazd   * trb(ji,jj,jk,jp13dia) / ( trb(ji,jj,jk,jpdia) + rtrn)
+               zgrazpo13  = zgrazpoc * trb(ji,jj,jk,jp13poc) / ( trb(ji,jj,jk,jppoc) + rtrn)
+#endif
+
+
                !  Mesozooplankton flux feeding on GOC
                !  ----------------------------------
@@ -150,4 +177,13 @@
                &           * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jppoc) * trb(ji,jj,jk,jpmes)
                zgrazfffp = zgrazffep * trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
+#ifdef key_n15
+               zgrazff15g = zgrazffeg * trb(ji,jj,jk,jp15goc) / (trb(ji,jj,jk,jpgoc) + rtrn)
+               zgrazff15p = zgrazffep * trb(ji,jj,jk,jp15poc) / (trb(ji,jj,jk,jppoc) + rtrn)
+#endif
+#ifdef key_c13
+               zgrazff13g = zgrazffeg * trb(ji,jj,jk,jp13goc) / (trb(ji,jj,jk,jpgoc) + rtrn)
+               zgrazff13p = zgrazffep * trb(ji,jj,jk,jp13poc) / (trb(ji,jj,jk,jppoc) + rtrn)
+#endif
+
               !
 # if ! defined key_kriest
@@ -173,4 +209,20 @@
               &   + zgrazpoc + zgrazffep + zgrazffeg
               zgraztotf = zgrazf + zgraznf + zgrazz * ferat3 + zgrazpof + zgrazfffp + zgrazfffg
+
+#ifdef key_n15
+              zfrac15   =  zfrac * trb(ji,jj,jk,jp15goc) / (trb(ji,jj,jk,jpgoc) + rtrn)
+              zgrazff15p = zproport * zgrazff15p
+              zgrazff15g = zproport * zgrazff15g
+              zgraztot15 = zgrazd15 + zgrazn15 + zgrazz * trb(ji,jj,jk,jp15zoo) / (trb(ji,jj,jk,jpzoo) + rtrn)  + zgrazpo15 + zgrazff15p + zgrazff15g
+#endif
+#ifdef key_c13
+              zfrac13   =  zfrac * trb(ji,jj,jk,jp13goc) / (trb(ji,jj,jk,jpgoc) + rtrn)
+              zgrazff13p = zproport * zgrazff13p
+              zgrazff13g = zproport * zgrazff13g
+              zgraztot13 = zgrazd13 + zgrazn13 + zgrazz * trb(ji,jj,jk,jp13zoo) / (trb(ji,jj,jk,jpzoo) + rtrn)  + zgrazpo13 + zgrazff13p + zgrazff13g
+#endif
+
+
+
 # else
               zgraztot  = zgrazd + zgrazz + zgrazn + zgrazpoc + zgrazffep
@@ -198,4 +250,25 @@
                 &       + ferat3 * ( ( 1. - epsher2 - unass2 ) /( 1. - epsher2 ) * ztortz2 )
                zgrapoc2  = zgraztot * unass2
+
+#ifdef key_n15
+               zgrarem2_15  = zgraztot15 * ( 1. - zepsherv - unass2 )            &
+                        + ( 1. - epsher2 - unass2 ) / ( 1. - epsher2 ) * ztortz2 &
+                        * trb(ji,jj,jk,jp15mes) / ( trb(ji,jj,jk,jpmes) + rtrn )
+               zgrarsig_15  = zgrarem2_15 * sigma2
+               zgrarsigex_15  =  zgraztot15 * ( 1. - epsher2 - unass2 ) * sigma2 * zgrasratn
+               zgrapoc2_15  = zgraztot15 * unass2
+               zmortzgoc_15  = ( unass2 / ( 1. - epsher2 ) * ztortz2 + zrespz2 )  &
+                        * trb(ji,jj,jk,jp15mes) / ( trb(ji,jj,jk,jpmes) + rtrn )
+#endif
+#ifdef key_c13
+               zgrarem2_13  = zgraztot13 * ( 1. - zepsherv - unass2 ) &
+                        + ( 1. - epsher2 - unass2 ) / ( 1. - epsher2 ) * ztortz2 &
+                        * trb(ji,jj,jk,jp13mes) / ( trb(ji,jj,jk,jpmes) + rtrn )
+               zgrarsig_13  = zgrarem2_13 * sigma2
+               zgrapoc2_13  = zgraztot13 * unass2
+               zmortzgoc_13  = ( unass2 / ( 1. - epsher2 ) * ztortz2 + zrespz2 ) &
+                        * trb(ji,jj,jk,jp13mes) / ( trb(ji,jj,jk,jpmes) + rtrn )
+#endif
+
 
                !   Update the arrays TRA which contain the biological sources and sinks
@@ -222,4 +295,26 @@
                tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) - zgrazf
 
+#if defined key_n15
+               tra(ji,jj,jk,jp15nh4) = tra(ji,jj,jk,jp15nh4) + (1. - e15n_ex2 / 1000.) * zgrarsigex_15   &
+               &  + (zgrarsig_15-zgrarsigex_15)
+               tra(ji,jj,jk,jp15doc) = tra(ji,jj,jk,jp15doc) + zgrarem2_15 - zgrarsig_15
+               tra(ji,jj,jk,jp15mes) = tra(ji,jj,jk,jp15mes) - zmortz2 * trb(ji,jj,jk,jp15mes) / (trb(ji,jj,jk,jpmes) + rtrn) &
+                                                             + zepsherv * zgraztot15 &
+                                                             + zgrarsigex_15 * (e15n_ex2 / 1000.) &
+               &                                             + zgrapoc2_15  * (e15n_ing2 / 1000.)
+               tra(ji,jj,jk,jp15dia) = tra(ji,jj,jk,jp15dia) - zgrazd15 
+               tra(ji,jj,jk,jp15zoo) = tra(ji,jj,jk,jp15zoo) - zgrazz * trb(ji,jj,jk,jp15zoo) / ( trb(ji,jj,jk,jpzoo) + rtrn )
+               tra(ji,jj,jk,jp15phy) = tra(ji,jj,jk,jp15phy) - zgrazn15 
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) + zgrarsig_13
+               tra(ji,jj,jk,jp13doc) = tra(ji,jj,jk,jp13doc) + zgrarem2_13 - zgrarsig_13
+               tra(ji,jj,jk,jp13mes) = tra(ji,jj,jk,jp13mes) - zmortz2 * trb(ji,jj,jk,jp13mes) / (trb(ji,jj,jk,jpmes) + rtrn) &
+                                                             + zepsherv * zgraztot13 
+               tra(ji,jj,jk,jp13dia) = tra(ji,jj,jk,jp13dia) - zgrazd13
+               tra(ji,jj,jk,jp13zoo) = tra(ji,jj,jk,jp13zoo) - zgrazz * trb(ji,jj,jk,jp13zoo) / ( trb(ji,jj,jk,jpzoo) + rtrn )
+               tra(ji,jj,jk,jp13phy) = tra(ji,jj,jk,jp13phy) - zgrazn13
+#endif
+
                ! calcite production
                zprcaca = xfracal(ji,jj,jk) * zgrazn
@@ -230,4 +325,11 @@
                tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca
                tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca
+#if defined key_c13
+               tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) - zprcaca *  ( 1.  - e13c_cal2/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) )
+               tra(ji,jj,jk,jp13cal) = tra(ji,jj,jk,jp13cal) + zprcaca *  ( 1.  - e13c_cal2/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) )
+#endif
+
 #if defined key_kriest
               znumpoc = trb(ji,jj,jk,jpnum) / ( trb(ji,jj,jk,jppoc) + rtrn )
@@ -243,4 +345,13 @@
               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ferat3 * zmortzgoc - zgrazfffg     &
                  &                + zgraztotf * unass2 - zfracfe
+#if defined key_n15
+              tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) - zgrazpo15 - zgrazff15p + zfrac15
+              tra(ji,jj,jk,jp15goc) = tra(ji,jj,jk,jp15goc) + zmortzgoc_15 - zgrazff15g       &
+                 &    + (1. - e15n_ing2/1000.) * zgrapoc2_15 - zfrac15
+#endif
+#if defined key_c13
+              tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) - zgrazpo13 - zgrazff13p + zfrac13
+              tra(ji,jj,jk,jp13goc) = tra(ji,jj,jk,jp13goc) + zmortzgoc_13 - zgrazff13g + zgrapoc2_13 - zfrac13
+#endif
 #endif
             END DO
@@ -289,5 +400,12 @@
       NAMELIST/nampismes/ part2, grazrat2, resrat2, mzrat2, xprefc, xprefp, xprefz,   &
          &                xprefpoc, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, &
+#if defined key_n15
+         &                e15n_ex2, e15n_ing2,   &
+#endif
+#if defined key_c13
+         &                e13c_cal2,  &
+#endif
          &                xthresh2, xkgraz2, epsher2, sigma2, unass2, grazflux
+
       INTEGER :: ios                 ! Local integer output status for namelist read
 
@@ -324,4 +442,11 @@
          WRITE(numout,*) '    Fraction of mesozoo excretion as DOM           sigma2       =', sigma2
          WRITE(numout,*) '    half sturation constant for grazing 2          xkgraz2      =', xkgraz2
+#if defined key_n15
+         WRITE(numout,*) '    N15 mesozoo excretion fractionation            e15n_ex2     =', e15n_ex2
+         WRITE(numout,*) '    N15 mesozoo ingestion fractionation            e15n_ing2    =', e15n_ing2
+#endif
+#if defined key_c13
+         WRITE(numout,*) '    C13 mesozoo calcification fractionation        e13c_cal2    =', e13c_cal2
+#endif
       ENDIF
 
@@ -340,3 +465,3 @@
 
    !!======================================================================
-END MODULE p4zmeso
+END MODULE  p4zmeso
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90	(revision 6534)
@@ -47,5 +47,11 @@
    REAL(wp), PUBLIC ::  sigma1      !: Fraction of microzoo excretion as DOM 
    REAL(wp), PUBLIC ::  epsher      !: half sturation constant for grazing 1 
-
+#if defined key_n15
+   REAL(wp), PUBLIC ::  e15n_ex     !: N15 microzoo excretion fractionation
+   REAL(wp), PUBLIC ::  e15n_ing     !: N15 microzoo ingestion fractionation
+#endif
+#if defined key_c13
+   REAL(wp), PUBLIC ::  e13c_calz   !: C13 microzoo calcification fractionation
+#endif
 
    !!* Substitution
@@ -79,4 +85,12 @@
       REAL(wp) :: zgrazp, zgrazm, zgrazsd
       REAL(wp) :: zgrazmf, zgrazsf, zgrazpf
+#if defined key_n15
+      REAL(wp) :: zgrazp15,zgrazm15,zgrazs15,zgraztot15
+      REAL(wp) :: zgrarem_15,zgrarsig_15,zgrarsigex_15,zmortz_15,zgrapoc_15
+#endif
+#if defined key_c13
+      REAL(wp) :: zgrazp13,zgrazm13,zgrazs13,zgraztot13
+      REAL(wp) :: zgrarem_13,zgrarsig_13,zmortz_13,zgrapoc_13
+#endif
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d
       CHARACTER (len=25) :: charout
@@ -126,8 +140,22 @@
                zgrazmf   = zgrazm  * trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
                zgrazsf   = zgrazsd * trb(ji,jj,jk,jpdfe) / (trb(ji,jj,jk,jpdia) + rtrn)
+
                !
                zgraztot  = zgrazp  + zgrazm  + zgrazsd 
                zgraztotf = zgrazpf + zgrazsf + zgrazmf 
                zgraztotn = zgrazp * quotan(ji,jj,jk) + zgrazm + zgrazsd * quotad(ji,jj,jk)
+               !
+#ifdef key_n15
+               zgrazp15   = zgrazp   * trb(ji,jj,jk,jp15phy) / ( trb(ji,jj,jk,jpphy) + rtrn)
+               zgrazm15   = zgrazm   * trb(ji,jj,jk,jp15poc) / ( trb(ji,jj,jk,jppoc) + rtrn)
+               zgrazs15  =  zgrazsd  * trb(ji,jj,jk,jp15dia) / ( trb(ji,jj,jk,jpdia) + rtrn)
+               zgraztot15 = zgrazp15 + zgrazm15 + zgrazs15
+#endif
+#ifdef key_c13
+               zgrazp13   = zgrazp   * trb(ji,jj,jk,jp13phy) / ( trb(ji,jj,jk,jpphy) + rtrn)
+               zgrazm13   = zgrazm   * trb(ji,jj,jk,jp13poc) / ( trb(ji,jj,jk,jppoc) + rtrn)
+               zgrazs13  =  zgrazsd  * trb(ji,jj,jk,jp13dia) / ( trb(ji,jj,jk,jpdia) + rtrn)
+               zgraztot13 = zgrazp13 + zgrazm13 + zgrazs13
+#endif
 
                ! Grazing by microzooplankton
@@ -143,4 +171,19 @@
                zgrarem   = zgraztot * ( 1. - zepsherv - unass )
                zgrapoc   = zgraztot * unass
+
+#ifdef key_n15
+               zgrarem_15   = zgraztot15 * ( 1. - zepsherv - unass )
+               zgrapoc_15   = zgraztot15 * unass
+               zgrarsig_15  = zgrarem_15 * sigma1
+               zgrarsigex_15 = (1. - epsher - unass ) * zgraztot15 * sigma1 * zgrasratn
+               zmortz_15 = (ztortz + zrespz) * trb(ji,jj,jk,jp15zoo) / ( trb(ji,jj,jk,jpzoo) + rtrn )
+#endif
+#ifdef key_c13
+               zgrarem_13   = zgraztot13 * ( 1. - zepsherv - unass )
+               zgrapoc_13   = zgraztot13 * unass
+               zgrarsig_13  = zgrarem_13 * sigma1
+               zmortz_13 = (ztortz + zrespz) * trb(ji,jj,jk,jp13zoo) / ( trb(ji,jj,jk,jpzoo) + rtrn )
+#endif
+
 
                !  Update of the TRA arrays
@@ -159,4 +202,18 @@
                tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + zgrapoc * xkr_dmicro
 #endif
+               !
+#if defined key_n15
+               tra(ji,jj,jk,jp15nh4) = tra(ji,jj,jk,jp15nh4) + (1. - e15n_ex / 1000.) * zgrarsigex_15   &
+               &   + (zgrarsig_15-zgrarsigex_15)
+               tra(ji,jj,jk,jp15doc) = tra(ji,jj,jk,jp15doc) + zgrarem_15 - zgrarsig_15
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) + (1. - e15n_ing / 1000.) * zgrapoc_15
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) + zgrarsig_13
+               tra(ji,jj,jk,jp13doc) = tra(ji,jj,jk,jp13doc) + zgrarem_13 - zgrarsig_13
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) + zgrapoc_13
+#endif
+
+               !
                !   Update the arrays TRA which contain the biological sources and sinks
                !   --------------------------------------------------------------------
@@ -186,4 +243,24 @@
                                                          - zgrazm * trb(ji,jj,jk,jpnum) / ( trb(ji,jj,jk,jppoc) + rtrn )
 #endif
+               !
+#if defined key_n15
+               tra(ji,jj,jk,jp15zoo) = tra(ji,jj,jk,jp15zoo) - zmortz_15 + zepsherv * zgraztot15    &
+               &                       + e15n_ex / 1000. * zgrarsigex_15 + e15n_ing / 1000. * zgrapoc_15
+               tra(ji,jj,jk,jp15phy) = tra(ji,jj,jk,jp15phy) - zgrazp15
+               tra(ji,jj,jk,jp15dia) = tra(ji,jj,jk,jp15dia) - zgrazs15
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) + zmortz_15 - zgrazm15
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13zoo) = tra(ji,jj,jk,jp13zoo) - zmortz_13 + zepsherv * zgraztot13
+               tra(ji,jj,jk,jp13phy) = tra(ji,jj,jk,jp13phy) - zgrazp13
+               tra(ji,jj,jk,jp13dia) = tra(ji,jj,jk,jp13dia) - zgrazs13
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) + zmortz_13 - zgrazm13
+
+               tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) - zprcaca *  ( 1.  - e13c_calz/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) )
+               tra(ji,jj,jk,jp13cal) = tra(ji,jj,jk,jp13cal) + zprcaca *  ( 1.  - e13c_calz/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) )
+#endif
+
             END DO
          END DO
@@ -228,4 +305,10 @@
       NAMELIST/nampiszoo/ part, grazrat, resrat, mzrat, xpref2c, xpref2p, &
          &                xpref2d,  xthreshdia,  xthreshphy,  xthreshpoc, &
+#ifdef key_n15
+         &                e15n_ex, e15n_ing,   &
+#endif
+#ifdef key_c13
+         &                e13c_calz, &
+#endif
          &                xthresh, xkgraz, epsher, sigma1, unass
       INTEGER :: ios                 ! Local integer output status for namelist read
@@ -259,4 +342,11 @@
          WRITE(numout,*) '    Fraction of microzoo excretion as DOM           sigma1      =', sigma1
          WRITE(numout,*) '    half sturation constant for grazing 1           xkgraz      =', xkgraz
+#if defined key_n15
+         WRITE(numout,*) '    N15 microzoo excretion fractionation            e15n_ex     =', e15n_ex
+         WRITE(numout,*) '    N15 microzoo ingestion fractionation            e15n_ing    =', e15n_ing
+#endif
+#if defined key_c13
+         WRITE(numout,*) '    C13 microzoo calcification fractionation        e13c_calz   =', e13c_calz
+#endif
       ENDIF
 
@@ -273,3 +363,3 @@
 
    !!======================================================================
-END MODULE p4zmicro
+END MODULE  p4zmicro
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90	(revision 6534)
@@ -33,4 +33,7 @@
    REAL(wp), PUBLIC :: mprat   !:
    REAL(wp), PUBLIC :: mprat2  !:
+#if defined key_c13
+   REAL(wp), PUBLIC :: e13c_cal !:
+#endif
 
 
@@ -121,4 +124,20 @@
                tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca
                tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca
+#if defined key_n15
+               tra(ji,jj,jk,jp15phy) = tra(ji,jj,jk,jp15phy) - zmortp   &
+                                                  * ( trb(ji,jj,jk,jp15phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+#endif
+
+#if defined key_c13
+               tra(ji,jj,jk,jp13phy) = tra(ji,jj,jk,jp13phy) - zmortp   &
+                                                  * ( trb(ji,jj,jk,jp13phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+
+               tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) - zprcaca *  ( 1. - e13c_cal/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) ) 
+               tra(ji,jj,jk,jp13cal) = tra(ji,jj,jk,jp13cal) + zprcaca *  ( 1. - e13c_cal/1000. ) &
+                                                  * ( trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) )
+#endif
+
+               !
 #if defined key_kriest
                tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortp
@@ -130,4 +149,18 @@
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ( 1. - zfracal ) * zmortp * zfactfe
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zfracal * zmortp * zfactfe
+#if defined key_n15
+               tra(ji,jj,jk,jp15goc) = tra(ji,jj,jk,jp15goc) + zfracal * zmortp & 
+                                                  * ( trb(ji,jj,jk,jp15phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) + ( 1. - zfracal ) * zmortp &
+                                                  * ( trb(ji,jj,jk,jp15phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+#endif
+
+#if defined key_c13
+               tra(ji,jj,jk,jp13goc) = tra(ji,jj,jk,jp13goc) + zfracal * zmortp &
+                                                  * ( trb(ji,jj,jk,jp13phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) + ( 1. - zfracal ) * zmortp &
+                                                  * ( trb(ji,jj,jk,jp13phy) / (trb(ji,jj,jk,jpphy) + rtrn ) )
+#endif
+
 #endif
             END DO
@@ -204,4 +237,14 @@
                tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zmortp2 * zfactsi
                tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zmortp2 * zfactsi
+#if defined key_n15
+               tra(ji,jj,jk,jp15dia) = tra(ji,jj,jk,jp15dia) - zmortp2 & 
+                                                  * ( trb(ji,jj,jk,jp15dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13dia) = tra(ji,jj,jk,jp13dia) - zmortp2 &
+                                                  * ( trb(ji,jj,jk,jp13dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+#endif
+
+
 #if defined key_kriest
                tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortp2  
@@ -213,4 +256,18 @@
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + 0.5 * ztortp2 * zfactfe
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ( zrespp2 + 0.5 * ztortp2 ) * zfactfe
+#if defined key_n15
+               tra(ji,jj,jk,jp15goc) = tra(ji,jj,jk,jp15goc) + ( zrespp2 + 0.5 * ztortp2 ) &
+                                                  * ( trb(ji,jj,jk,jp15dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) + 0.5 * ztortp2 &
+                                                  * ( trb(ji,jj,jk,jp15dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+#endif
+
+#if defined key_n15
+               tra(ji,jj,jk,jp13goc) = tra(ji,jj,jk,jp13goc) + ( zrespp2 + 0.5 * ztortp2 ) &
+                                                  * ( trb(ji,jj,jk,jp13dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) + 0.5 * ztortp2 &
+                                                  * ( trb(ji,jj,jk,jp13dia) / (trb(ji,jj,jk,jpdia) + rtrn ) )
+#endif
+
 #endif
             END DO
@@ -242,5 +299,9 @@
       !!----------------------------------------------------------------------
 
-      NAMELIST/nampismort/ wchl, wchld, wchldm, mprat, mprat2
+      NAMELIST/nampismort/ wchl, wchld, wchldm, mprat, &
+#ifdef key_c13
+                &     e13c_cal,                        &
+#endif
+                &     mprat2
       INTEGER :: ios                 ! Local integer output status for namelist read
 
@@ -263,4 +324,7 @@
          WRITE(numout,*) '    phytoplankton mortality rate              mprat     =', mprat
          WRITE(numout,*) '    Diatoms mortality rate                    mprat2    =', mprat2
+#ifdef key_c13
+          WRITE(numout,*) '    C13 Calcification Fractionation          e13c_cal  =', e13c_cal
+#endif
       ENDIF
 
@@ -277,3 +341,3 @@
 
    !!======================================================================
-END MODULE p4zmort
+END MODULE  p4zmort
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90	(revision 6534)
@@ -45,4 +45,11 @@
    REAL(wp), PUBLIC ::  fecdm           !:
    REAL(wp), PUBLIC ::  grosip          !:
+#if defined key_n15
+   REAL(wp), PUBLIC ::  e15n_prod       !:
+#endif
+#if defined key_c13
+   REAL(wp), PUBLIC ::  e13c_min       !:
+   REAL(wp), PUBLIC ::  e13c_max      !:
+#endif
 
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   prmax    !: optimal production = f(temperature)
@@ -76,8 +83,15 @@
       INTEGER, INTENT(in) :: kt, knt
       !
-      INTEGER  ::   ji, jj, jk
+      INTEGER  ::   ji, jj, jk, jm
       REAL(wp) ::   zsilfac, znanotot, zdiattot, zconctemp, zconctemp2
       REAL(wp) ::   zratio, zmax, zsilim, ztn, zadap
       REAL(wp) ::   zlim, zsilfac2, zsiborn, zprod, zproreg, zproreg2
+#ifdef key_n15
+      REAL(wp) ::   zr15_reg, zr15_new, zu_15, zun_15
+#endif
+#ifdef key_c13
+      REAL(wp) ::   zr13, zr13_2, z_e13c_prod, z_e13c_prod2
+      REAL(wp) ::   ztc, zft, zrhop, zfco3, zbot, zdic, zph, zalka, zalk, zah2
+#endif
       REAL(wp) ::   zmxltst, zmxlday, zmaxday
       REAL(wp) ::   zpislopen  , zpislope2n
@@ -88,4 +102,7 @@
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zpislopead, zpislopead2, zprdia, zprbio, zprdch, zprnch, zysopt, zw3d   
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zprorca, zprorcad, zprofed, zprofen, zprochln, zprochld, zpronew, zpronewd
+#ifdef key_c13
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: za_g, za_dic, zh2co3, z_co3
+#endif
       !!---------------------------------------------------------------------
       !
@@ -96,4 +113,7 @@
       CALL wrk_alloc( jpi, jpj, jpk, zpislopead, zpislopead2, zprdia, zprbio, zprdch, zprnch, zysopt            ) 
       CALL wrk_alloc( jpi, jpj, jpk, zprorca, zprorcad, zprofed, zprofen, zprochln, zprochld, zpronew, zpronewd )
+#ifdef key_c13
+      CALL wrk_alloc( jpi, jpj, jpk, za_g, za_dic, zh2co3, z_co3)
+#endif
       !
       zprorca (:,:,:) = 0._wp
@@ -294,4 +314,60 @@
          END DO
       END DO
+
+      ! Computation of the C13 fractionation factor from Laws et al. 1995
+#if defined key_c13
+      DO jm = 1, 10
+!CDIR NOVERRCHK
+         DO jk = 1, jpkm1
+!CDIR NOVERRCHK
+           DO jj = 1, jpj
+!CDIR NOVERRCHK
+             DO ji = 1, jpi
+
+               ! DUMMY VARIABLES FOR DIC, H+, AND BORATE
+               zbot  = borat(ji,jj,jk)
+               zrhop = rhop(ji,jj,jk) / 1000. + rtrn
+               zdic  = trb(ji,jj,jk,jpdic) / zrhop
+               zph   = MAX( hi(ji,jj,jk), 1.e-10 ) / zrhop
+               zalka = trb(ji,jj,jk,jptal) / zrhop
+
+               ! CALCULATE [ALK]([CO3--], [HCO3-])
+               zalk  = zalka - (  akw3(ji,jj,jk) / zph - zph + zbot / ( 1.+ zph / akb3(ji,jj,jk) )  )
+
+               ! CALCULATE [H+] AND [H2CO3]
+               zah2   = SQRT(  (zdic-zalk)**2 + 4.* ( zalk * ak23(ji,jj,jk)   &
+                  &                               / ak13(ji,jj,jk) ) * ( 2.* zdic - zalk )  )
+               zah2   = 0.5 * ak13(ji,jj,jk) / zalk * ( ( zdic - zalk ) + zah2 )
+               zh2co3(ji,jj,jk) = ( 2.* zdic - zalk ) / ( 2.+ ak13(ji,jj,jk) / zah2) * zrhop
+               z_co3(ji,jj,jk) = zalk / ( 2. + zah2 / ak23(ji,jj,jk) ) * zrhop
+               hi(ji,jj,jk)   = zah2 * zrhop
+             END DO
+           END DO
+         END DO
+      END DO
+
+!CDIR NOVERRCHK
+         DO jk = 1, jpkm1
+!CDIR NOVERRCHK
+           DO jj = 1, jpj
+!CDIR NOVERRCHK
+             DO ji = 1, jpi
+
+
+!           ! Compute fractionation factors for C13 from Zhang et al. 1995
+
+               zfco3 = MAX(0.05,(z_co3(ji,jj,jk)/trb(ji,jj,jk,jpdic)+rtrn))
+               zfco3 = MIN(0.2 , zfco3)
+               ztc = MIN( 35., tsn(ji,jj,jk,jp_tem) )
+               zft = MIN(25.,ztc)
+               zft = MAX( 5.,zft)
+               za_g(ji,jj,jk) = 1.   + ( -0.0049 * zft - 1.31 ) / 1000.
+               za_dic(ji,jj,jk) = 1. + ( 0.014 * zft * zfco3 &
+                   - 0.105 * zft + 10.53 ) / 1000.
+            END DO
+         END DO
+      END DO
+
+#endif
 
       ! Computation of the various production terms 
@@ -390,4 +466,34 @@
               zproreg  = zprorca(ji,jj,jk) - zpronew(ji,jj,jk)
               zproreg2 = zprorcad(ji,jj,jk) - zpronewd(ji,jj,jk)
+#ifdef key_n15
+              zu_15 = MIN(1.,MAX(0., 1. - (zpronew(ji,jj,jk)+zpronewd(ji,jj,jk))/(rtrn + trb(ji,jj,jk,jp15no3))))
+              zun_15 = MIN(1.,MAX(0., 1. - (zproreg+zproreg2) / (rtrn + trb(ji,jj,jk,jp15nh4))))
+
+              zr15_new = ( 1. - e15n_prod/1000. * zu_15) * trb(ji,jj,jk,jp15no3) / (trb(ji,jj,jk,jpno3) + rtrn)
+              zr15_reg = ( 1. - e15n_prod/1000. * zun_15) * trb(ji,jj,jk,jp15nh4) / (trb(ji,jj,jk,jpnh4) + rtrn)
+#endif
+
+#ifdef key_c13
+
+              z_e13c_prod  = MAX(e13c_min,MIN(e13c_max,( (86400 * zprorca(ji,jj,jk) / (rtrn + trb(ji,jj,jk,jpphy) ) &
+              &                               / ( rfact2 * (zh2co3(ji,jj,jk)/1025*1.e9 + rtrn)  ) - 0.371 ) / (-0.015))))
+              
+              z_e13c_prod2 = MAX(e13c_min,MIN(e13c_max,( ( 86400 * zprorcad(ji,jj,jk) / (rtrn + trb(ji,jj,jk,jpdia) ) &
+              &                               / ( rfact2 * (zh2co3(ji,jj,jk)/1025*1.e9 + rtrn)  ) - 0.371 ) / (-0.015))))
+
+              !
+              if ((mig(ji).eq.74).and.(mig(jj).eq.74).and.(jk.eq.1)) THEN
+                      write(numout,*) z_e13c_prod, 86400*(zprorca(ji,jj,jk)/(rtrn + trb(ji,jj,jk,jpphy))/rfact2) , (zh2co3(ji,jj,jk)) , &
+               &                      ((86400 * zprorca(ji,jj,jk) / (rtrn + trb(ji,jj,jk,jpphy) ) &
+               &                             / ( rfact2 * (zh2co3(ji,jj,jk)/1025*1.e9 + rtrn)  ) - 0.371 ) / (-0.015))
+              endif
+
+
+              zr13 = ( 1. - z_e13c_prod/1000. ) * trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn )  !&
+              !&                                 * za_g(ji,jj,jk) / (za_dic(ji,jj,jk) + rtrn)
+              zr13_2 = ( 1. - z_e13c_prod2/1000. ) * trb(ji,jj,jk,jp13dic) / (trb(ji,jj,jk,jpdic) + rtrn ) !&
+              !&                                 * za_g(ji,jj,jk) / (za_dic(ji,jj,jk) + rtrn)
+#endif
+
               tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) - zprorca(ji,jj,jk) - zprorcad(ji,jj,jk)
               tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) - zpronew(ji,jj,jk) - zpronewd(ji,jj,jk)
@@ -408,4 +514,25 @@
               tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * ( zpronew(ji,jj,jk) + zpronewd(ji,jj,jk) ) &
                  &                                      - rno3 * ( zproreg + zproreg2 )
+#ifdef key_n15
+
+              tra(ji,jj,jk,jp15no3) = tra(ji,jj,jk,jp15no3) - zr15_new * ( zpronew(ji,jj,jk) + zpronewd(ji,jj,jk) )
+              tra(ji,jj,jk,jp15nh4) = tra(ji,jj,jk,jp15nh4) - zr15_reg * ( zproreg + zproreg2 )
+              tra(ji,jj,jk,jp15phy) = tra(ji,jj,jk,jp15phy) + zr15_new * zpronew(ji,jj,jk) * texcret  &
+                                                            + zr15_reg * zproreg * texcret    
+              tra(ji,jj,jk,jp15dia) = tra(ji,jj,jk,jp15dia) + zr15_new * zpronewd(ji,jj,jk) * texcret2  &
+                                                            + zr15_reg * zproreg2 * texcret2 
+              tra(ji,jj,jk,jp15doc) = tra(ji,jj,jk,jp15doc) + zr15_new * zpronew(ji,jj,jk) * excret     &
+                                                            + zr15_reg * zproreg * excret     &
+                                                            + zr15_new * zpronewd(ji,jj,jk) * excret2   &
+                                                            + zr15_reg * zproreg2 * excret2
+#endif
+
+#ifdef key_c13
+              tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) - zr13 * zprorca(ji,jj,jk) - zr13_2 * zprorcad(ji,jj,jk) 
+              tra(ji,jj,jk,jp13phy) = tra(ji,jj,jk,jp13phy) + zr13 * zprorca(ji,jj,jk) * texcret  
+              tra(ji,jj,jk,jp13dia) = tra(ji,jj,jk,jp13dia) + zr13_2 * zprorcad(ji,jj,jk) * texcret2 
+              tra(ji,jj,jk,jp13doc) = tra(ji,jj,jk,jp13doc) + zr13 * zprorca(ji,jj,jk) * excret       &   
+                                                            + zr13_2 * zprorcad(ji,jj,jk) * excret2
+#endif
           END DO
         END DO
@@ -548,4 +675,7 @@
      CALL wrk_dealloc( jpi, jpj, jpk, zpislopead, zpislopead2, zprdia, zprbio, zprdch, zprnch, zysopt            ) 
      CALL wrk_dealloc( jpi, jpj, jpk, zprorca, zprorcad, zprofed, zprofen, zprochln, zprochld, zpronew, zpronewd )
+#ifdef key_c13
+     CALL wrk_dealloc( jpi, jpj, jpk, za_g, za_dic, zh2co3, z_co3)
+#endif
      !
      IF( nn_timing == 1 )  CALL timing_stop('p4z_prod')
@@ -567,5 +697,13 @@
       !
       NAMELIST/nampisprod/ pislope, pislope2, xadap, ln_newprod, bresp, excret, excret2,  &
+#if defined key_n15
+         &                 e15n_prod, &
+#endif
+#if defined key_c13
+         &                 e13c_min, e13c_max, &
+#endif
+
          &                 chlcnm, chlcdm, chlcmin, fecnm, fecdm, grosip
+
       INTEGER :: ios                 ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
@@ -599,4 +737,12 @@
          WRITE(numout,*) '    Maximum Fe/C in nanophytoplankton         fecnm        =', fecnm
          WRITE(numout,*) '    Minimum Fe/C in diatoms                   fecdm        =', fecdm
+#if defined key_n15
+         WRITE(numout,*) '    N15 assimilation fractionation            e15n_prod    =', e15n_prod
+#endif
+#if defined key_c13
+         WRITE(numout,*) '    C13 assimilation fractionation min       e13c_min    =', e13c_min
+         WRITE(numout,*) '    C13 assimilation fractionation max       e13c_max   =', e13c_max
+#endif
+
       ENDIF
       !
@@ -629,3 +775,3 @@
 
    !!======================================================================
-END MODULE p4zprod
+END MODULE  p4zprod
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90	(revision 6534)
@@ -45,5 +45,8 @@
    REAL(wp), PUBLIC ::  xsilab     !: fraction of labile biogenic silica 
    REAL(wp), PUBLIC ::  oxymin     !: halk saturation constant for anoxia 
-
+#if defined key_n15
+   REAL(wp), PUBLIC ::  e15n_denit !: N15 denitrification fractionation 
+   REAL(wp), PUBLIC ::  e15n_nit !: N15 nitrification fractionation 
+#endif
 
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   denitr     !: denitrification array
@@ -141,8 +144,8 @@
                ! -----------------------------------------------------
                zolimit = zremik * ( 1.- nitrfac(ji,jj,jk) ) * trb(ji,jj,jk,jpdoc) 
-               zolimi(ji,jj,jk) = MIN( ( trb(ji,jj,jk,jpoxy) - rtrn ) / o2ut, zolimit ) 
+               zolimi(ji,jj,jk) = MIN( ( trb(ji,jj,jk,jpoxy) - 1.E-8 ) / o2ut, zolimit ) 
                ! Ammonification in suboxic waters with denitrification
                ! -------------------------------------------------------
-               denitr(ji,jj,jk)  = MIN(  ( trb(ji,jj,jk,jpno3) - rtrn ) / rdenit,   &
+               denitr(ji,jj,jk)  = MIN(  ( trb(ji,jj,jk,jpno3) - 1.E-8 ) / rdenit,   &
                   &                     zremik * nitrfac(ji,jj,jk) * trb(ji,jj,jk,jpdoc)  )
                !
@@ -169,8 +172,19 @@
                ! Update of the tracers trends
                ! ----------------------------
+               denitnh4(ji,jj,jk)  = MIN(  ( trb(ji,jj,jk,jpno3) - 1.E-8 ) / rdenita,   &
+                  &                     denitnh4(ji,jj,jk)  )
+               denitnh4(ji,jj,jk) = MAX(0., denitnh4(ji,jj,jk) )
+
                tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) - zonitr - denitnh4(ji,jj,jk)
                tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zonitr - rdenita * denitnh4(ji,jj,jk)
                tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2nit * zonitr
                tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2 * rno3 * zonitr + rno3 * ( rdenita - 1. ) * denitnh4(ji,jj,jk)
+# if defined key_n15
+               tra(ji,jj,jk,jp15nh4) = tra(ji,jj,jk,jp15nh4) - ( zonitr * (1. - e15n_nit / 1000.) + denitnh4(ji,jj,jk) ) &
+                                 &     *   trb(ji,jj,jk,jp15nh4) / (trb(ji,jj,jk,jpnh4) + rtrn )
+               tra(ji,jj,jk,jp15no3) = tra(ji,jj,jk,jp15no3) + zonitr * (1. - e15n_nit / 1000.) &  
+                                 &     *  trb(ji,jj,jk,jp15nh4) / (trb(ji,jj,jk,jpnh4) + rtrn )  &
+                                 &    - rdenita * denitnh4(ji,jj,jk) *   trb(ji,jj,jk,jp15no3) / (trb(ji,jj,jk,jpno3) + rtrn )
+# endif
             END DO
          END DO
@@ -242,4 +256,11 @@
                tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zorem
                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer
+#if defined key_n15
+               tra(ji,jj,jk,jp15doc) = tra(ji,jj,jk,jp15doc) + zorem * trb(ji,jj,jk,jp15poc) / (trb(ji,jj,jk,jppoc) + rtrn )
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13doc) = tra(ji,jj,jk,jp13doc) + zorem * trb(ji,jj,jk,jp13poc) / (trb(ji,jj,jk,jppoc) + rtrn )
+#endif
+
 #if defined key_kriest
                tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zorem
@@ -247,8 +268,18 @@
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zofer
 #else
-               tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zorem2 - zorem
-               tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) - zorem2
+               tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zorem2 - zorem 
+               tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) - zorem2 
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zofer2 - zofer
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zofer2
+#if defined key_n15
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) + zorem2 * trb(ji,jj,jk,jp15goc) / (trb(ji,jj,jk,jpgoc) + rtrn ) &
+                                      - zorem * trb(ji,jj,jk,jp15poc) / (trb(ji,jj,jk,jppoc) + rtrn )
+               tra(ji,jj,jk,jp15goc) = tra(ji,jj,jk,jp15goc) - zorem2 * trb(ji,jj,jk,jp15goc) / (trb(ji,jj,jk,jpgoc) + rtrn )
+#endif
+#if defined key_c13
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) + zorem2 * trb(ji,jj,jk,jp13goc) / (trb(ji,jj,jk,jpgoc) + rtrn ) &
+                                      - zorem * trb(ji,jj,jk,jp13poc) / (trb(ji,jj,jk,jppoc) + rtrn )
+               tra(ji,jj,jk,jp13goc) = tra(ji,jj,jk,jp13goc) - zorem2 * trb(ji,jj,jk,jp13goc) / (trb(ji,jj,jk,jpgoc) + rtrn )
+#endif
 #endif
 
@@ -313,4 +344,19 @@
          tra(:,:,jk,jpdic) = tra(:,:,jk,jpdic) + zolimi (:,:,jk) + denitr(:,:,jk)
          tra(:,:,jk,jptal) = tra(:,:,jk,jptal) + rno3 * ( zolimi(:,:,jk) + ( rdenit + 1.) * denitr(:,:,jk) )
+#if defined key_n15
+         tra(:,:,jk,jp15nh4) = tra(:,:,jk,jp15nh4) + ( zolimi (:,:,jk) + denitr(:,:,jk) ) &
+                                                   * trb(:,:,jk,jp15doc) / ( trb(:,:,jk,jpdoc) + rtrn )
+         tra(:,:,jk,jp15no3) = tra(:,:,jk,jp15no3) - (1. - e15n_denit / 1000.) * denitr (:,:,jk) * rdenit &
+                                                   * trb(:,:,jk,jp15no3) / ( trb(:,:,jk,jpno3) + rtrn )
+         tra(:,:,jk,jp15doc) = tra(:,:,jk,jp15doc) - ( zolimi (:,:,jk) + denitr(:,:,jk) ) &
+                                                   * trb(:,:,jk,jp15doc) / ( trb(:,:,jk,jpdoc) + rtrn )
+#endif
+#if defined key_c13
+         tra(:,:,jk,jp13dic) = tra(:,:,jk,jp13dic) + ( zolimi (:,:,jk) + denitr(:,:,jk) ) &
+                                                   * trb(:,:,jk,jp13doc) / ( trb(:,:,jk,jpdoc) + rtrn )
+         tra(:,:,jk,jp13doc) = tra(:,:,jk,jp13doc) - ( zolimi (:,:,jk) + denitr(:,:,jk) ) &
+                                                   * trb(:,:,jk,jp13doc) / ( trb(:,:,jk,jpdoc) + rtrn )
+#endif
+
       END DO
 
@@ -357,6 +403,11 @@
       !!
       !!----------------------------------------------------------------------
+#if ! defined key_n15
       NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab,   &
       &                   oxymin
+#else
+      NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab,   &
+      &                   oxymin, e15n_denit, e15n_nit
+#endif
       INTEGER :: ios                 ! Local integer output status for namelist read
 
@@ -381,4 +432,8 @@
          WRITE(numout,*) '    NH4 nitrification rate                    nitrif    =', nitrif
          WRITE(numout,*) '    halk saturation constant for anoxia       oxymin    =', oxymin
+#if defined key_n15
+         WRITE(numout,*) '    N15 denitrification fractionation        e15n_denit =', e15n_denit
+         WRITE(numout,*) '    N15 nitrification fractionation          e15n_nit =', e15n_nit
+#endif
       ENDIF
       !
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90	(revision 6534)
@@ -42,4 +42,15 @@
    REAL(wp), PUBLIC  :: concfediaz  !: Fe half-saturation Cste for diazotrophs 
    REAL(wp)          :: hratio      !: Fe:3He ratio assumed for vent iron supply
+#ifdef key_n15
+   REAL(wp), PUBLIC  :: e15n_bdenit  !: N15 Benthic denit fractionation
+   REAL(wp), PUBLIC  :: d15n_nfix    !: dN15 Nfix
+   REAL(wp), PUBLIC  :: d15n_river   !: dN15 river input
+   REAL(wp), PUBLIC  :: d15n_depo    !: dN15 depo
+#endif
+#ifdef key_c13
+   REAL(wp), PUBLIC  :: d13c_rivdoc   !: dC13 river doc input
+   REAL(wp), PUBLIC  :: d13c_rivdic   !: dC13 river dic input
+#endif
+
 
    LOGICAL , PUBLIC  :: ll_sbc
@@ -204,9 +215,16 @@
       TYPE(FLD_N) ::   sn_riverdoc, sn_riverdic, sn_riverdsi   ! informations about the fields to be read
       TYPE(FLD_N) ::   sn_riverdin, sn_riverdon, sn_riverdip, sn_riverdop
-      !
+      
       NAMELIST/nampissbc/cn_dir, sn_dust, sn_solub, sn_riverdic, sn_riverdoc, sn_riverdin, sn_riverdon,     &
         &                sn_riverdip, sn_riverdop, sn_riverdsi, sn_ndepo, sn_ironsed, sn_hydrofe, &
         &                ln_dust, ln_solub, ln_river, ln_ndepo, ln_ironsed, ln_ironice, ln_hydrofe,    &
+#ifdef key_n15
+        &                e15n_bdenit, d15n_nfix, d15n_river, d15n_depo,   &
+#endif
+#ifdef key_c13
+        &                d13c_rivdic, d13c_rivdoc,     &
+#endif
         &                sedfeinput, dustsolub, icefeinput, wdust, mfrac, nitrfix, diazolight, concfediaz, hratio
+
       !!----------------------------------------------------------------------
       !
@@ -252,4 +270,14 @@
          WRITE(numout,*) '    fe half-saturation cste for diazotrophs  concfediaz  = ', concfediaz
          WRITE(numout,*) '    Fe to 3He ratio assumed for vent iron supply hratio  = ', hratio
+#ifdef key_n15
+         WRITE(numout,*) '    N15 Benthic denit fractionation          e15n_bdenit = ', e15n_bdenit
+         WRITE(numout,*) '    dN15 nfix                                d15n_nfix   = ', d15n_nfix
+         WRITE(numout,*) '    dN15 river input                         d15n_river  = ', d15n_river
+         WRITE(numout,*) '    dN15 depo                                d15n_depo   = ', d15n_depo
+#endif
+#ifdef key_c13
+         WRITE(numout,*) '    dC13 river dic input                     d13c_rivdic = ', d13c_rivdic
+         WRITE(numout,*) '    dC13 river doc input                     d13c_rivdoc = ', d13c_rivdoc
+#endif
       END IF
 
@@ -519,3 +547,3 @@
 
    !!======================================================================
-END MODULE p4zsbc
+END MODULE  p4zsbc
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90	(revision 6534)
@@ -69,4 +69,6 @@
       REAL(wp) ::  zsiloss, zcaloss, zws3, zws4, zwsc, zdep, zwstpoc
       REAL(wp) ::  ztrfer, ztrpo4, zwdust, zlight
+      REAL(wp) :: zwstpoc15, zwstpoc13
+
       !
       CHARACTER (len=25) :: charout
@@ -175,4 +177,11 @@
                   tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) +  rivdic(ji,jj) * rfact2
                   tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) +  ( rivalk(ji,jj) - rno3 * rivdin(ji,jj) ) * rfact2
+#if defined key_n15
+                  tra(ji,jj,jk,jp15no3) = tra(ji,jj,jk,jp15no3) + ( 1. + d15n_river / 1000. ) * rivdin(ji,jj) * rfact2
+#endif
+#if defined key_c13
+                  tra(ji,jj,jk,jp13dic) = tra(ji,jj,jk,jp13dic) + ( 1. + d13c_rivdic / 1000. ) * rivdic(ji,jj) * rfact2
+#endif
+
                ENDDO
             ENDDO
@@ -185,4 +194,7 @@
          tra(:,:,1,jpno3) = tra(:,:,1,jpno3) + nitdep(:,:) * rfact2
          tra(:,:,1,jptal) = tra(:,:,1,jptal) - rno3 * nitdep(:,:) * rfact2
+#if defined key_n15
+         tra(:,:,1,jp15no3) = tra(:,:,1,jp15no3) + (1. + d15n_depo / 1000.) * nitdep(:,:) * rfact2
+#endif
       ENDIF
 
@@ -305,4 +317,11 @@
             zsedcal(ji,jj) = (1.0 - zrivalk) * zcaloss / zdep
             zsedsi (ji,jj) = (1.0 - zrivsil) * zsiloss / zdep
+#if defined key_c13
+            tra(ji,jj,ikt,jp13dic) =  tra(ji,jj,ikt,jp13dic) + zcaloss * zrivalk &
+              &                     * trb(ji,jj,ikt,jp13cal) / ( trb(ji,jj,ikt,jpcal) + rtrn)
+            tra(ji,jj,ikt,jp13cal) =  tra(ji,jj,ikt,jp13cal) - zcaloss * trb(ji,jj,ikt,jp13cal) &
+              &                     / ( trb(ji,jj,ikt,jpcal) + rtrn)
+#endif
+
 #endif
          END DO
@@ -322,4 +341,15 @@
             tra(ji,jj,ikt,jpsfe) = tra(ji,jj,ikt,jpsfe) - trb(ji,jj,ikt,jpsfe) * zws3
             zwstpoc              = trb(ji,jj,ikt,jpgoc) * zws4 + trb(ji,jj,ikt,jppoc) * zws3
+#if defined key_n15
+            tra(ji,jj,ikt,jp15goc) = tra(ji,jj,ikt,jp15goc) - trb(ji,jj,ikt,jp15goc) * zws4
+            tra(ji,jj,ikt,jp15poc) = tra(ji,jj,ikt,jp15poc) - trb(ji,jj,ikt,jp15poc) * zws3
+            zwstpoc15              = trb(ji,jj,ikt,jp15goc) * zws4 + trb(ji,jj,ikt,jp15poc) * zws3
+#endif
+#if defined key_c13
+            tra(ji,jj,ikt,jp13goc) = tra(ji,jj,ikt,jp13goc) - trb(ji,jj,ikt,jp13goc) * zws4
+            tra(ji,jj,ikt,jp13poc) = tra(ji,jj,ikt,jp13poc) - trb(ji,jj,ikt,jp13poc) * zws3
+            zwstpoc13              = trb(ji,jj,ikt,jp13goc) * zws4 + trb(ji,jj,ikt,jp13poc) * zws3
+#endif
+
 # else
             tra(ji,jj,ikt,jpnum) = tra(ji,jj,ikt,jpnum) - trb(ji,jj,ikt,jpnum) * zws4 
@@ -345,4 +375,20 @@
             sdenit(ji,jj) = rdenit * zpdenit / zdep
             zsedc(ji,jj)   = (1. - zrivno3) * zwstpoc / zdep
+#if defined key_n15
+            tra(ji,jj,ikt,jp15doc) = tra(ji,jj,ikt,jp15doc) + zwstpoc15 * zrivno3  &
+                                   - zolimit * zwstpoc15 / ( zwstpoc + rtrn ) &
+                          &        - (zpdenit + zdenitt) * zwstpoc15 / ( zwstpoc + rtrn )
+            tra(ji,jj,ikt,jp15nh4) = tra(ji,jj,ikt,jp15nh4) + zolimit * zwstpoc15 / ( zwstpoc + rtrn ) &
+                          &        + (zpdenit + zdenitt) * zwstpoc15 / ( zwstpoc + rtrn )
+            tra(ji,jj,ikt,jp15no3) = tra(ji,jj,ikt,jp15no3) - (1. - e15n_bdenit / 1000.) * rdenit * (zpdenit + zdenitt) & 
+                          &        * trb(ji,jj,ikt,jp15no3) / ( trb(ji,jj,ikt,jpno3) + rtrn )
+#endif
+#if defined key_c13
+            tra(ji,jj,ikt,jp13doc) = tra(ji,jj,ikt,jp13doc) + zwstpoc13 * zrivno3  &
+                                   - zolimit * zwstpoc13 / ( zwstpoc + rtrn ) &
+                          &        - (zpdenit + zdenitt) * zwstpoc13 / ( zwstpoc + rtrn )
+            tra(ji,jj,ikt,jp13dic) = tra(ji,jj,ikt,jp13dic) + zolimit * zwstpoc13 / ( zwstpoc + rtrn ) &
+                          &        + (zpdenit + zdenitt) * zwstpoc13 / ( zwstpoc + rtrn )
+#endif
 #endif
          END DO
@@ -385,4 +431,7 @@
                &                     * 0.002 * trb(ji,jj,jk,jpdoc) * xstep
                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + 0.002 * 4E-10 * zsoufer(ji,jj,jk) * xstep
+#if defined key_n15
+               tra(ji,jj,jk,jp15nh4) = tra(ji,jj,jk,jp15nh4) +  (1. + d15n_nfix / 1000.)  * zfact
+#endif
            END DO
          END DO 
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90	(revision 6534)
@@ -40,5 +40,10 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sinkfer2           !: Big iron sinking fluxes
 #endif
-
+#if defined key_n15
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sink15n, sink15n2   !: N15 POC sinking fluxes
+#endif
+#if defined key_c13
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sink13c, sink13c2, sink13cal   !: C13 POC sinking fluxes
+#endif
    INTEGER  :: ik100
 
@@ -94,4 +99,10 @@
       REAL(wp) ::   zagg , zaggfe, zaggdoc, zaggdoc2, zaggdoc3
       REAL(wp) ::   zfact, zwsmax, zmax, zstep
+#if defined key_n15 
+      REAL(wp) :: zagg_15, zratdoc_15
+#endif
+#if defined key_c13
+      REAL(wp) :: zagg_13, zratdoc_13
+#endif
       CHARACTER (len=25) :: charout
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zw3d
@@ -172,5 +183,4 @@
       sinksil (:,:,:) = 0.e0
       sinkfer2(:,:,:) = 0.e0
-
       !   Compute the sedimentation term using p4zsink2 for all the sinking particles
       !   -----------------------------------------------------
@@ -187,4 +197,26 @@
       END DO
 
+#if defined key_n15
+      sink15n (:,:,:) = 0.e0
+      sink15n2(:,:,:) = 0.e0
+      DO jit = 1, iiter1
+        CALL p4z_sink2( wsbio3, sink15n , jp15poc, iiter1 )
+      END DO
+      DO jit = 1, iiter2
+        CALL p4z_sink2( wsbio4, sink15n2, jp15goc, iiter2 )
+      END DO
+#endif
+#ifdef key_c13
+      sink13c (:,:,:) = 0.e0
+      sink13c2(:,:,:) = 0.e0
+      sink13cal(:,:,:) = 0.e0
+      DO jit = 1, iiter1
+        CALL p4z_sink2( wsbio3, sink13c , jp13poc, iiter1 )
+      END DO
+      DO jit = 1, iiter2
+        CALL p4z_sink2( wsbio4, sink13c2, jp13goc, iiter2 )
+        CALL p4z_sink2( wscal , sink13cal,jp13cal, iiter2 )
+      END DO
+#endif
       !  Exchange between organic matter compartments due to coagulation/disaggregation
       !  ---------------------------------------------------
@@ -230,4 +262,21 @@
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zaggfe
                tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) - zaggdoc - zaggdoc2 - zaggdoc3
+#if defined key_n15
+               zagg_15 = zagg * trb(ji,jj,jk,jp15poc) / ( trb(ji,jj,jk,jppoc) + rtrn )
+               zratdoc_15 = trb(ji,jj,jk,jp15doc)/ ( trb(ji,jj,jk,jpdoc) + rtrn )
+               tra(ji,jj,jk,jp15poc) = tra(ji,jj,jk,jp15poc) - zagg_15 + ( zaggdoc + zaggdoc3 ) * zratdoc_15
+               tra(ji,jj,jk,jp15goc) = tra(ji,jj,jk,jp15goc) + zagg_15 + zaggdoc2 * zratdoc_15
+               tra(ji,jj,jk,jp15doc) = tra(ji,jj,jk,jp15doc) - ( zaggdoc + zaggdoc2 + zaggdoc3 )   &
+               &                       * zratdoc_15
+#endif
+#if defined key_c13
+               zagg_13 = zagg * trb(ji,jj,jk,jp13poc) / ( trb(ji,jj,jk,jppoc) + rtrn )
+               zratdoc_13 = trb(ji,jj,jk,jp13doc)/ ( trb(ji,jj,jk,jpdoc) + rtrn)
+               tra(ji,jj,jk,jp13poc) = tra(ji,jj,jk,jp13poc) - zagg_13 + ( zaggdoc + zaggdoc3 ) * zratdoc_13
+               tra(ji,jj,jk,jp13goc) = tra(ji,jj,jk,jp13goc) + zagg_13 + zaggdoc2 * zratdoc_13
+               tra(ji,jj,jk,jp13doc) = tra(ji,jj,jk,jp13doc) - ( zaggdoc + zaggdoc2 + zaggdoc3 )   &
+               &                       * zratdoc_13
+#endif
+
                !
             END DO
@@ -891,5 +940,11 @@
       ALLOCATE( wsbio3 (jpi,jpj,jpk) , wsbio4  (jpi,jpj,jpk) , wscal(jpi,jpj,jpk) ,     &
          &      sinking(jpi,jpj,jpk) , sinking2(jpi,jpj,jpk)                      ,     &                
-         &      sinkcal(jpi,jpj,jpk) , sinksil (jpi,jpj,jpk)                      ,     &                
+         &      sinkcal(jpi,jpj,jpk) , sinksil (jpi,jpj,jpk)                      ,     &  
+#if defined key_n15
+         &      sink15n(jpi,jpj,jpk) , sink15n2(jpi,jpj,jpk)                      ,     &
+#endif       
+#if defined key_c13
+         &      sink13c(jpi,jpj,jpk) , sink13c2(jpi,jpj,jpk), sink13cal(jpi,jpj,jpk) ,     &
+#endif        
 #if defined key_kriest
          &      xnumm(jpk)                                                        ,     &                
@@ -913,3 +968,3 @@
 
    !!======================================================================
-END MODULE p4zsink
+END MODULE  p4zsink
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90	(revision 6534)
@@ -71,4 +71,7 @@
       INTEGER ::   ji, jj, jk, jnt, jn, jl
       REAL(wp) ::  ztra
+#if defined key_n15
+      REAL(wp) ::  zn15min, zn15max
+#endif
 #if defined key_kriest
       REAL(wp) ::  zcoef1, zcoef2
@@ -133,8 +136,7 @@
          !
          CALL p4z_bio( kt, jnt )   ! Biology
+         CALL p4z_sed( kt, jnt )   ! Sedimentation
          CALL p4z_lys( kt, jnt )   ! Compute CaCO3 saturation
-         CALL p4z_sed( kt, jnt )   ! Surface and Bottom boundary conditions
          CALL p4z_flx( kt, jnt )   ! Compute surface fluxes
-         !
          xnegtr(:,:,:) = 1.e0
          DO jn = jp_pcs0, jp_pcs1
@@ -153,5 +155,5 @@
          !                                ! 
          DO jn = jp_pcs0, jp_pcs1
-           trb(:,:,:,jn) = trb(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn)
+           trb(:,:,:,jn) = trb(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn) * tmask(:,:,:)
          END DO
         !
@@ -159,5 +161,13 @@
             tra(:,:,:,jn) = 0._wp
          END DO
-         !
+        !
+#if defined key_n15
+         zn15min = -30.
+         zn15max = 30.
+         !
+         trb(:,:,:,jp15no3) = MIN(trb(:,:,:,jp15no3) , ( (zn15max/1000.)  + 1.)*trb(:,:,:,jpno3) )
+         trb(:,:,:,jp15no3) = MAX(trb(:,:,:,jp15no3) , ( (zn15min/1000.)  + 1.)*trb(:,:,:,jpno3) )
+#endif
+      !
          IF( ln_top_euler ) THEN
             DO jn = jp_pcs0, jp_pcs1
@@ -166,5 +176,5 @@
          ENDIF
       END DO
-
+      !
 #if defined key_kriest
       ! 
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90	(revision 6534)
@@ -75,4 +75,5 @@
    INTEGER, PUBLIC, PARAMETER ::   jpnh4 = 23    !: Ammonium Concentration
 
+
 #elif defined key_pisces
    !!---------------------------------------------------------------------
@@ -82,12 +83,16 @@
    LOGICAL, PUBLIC, PARAMETER ::   lk_p4z        = .TRUE.  !: p4z flag 
    LOGICAL, PUBLIC, PARAMETER ::   lk_kriest     = .FALSE. !: Kriest flag 
-   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 24      !: number of PISCES passive tracers
-   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 13      !: additional 2d output 
-   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 11      !: additional 3d output 
-   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
 
    ! assign an index in trc arrays for each LOBSTER prognostic variables
    !    WARNING: be carefull about the order when reading the restart
         !   !!gm  this warning should be obsolet with IOM
+#if ! defined key_n15  && ! defined key_c13
+   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .FALSE.  !: N15 flag 
+   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .FALSE.  !: C13 flag 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 24      !: number of PISCES passive tracers
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
+#endif
    INTEGER, PUBLIC, PARAMETER ::   jpdic =  1    !: dissolved inoganic carbon concentration 
    INTEGER, PUBLIC, PARAMETER ::   jptal =  2    !: total alkalinity 
@@ -115,4 +120,71 @@
    INTEGER, PUBLIC, PARAMETER ::   jpnh4 = 24    !: Ammonium Concentration
 
+#if defined key_n15  && defined key_c13
+   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .TRUE.  !: N15 flag 
+   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .TRUE.  !: C13 flag 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 42      !: number of PISCES passive tracers
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
+
+   INTEGER, PUBLIC, PARAMETER ::   jp15poc = 25  !: 15N small particulate organic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15phy = 26  !: 15N phytoplancton concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp15zoo = 27  !: 15N zooplancton concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15doc = 28  !: 15N dissolved organic carbon concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp15dia = 29  !: 15N Diatoms Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15mes = 30  !: 15N Mesozooplankton Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15goc = 31  !: 15N big particulate organic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15no3 = 32  !: 15N Nitrates Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15nh4 = 33  !: 15N Ammonium Concentration
+
+   INTEGER, PUBLIC, PARAMETER ::   jp13dic = 34  !: 13C dissolved inorganic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13cal = 35  !: 13C calcite Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13poc = 36  !: 13C small particulate organic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13phy = 37  !: 13C phytoplancton concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp13zoo = 38  !: 13C zooplancton concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13doc = 39  !: 13C dissolved organic carbon concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp13dia = 40  !: 13C Diatoms Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13mes = 41  !: 13C Mesozooplankton Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13goc = 42  !: 13C big particulate organic concentration
+#endif
+#if defined key_n15  && ! defined key_c13
+   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .TRUE.  !: N15 flag 
+   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .FALSE.  !: C13 flag 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 33      !: number of PISCES passive tracers
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
+
+   INTEGER, PUBLIC, PARAMETER ::   jp15poc = 25  !: 15N small particulate organic phosphate concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15phy = 26  !: 15N phytoplancton concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp15zoo = 27  !: 15N zooplancton concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15doc = 28  !: 15N dissolved organic carbon concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp15dia = 29  !: 15N Diatoms Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15mes = 30  !: 15N Mesozooplankton Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15goc = 31  !: 15N big particulate organic phosphate concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15no3 = 32  !: 15N Nitrates Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp15nh4 = 33  !: 15N Ammonium Concentration
+#endif
+#if ! defined key_n15  &&  defined key_c13
+   LOGICAL, PUBLIC, PARAMETER ::   lk_n15        = .FALSE.  !: N15 flag 
+   LOGICAL, PUBLIC, PARAMETER ::   lk_c13        = .TRUE.  !: C13 flag 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces     = 33      !: number of PISCES passive tracers
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_2d  = 20      !: additional 2d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_3d  = 14      !: additional 3d output 
+   INTEGER, PUBLIC, PARAMETER ::   jp_pisces_trd =  1      !: number of sms trends for PISCES
+
+   INTEGER, PUBLIC, PARAMETER ::   jp13dic = 25  !: 13C dissolved inorganic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13cal = 26  !: 13C calcite Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13poc = 27  !: 13C small particulate organic concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13phy = 28  !: 13C phytoplancton concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp13zoo = 29  !: 13C zooplancton concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13doc = 30  !: 13C dissolved organic carbon concentration 
+   INTEGER, PUBLIC, PARAMETER ::   jp13dia = 31  !: 13C Diatoms Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13mes = 32  !: 13C Mesozooplankton Concentration
+   INTEGER, PUBLIC, PARAMETER ::   jp13goc = 33  !: 13C big particulate organic concentration
+#endif
+
+#endif
+
 #else
    !!---------------------------------------------------------------------
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90	(revision 6534)
@@ -159,4 +159,27 @@
       END IF
 
+#if defined key_n15
+      trn(:,:,:,jp15doc) = trn(:,:,:,jpdoc)
+      trn(:,:,:,jp15poc) = trn(:,:,:,jppoc)
+      trn(:,:,:,jp15goc) = trn(:,:,:,jpgoc)
+      trn(:,:,:,jp15phy) = trn(:,:,:,jpphy)
+      trn(:,:,:,jp15dia) = trn(:,:,:,jpdia)
+      trn(:,:,:,jp15zoo) = trn(:,:,:,jpzoo)
+      trn(:,:,:,jp15mes) = trn(:,:,:,jpmes)
+      trn(:,:,:,jp15no3) = trn(:,:,:,jpno3)
+      trn(:,:,:,jp15nh4) = trn(:,:,:,jpnh4)
+#endif
+#if defined key_c13
+      trn(:,:,:,jp13doc) = trn(:,:,:,jpdoc)
+      trn(:,:,:,jp13poc) = trn(:,:,:,jppoc)
+      trn(:,:,:,jp13goc) = trn(:,:,:,jpgoc)
+      trn(:,:,:,jp13phy) = trn(:,:,:,jpphy)
+      trn(:,:,:,jp13dia) = trn(:,:,:,jpdia)
+      trn(:,:,:,jp13zoo) = trn(:,:,:,jpzoo)
+      trn(:,:,:,jp13mes) = trn(:,:,:,jpmes)
+      trn(:,:,:,jp13dic) = trn(:,:,:,jpdic)
+      trn(:,:,:,jp13cal) = trn(:,:,:,jpcal)
+#endif
+
 
       CALL p4z_sink_init      !  vertical flux of particulate organic matter
Index: /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90
===================================================================
--- /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90	(revision 6533)
+++ /branches/CNRS/dev_r6526_PISCES_ISO/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90	(revision 6534)
@@ -47,4 +47,7 @@
          zfact = 1.0e+6 
          IF( jn == jpno3 .OR. jn == jpnh4 ) zfact = rno3 * 1.0e+6 
+#if defined key_n15
+         IF( jn == jp15no3 .OR. jn == jp15nh4 ) zfact = rno3 * 1.0e+6
+#endif
          IF( jn == jppo4  )                 zfact = po4r * 1.0e+6
          cltra = TRIM( ctrcnm(jn) )                  ! short title for tracer
