Changeset 2451


Ignore:
Timestamp:
03/27/15 13:27:43 (7 years ago)
Author:
aclsce
Message:

Modifications to use new NEMO revision in IPSLCM6_rc1 configuration.

  • Removed any SOURCES/NEMO not needed anymore (modifs commited in NEMO reference version).
  • Adapted *.xml, namelist, AA_make, *.card, *.driver related to NEMO to fit with new NEMO reference version.
  • Removed files not needed anymore : files related to L19 LMDZ vertical levels
Location:
CONFIG/UNIFORM/v6/IPSLCM6
Files:
29 deleted
20 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/IPSLCM6/AA_make

    r2390 r2451  
    121121 
    122122orca2lim2:  ../../modeles/NEMOGCM/CONFIG 
    123         (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM_PISCES -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim2_vp key_coupled key_oasis3 key_oasis_mct key_diahth key_diaar5" del_key "key_nosignedzero " ; cp ORCA2_LIM_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
     123        (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM_PISCES -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim2_vp key_oasis3 key_diahth key_diaar5" del_key "key_nosignedzero " ; cp ORCA2_LIM_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
    124124 
    125125orca2lim2nopisces:  ../../modeles/NEMOGCM/CONFIG 
    126         (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim2_vp key_coupled key_oasis3 key_oasis_mct key_diahth key_diaar5" del_key "key_nosignedzero" ; cp ORCA2_LIM/BLD/bin/nemo.exe ../../../bin/opa ; ) 
     126        (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim2_vp key_oasis3 key_diahth key_diaar5" del_key "key_nosignedzero" ; cp ORCA2_LIM/BLD/bin/nemo.exe ../../../bin/opa ; ) 
    127127 
    128128orca2lim3:  ../../modeles/NEMOGCM/CONFIG 
    129         (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM3_PISCES -d "OPA_SRC LIM_SRC_3 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim3 key_coupled key_oasis3 key_oasis_mct key_top key_pisces key_diahth key_diaar5" del_key "key_lim2 key_nosignedzero" ; cp ORCA2_LIM3_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
     129        (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA2_LIM3_PISCES -d "OPA_SRC LIM_SRC_3 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim3 key_oasis3 key_top key_pisces key_diahth key_diaar5" del_key "key_lim2 key_nosignedzero" ; cp ORCA2_LIM3_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
    130130 
    131131orca1lim3: ../../modeles/NEMOGCM/CONFIG 
    132         (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA1_LIM3_PISCES -d "OPA_SRC LIM_SRC_3 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim3 key_oasis3 key_coupled key_oasis_mct  key_top key_pisces key_diahth key_diaar5 key_vvl" del_key "key_lim2 key_lim2_vp key_nosignedzero" ; cp ORCA1_LIM3_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
     132        (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA1_LIM3_PISCES -d "OPA_SRC LIM_SRC_3 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_lim3 key_oasis3 key_top key_pisces key_diahth key_diaar5 key_vvl" del_key "key_lim2 key_lim2_vp key_nosignedzero" ; cp ORCA1_LIM3_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
    133133 
    134134orca1lim2:  ../../modeles/NEMOGCM/CONFIG 
    135         (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA1_LIM2_PISCES -d "OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_top key_pisces key_lim2 key_lim2_vp key_coupled key_oasis3 key_oasis_mct key_diahth key_diaar5 key_vvl" del_key "key_lim3 key_nosignedzero " ; cp ORCA1_LIM2_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
     135        (cd ../../modeles/NEMOGCM/CONFIG; cp ../../../config/IPSLCM6/SOURCES/NEMO/arch-${FCM_ARCH}.fcm ../ARCH/. ; cp ../../../config/IPSLCM6/SOURCES/NEMO/*.cfg ../TOOLS/COMPILE/. ; ./makenemo -m ${FCM_ARCH} -n ORCA1_LIM2_PISCES -d "OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC" -e "$(MYSRC_DIR)/SOURCES/NEMO" -j16 add_key "key_top key_pisces key_lim2 key_lim2_vp key_oasis3 key_diahth key_diaar5 key_vvl" del_key "key_lim3 key_nosignedzero " ; cp ORCA1_LIM2_PISCES/BLD/bin/nemo.exe ../../../bin/opa ; ) 
    136136 
    137137verif: ../../bin 
  • CONFIG/UNIFORM/v6/IPSLCM6/EXPERIMENTS/IPSLCM/piControl/COMP/opa9.card

    r2446 r2451  
    5555        (${config_UserChoices_JobName}_1m_grid_U.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_U.nc, Post_1M_grid_U),\ 
    5656        (${config_UserChoices_JobName}_1m_grid_V.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_V.nc, Post_1M_grid_V),\ 
    57         (${config_UserChoices_JobName}_1m_grid_W.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_W.nc, NONE),\ 
     57        (${config_UserChoices_JobName}_1m_grid_W.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_W.nc, Post_1M_grid_W),\ 
    5858        (${config_UserChoices_JobName}_1m_diaptr.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_diaptr.nc, Post_1M_diaptr),\ 
    5959        (${config_UserChoices_JobName}_1d_grid_T.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_grid_T.nc, Post_1D_grid_T),\ 
     
    109109TimeSeriesVars3D = (difvmo, difvmto, difvtrto, voddmavs, voveeivw, votkeavt, vovecrtz, wmo, wmosq) 
    110110ChunckJob3D = 50Y 
    111 Seasonal=OFF 
     111Seasonal=ON 
    112112 
    113113[Post_1M_diaptr] 
    114114Patches = () 
    115 GatherWithInternal = (nav_lat, time_counter, deptht, time_counter_bounds) 
     115GatherWithInternal = (nav_lat, time_counter, deptht, depthw, time_counter_bounds) 
    116116TimeSeriesVars2D = (zotemglo, zosalglo, zomsfglo, zotematl, zosalatl, zomsfatl, zotempac, zosalpac, zomsfpac, zotemind, zosalind, zomsfind, zotemipc, zosalipc, zomsfipc, sophtadv, sophtldf, sopstadv) 
    117117ChunckJob2D = NONE 
     
    198198[Post_1Y_diaptr] 
    199199Patches = () 
    200 GatherWithInternal = (nav_lat, time_counter, deptht, time_counter_bounds) 
     200GatherWithInternal = (nav_lat, time_counter, deptht, depthw, time_counter_bounds) 
    201201TimeSeriesVars2D = (zotemglo, zosalglo, zomsfglo, zotematl, zosalatl, zomsfatl, zotempac, zosalpac, zomsfpac, zotemind, zosalind, zomsfind, zotemipc, zosalipc, zomsfipc, sophtadv, sophtldf, sopstadv) 
    202202ChunckJob2D = NONE 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/DRIVER/opa9.driver

    r2384 r2451  
    144144    fi 
    145145     
    146     ##--Write Frequency Purpose .... 
    147     frequency=$( echo ${config_OCE_WriteFrequency} | awk "-F " '{print $1}' ) 
    148     factor=$( echo ${frequency} | sed -e "s/[yYmMdD]//" ) 
    149     case ${frequency} in 
    150         1Y|1y) 
    151             (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInYear  ${year} ) ))          ; 
    152             PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}y ; 
    153             R_OUT_OCE_NWRITE=${R_OUT_OCE_O_Y} ; 
    154             WF1=${factor}Y ;; 
    155         1M|1m) 
    156             (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInMonth ${year} ${month} ) )) ; 
    157             PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}m ; 
    158             R_OUT_OCE_NWRITE=${R_OUT_OCE_O_M} ; 
    159             WF1=${factor}M ;; 
    160         *D|*d) 
    161             (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor  ))  ; 
    162             PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}d ; 
    163             R_OUT_OCE_NWRITE=${R_OUT_OCE_O_D} ; 
    164             WF1=${factor}D ;; 
    165         *) 
    166             (( ORCA_NWRITE = 0 ))                                                            ;; 
    167     esac 
    168  
    169146    ##--Write Frequency for iomput 
    170147    ## Differents frequencies are allowed for grid_[TUVW] and icemod files 
     
    266243        -e "s%${PAT_NLEAPY}%       nn_leapy=${ORCA_NLEAPY}%"                    \ 
    267244        -e "s/${PAT_NSTOCK}/       nn_stock=${ORCA_NSTOCK}/"                    \ 
    268         -e "s/${PAT_NWRITE}/       nn_write=${ORCA_NWRITE}/"                    \ 
    269245        -e "s/${PAT_RESTAR}/       ln_rstart=${ORCA_LRSTAR}/"                   \ 
    270246        -e "s/${PAT_NRSTAR}/       nn_rstctl=${ORCA_NRSTDT}/"                   \ 
    271247        -e "s/${PAT_NMSH}/         nn_msh=${ORCA_NMSH}/"                        \ 
    272         -e "s/${PAT_NN_WRITE}/     nn_prg=${ORCA_NWRITE}/"                      \ 
    273         -e "s/${PAT_NN_FWRI}/      nn_fwri=${ORCA_NWRITE}/"                     \ 
    274248        -e "s/${PAT_NN_FSBC}/      nn_fsbc=${ORCA_NN_FSBC}/"                    \ 
    275249        -e "s/${PAT_ICE_EMBD}/     nn_ice_embd=${NN_ICE_EMBD}/"                 \ 
     
    323297    IGCM_debug_PushStack "OCE_Finalize" 
    324298 
    325     if [ -f date.file ] ; then 
    326         # Prefix use in opa9.card AND in lim2.card : 
    327         DATE_OPA=$( cat date.file | \ 
    328             sed "s/\ ${config_UserChoices_JobName}_[0-9]*[a-z]_\([0-9]*_[0-9]*\)_\ */\1/g" ) 
    329         MainPrefix=${config_UserChoices_JobName}_1d_${DATE_OPA} 
    330         SecondPrefix=${config_UserChoices_JobName}_1m_${DATE_OPA} 
    331     fi 
    332  
    333299    IGCM_debug_Print 1 FINALIZE OCE !!! 
    334300 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/DRIVER/pisces.driver

    r2372 r2451  
    3030    PAT_TOP_NRSTTR=$( supergrep nn_rsttr          ${NAMELIST_TOP_CFG} ) 
    3131    PAT_PIS_ATCCO2=$( supergrep atcco2            ${NAMELIST_PISCES_CFG} ) 
    32     PAT_PIS_LRIVER=$( supergrep ln_river          ${NAMELIST_PISCES_CFG} ) 
    3332 
    3433    IGCM_debug_PopStack "MBG_Initialize" 
     
    4746        TOP_NRSTTR=0 
    4847 
    49         PISCES_LRIVER=.FALSE. 
    50  
    5148    elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_MBG_Restart}" = "y" ] ) ; then 
    5249 
     
    5451        TOP_LRSTTR=.TRUE. 
    5552        TOP_NRSTTR=0 
    56  
    57         PISCES_LRIVER=.TRUE. 
    5853 
    5954        # If we start from IPSLCM5* restart files. 
     
    6762        TOP_LRSTTR=.TRUE. 
    6863        TOP_NRSTTR=2 
    69  
    70         PISCES_LRIVER=.FALSE. 
    7164 
    7265    fi 
     
    118111    IGCM_sys_Mv namelist_top_cfg.tmp namelist_top_cfg 
    119112 
    120     sed -e "s%${PAT_PIS_LRIVER}%          ln_river=${PISCES_LRIVER}%"                \ 
    121         namelist_pisces_cfg > namelist_pisces_cfg.tmp 
    122  
    123     IGCM_sys_Mv namelist_pisces_cfg.tmp namelist_pisces_cfg 
    124  
    125  
    126113    IGCM_debug_PopStack "MBG_Update" 
    127114} 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/config.def_preind

    r2276 r2451  
    1212 
    1313##  Facteur additif pour l'albedo 
    14 pmagic=0.0 
     14pmagic=-0.01 
    1515 
    1616# 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/context_nemo.xml

    r2436 r2451  
    9191        </file> 
    9292 
    93         <file id="file7" name_suffix="_icemod" description="ice variables" > 
     93        <file id="file7" name_suffix="_grid_W" description="ocean W grid variables" > 
     94          <field field_ref="woce"         name="vovecrtz" /> 
     95          <field field_ref="avt"          name="votkeavt" /> 
     96          <field field_ref="avs"          name="voddmavs" /> 
     97        </file> 
     98 
     99        <file id="file8" name_suffix="_icemod" description="ice variables" > 
    94100          <field field_ref="ice_pres"                     /> 
    95101          <field field_ref="snowthic_cea" name="isnowthi" /> 
     
    106112        </file> 
    107113 
    108     <file id="file8" name_suffix="_diaptr" description="diaptr variables" > 
     114    <file id="file9" name_suffix="_diaptr" description="diaptr variables" > 
    109115       <field field_ref="zotemglo"   /> 
    110116       <field field_ref="zosalglo"     /> 
    111117       <field field_ref="zosrfglo"    /> 
    112118       <field field_ref="zomsfglo"    /> 
    113        <field field_ref="zomsfeiv"    /> 
    114119       <field field_ref="zotematl"   /> 
    115120       <field field_ref="zosalatl"     /> 
     
    129134       <field field_ref="zomsfipc"     /> 
    130135       <field field_ref="sophtadv"      /> 
    131        <field field_ref="sophteiv"      /> 
    132136       <field field_ref="sophtldf"      /> 
    133        <field field_ref="sophtove"      /> 
    134        <field field_ref="sohtatl"       /> 
    135        <field field_ref="sohtpac"       /> 
    136        <field field_ref="sohtind"       /> 
    137        <field field_ref="sohtipc"       /> 
    138137       <field field_ref="sopstadv"      /> 
    139        <field field_ref="sopsteiv"      /> 
    140138       <field field_ref="sopstldf"      /> 
    141        <field field_ref="sopstove"      /> 
    142        <field field_ref="sostatl"       /> 
    143        <field field_ref="sostpac"       /> 
    144        <field field_ref="sostind"       /> 
    145        <field field_ref="sostipc"       /> 
    146139    </file> 
    147140 
    148         <file id="file9" name_suffix="_ptrc_T" description="pisces sms variables" > 
     141        <file id="file10" name_suffix="_ptrc_T" description="pisces sms variables" > 
    149142       <field field_ref="DIC"      /> 
    150143       <field field_ref="Alkalini" /> 
     
    173166        </file> 
    174167         
    175         <file id="file10" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     168        <file id="file11" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    176169       <field field_ref="PH"       /> 
    177170       <field field_ref="CO3"      /> 
     
    224217      <file_group id="1y"  output_freq="1y" output_level="10" enabled="_1Y_ENABLE_"> <!-- real yearly files --> 
    225218 
    226         <file id="file11" name_suffix="_ptrc_T" description="pisces sms variables" > 
     219        <file id="file12" name_suffix="_ptrc_T" description="pisces sms variables" > 
    227220       <field field_ref="DIC"      /> 
    228221       <field field_ref="Alkalini" /> 
     
    251244        </file> 
    252245 
    253         <file id="file12" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     246        <file id="file13" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    254247       <field field_ref="PH"       /> 
    255248       <field field_ref="CO3"      /> 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/field_def_nemo.xml

    r2372 r2451  
    1616       
    1717     <field_group id="grid_T" grid_ref="grid_T_2D" > 
    18          <field id="toce"         long_name="temperature"                               unit="degC" grid_ref="grid_T_3D"/> 
    19          <field id="soce"         long_name="salinity"                                  unit="psu"  grid_ref="grid_T_3D"/> 
     18         <field id="e3t"          long_name="T-cell thickness"                          unit="m"      grid_ref="grid_T_3D"/> 
     19 
     20         <field id="toce"         long_name="temperature"                               unit="degC"   grid_ref="grid_T_3D"/> 
     21         <field id="toce_e3t"     long_name="temperature * e3t"                         unit="degC*m" grid_ref="grid_T_3D" > toce * e3t </field > 
     22         <field id="soce"         long_name="salinity"                                  unit="psu"    grid_ref="grid_T_3D"/> 
     23         <field id="soce_e3t"     long_name="salinity * e3t"                            unit="psu*m"  grid_ref="grid_T_3D" > soce * e3t </field > 
     24 
    2025         <field id="sst"          long_name="sea surface temperature"                   unit="degC"                     /> 
    21          <field id="sst2"         long_name="square of sea surface temperature"         unit="degC2"                    /> 
     26         <field id="sst2"         long_name="square of sea surface temperature"         unit="degC2"                     >  sst * sst </field > 
     27         <field id="sstmax"       long_name="max of sea surface temperature"       field_ref="sst"  operation="maximum" /> 
     28         <field id="sstmin"       long_name="min of sea surface temperature"       field_ref="sst"  operation="minimum" /> 
    2229         <field id="sstgrad"      long_name="module of sst gradient"                    unit="degC/m"                   /> 
    2330         <field id="sstgrad2"     long_name="square of module of sst gradient"          unit="degC2/m2"                 /> 
     31         <field id="sbt"          long_name="sea bottom temperature"                    unit="degC"                     /> 
     32          
    2433         <field id="sss"          long_name="sea surface salinity"                      unit="psu"                      /> 
    25          <field id="sss2"         long_name="square of sea surface salinity"            unit="psu2"                     /> 
     34         <field id="sss2"         long_name="square of sea surface salinity"            unit="psu2"                      >  sss * sss </field > 
     35         <field id="sssmax"       long_name="max of sea surface salinity"          field_ref="sss"  operation="maximum" /> 
     36         <field id="sssmin"       long_name="min of sea surface salinity"          field_ref="sss"  operation="minimum" /> 
     37         <field id="sbs"          long_name="sea bottom salinity"                       unit="psu"                      /> 
     38 
    2639         <field id="ssh"          long_name="sea surface height"                        unit="m"                        /> 
    27          <field id="ssh2"         long_name="square of sea surface height"              unit="m2"                       /> 
     40         <field id="ssh2"         long_name="square of sea surface height"              unit="m2"                        >  ssh * ssh </field > 
     41         <field id="sshmax"       long_name="max of sea surface height"            field_ref="ssh"  operation="maximum" /> 
     42 
    2843         <field id="mldkz5"       long_name="mixing layer depth (Turbocline)"           unit="m"                        /> 
    2944         <field id="mldr10_1"     long_name="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        /> 
    30          <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3" grid_ref="grid_T_3D"/> 
    31          <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                      /> 
     45         <field id="mldr10_1max"  long_name="max of Mixed Layer Depth 0.01 ref.10m" field_ref="mldr10_1"  operation="maximum" /> 
     46         <field id="mldr10_1min"  long_name="min of Mixed Layer Depth 0.01 ref.10m" field_ref="mldr10_1"  operation="minimum" /> 
     47         <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                     /> 
    3248         <field id="saltc"        long_name="Salt content vertically integrated"        unit="PSU*kg/m2"                /> 
    33          <field id="eken"         long_name="kinetic energy"                            unit="m2/s2" grid_ref="grid_T_3D"/> 
    34          <field id="hdiv"         long_name="horizontal divergence"                     unit="s-1"   grid_ref="grid_T_3D"/> 
     49         <!-- EOS --> 
     50         <field id="alpha"        long_name="thermal expansion"                         unit="1/degC" grid_ref="grid_T_3D"/> 
     51         <field id="beta"         long_name="haline contraction"                        unit="1/psu"  grid_ref="grid_T_3D"/> 
     52         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"           unit="1/s"    grid_ref="grid_T_3D"/> 
     53         <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3"  grid_ref="grid_T_3D"/> 
     54         <!-- Energy - horizontal divergence --> 
     55         <field id="eken"         long_name="kinetic energy"                            unit="m2/s2"  grid_ref="grid_T_3D"/> 
     56         <field id="hdiv"         long_name="horizontal divergence"                     unit="s-1"    grid_ref="grid_T_3D"/> 
    3557         <!-- variables available with MLE --> 
    3658         <field id="Lf_NHpf"      long_name="MLE: Lf = N H / f"                         unit="m"                        /> 
     
    5274         <field id="botpres"      long_name="Pressure at sea floor"                     unit="dbar"                     /> 
    5375         <!-- variables available with key_diaar5/key_vvl --> 
    54          <field id="cellthc"      long_name="Cell thickness"                            unit="m"    grid_ref="grid_T_3D"/> 
    5576         <!-- variables available with key_vvl --> 
    5677         <field id="tpt_dep"      long_name="T-point depth"                             unit="m"    grid_ref="grid_T_3D"/> 
     
    140161 
    141162         <field id="empmr"        long_name="Net Upward Water Flux"                                        unit="kg/m2/s"  /> 
    142          <field id="saltflx"      long_name="Downward salt flux"                                           unit="PSU/m2/s"  /> 
     163         <field id="saltflx"      long_name="Downward salt flux"                                           unit="PSU/m2/s" /> 
    143164         <field id="fmmflx"       long_name="Water flux due to freezing/melting"                           unit="kg/m2/s"  /> 
    144165         <field id="snowpre"      long_name="Snow precipitation"                                           unit="kg/m2/s"  /> 
    145166         <field id="runoffs"      long_name="River Runoffs"                                                unit="Kg/m2/s"  /> 
    146          <field id="precip"       long_name="Total precipitation"                                           unit="kg/m2/s"  /> 
     167         <field id="precip"       long_name="Total precipitation"                                          unit="kg/m2/s"  /> 
    147168 
    148169 
     
    158179         <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    159180         <field id="ssh_ib"       long_name="Inverse barometer sea surface height"                         unit="m"        /> 
    160          
     181 
     182         <!-- * variable related to ice shelf forcing * --> 
     183         <field id="fwfisf"       long_name="Ice shelf melting"                                            unit="Kg/m2/s"  /> 
     184         <field id="qisf"         long_name="Ice Shelf Heat Flux"                                          unit="W/m2"     /> 
     185         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "                 unit="m/s"      /> 
     186         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "                 unit="m/s"      /> 
     187             <field id="stbl"         long_name="salinity in the Losh tbl                    "                 unit="PSU"      /> 
     188             <field id="ttbl"         long_name="temperature in the Losh tbl                 "                 unit="C"        /> 
     189 
    161190         <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
    162191         <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
     
    189218         <field id="snowmel_cea"   long_name="Snow Melt Rate (cell average)"                               unit="kg/m2/s"  /> 
    190219         <field id="sntoice_cea"   long_name="Snow-Ice Formation Rate (cell average)"                      unit="kg/m2/s"  /> 
    191          <field id="ticemel_cea"      long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"     unit="kg/m2/s"  /> 
     220         <field id="ticemel_cea"   long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"     unit="kg/m2/s"  /> 
    192221 
    193222         <!-- ice fields --> 
     
    205234         <field id="icethic_cea"  long_name="Ice thickness (cell average)"                                 unit="m"        /> 
    206235         <field id="iceprod_cea"  long_name="Ice production (cell average)"                                unit="m/s"      /> 
     236         <field id="iiceconc"     long_name="Ice concentration"                                            unit=""         /> 
    207237          
    208238         <field id="ice_pres"     long_name="Ice presence"                                                 unit="-"        /> 
     
    219249         <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                           unit="kgC/m2/s" /> 
    220250         <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                         unit="kgPSU/m2/s" />         
    221           
    222           
     251        
    223252         <field id="iceconc"      long_name="ice concentration"                                            unit="%"        /> 
    224253         <field id="uice_ipa"     long_name="Ice velocity along i-axis at I-point (ice presence average)"  unit="m/s"      /> 
     
    246275 
    247276         <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     /> 
    248          <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"   />  
     277         <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"    />  
    249278         <field id="isnowhc"      long_name="snow total heat content"                                      unit="10^9J"    /> 
    250279         <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     /> 
     
    255284         <field id="idive"        long_name="divergence"                                                   unit="10-8s-1"  /> 
    256285         <field id="ishear"       long_name="shear"                                                        unit="10-8s-1"  /> 
    257          <field id="icevolu"      long_name="ice volume"                                                   unit="m"      /> 
    258          <field id="snowvol"      long_name="snow volume"                                                  unit="m"      /> 
     286         <field id="icevolu"      long_name="ice volume"                                                   unit="m"        /> 
     287         <field id="snowvol"      long_name="snow volume"                                                  unit="m"        /> 
    259288 
    260289         <field id="icetrp"       long_name="ice volume transport"                                         unit="m/day"   /> 
    261290         <field id="snwtrp"       long_name="snw volume transport"                                         unit="m/day"   /> 
    262          <field id="deitrp"       long_name="advected ice enhalpy"                                         unit="W/m2"   /> 
    263          <field id="destrp"       long_name="advected snw enhalpy"                                         unit="W/m2"   /> 
     291         <field id="saltrp"       long_name="salt content transport"                                       unit="psu*kg/m2/day" /> 
     292         <field id="deitrp"       long_name="advected ice enthalpy"                                         unit="W/m2"    /> 
     293         <field id="destrp"       long_name="advected snw enthalpy"                                         unit="W/m2"    /> 
    264294 
    265295         <field id="sfxbri"       long_name="brine salt flux"                                              unit="psu*kg/m2/day" /> 
     
    285315         <field id="vfxspr"       long_name="snw precipitation on ice"                                     unit="m/day"   /> 
    286316 
     317         <field id="afxtot"       long_name="area tendency (total)"                                        unit="day-1"   /> 
     318         <field id="afxdyn"       long_name="area tendency (dynamics)"                                     unit="day-1"   /> 
     319         <field id="afxthd"       long_name="area tendency (thermo)"                                       unit="day-1"   /> 
     320 
    287321         <field id="hfxsum"   long_name="heat fluxes causing surface ice melt"            unit="W/m2"  /> 
    288322         <field id="hfxbom"   long_name="heat fluxes causing bottom ice melt"             unit="W/m2"  /> 
     
    301335         <field id="hfxres"   long_name="heat fluxes from ice-ocean mass exchange during resultant"           unit="W/m2"  /> 
    302336         <field id="hfxsub"   long_name="heat fluxes from ice-atm. mass exchange during sublimation"          unit="W/m2"  /> 
    303          <field id="hfxspr"   long_name="heat fluxes from ice-atm. mass exchange during snow precip"          unit="W/m2" /> 
     337         <field id="hfxspr"   long_name="heat fluxes from ice-atm. mass exchange during snow precip"          unit="W/m2"  /> 
    304338 
    305339         <!-- diags --> 
    306340         <field id="hfxdhc"    long_name="Heat content variation in snow and ice"   unit="W/m2" /> 
    307          <field id="hfxtur"    long_name="turbulent heat flux at the ice base"      unit="W/m2"  /> 
     341         <field id="hfxtur"    long_name="turbulent heat flux at the ice base"      unit="W/m2" /> 
    308342 
    309343      </field_group> 
     
    312346       
    313347      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    314          <field id="utau"         long_name="Wind Stress along i-axis"                    unit="N/m2" /> 
    315          <field id="suoce"        long_name="ocean surface current along i-axis"          unit="m/s"  /> 
    316          <field id="uoce"         long_name="ocean current along i-axis"                  unit="m/s"  grid_ref="grid_U_3D" /> 
    317          <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"      unit="m3/s" grid_ref="grid_U_3D" /> 
     348         <field id="e3u"          long_name="U-cell thickness"                               unit="m"        grid_ref="grid_U_3D" /> 
     349         <field id="utau"         long_name="Wind Stress along i-axis"                       unit="N/m2"                          /> 
     350         <field id="uoce"         long_name="ocean current along i-axis"                     unit="m/s"      grid_ref="grid_U_3D" /> 
     351         <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"               unit="m2/s"     grid_ref="grid_U_3D"  > uoce * e3u </field> 
     352         <field id="ssu"          long_name="ocean surface current along i-axis"             unit="m/s"                           /> 
     353         <field id="sbu"          long_name="ocean bottom current along i-axis"              unit="m/s"                           /> 
     354         <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"         unit="m3/s"     grid_ref="grid_U_3D" /> 
    318355         <field id="uocet"        long_name="ocean transport along i-axis times temperature" unit="degC.m/s" grid_ref="grid_U_3D" /> 
    319356         <field id="uoces"        long_name="ocean transport along i-axis times salinity"    unit="psu.m/s"  grid_ref="grid_U_3D" /> 
    320357         <!-- variables available with MLE --> 
    321          <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"             unit="m3/s"  grid_ref="grid_U_3D" /> 
     358         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"                unit="m3/s"     grid_ref="grid_U_3D" /> 
    322359         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    323          <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
     360         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"                 unit="m/s"      grid_ref="grid_U_3D" /> 
    324361         <!-- uoce_eiv: available with key_trabbl --> 
    325          <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
    326          <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
     362         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"                 unit="m/s"      grid_ref="grid_U_3D" /> 
     363             <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"                unit="m3/s"                          /> 
     364         <!-- variable for ice shelves --> 
     365         <field id="utbl"         long_name="zonal current in the Losh tbl"                  unit="m/s"                           /> 
    327366         <!-- variables available with key_diaar5 --> 
    328          <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"  unit="kg/s" grid_ref="grid_U_3D" /> 
    329          <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"  unit="W"    /> 
    330          <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    331          <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"     unit="W"    /> 
    332          <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" unit="W"    /> 
     367         <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"     unit="kg/s"    grid_ref="grid_U_3D" /> 
     368         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"     unit="W"                             /> 
     369         <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"     unit="PSU*kg/s"                      /> 
     370         <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"        unit="W"                             /> 
     371         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis"    unit="W"                             /> 
    333372     </field_group> 
    334373       
     
    336375       
    337376      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    338          <field id="vtau"         long_name="Wind Stress along j-axis"                    unit="N/m2" /> 
    339          <field id="svoce"        long_name="ocean surface current along j-axis"          unit="m/s"  /> 
    340          <field id="voce"         long_name="ocean current along j-axis"                  unit="m/s"  grid_ref="grid_V_3D" /> 
    341          <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"      unit="m3/s" grid_ref="grid_V_3D" /> 
     377         <field id="e3v"          long_name="V-cell thickness"                               unit="m"        grid_ref="grid_V_3D" /> 
     378         <field id="vtau"         long_name="Wind Stress along j-axis"                       unit="N/m2"                          /> 
     379         <field id="voce"         long_name="ocean current along j-axis"                     unit="m/s"      grid_ref="grid_V_3D" /> 
     380         <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"               unit="m2/s"     grid_ref="grid_V_3D"  > voce * e3v </field> 
     381         <field id="ssv"          long_name="ocean surface current along j-axis"             unit="m/s"                           /> 
     382         <field id="sbv"          long_name="ocean bottom current along j-axis"              unit="m/s"                           /> 
     383         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"         unit="m3/s"     grid_ref="grid_V_3D" /> 
    342384         <field id="vocet"        long_name="ocean transport along j-axis times temperature" unit="degC.m/s" grid_ref="grid_V_3D" /> 
    343385         <field id="voces"        long_name="ocean transport along j-axis times salinity"    unit="psu.m/s"  grid_ref="grid_V_3D" /> 
    344386         <!-- variables available with MLE --> 
    345          <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_V_3D" /> 
     387         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"                unit="m3/s"     grid_ref="grid_V_3D" /> 
    346388         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    347          <field id="voce_eiv"     long_name="EIV ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
     389         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"                 unit="m/s"      grid_ref="grid_V_3D" /> 
    348390         <!-- voce_eiv: available with key_trabbl --> 
    349          <field id="voce_bbl"     long_name="BBL ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
    350          <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"             unit="m3/s"                   /> 
     391         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"                 unit="m/s"      grid_ref="grid_V_3D" /> 
     392             <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"                unit="m3/s"                          /> 
     393         <!-- variable for ice shelves --> 
     394             <field id="vtbl"         long_name="meridional current in the Losh tbl"             unit="m/s"                           /> 
    351395         <!-- variables available with key_diaar5 --> 
    352          <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"  unit="kg/s" grid_ref="grid_V_3D" /> 
    353          <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"  unit="W"     /> 
    354          <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    355          <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"     unit="W"     /> 
    356          <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" unit="W"     /> 
     396         <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"     unit="kg/s"    grid_ref="grid_V_3D" /> 
     397         <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"     unit="W"                             /> 
     398         <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"     unit="PSU*kg/s"                      /> 
     399         <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"        unit="W"                             /> 
     400         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis"    unit="W"                             /> 
    357401      </field_group> 
    358402       
     
    360404       
    361405      <field_group id="grid_W" grid_ref="grid_W_3D"> 
     406        <field id="e3w"          long_name="W-cell thickness"                            unit="m"   /> 
    362407        <field id="woce"         long_name="ocean vertical velocity"                  unit="m/s"    /> 
     408        <field id="woce_e3w"     long_name="ocean vertical velocity * e3v"            unit="m2/s"    > woce * e3w </field> 
    363409        <field id="wocetr_eff"   long_name="effective ocean vertical transport"       unit="m3/s"   /> 
    364410        <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     
    386432 
    387433      <field_group id="scalar"  domain_ref="1point" > 
    388      <field id="voltot"     long_name="global mean volume"                         unit="m3"   /> 
    389      <field id="sshtot"     long_name="global mean ssh"                            unit="m"    /> 
    390      <field id="sshsteric"  long_name="global mean ssh steric"                     unit="m"    /> 
    391      <field id="sshthster"  long_name="global mean ssh thermosteric"               unit="m"    /> 
    392      <field id="masstot"    long_name="global mean mass"                           unit="kg"   /> 
    393      <field id="temptot"    long_name="global mean temperature"                    unit="degC" /> 
    394      <field id="saltot"     long_name="global mean salinity"                       unit="psu"  /> 
    395      <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait" unit="kg/s" /> 
     434         <field id="voltot"     long_name="global mean volume"                         unit="m3"   /> 
     435         <field id="sshtot"     long_name="global mean ssh"                            unit="m"    /> 
     436         <field id="sshsteric"  long_name="global mean ssh steric"                     unit="m"    /> 
     437         <field id="sshthster"  long_name="global mean ssh thermosteric"               unit="m"    /> 
     438         <field id="masstot"    long_name="global mean mass"                           unit="kg"   /> 
     439         <field id="temptot"    long_name="global mean temperature"                    unit="degC" /> 
     440         <field id="saltot"     long_name="global mean salinity"                       unit="psu"  /> 
     441         <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait" unit="kg/s" /> 
     442 
    396443         <!-- available with ln_diahsb --> 
    397          <field id="bgtemper" long_name="global mean temperature variation"            unit="degC"/> 
    398          <field id="bgsaline" long_name="global mean salinity variation"               unit="psu"/> 
    399          <field id="bgheatco" long_name="global mean heat content variation"           unit="10^20J"/> 
    400          <field id="bgsaltco" long_name="global mean salt content variation"           unit="psu*km3" /> 
    401          <field id="bgvolssh" long_name="global mean volume variation (ssh)"           unit="km3"/> 
    402          <field id="bgvole3t" long_name="global mean volume variation (e3t)"           unit="km3"/> 
    403          <field id="bgfrcvol" long_name="global mean volume variation from forcing"    unit="km3"/> 
    404          <field id="bgfrctem" long_name="global mean forcing from heat content variation"   unit="degC"/> 
    405          <field id="bgfrcsal" long_name="global mean forcing salt content variation"        unit="psu"/> 
    406          <field id="bgmistem" long_name="global mean temperature error due to free surface" unit="degC"/> 
    407          <field id="bgmissal" long_name="global mean salinity error due to free surface"    unit="psu"/> 
     444                 <field id="bgtemper"     long_name="global mean temperature"                  unit="degC"   /> 
     445                 <field id="bgsaline"     long_name="global mean salinity"                     unit="psu"    /> 
     446                 <field id="bgheatco"     long_name="global mean heat content"                 unit="10^9J"  /> 
     447                 <field id="bgsaltco"     long_name="global mean salt content"                 unit="psu*m3" /> 
     448                 <field id="bgvolssh"     long_name="global mean ssh volume"                   unit="km3"     /> 
     449             <field id="bgvole3t"     long_name="global mean volume variation (e3t)"           unit="km3"/> 
     450                 <field id="bgvoltot"     long_name="global mean volume"                       unit="km3"     /> 
     451                 <field id="bgsshtot"     long_name="global mean ssh"                          unit="m"      /> 
     452                 <field id="bgfrcvol"     long_name="global mean volume from forcing"          unit="km3"     /> 
     453                 <field id="bgfrctem"     long_name="global mean heat content from forcing"    unit="10^9J"  /> 
     454                 <field id="bgfrcsal"     long_name="global mean salt content from forcing"    unit="psu*km3" /> 
     455                 <field id="bgmistem"     long_name="global mean temperature error due to free surface" unit="degC" /> 
     456                 <field id="bgmissal"     long_name="global mean salinity error due to free surface"    unit="psu"  /> 
    408457      </field_group> 
    409458 
    410459      <field_group id="SBC_scalar"  domain_ref="1point" > 
    411460         <!-- available with ln_limdiaout --> 
    412          <field id="ibgvoltot"    long_name="global mean ice volume"                   unit="km3"   /> 
    413          <field id="sbgvoltot"    long_name="global mean snow volume"                  unit="km3"   /> 
    414          <field id="ibgarea"      long_name="global mean ice area"                     unit="km2"   /> 
    415          <field id="ibgsaline"    long_name="global mean ice salinity"                 unit="psu"   /> 
    416          <field id="ibgtemper"    long_name="global mean ice temperature"              unit="degC"   /> 
     461         <field id="ibgvoltot"    long_name="global mean ice volume"                   unit="km3"      /> 
     462         <field id="sbgvoltot"    long_name="global mean snow volume"                  unit="km3"      /> 
     463         <field id="ibgarea"      long_name="global mean ice area"                     unit="km2"      /> 
     464         <field id="ibgsaline"    long_name="global mean ice salinity"                 unit="psu"      /> 
     465         <field id="ibgtemper"    long_name="global mean ice temperature"              unit="degC"     /> 
    417466         <field id="ibgheatco"    long_name="global mean ice heat content"             unit="10^20J"   /> 
    418467         <field id="sbgheatco"    long_name="global mean snow heat content"            unit="10^20J"   /> 
    419          <field id="ibgsaltco"    long_name="global mean ice salt content"             unit="psu*km3"   /> 
     468         <field id="ibgsaltco"    long_name="global mean ice salt content"             unit="psu*km3"  /> 
    420469 
    421470         <field id="ibgvfx"       long_name="global mean volume flux (emp)"               unit="m/day"   /> 
     
    446495 
    447496        <field id="ibghfxthd"   long_name="heat fluxes from ice-ocean exchange during thermo"              unit="W"  /> 
    448         <field id="ibghfxsum"   long_name="heat fluxes causing surface ice melt"              unit="W"  /> 
    449         <field id="ibghfxbom"   long_name="heat fluxes causing bottom ice melt"              unit="W"  /> 
    450         <field id="ibghfxbog"   long_name="heat fluxes causing bottom ice growth"              unit="W"  /> 
    451         <field id="ibghfxdif"   long_name="heat fluxes causing ice temperature change"              unit="W"  /> 
    452         <field id="ibghfxopw"   long_name="heat fluxes causing open water ice formation"              unit="W"  /> 
     497        <field id="ibghfxsum"   long_name="heat fluxes causing surface ice melt"                           unit="W"  /> 
     498        <field id="ibghfxbom"   long_name="heat fluxes causing bottom ice melt"                            unit="W"  /> 
     499        <field id="ibghfxbog"   long_name="heat fluxes causing bottom ice growth"                          unit="W"  /> 
     500        <field id="ibghfxdif"   long_name="heat fluxes causing ice temperature change"                     unit="W"  /> 
     501        <field id="ibghfxopw"   long_name="heat fluxes causing open water ice formation"                   unit="W"  /> 
    453502        <field id="ibghfxdyn"   long_name="heat fluxes from ice-ocean exchange during dynamic"             unit="W"  /> 
    454503        <field id="ibghfxres"   long_name="heat fluxes from ice-ocean exchange during resultant"           unit="W"  /> 
     
    459508 
    460509         <field id="ibgfrcvol"    long_name="global mean forcing volume (emp)"         unit="km3"      /> 
    461          <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"         unit="psu*km3"   /> 
     510         <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"         unit="psu*km3"  /> 
    462511         <field id="ibgvolgrm"    long_name="global mean ice growth+melt volume"       unit="km3"      /> 
    463512      </field_group> 
     
    491540      </field_group> 
    492541 
    493       <!-- PTR --> 
    494       <field_group id="diaptr" domain_ref="ptr"   grid_ref="grid_T_3D"  enabled=".FALSE."> <!-- time step automaticaly defined based on nn_fptr --> 
    495         <field id="zotemglo"          long_name="Zonal Mean Temperature"                       unit="C"   /> 
    496         <field id="zosalglo"          long_name="Zonal Mean Salinity"                          unit="PSU"   /> 
    497         <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"   /> 
    498         <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"   /> 
    499         <field id="zomsfeiv"          long_name="Bolus Meridional Stream-Function: Global"     unit="Sv"   /> 
    500         <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="C"   /> 
    501         <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="PSU"   /> 
    502         <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"   /> 
    503         <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"   /> 
    504         <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="C"   /> 
    505         <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="PSU"   /> 
    506         <field id="zosrfpac"          long_name="Zonal Mean Surface : Pacific"                 unit="m2"   /> 
    507         <field id="zomsfpac"          long_name="Meridional Stream-Function: Pacific"          unit="Sv"   /> 
    508         <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="C"   /> 
    509         <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="PSU"   /> 
    510         <field id="zosrfind"          long_name="Zonal Mean Surface : Indian"                  unit="m2"   /> 
    511         <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"   /> 
    512         <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="C"   /> 
    513         <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="PSU"   /> 
    514         <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"   /> 
    515         <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"    /> 
    516         <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"       grid_ref="grid_T_2D"   /> 
    517         <field id="sophteiv"          long_name="Bolus Advective Heat Transport"               unit="PW"       grid_ref="grid_T_2D"   /> 
    518         <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"       grid_ref="grid_T_2D"   /> 
    519         <field id="sophtove"          long_name="Overtuning Heat Transport"                    unit="PW"       grid_ref="grid_T_2D"   /> 
    520         <field id="sohtatl"           long_name="Heat Transport Atlantic"                      unit="PW"       grid_ref="grid_T_2D"   /> 
    521         <field id="sohtpac"           long_name="Heat Transport Pacific"                       unit="PW"       grid_ref="grid_T_2D"   /> 
    522         <field id="sohtind"           long_name="Heat Transport Indian"                        unit="PW"       grid_ref="grid_T_2D"   /> 
    523         <field id="sohtipc"           long_name="Heat Transport Pacific+Indian"                unit="PW"       grid_ref="grid_T_2D"   /> 
    524         <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    525         <field id="sopsteiv"          long_name="Bolus Advective Salt Transport"               unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    526         <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    527         <field id="sopstove"          long_name="Overtuning Salt Transport"                    unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    528         <field id="sostatl"           long_name="Salt Transport Atlantic"                      unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    529         <field id="sostpac"           long_name="Salt Transport Pacific"                       unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    530         <field id="sostind"           long_name="Salt Transport Indian"                        unit="Giga g/s" grid_ref="grid_T_2D"   /> 
    531         <field id="sostipc"           long_name="Salt Transport Pacific+Indian"                unit="Giga g/s" grid_ref="grid_T_2D"   /> 
     542      <!-- Poleward transport : ptr -->      
     543      <field_group id="diaptr" domain_ref="ptr"  >  
     544        <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"       grid_ref="grid_W_3D"  /> 
     545        <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"       grid_ref="grid_W_3D"  /> 
     546        <field id="zomsfpac"          long_name="Meridional Stream-Function: Pacific"          unit="Sv"       grid_ref="grid_W_3D"  /> 
     547        <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"       grid_ref="grid_W_3D"  /> 
     548        <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"       grid_ref="grid_W_3D"   /> 
     549        <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="C"        grid_ref="grid_T_3D" /> 
     550        <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="C"        grid_ref="grid_T_3D" /> 
     551        <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="C"        grid_ref="grid_T_3D" /> 
     552        <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="C"        grid_ref="grid_T_3D" /> 
     553        <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="C"        grid_ref="grid_T_3D" /> 
     554        <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="PSU"      grid_ref="grid_T_3D"   /> 
     555        <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="PSU"      grid_ref="grid_T_3D"   /> 
     556        <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="PSU"      grid_ref="grid_T_3D"   /> 
     557        <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="PSU"      grid_ref="grid_T_3D"   /> 
     558        <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="PSU"      grid_ref="grid_T_3D"   /> 
     559        <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"       grid_ref="grid_T_3D"  /> 
     560        <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"       grid_ref="grid_T_3D"  /> 
     561        <field id="zosrfpac"          long_name="Zonal Mean Surface : Pacific"                 unit="m2"       grid_ref="grid_T_3D"  /> 
     562        <field id="zosrfind"          long_name="Zonal Mean Surface : Indian"                  unit="m2"       grid_ref="grid_T_3D"  /> 
     563        <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"       grid_ref="grid_T_3D"  /> 
     564        <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"       grid_ref="grid_T_2D" /> 
     565        <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"       grid_ref="grid_T_2D" /> 
     566        <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s" grid_ref="grid_T_2D" /> 
     567        <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s" grid_ref="grid_T_2D" /> 
    532568      </field_group> 
     569 
    533570 
    534571      <!-- ptrc on T grid --> 
     
    620657       <field id="EPSI100"     long_name="Export of Silicate at 100 m"             unit="mol/m2/s"                     /> 
    621658       <field id="EPCAL100"    long_name="Export of Calcite at 100 m"              unit="mol/m2/s"                     /> 
     659       <field id="EXPC"        long_name="Export of carbon"                        unit="molC/m2/s"  grid_ref="grid_T_3D" /> 
     660       <field id="EXPFE"       long_name="Export of biogenic iron"                 unit="molFe/m2/s" grid_ref="grid_T_3D" /> 
     661       <field id="EXPSI"       long_name="Export of Silicate"                      unit="molSi/m2/s" grid_ref="grid_T_3D" /> 
     662       <field id="EXPCAL"      long_name="Export of Calcite"                       unit="molC/m2/s"  grid_ref="grid_T_3D" /> 
    622663       <field id="Cflx"        long_name="DIC flux"                                unit="mol/m2/s"                     /> 
    623664       <field id="Oflx"        long_name="Oxygen flux"                             unit="mol/m2/s"                     /> 
     
    629670       <field id="Ironsed"     long_name="Iron deposition from sediment"           unit="mol/m2/s"  grid_ref="grid_T_3D"/> 
    630671 
     672 
    631673       <!-- PISCES with Kriest parametisation : variables available with key_kriest --> 
    632        <field id="POCFlx"      long_name="Particulate organic C flux"              unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    633        <field id="NumFlx"      long_name="Particle number flux"                    unit="nbr/m2/s"   grid_ref="grid_T_3D" /> 
    634        <field id="SiFlx"       long_name="Biogenic Si flux"                        unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    635        <field id="CaCO3Flx"    long_name="CaCO3 flux"                              unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    636        <field id="xnum"        long_name="Number of particles in aggregats"        unit="-"          grid_ref="grid_T_3D" /> 
    637        <field id="W1"          long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" /> 
    638        <field id="W2"          long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" /> 
     674       <field id="EPN100"      long_name="Particulate number flux at 100 m"        unit="mol/m2/s"                     /> 
     675       <field id="EXPN"        long_name="Particulate number flux"                 unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
     676       <field id="XNUM"        long_name="Number of particles in aggregats"        unit="-"          grid_ref="grid_T_3D" /> 
     677       <field id="WSC"         long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" /> 
     678       <field id="WSN"         long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" /> 
     679 
     680      <!-- dbio_T on T grid : variables available with key_diaar5 --> 
     681       <field id="TPP"         long_name="Total Primary production of phyto"                   unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     682       <field id="TPNEW"       long_name="New Primary production of phyto"                     unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     683       <field id="TPBFE"       long_name="Total biogenic iron production"                      unit="molFe/m3/s" grid_ref="grid_T_3D" /> 
     684       <field id="INTDIC"      long_name="DIC content"                                         unit="kg/m2"   /> 
     685       <field id="O2MIN"       long_name="Oxygen minimum concentration"                        unit="molC/L" /> 
     686       <field id="ZO2MIN"      long_name="Depth of oxygen minimum concentration"               unit="m"       /> 
     687       <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="molN/m2/s"  /> 
     688       <field id="INTPPPHY"    long_name="Vertically integrated primary production by nanophy" unit="molC/m2/s"  /> 
     689       <field id="INTPPPHY2"   long_name="Vertically integrated primary production by diatom"  unit="molC/m2/s"  /> 
     690       <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="molC/m2/s"  /> 
     691       <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="molC/m2/s"  /> 
     692       <field id="INTPBFE"     long_name="Vertically integrated of biogenic iron production"   unit="molFe/m2/s" /> 
     693       <field id="INTPBSI"     long_name="Vertically integrated of biogenic Si production"     unit="molSi/m2/s" /> 
     694       <field id="INTPCAL"     long_name="Vertically integrated of calcite production"         unit="molC/m2/s" /> 
    639695 
    640696       <!-- PISCES light : variables available with key_pisces_reduced --> 
     
    668724       <field id="qintC14b"    long_name="Cumulative air-sea flux of Bomb C14"      unit="mol/m2"     /> 
    669725       <field id="fdecay"      long_name="Radiactive decay of Bomb C14"             unit="mol/m3"  grid_ref="grid_T_3D"  /> 
     726     </field_group> 
     727 
     728    <field_group id="PISCES_scalar"  domain_ref="1point" > 
     729       <field id="pno3tot"         long_name="global mean nitrate concentration"                  unit="mol/l"   /> 
     730       <field id="psiltot"         long_name="global mean silicate concentration"                 unit="mol/l"   /> 
     731       <field id="palktot"         long_name="global mean alkalinity concentration"               unit="mol/l"   /> 
     732       <field id="pfertot"         long_name="global mean iron concentration"                     unit="mol/l"   /> 
     733       <field id="tcflx"           long_name="total Flux of Carbon out of the ocean"              unit="molC/s" /> 
     734       <field id="tcflxcum"        long_name="cumulative total Flux of Carbon out of the ocean"   unit="molC/s" /> 
     735       <field id="tcexp"           long_name="total Carbon export at 100m"                        unit="molC/s" /> 
     736       <field id="tintpp"          long_name="global total integrated primary production"         unit="molC/s" /> 
     737       <field id="tnfix"           long_name="global total nitrogen fixation"                     unit="molN/s" /> 
     738       <field id="tdenit"          long_name="Total denitrification"                              unit="molN/s" /> 
    670739     </field_group> 
    671740 
     
    734803    <field_group id="groupU" > 
    735804      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    736       <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     805      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    737806      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    738807    </field_group> 
     
    740809    <field_group id="groupV" > 
    741810      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    742       <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     811      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    743812      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    744813    </field_group> 
     
    748817    </field_group> 
    749818     
     819    <!--  
     820============================================================================================================ 
     821                     Trend diagnostics : temperature, KE, PE, momentum 
     822============================================================================================================ 
     823    --> 
     824 
     825    <field_group id="trendT" grid_ref="grid_T_3D"> 
     826      <!-- variables available with ln_tra_trd --> 
     827      <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                 unit="degC/s" /> 
     828      <field id="strd_xad"      long_name="salinity   -trend: i-advection"                 unit="psu/s"  /> 
     829      <field id="ttrd_yad"      long_name="temperature-trend: j-advection"                 unit="degC/s" /> 
     830      <field id="strd_yad"      long_name="salinity   -trend: j-advection"                 unit="psu/s"  /> 
     831      <field id="ttrd_zad"      long_name="temperature-trend: k-advection"                 unit="degC/s" /> 
     832      <field id="strd_zad"      long_name="salinity   -trend: k-advection"                 unit="psu/s"  /> 
     833      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (no-vvl)"       unit="degC/s" grid_ref="grid_T_2D" /> 
     834      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (no-vvl)"       unit="psu/s"  grid_ref="grid_T_2D" /> 
     835      <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"          unit="degC/s" /> 
     836      <field id="strd_ldf"      long_name="salinity   -trend: lateral  diffusion"          unit="psu/s"  /> 
     837      <field id="ttrd_zdf"      long_name="temperature-trend: vertical diffusion"          unit="degC/s" /> 
     838      <field id="strd_zdf"      long_name="salinity   -trend: vertical diffusion"          unit="psu/s"  /> 
     839      <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 
     840      <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"        unit="degC/s" /> 
     841      <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"        unit="psu/s"  /> 
     842      <!-- --> 
     843      <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"          unit="degC/s" /> 
     844      <field id="strd_dmp"      long_name="salinity   -trend: interior restoring"          unit="psu/s"  /> 
     845      <field id="ttrd_bbl"      long_name="temperature-trend: bottom boundary layer"       unit="degC/s" /> 
     846      <field id="strd_bbl"      long_name="salinity   -trend: bottom boundary layer"       unit="psu/s"  /> 
     847      <field id="ttrd_npc"      long_name="temperature-trend: non-penetrative conv."       unit="degC/s" /> 
     848      <field id="strd_npc"      long_name="salinity   -trend: non-penetrative conv."       unit="psu/s"  /> 
     849      <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"     unit="degC/s" /> 
     850      <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"     unit="degC/s" /> 
     851      <field id="ttrd_qsr"      long_name="temperature-trend: solar penetr. heating"       unit="degC/s" /> 
     852      <field id="ttrd_bbc"      long_name="temperature-trend: geothermal heating"          unit="degC/s" /> 
     853      <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"         unit="degC/s" /> 
     854      <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"         unit="psu/s"  /> 
     855      <!-- variables available with ln_KE_trd --> 
     856      <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"        unit="W/s^3"  /> 
     857      <field id="ketrd_spg"     long_name="ke-trend: surface     pressure gradient"        unit="W/s^3"  /> 
     858      <field id="ketrd_spgexp"  long_name="ke-trend: surface pressure gradient (explicit)" unit="W/s^3"  /> 
     859      <field id="ketrd_spgflt"  long_name="ke-trend: surface pressure gradient (filter)"   unit="W/s^3"  /> 
     860      <field id="ssh_flt"       long_name="filtered contribution to ssh (dynspg_flt)"      unit="m"     grid_ref="grid_T_2D"   /> 
     861      <field id="w0"            long_name="surface vertical velocity"                      unit="m/s"   grid_ref="grid_T_2D"   /> 
     862      <field id="pw0_exp"       long_name="surface pressure flux due to ssh"               unit="W/s^2" grid_ref="grid_T_2D"   /> 
     863      <field id="pw0_flt"       long_name="surface pressure flux due to filtered ssh"      unit="W/s^2" grid_ref="grid_T_2D"   /> 
     864      <field id="ketrd_keg"     long_name="ke-trend: KE gradient         or hor. adv."     unit="W/s^3"  /> 
     865      <field id="ketrd_rvo"     long_name="ke-trend: relative  vorticity or metric term"   unit="W/s^3"  /> 
     866      <field id="ketrd_pvo"     long_name="ke-trend: planetary vorticity"                  unit="W/s^3"  /> 
     867      <field id="ketrd_zad"     long_name="ke-trend: vertical  advection"                  unit="W/s^3"  /> 
     868      <field id="ketrd_udx"     long_name="ke-trend: U.dx[U]"                              unit="W/s^3"  /> 
     869      <field id="ketrd_ldf"     long_name="ke-trend: lateral   diffusion"                  unit="W/s^3"  /> 
     870      <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                  unit="W/s^3"  /> 
     871      <field id="ketrd_tau"     long_name="ke-trend: wind stress "                         unit="W/s^3" grid_ref="grid_T_2D"   /> 
     872      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"           unit="W/s^3"  />    
     873      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"           unit="W/s^3"  />    
     874      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"            unit="W/s^3"  />   
     875      <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"    unit="W/s^3"  /> 
     876      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                  unit="W/s^2"  />    
     877      <!-- variables available with ln_PE_trd --> 
     878      <field id="petrd_xad"     long_name="pe-trend: i-advection"                          unit="W/m^3"  /> 
     879      <field id="petrd_yad"     long_name="pe-trend: j-advection"                          unit="W/m^3"  /> 
     880      <field id="petrd_zad"     long_name="pe-trend: k-advection"                          unit="W/m^3"  /> 
     881      <field id="petrd_sad"     long_name="pe-trend: surface adv. (no-vvl)"                unit="W/m^3"  grid_ref="grid_T_2D" /> 
     882      <field id="petrd_ldf"     long_name="pe-trend: lateral  diffusion"                   unit="W/m^3"  /> 
     883      <field id="petrd_zdf"     long_name="pe-trend: vertical diffusion"                   unit="W/m^3"  /> 
     884      <field id="petrd_zdfp"    long_name="pe-trend: pure vert. diffusion"                 unit="W/m^3"  /> 
     885      <field id="petrd_dmp"     long_name="pe-trend: interior restoring"                   unit="W/m^3"  /> 
     886      <field id="petrd_bbl"     long_name="pe-trend: bottom boundary layer"                unit="W/m^3"  /> 
     887      <field id="petrd_npc"     long_name="pe-trend: non-penetrative conv."                unit="W/m^3"  /> 
     888      <field id="petrd_nsr"     long_name="pe-trend: surface forcing + runoff"             unit="W/m^3"  /> 
     889      <field id="petrd_qsr"     long_name="pe-trend: solar penetr. heating"                unit="W/m^3"  /> 
     890      <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"                   unit="W/m^3"  /> 
     891      <field id="petrd_atf"     long_name="pe-trend: asselin time filter"                  unit="W/m^3"  /> 
     892      <field id="PEanom"        long_name="potential energy anomaly"                       unit="SI"     />    
     893      <field id="alphaPE"       long_name="- partial deriv. of PEanom wrt T"               unit="/degC"  />    
     894      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"                 unit="/psu"   />    
     895    </field_group> 
     896 
     897    <field_group id="trendU" grid_ref="grid_U_3D"> 
     898     <!-- variables available with ln_dyn_trd --> 
     899     <field id="utrd_hpg"       long_name="i-trend: hydrostatic pressure gradient"         unit="m/s^2"                      /> 
     900     <field id="utrd_spg"       long_name="i-trend: surface     pressure gradient"         unit="m/s^2"                      /> 
     901     <field id="utrd_spgexp"    long_name="i-trend: surface pressure gradient (explicit)"  unit="m/s^2"                      /> 
     902     <field id="utrd_spgflt"    long_name="i-trend: surface pressure gradient (filtered)"  unit="m/s^2"                      /> 
     903     <field id="utrd_keg"       long_name="i-trend: KE gradient         or hor. adv."      unit="m/s^2"                      /> 
     904     <field id="utrd_rvo"       long_name="i-trend: relative  vorticity or metric term"    unit="m/s^2"                      /> 
     905     <field id="utrd_pvo"       long_name="i-trend: planetary vorticity"                   unit="m/s^2"                      /> 
     906     <field id="utrd_zad"       long_name="i-trend: vertical  advection"                   unit="m/s^2"                      /> 
     907     <field id="utrd_udx"       long_name="i-trend: U.dx[U]"                               unit="m/s^2"                      /> 
     908     <field id="utrd_ldf"       long_name="i-trend: lateral   diffusion"                   unit="m/s^2"                      /> 
     909     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                   unit="m/s^2"                      /> 
     910     <field id="utrd_tau"       long_name="i-trend: wind stress "                          unit="m/s^2" grid_ref="grid_U_2D" /> 
     911     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"            unit="m/s^2"                      />    
     912     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"            unit="m/s^2"                      />    
     913     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"       unit="m/s^2"                      />    
     914     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"             unit="m/s^2"                      />    
     915    </field_group> 
     916 
     917    <field_group id="trendV" grid_ref="grid_V_3D"> 
     918     <!-- variables available with ln_dyn_trd --> 
     919     <field id="vtrd_hpg"       long_name="j-trend: hydrostatic pressure gradient"         unit="m/s^2"                      /> 
     920     <field id="vtrd_spg"       long_name="j-trend: surface     pressure gradient"         unit="m/s^2"                      /> 
     921     <field id="vtrd_spgexp"    long_name="j-trend: surface pressure gradient (explicit)"  unit="m/s^2"                      /> 
     922     <field id="vtrd_spgflt"    long_name="j-trend: surface pressure gradient (filtered)"  unit="m/s^2"                      /> 
     923     <field id="vtrd_keg"       long_name="j-trend: KE gradient         or hor. adv."      unit="m/s^2"                      /> 
     924     <field id="vtrd_rvo"       long_name="j-trend: relative  vorticity or metric term"    unit="m/s^2"                      /> 
     925     <field id="vtrd_pvo"       long_name="j-trend: planetary vorticity"                   unit="m/s^2"                      /> 
     926     <field id="vtrd_zad"       long_name="j-trend: vertical  advection"                   unit="m/s^2"                      /> 
     927     <field id="vtrd_vdy"       long_name="i-trend: V.dx[V]"                               unit="m/s^2"                      /> 
     928     <field id="vtrd_ldf"       long_name="j-trend: lateral   diffusion"                   unit="m/s^2"                      /> 
     929     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                   unit="m/s^2"                      /> 
     930     <field id="vtrd_tau"       long_name="j-trend: wind stress "                          unit="m/s^2" grid_ref="grid_V_2D" /> 
     931     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"            unit="m/s^2"                      />    
     932     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"            unit="m/s^2"                      />    
     933     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"       unit="m/s^2"                      />    
     934     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"            unit="m/s^2"                       />    
     935    </field_group> 
     936 
    750937    </field_definition> 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2_cfg

    r2295 r2451  
    8888                           !     (also = the frequency of sea-ice model call) 
    8989   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
    90    ln_cpl      = .true.    !  Coupled formulation                       (T => fill namsbc_cpl ) 
    9190   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    9291   nn_ice      = 2         !  =0 no ice boundary condition   , 
    9392                           !  =1 use observed ice-cover      , 
    9493                           !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    95    cn_iceflx   = 'linear'  !  AUTO -  
    96                            !  ="ave"    - Average flux over all ice categories ("key_lim3") 
    97                            !  ="linear" - Linear repartition of fluxes over ice categories ("key_lim3") 
    9894   nn_ice_embd = 0         !  AUTO -  
    9995                           !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
     
    279275!----------------------------------------------------------------------- 
    280276   ln_diaptr  = .true.     !  Poleward heat and salt transport (T) or not (F) 
    281    ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    282277   ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    283278                           !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    284    ln_ptrcomp = .true.     !  Add decomposition : overturning 
    285    nn_fptr    =  1         !  Frequency of ptr computation [time step] 
    286    nn_fwri    =  15        !  AUTO - Frequency of ptr outputs [time step] 
    287279/ 
    288280!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2_ref

    r2372 r2451  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun, namcfg) 
    3 !! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
     3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
     
    1010!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    1111!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
    12 !!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namptr, namhsb) 
     12!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb) 
    1313!!             10 - miscellaneous    (namsol, nammpp, namctl) 
    1414!!             11 - Obs & Assim      (namobs, nam_asminc) 
     
    1818!!                   ***  Run management namelists  *** 
    1919!!====================================================================== 
    20 !!   namrun        parameters of the run 
     20!!   namrun       parameters of the run 
    2121!!====================================================================== 
    2222! 
     
    3131   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3232   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    33    nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=.true. 
    34    nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
     33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
     34   nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    3535                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    3636                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     
    4747/ 
    4848! 
    49 !----------------------------------------------------------------------- 
    50 &namcfg     !   default parameters of the configuration       
     49!!====================================================================== 
     50!!                      ***  Domain namelists  *** 
     51!!====================================================================== 
     52!!   namcfg       parameters of the configuration       
     53!!   namzgr       vertical coordinate 
     54!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     55!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     56!!   namtsd       data: temperature & salinity 
     57!!====================================================================== 
     58! 
     59!----------------------------------------------------------------------- 
     60&namcfg     !   parameters of the configuration       
    5161!----------------------------------------------------------------------- 
    5262   cp_cfg      =  "default"            !  name of the configuration 
    53    cp_cfz      =         ''            !  name of the zoom of configuration 
     63   cp_cfz      =  "no zoom"            !  name of the zoom of configuration 
    5464   jp_cfg      =       0               !  resolution of the configuration 
    5565   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
     
    5767   jpkdta      =      31               !  number of levels      ( >= jpk ) 
    5868   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
    59    jpjglo      =      12               !  2nd    -                  -    --> j  =jpjdta 
     69   jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta 
    6070   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    6171   jpjzoom     =       1               !  in data domain indices 
     
    6676                                       !  = 5 North fold F-point pivot 
    6777                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    68 / 
    69 !!====================================================================== 
    70 !!                      ***  Domain namelists  *** 
    71 !!====================================================================== 
    72 !!   namzgr       vertical coordinate 
    73 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    74 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    75 !!   namtsd       data: temperature & salinity 
    76 !!====================================================================== 
    77 ! 
     78   ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present  
     79                                       !  in netcdf input files, as the start j-row for reading 
     80/ 
    7881!----------------------------------------------------------------------- 
    7982&namzgr        !   vertical coordinate 
     
    8285   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    8386   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     87   ln_isfcav   = .false.   !  ice shelf cavity               (T/F) 
    8488/ 
    8589!----------------------------------------------------------------------- 
     
    215219!!   namtra_qsr      penetrative solar radiation 
    216220!!   namsbc_rnf      river runoffs 
     221!!   namsbc_isf      ice shelf melting/freezing 
    217222!!   namsbc_apr      Atmospheric Pressure 
    218223!!   namsbc_ssr      sea surface restoring term (for T and/or S) 
     
    230235   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    231236   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    232    ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
    233237   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    234238   nn_ice      = 2         !  =0 no ice boundary condition   , 
    235239                           !  =1 use observed ice-cover      , 
    236                            !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     240                           !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    237241   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    238242                           !  =1 levitating ice with mass and salt exchange but no presure effect 
    239243                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    240244   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    241    ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     245   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
     246   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
     247                           !  0 =no isf                  1 = presence of ISF  
     248                           !  2 = bg03 parametrisation   3 = rnf file for isf    
     249                           !  4 = ISF fwf specified 
     250                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    242251   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    243252   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     
    249258   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    250259                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    251    cn_iceflx = 'linear'    !  redistribution of solar input into ice categories during coupling ice/atm. 
     260   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
     261                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
     262                           !  = 0  Average per-category fluxes (forced and coupled mode) 
     263                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
     264                           !  = 2  Redistribute a single flux over categories (coupled mode only) 
    252265/ 
    253266!----------------------------------------------------------------------- 
     
    305318 
    306319   cn_dir      = './'      !  root directory for the location of the bulk files 
    307    ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    308320   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    309    ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
    310    rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
    311    rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
     321   rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
     322   rn_zu       = 10.        !  Wind vector reference height (m)                  
    312323   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    313324   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     
    336347!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    337348! send 
    338 sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    339 sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    340 sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
    341 sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    342 sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     349   sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     350   sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     351   sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     352   sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     353   sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    343354! receive 
    344 sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    345 sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    346 sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    347 sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    348 sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    349 sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    350 sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    351 sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    352 sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    353 sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     355   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     356   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     357   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     358   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     359   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     360   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     361   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     362   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     363   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     364   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     365! 
     366   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     367   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models 
     368                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    354369/ 
    355370!----------------------------------------------------------------------- 
     
    405420   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
    406421   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
     422/ 
     423!----------------------------------------------------------------------- 
     424&namsbc_isf    !  Top boundary layer (ISF)  
     425!----------------------------------------------------------------------- 
     426!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     427!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     428! nn_isf == 4 
     429   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     430   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     431! nn_isf == 3 
     432   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     433! nn_isf == 2 and 3 
     434   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     435   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     436! nn_isf == 2 
     437   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     438! for all case 
     439   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 
     440! only for nn_isf = 1 or 2 
     441   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula 
     442   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula 
     443! only for nn_isf = 1 
     444   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 
     445   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008) 
     446                          ! 0 => thickness of the tbl = thickness of the first wet cell 
     447   ln_conserve = .true.   ! conservative case (take into account meltwater advection) 
     448   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s) 
     449                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     450                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 
     451                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999 
    407452/ 
    408453!----------------------------------------------------------------------- 
     
    631676   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    632677   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
    633    ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    634678   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    635679   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     680   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
     681   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
     682   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
     683   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
     684   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
     685   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
     686   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T) 
     687 
    636688   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
     689   ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    637690/ 
    638691!----------------------------------------------------------------------- 
     
    659712!!   nameos        equation of state 
    660713!!   namtra_adv    advection scheme 
     714!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    661715!!   namtra_ldf    lateral diffusion scheme 
    662716!!   namtra_dmp    T & S newtonian damping 
     
    666720&nameos        !   ocean physical parameters 
    667721!----------------------------------------------------------------------- 
    668    nn_eos      =   0       !  type of equation of state and Brunt-Vaisala frequency 
    669                            !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
    670                            !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    671                            !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    672    rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
    673    rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     722   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency 
     723                                 !  =-1, TEOS-10  
     724                                 !  = 0, EOS-80  
     725                                 !  = 1, S-EOS   (simplified eos) 
     726   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     727   !                             ! 
     728   !                     ! S-EOS coefficients : 
     729   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     730   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     731   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     732   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     733   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     734   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     735   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     736   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
    674737/ 
    675738!----------------------------------------------------------------------- 
    676739&namtra_adv    !   advection scheme for tracer 
    677740!----------------------------------------------------------------------- 
    678    ln_traadv_cen2   =  .false.  !  2nd order centered scheme 
    679    ln_traadv_tvd    =  .true.   !  TVD scheme 
    680    ln_traadv_muscl  =  .false.  !  MUSCL scheme 
    681    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    682    ln_traadv_ubs    =  .false.  !  UBS scheme 
    683    ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    684    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     741   ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
     742   ln_traadv_tvd    =  .true.    !  TVD scheme 
     743   ln_traadv_muscl  =  .false.   !  MUSCL scheme 
     744   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
     745   ln_traadv_ubs    =  .false.   !  UBS scheme 
     746   ln_traadv_qck    =  .false.   !  QUICKEST scheme 
     747   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     748   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
    685749/ 
    686750!----------------------------------------------------------------------- 
     
    729793!----------------------------------------------------------------------- 
    730794   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    731    nn_hdmp     =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
    732                            !                   =XX, damping poleward of XX degrees (XX>0) 
    733                            !                      + F(distance-to-coast) + Red and Med Seas 
    734795   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column 
    735796                           !                   =1 no damping in the mixing layer (kz  criteria) 
    736797                           !                   =2 no damping in the mixed  layer (rho crieria) 
    737    rn_surf     =   50.     !  surface time scale of damping   [days] 
    738    rn_bot      =  360.     !  bottom  time scale of damping   [days] 
    739    rn_dep      =  800.     !  depth of transition between rn_surf and rn_bot [meters] 
    740    nn_file     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     798   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this) 
    741799/ 
    742800 
     
    757815   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    758816   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     817   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection 
    759818/ 
    760819!----------------------------------------------------------------------- 
     
    779838   ln_dynvor_mix = .false. !  mixed scheme 
    780839   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     840   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
    781841/ 
    782842!----------------------------------------------------------------------- 
     
    786846   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    787847   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
     848   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf 
    788849   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    789850   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
     
    903964&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    904965!----------------------------------------------------------------------- 
    905    rn_emin       = 1.e-6   !  minimum value of e   [m2/s2] 
     966   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
    906967   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3] 
    907968   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988) 
    908    rn_clim_galp  = 0.53    !  galperin limit 
    909    ln_crban      = .true.  !  Use Craig & Banner (1994) surface wave mixing parametrisation 
     969   rn_clim_galp  = 0.267   !  galperin limit 
    910970   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case 
    911971   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux 
    912972   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    913    nn_tkebc_surf =   1     !  surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
    914    nn_tkebc_bot  =   1     !  bottom tke condition (0/1=Dir/Neum) 
    915    nn_psibc_surf =   1     !  surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
    916    nn_psibc_bot  =   1     !  bottom psi condition (0/1=Dir/Neum) 
    917    nn_stab_func  =   2     !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
    918    nn_clos       =   1     !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
     973   rn_hsro       =  0.02   !  Minimum surface roughness 
     974   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2) 
     975   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2) 
     976   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
     977   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
     978   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
     979   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
    919980/ 
    920981!----------------------------------------------------------------------- 
     
    938999!!                  ***  Miscellaneous namelists  *** 
    9391000!!====================================================================== 
     1001!!   namsol            elliptic solver / island / free surface 
    9401002!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    9411003!!   namctl            Control prints & Benchmark 
    942 !!   namsol            elliptic solver / island / free surface 
     1004!!   namc1d            1D configuration options                         ("key_c1d") 
     1005!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
     1006!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
    9431007!!====================================================================== 
    9441008! 
     
    9991063   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    10001064/ 
     1065 
    10011066!!====================================================================== 
    10021067!!                  ***  Diagnostics namelists  *** 
    10031068!!====================================================================== 
    10041069!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    1005 !!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     1070!!   namtrd       dynamics and/or tracer trends 
     1071!!   namptr       Poleward Transport Diagnostics 
    10061072!!   namflo       float parameters                                      ("key_float") 
    1007 !!   namptr       Poleward Transport Diagnostics 
    10081073!!   namhsb       Heat and salt budgets 
    10091074!!====================================================================== 
     
    10211086/ 
    10221087!----------------------------------------------------------------------- 
    1023 &namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra") 
    1024 !              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor") 
    1025 !----------------------------------------------------------------------- 
    1026    nn_trd      = 365       !  time step frequency dynamics and tracers trends 
    1027    nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    1028    rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    1029    cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
    1030    cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
    1031    ln_trdmld_restart = .false.         !  restart for ML diagnostics 
    1032    ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    1033 / 
     1088&namtrd        !   diagnostics on dynamics and/or tracer trends 
     1089!              !       and/or mixed-layer trends and/or barotropic vorticity 
     1090!----------------------------------------------------------------------- 
     1091   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
     1092   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     1093   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     1094   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     1095   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
     1096   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
     1097   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output 
     1098   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
     1099   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
     1100/ 
     1101!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     1102!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     1103!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
     1104!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
     1105!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
     1106!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
     1107!!gm 
    10341108!----------------------------------------------------------------------- 
    10351109&namflo       !   float parameters                                      ("key_float") 
     
    10501124!----------------------------------------------------------------------- 
    10511125   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    1052    ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    1053    ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    1054                            !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    1055    ln_ptrcomp = .true.     !  Add decomposition : overturning 
    1056    nn_fptr    =  1         !  Frequency of ptr computation [time step] 
    1057    nn_fwri    =  15        !  Frequency of ptr outputs [time step] 
     1126   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    10581127/ 
    10591128!----------------------------------------------------------------------- 
     
    11031172                           !     ln_ssh                  Logical switch for SSH observations 
    11041173 
    1105    ln_sst     = .true.     ! Logical switch for SST observations 
    1106    ln_reysst  = .true.     !     ln_reysst               Logical switch for Reynolds observations 
     1174   ln_sst     = .false.     ! Logical switch for SST observations 
     1175   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    11071176   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
    11081177 
    11091178   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
    11101179                           !     ln_sss                  Logical switch for SSS observations 
    1111                            !     ln_seaice              Logical switch for Sea Ice observations 
     1180   ln_seaice  = .false.    ! Logical switch for Sea Ice observations 
    11121181                           !     ln_vel3d                Logical switch for velocity observations 
    11131182                           !     ln_velavcur             Logical switch for velocity daily av. cur. 
     
    11301199                           !     sstfiles                GHRSST input observation file name 
    11311200   !                       ! sstfbfiles: Feedback SST input observation file name 
    1132    sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    1133                            !     seaicefiles             Sea Ice input observation file name 
     1201   sstfbfiles = 'sst_01.nc' 
     1202                           !     seaicefiles             Sea Ice input observation file names 
     1203   seaicefiles = 'seaice_01.nc'   
    11341204                           !     velavcurfiles           Vel. cur. daily av. input file name 
    11351205                           !     velhvcurfiles           Vel. cur. high freq. input file name 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_lim2_ORCA2_cfg

    r2221 r2451  
    1414&namicedyn     !   ice dynamic 
    1515!----------------------------------------------------------------------- 
    16       creepl      =   1.0e-08 !  creep limit 
    1716/ 
    1817!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_lim2_ORCA2_ref

    r2221 r2451  
    4848   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    4949   etamn       =   0.0e+07 !  minimun value for viscosity 
    50    creepl      =   1.0e-08 !  creep limit 
    51    ecc         =   2.0     !  eccentricity of the elliptical yield curve 
     50   rn_creepl   =   1.0e-08 !  creep limit 
     51   rn_ecc      =   2.0     !  eccentricity of the elliptical yield curve 
    5252   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    53    nevp        =   120     !  number of EVP subcycling iterations 
     53   nn_nevp     =   120     !  number of EVP subcycling iterations 
    5454   telast      =   9600    !  timescale for EVP elastic waves 
    5555   alphaevp    =   1.0     !  coefficient for the solution of EVP int. stresses 
    56    hminrhg     =   0.05     !  ice thickness (m) below which ice velocity equal ocean velocity 
    5756/ 
    5857!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_lim3_ORCA2_cfg

    r2372 r2451  
    2222&namicedyn     !   ice dynamic 
    2323!----------------------------------------------------------------------- 
    24    pstar       =   4.0e+04 !  1st bulk-rheology parameter 
    2524/ 
    2625!----------------------------------------------------------------------- 
     
    3130&namicethd     !   ice thermodynamic 
    3231!----------------------------------------------------------------------- 
    33    parsub      = 0.0       !  switch for snow sublimation or not - O.M. attention : 0 en couple 
    3432/ 
    3533!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_lim3_ORCA2_ref

    r2372 r2451  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/LIM3 :  1 - dynamics/advection/thermo          (namicerun) 
    3 !! namelists    2 - ice intialisation                  (namiceini) 
    4 !!              3 - ice dynamic                        (namicedyn) 
    5 !!              4 - ice advection                      (namicetrp) 
    6 !!              5 - thermodynamic                      (namicethd) 
    7 !!              6 - ice salinity                       (namicesal) 
    8 !!              7 - mechanical redistribution of ice   (namiceitdme) 
    9 !!              8 - ice diagnostics                    (namicedia) 
    10 !!              9 - ice outputs                        (namiceout) 
     2!! LIM3 namelist :   
     3!!              1 - Generic parameters                 (namicerun) 
     4!!              2 - Ice initialization                 (namiceini) 
     5!!              3 - Ice discretization                 (namiceitd) 
     6!!              4 - Ice dynamics and transport         (namicedyn) 
     7!!              5 - Ice thermodynamics                 (namicethd) 
     8!!              6 - Ice salinity                       (namicesal) 
     9!!              7 - Ice mechanical redistribution      (namiceitdme) 
    1110!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    12  
    13 !----------------------------------------------------------------------- 
    14 &namicerun     !   Share parameters for dynamics/advection/thermo 
    15 !----------------------------------------------------------------------- 
     11! 
     12!------------------------------------------------------------------------------ 
     13&namicerun     !   Generic parameters 
     14!------------------------------------------------------------------------------ 
     15   jpl            =    5           !  number of ice  categories 
     16   nlay_i         =    2           !  number of ice  layers 
     17   nlay_s         =    1           !  number of snow layers (only 1 is working) 
    1618   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input) 
    1719   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output) 
    1820   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F) 
    19    amax          = 0.999           !  maximum ice concentration 
    20    cai           = 1.40e-3         !  atmospheric drag over sea ice (clio) 
    21    cao           = 1.00e-3         !  atmospheric drag over ocean   (clio) 
    22    ln_nicep      = .false.         !  Ice points output for debug (yes or no) 
    23    ln_limdiahsb  = .false.          !  check the heat and salt budgets (T) or not (F) 
     21   rn_amax       = 0.999           !  maximum tolerated ice concentration  
     22   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F) 
    2423   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F) 
     24   ln_icectl     = .false.         !  ice points output for debug (T or F) 
     25   iiceprt       = 10              !  i-index for debug 
     26   jiceprt       = 10              !  j-index for debug 
    2527/ 
    26 !----------------------------------------------------------------------- 
    27 &namiceini     !   ice initialisation 
    28 !----------------------------------------------------------------------- 
    29    ln_limini   = .false.   !  activate ice initialization (T) or not (F) 
    30    thres_sst   =  0.0      !  threshold water temperature for initial sea ice 
    31    hts_ini_n   =  0.3      !  initial snow thickness in the north 
    32    hts_ini_s   =  0.3      !        "            "          south 
    33    hti_ini_n   =  1.0      !  initial ice thickness in the north 
    34    hti_ini_s   =  1.0      !        "            "         south 
    35    ati_ini_n   =  0.9      !  initial ice concentration in the north 
    36    ati_ini_s   =  0.9      !        "            "             south 
    37    smi_ini_n   =  6.301    !  initial ice salinity in the north 
    38    smi_ini_s   =  6.301    !        "            "    south 
    39    tmi_ini_n   =  270.     !  initial ice/snw temp in the north 
    40    tmi_ini_s   =  270.     !  initial ice/snw temp in the south 
     28!------------------------------------------------------------------------------ 
     29&namiceini     !   Ice initialization 
     30!------------------------------------------------------------------------------ 
     31   ln_iceini      = .true.         !  activate ice initialization (T) or not (F) 
     32   rn_thres_sst   =  2.0           !  maximum water temperature with initial ice (degC) 
     33   rn_hts_ini_n   =  0.3           !  initial real snow thickness (m), North 
     34   rn_hts_ini_s   =  0.3           !        "            "             South 
     35   rn_hti_ini_n   =  3.0           !  initial real ice thickness  (m), North 
     36   rn_hti_ini_s   =  1.0           !        "            "             South 
     37   rn_ati_ini_n   =  0.9           !  initial ice concentration   (-), North 
     38   rn_ati_ini_s   =  0.9           !        "            "             South 
     39   rn_smi_ini_n   =  6.3           !  initial ice salinity     (g/kg), North 
     40   rn_smi_ini_s   =  6.3           !        "            "             South 
     41   rn_tmi_ini_n   =  270.          !  initial ice/snw temperature (K), North 
     42   rn_tmi_ini_s   =  270.          !        "            "             South 
    4143/ 
    42 !----------------------------------------------------------------------- 
    43 &namicedyn     !   ice dynamic 
    44 !----------------------------------------------------------------------- 
    45    epsd        =   1.0e-20 !  tolerance parameter 
    46    om          =   0.5     !  relaxation constant  
    47    cw          =   5.0e-03 !  drag coefficient for oceanic stress 
    48    angvg       =   0.0     !  turning angle for oceanic stress 
    49    pstar       =   2.0e+04 !  1st bulk-rheology parameter 
    50    c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    51    creepl      =   1.0e-12 !  creep limit 
    52    ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    53    ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    54    nevp        = 300       !  number of iterations for subcycling in EVP, SB, 400 
    55    telast      =9600.0     !  timescale for elastic waves, SB, 720.0 
    56    alphaevp    =   1.0     !  coefficient for the solution of internal ice stresses 
    57    hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equal ocean velocity 
     44!------------------------------------------------------------------------------ 
     45&namiceitd     !   Ice discretization 
     46!------------------------------------------------------------------------------ 
     47   nn_catbnd      =    1           !  computation of ice category boundaries based on 
     48                                   !      1: tanh function 
     49                                   !      2: h^(-alpha), function of rn_himean 
     50   rn_himean      =    2.5         !  expected domain-average ice thickness (m), nn_catbnd = 2 only 
    5851/ 
    59 !----------------------------------------------------------------------- 
    60 &namicethd     !   ice thermodynamic 
    61 !----------------------------------------------------------------------- 
    62    hmelt       = -0.15     !  maximum melting at the bottom 
    63    hiccrit     = 0.1       !  ice thickness for lateral accretion  
    64                            !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ???? 
    65    fraz_swi    = 0         !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
    66    maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom 
    67    vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom 
    68    Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom 
    69    hiclim      = 0.10      !  minimum ice thickness 
    70    hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
    71    parsub      = 1.0       !  switch for snow sublimation or not 
    72    betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66 
    73                            !        betas = 1 -> equipartition, betas < 1 -> more on leads 
    74    kappa_i     = 1.0       !  extinction radiation parameter in sea ice (1.0) 
    75    nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation 
    76    maxer_i_thd = 0.0001    !  maximal error in temperature for heat diffusion computation 
    77    thcon_i_swi = 1         !  switch for computation of thermal conductivity in the ice 
    78                            !        (0) Untersteiner (1964), (1) Pringle et al. (2007) 
     52!------------------------------------------------------------------------------ 
     53&namicedyn     !   Ice dynamics and transport 
     54!------------------------------------------------------------------------------ 
     55   nn_icestr      =    0           !  ice strength parameteriztaion                       
     56                                   !     0: Hibler_79     P = pstar*<h>*exp(-c_rhg*A) 
     57                                   !     1: Rothrock_75   P = Cf*coeff*integral(wr.h^2)     
     58   ln_icestr_bvf  =    .false.     !  ice strength function brine volume (T) or not (F)      
     59   rn_pe_rdg      =   17.0         !  ridging work divided by pot. energy change in ridging, if nn_icestr = 1 
     60   rn_pstar       =    2.0e+04     !  ice strength thickness parameter (N/m2), nn_icestr = 0  
     61   rn_crhg        =   20.0         !  ice strength conc. parameter (-), nn_icestr = 0        
     62   rn_cio         =    5.0e-03     !  ice-ocean drag coefficient           (-)              
     63   rn_creepl      =    1.0e-12     !  creep limit (s-1)                                    
     64   rn_ecc         =    2.0         !  eccentricity of the elliptical yield curve           
     65   nn_nevp        =  120           !  number of EVP subcycles                              
     66   rn_relast      =    0.333       !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
     67                                   !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     68   nn_ahi0        =    2           !  horizontal diffusivity computation 
     69                                   !     0: use rn_ahi0_ref 
     70                                   !     1: use rn_ahi0_ref x mean grid cell length / ( 2deg mean grid cell length ) 
     71                                   !     2: use rn_ahi0_ref x grid cell length      / ( 2deg mean grid cell length ) 
     72   rn_ahi0_ref    = 350.0          !  horizontal sea ice diffusivity (m2/s)  
     73                                   !     if nn_ahi0 > 0, rn_ahi0_ref is the reference value at a nominal 2 deg resolution 
    7974/ 
    80 !----------------------------------------------------------------------- 
    81 &namicesal     !   ice salinity 
    82 !----------------------------------------------------------------------- 
    83    num_sal     =  2        !  salinity option: 1 -> S = bulk_sal 
    84                            !                   2 -> S = S(z,t) with a simple parameterization 
    85                            !                   3 -> S = S(z) profile of Scwharzacher [1959] 
    86                            !                   4 -> S = S(h) Cox and Weeks [1974] 
    87    bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity 
    88    sal_G       =  5.00     !  restoring salinity for GD 
    89    time_G      =  1.728e+6 !  restoring time for GD 
    90    sal_F       =  2.00     !  restoring salinity for flushing 
    91    time_F      =  8.640e+5 !  restoring time for flushing 
    92    s_i_max     = 20.0      !  Maximum salinity  
    93    s_i_min     =  0.1      !  Minimum tolerated ice salinity 
    94    s_i_0       =  3.5      !  1st salinity for salinity profile 
    95    s_i_1       =  4.5      !  2nd salinity for salinity profile 
     75!------------------------------------------------------------------------------ 
     76&namicethd     !   Ice thermodynamics 
     77!------------------------------------------------------------------------------ 
     78   rn_hnewice  = 0.1               !  thickness for new ice formation in open water (m) 
     79   ln_frazil   = .false.           !  use frazil ice collection thickness as a function of wind (T) or not (F) 
     80   rn_maxfrazb = 0.0               !  maximum fraction of frazil ice collecting at the ice base 
     81   rn_vfrazb   = 0.417             !  thresold drift speed for frazil ice collecting at the ice bottom (m/s) 
     82   rn_Cfrazb   = 5.0               !  squeezing coefficient for frazil ice collecting at the ice bottom 
     83   rn_himin    = 0.10              !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
     84   rn_betas    = 0.66              !  exponent in lead-ice repratition of snow precipitation 
     85                                   !     betas = 1 -> equipartition, betas < 1 -> more on leads 
     86   rn_kappa_i  = 1.0               !  radiation attenuation coefficient in sea ice (m-1) 
     87   nn_conv_dif = 50                !  maximal number of iterations for heat diffusion computation 
     88   rn_terr_dif = 0.0001            !  maximum temperature after heat diffusion (degC) 
     89   nn_ice_thcon= 1                 !  sea ice thermal conductivity 
     90                                   !     0: k = k0 + beta.S/T (Untersteiner, 1964) 
     91                                   !     1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 
     92   nn_monocat  = 0                 !  virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0) 
     93                                   !     2: simple piling instead of ridging --- temporary option 
     94                                   !     3: activate G(he) only              --- temporary option 
     95                                   !     4: activate lateral melting only    --- temporary option 
     96  ln_it_qnsice = .true.            !  iterate the surface non-solar flux with surface temperature (T) or not (F) 
    9697/ 
    97 !----------------------------------------------------------------------- 
    98 &namiceitdme   !   parameters for mechanical redistribution of ice  
    99 !----------------------------------------------------------------------- 
    100    ridge_scheme_swi = 0      !  which ridging scheme using (1=Rothrock,else=Hibler79) 
    101    Cs               =   0.50 !  shearing energy contribution to ridging 
    102    Cf               =  17.0  !  ratio of ridging work to PE change in ridging 
    103    fsnowrdg         =   0.5  !  snow fraction that survives in ridging 
    104    fsnowrft         =   0.5  !  snow fraction that survives in rafting 
    105    Gstar            =   0.15 !  fractional area of thin ice being ridged 
    106    astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
    107    Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice 
    108    raft_swi         =   1    !  rafting or not 
    109    hparmeter        =   0.75 !  threshold thickness for rafting or not 
    110    Craft            =   5.0  !  coefficient used in the rafting function 
    111    ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30) 
    112    partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1) 
    113    brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength 
     98!------------------------------------------------------------------------------ 
     99&namicesal     !   Ice salinity 
     100!------------------------------------------------------------------------------ 
     101   nn_icesal   =  2                !  ice salinity option 
     102                                   !     1: constant ice salinity (S=rn_icesal) 
     103                                   !     2: varying salinity parameterization S(z,t) 
     104                                   !     3: prescribed salinity profile S(z), Schwarzacher, 1959 
     105   rn_icesal   =  4.               !  ice salinity (g/kg, nn_icesal = 1 only) 
     106   rn_sal_gd   =  5.               !  restoring ice salinity, gravity drainage (g/kg) 
     107   rn_time_gd  =  1.73e+6          !  restoring time scale, gravity drainage  (s) 
     108   rn_sal_fl   =  2.               !  restoring ice salinity, flushing (g/kg) 
     109   rn_time_fl  =  8.64e+5          !  restoring time scale, flushing (s) 
     110   rn_simax    = 20.               !  maximum tolerated ice salinity (g/kg) 
     111   rn_simin    =  0.1              !  minimum tolerated ice salinity (g/kg) 
    114112/ 
    115 !----------------------------------------------------------------------- 
    116 &namicedia     !   ice diagnostics 
    117 !----------------------------------------------------------------------- 
    118    fmtinf      ='1PE13.5 ' !  format of the output values 
    119    nfrinf      = 4         !  number of variables written in one line 
    120    ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
    121    ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     113!------------------------------------------------------------------------------ 
     114&namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
     115!------------------------------------------------------------------------------ 
     116   rn_Cs       =   0.5             !  fraction of shearing energy contributing to ridging 
     117   rn_fsnowrdg =   0.5             !  snow volume fraction that survives in ridging 
     118   rn_fsnowrft =   0.5             !  snow volume fraction that survives in rafting 
     119   nn_partfun  =   1               !  type of ridging participation function 
     120                                   !     0: linear (Thorndike et al, 1975) 
     121                                   !     1: exponential (Lipscomb, 2007 
     122   rn_gstar    =   0.15            !  fractional area of thin ice being ridged (nn_partfun = 0) 
     123   rn_astar    =   0.05            !  exponential measure of ridging ice fraction (nn_partfun = 1) 
     124   rn_hstar    = 100.0             !  determines the maximum thickness of ridged ice (m) (Hibler, 1980) 
     125   ln_rafting  =   .true.          !  rafting activated (T) or not (F) 
     126   rn_hraft    =   0.75            !  threshold thickness for rafting (m) 
     127   rn_craft    =   5.0             !  squeezing coefficient used in the rafting function 
     128   rn_por_rdg  =   0.3             !  porosity of newly ridged ice (Lepparanta et al., 1995) 
    122129/ 
    123 !!----------------------------------------------------------------------- 
    124 !&namicehsb       !  Heat and salt budgets  
    125 !!----------------------------------------------------------------------- 
    126 !/ 
    127  
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_pisces_ORCA2_cfg

    r2221 r2451  
    1111&nampisatm     !  Atmospheric pressure  
    1212!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    13   ln_presatm  = .false.   ! constant atmopsheric pressure (F) or from a file (T) 
    1413/ 
    1514!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     
    5655&nampissbc     !   parameters for inputs deposition 
    5756!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    58    ln_river    =  .true.   ! AUTO - boolean for river input of nutrients 
    59    ln_hydrofe  =  .false.  ! AUTO - boolean for hydrothermal vents 
    6057/ 
    6158!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_pisces_ORCA2_ref

    r2221 r2451  
    4848!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    4949   concnno3   =  1.e-6    ! Nitrate half saturation of nanophytoplankton 
    50    concdno3   =  3.E-6    ! Phosphate half saturation for diatoms 
     50   concdno3   =  3.E-6    ! Nitrate half saturation for diatoms 
    5151   concnnh4   =  1.E-7    ! NH4 half saturation for phyto 
    5252   concdnh4   =  3.E-7    ! NH4 half saturation for diatoms 
    5353   concnfer   =  1.E-9    ! Iron half saturation for phyto 
    5454   concdfer   =  3.E-9    ! Iron half saturation for diatoms 
    55    concbfe    =  1.E-11   ! Half-saturation for Fe limitation of Bacteria 
    56    concbnh4   =  2.E-8    ! NH4 half saturation for phyto 
    57    concbno3   =  2.E-7    ! Phosphate half saturation for diatoms 
     55   concbfe    =  1.E-11   ! Iron half-saturation for DOC remin. 
     56   concbnh4   =  2.E-8    ! NH4 half saturation for DOC remin. 
     57   concbno3   =  2.E-7    ! Nitrate half saturation for DOC remin. 
    5858   xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
    5959   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     
    6161   xsizerd    =  3.0      ! Size ratio for diatoms 
    6262   xksi1      =  2.E-6    ! half saturation constant for Si uptake 
    63    xksi2      =  20E-6  ! half saturation constant for Si/C 
     63   xksi2      =  20E-6    ! half saturation constant for Si/C 
    6464   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
    6565   qnfelim    =  7.E-6    ! Optimal quota of phyto 
     
    8686   excret2    =  0.05     ! excretion ratio of diatoms 
    8787   ln_newprod =  .true.   ! Enable new parame. of production (T/F)  
    88    bresp      =  0.00333  ! Basal respiration rate 
    89    chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton 
    90    chlcdm     =  0.05     ! Minimum Chl/C in diatoms 
    91    chlcmin    =  0.004    ! Maximum Chl/c in phytoplankton 
     88   bresp      =  0.033    ! Basal respiration rate 
     89   chlcnm     =  0.033    ! Maximum Chl/C in nanophytoplankton 
     90   chlcdm     =  0.05     ! Maximum Chl/C in diatoms 
     91   chlcmin    =  0.004    ! Minimum Chl/c in phytoplankton 
    9292   fecnm      =  40E-6    ! Maximum Fe/C in nanophytoplankton 
    93    fecdm      =  40E-6    ! Minimum Fe/C in diatoms 
     93   fecdm      =  40E-6    ! Maximum Fe/C in diatoms 
    9494   grosip     =  0.159    ! mean Si/C ratio 
    9595/ 
     
    110110   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    111111   mzrat2     =  0.03     ! mesozooplankton mortality rate 
    112    xprefc     =  1.       ! zoo preference for phyto 
    113    xprefp     =  0.3      ! zoo preference for POC 
    114    xprefz     =  1.       ! zoo preference for zoo 
    115    xprefpoc   =  0.3      ! zoo preference for poc 
     112   xprefc     =  1.       ! mesozoo preference for diatoms 
     113   xprefp     =  0.3      ! mesozoo preference for nanophyto. 
     114   xprefz     =  1.       ! mesozoo preference for microzoo. 
     115   xprefpoc   =  0.3      ! mesozoo preference for poc 
    116116   xthresh2zoo = 1E-8     ! zoo feeding threshold for mesozooplankton  
    117117   xthresh2dia = 1E-8     ! diatoms feeding threshold for mesozooplankton  
     
    119119   xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton  
    120120   xthresh2   =  3E-7     ! Food threshold for grazing 
    121    xkgraz2    =  20.E-6   ! half sturation constant for meso grazing 
     121   xkgraz2    =  20.E-6   ! half saturation constant for meso grazing 
    122122   epsher2    =  0.35     ! Efficicency of Mesozoo growth 
    123123   sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
     
    156156&nampisrem     !   parameters for remineralization 
    157157!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    158    xremik    =  0.35      ! remineralization rate of DOC 
     158   xremik    =  0.3       ! remineralization rate of DOC 
    159159   xremip    =  0.025     ! remineralisation rate of POC 
    160160   nitrif    =  0.05      ! NH4 nitrification rate 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_top_ORCA2_cfg

    r2186 r2451  
    8787/ 
    8888!----------------------------------------------------------------------- 
     89&namtrc_dmp    !   passive tracer newtonian damping 
     90!----------------------------------------------------------------------- 
     91/ 
     92!----------------------------------------------------------------------- 
     93&namtrc_trd       !   diagnostics on tracer trends        ('key_trdtrc') 
     94!                          or mixed-layer trends          ('key_trdmld_trc') 
     95!---------------------------------------------------------------------- 
     96/ 
     97!----------------------------------------------------------------------- 
    8998&namtrc_dia       !   parameters for passive tracer additional diagnostics 
    9099!---------------------------------------------------------------------- 
    91100/ 
     101!---------------------------------------------------------------------- 
     102! namtrc_bc       !   data for boundary conditions 
     103!----------------------------------------------------------------------- 
     104&namtrc_bc 
     105! 
     106/ 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_top_ORCA2_ref

    r2221 r2451  
    7676&namtrc_dmp    !   passive tracer newtonian damping    
    7777!----------------------------------------------------------------------- 
    78    nn_hdmp_tr  =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
    79                            !                   =XX, damping poleward of XX degrees (XX>0) 
    80                            !                      + F(distance-to-coast) + Red and Med Seas 
    8178   nn_zdmp_tr  =    1      !  vertical   shape =0    damping throughout the water column 
    8279                           !                   =1 no damping in the mixing layer (kz  criteria) 
    8380                           !                   =2 no damping in the mixed  layer (rho crieria) 
    84    rn_surf_tr  =   50.     !  surface time scale of damping   [days] 
    85    rn_bot_tr   =  360.     !  bottom  time scale of damping   [days] 
    86    rn_dep_tr   =  800.     !  depth of transition between rn_surf and rn_bot [meters] 
    87    nn_file_tr  =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     81   cn_resto_tr  = 'resto_tr.nc'    !  create a damping.coeff NetCDF file (=1) or not (=0) 
    8882/ 
    8983!----------------------------------------------------------------------- 
  • CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/arch-X64_ADA.fcm

    r2214 r2451  
    3939%OASIS_LIB           -L%OASIS_HOME/lib -lpsmile.MPI1 -lscrip -lmct -lmpeu 
    4040 
     41%CPP               cpp 
    4142%FC                  mpiifort -c -cpp 
    4243# for Ada 
    43 %FCFLAGS             -DCPP_PARA -i4 -r8 -O3 
     44%FCFLAGS             -DCPP_PARA -i4 -r8 -O3 -xAVX -fp-model precise 
    4445%FFLAGS              %FCFLAGS 
    4546%LD                  mpiifort 
  • CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/arch-X64_CURIE.fcm

    r2194 r2451  
    3939%OASIS_LIB           -L%OASIS_HOME/lib -lpsmile.MPI1 -lscrip -lmct -lmpeu 
    4040 
     41%CPP                 cpp 
    4142%FC                  mpif90 -c -cpp 
    4243# for Curie Fat Node 
  • CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/p4zopt.F90

    r2268 r2451  
    4949   REAL(wp), DIMENSION(3,jp_rgb), PUBLIC ::   xkrgb   !: tabulated attenuation coefficients for RGB absorption 
    5050 
    51       
     51 
     52    
    5253   !!* Substitution 
    5354#  include "top_substitute.h90" 
     
    101102               zchl = ( trn(ji,jj,jk,jpnch) + trn(ji,jj,jk,jpdch) + rtrn ) * 1.e6 
    102103               zchl = MIN(  10. , MAX( 0.05, zchl )  ) 
    103                irgb = MAX ( 1, MIN ( jp_rgb, NINT( 41 + 20.* LOG10( zchl ) + rtrn ))) 
    104                !                                                          
     104               irgb = MAX ( 1, MIN ( jp_rgb, NINT( 41 + 20.* LOG10( zchl ) + rtrn )))               !                                                          
    105105               zekb(ji,jj,jk) = xkrgb(1,irgb) * fse3t(ji,jj,jk) 
    106106               zekg(ji,jj,jk) = xkrgb(2,irgb) * fse3t(ji,jj,jk) 
     
    207207              IF( etot(ji,jj,jk) * tmask(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )  THEN 
    208208                 neln(ji,jj) = jk+1                    ! Euphotic level : 1rst T-level strictly below Euphotic layer 
    209                  !                                     ! nb: ensure the compatibility with nmld_trc definition in trd_mld_trc_zint 
     209                 !                                     ! nb: ensure the compatibility with nmld_trc definition in trd_mxl_trc_zint 
    210210                 heup(ji,jj) = fsdepw(ji,jj,jk+1)      ! Euphotic layer depth 
    211211              ENDIF 
     
    254254      END DO 
    255255 
    256       IF( ln_diatrc ) THEN        ! save output diagnostics 
    257         ! 
    258         IF( lk_iomput ) THEN 
    259            IF( jnt == nrdttrc ) THEN 
    260               CALL iom_put( "Heup", heup(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht 
    261               CALL iom_put( "PAR" , emoy(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation 
    262            ENDIF 
    263         ELSE 
    264            trc2d(:,:,  jp_pcs0_2d + 10) = heup(:,:  ) * tmask(:,:,1)   
    265            trc3d(:,:,:,jp_pcs0_3d + 3)  = etot(:,:,:) * tmask(:,:,:) 
     256      IF( lk_iomput ) THEN 
     257        IF( jnt == nrdttrc  ) THEN 
     258           IF( iom_use( "Heup" ) ) CALL iom_put( "Heup", heup(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht 
     259           IF( iom_use( "PAR"  ) ) CALL iom_put( "PAR" , emoy(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation 
    266260        ENDIF 
    267         ! 
     261      ELSE 
     262         IF( ln_diatrc ) THEN        ! save output diagnostics 
     263            trc2d(:,:,  jp_pcs0_2d + 10) = heup(:,:  ) * tmask(:,:,1)   
     264            trc3d(:,:,:,jp_pcs0_3d + 3)  = etot(:,:,:) * tmask(:,:,:) 
     265         ENDIF 
    268266      ENDIF 
    269267      ! 
Note: See TracChangeset for help on using the changeset viewer.