Changeset 2061


Ignore:
Timestamp:
04/26/13 13:01:38 (11 years ago)
Author:
cetlod
Message:

IPSLCM6_rc0 : update ORCA2 namelists

Location:
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM
Files:
3 added
5 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/iodef.xml

    r1437 r2061  
    11<?xml version="1.0"?> 
    22<simulation>  
    3    
    4   <context id="nemo"> 
    5  
     3 
     4 <context id="nemo" time_origin="1900-01-01 00:00:00" > 
     5     
    66    <!-- $id$ --> 
    77     
     
    1212============================================================================================================ 
    1313    --> 
    14      
    15     <field_definition level="0" prec="4" operation="ave(X)" enabled=".TRUE."> <!-- time step automaticaly defined --> 
    16  
    17       <!-- T grid --> 
    18        
    19       <group id="grid_T" axis_ref="none" grid_ref="grid_T"> 
    20         <field id="toce"         description="temperature"                               unit="degC" axis_ref="deptht"   /> 
    21         <field id="soce"         description="salinity"                                  unit="psu"  axis_ref="deptht"   /> 
    22         <field id="sst"          description="sea surface temperature"                   unit="degC"                     /> 
    23         <field id="sst2"         description="square of sea surface temperature"         unit="degC2"                    /> 
    24         <field id="|sstgrad|"    description="module of sst gradient"                    unit="degC/m"                   /> 
    25         <field id="|sstgrad|2"   description="square of module of sst gradient"          unit="degC2/m2"                 /> 
    26         <field id="sss"          description="sea surface salinity"                      unit="psu"                      /> 
    27         <field id="sss2"         description="square of sea surface salinity"            unit="psu2"                     /> 
    28         <field id="ssh"          description="sea surface height"                        unit="m"                        /> 
    29         <field id="ssh2"         description="square of sea surface height"              unit="m2"                       /> 
    30         <field id="mldkz5"       description="mixing layer depth (Turbocline)"           unit="m"                        /> 
    31         <field id="mldr10_1"     description="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        /> 
    32         <field id="rhop"         description="potential density (sigma0)"                unit="kg/m3" axis_ref="deptht"  /> 
    33         <!-- next variables available with key_diahth --> 
    34         <field id="mlddzt"       description="Thermocline Depth (max dT/dz)"             unit="m"                        /> 
    35         <field id="mldr10_3"     description="Mixed Layer Depth dr=0.03 (ref.10m)"       unit="m"                        /> 
    36         <field id="mldr0_1"      description="Mixed Layer Depth dr=0.01 (ref.surf)"      unit="m"                        /> 
    37         <field id="mldr0_3"      description="Mixed Layer Depth dr=0.03 (ref.surf)"      unit="m"                        /> 
    38         <field id="mld|dt|"      description="Mixed Layer Depth |dt|=0.2 (ref.10m)"      unit="m"                        /> 
    39         <field id="topthdep"     description="Top of the thermocline dt=-0.2 (ref.10m)"  unit="m"                        /> 
    40         <field id="pycndep"      description="Pycnocline depth dr~dt=-0.2 (ref.10m)"     unit="m"                        /> 
    41         <field id="BLT"          description="Barrier Layer Thickness"                   unit="m"                        /> 
    42         <field id="tinv"         description="Max of vertical invertion of temperature"  unit="degC"                     /> 
    43         <field id="depti"        description="Depth of max. vert. inv. of temperature"   unit="m"                        /> 
    44         <field id="20d"          description="Depth of 20C isotherm"                     unit="m"                        /> 
    45         <field id="28d"          description="Depth of 28C isotherm"                     unit="m"                        /> 
    46         <field id="hc300"        description="Heat content 300 m"                        unit="W"                        /> 
    47         <!-- variables available with key_diaar5 --> 
    48         <field id="botpres"      description="Pressure at sea floor"                     unit="dbar"                     /> 
    49         <field id="cellthc"      description="Cell thickness"                            unit="m"     axis_ref="deptht"  /> 
    50      </group> 
    51  
    52       <!-- SBC --> 
    53        
    54       <group id="SBC" axis_ref="none" grid_ref="grid_T" > <!-- time step automaticaly defined based on nn_fsbc --> 
    55  
    56         <field id="emp"          description="Net Upward Water Flux"                                        unit="kg/m2/s"  /> 
    57         <field id="emps"         description="concentration/dilution water flux"                            unit="kg/m2/s"  /> 
    58         <field id="snowpre"      description="Snow precipitation"                                           unit="kg/m2/s"  /> 
    59         <field id="runoffs"      description="River Runoffs"                                                unit="kg/m2/s"  /> 
    60  
    61         <field id="qns+qsr"      description="Net Downward Heat Flux"                                       unit="W/m2"     /> 
    62         <field id="qns"          description="non solar Downward Heat Flux"                                 unit="W/m2"     /> 
    63         <field id="qsr"          description="Shortwave Radiation"                                          unit="W/m2"     /> 
    64         <field id="qsr3d"        description="Shortwave Radiation 3D distribution"        axis_ref="deptht" unit="W/m2"     /> 
    65         <field id="qrp"          description="Surface Heat Flux: Damping"                                   unit="W/m2"     /> 
    66         <field id="erp"          description="Surface Water Flux: Damping"                                  unit="kg/m2/s"  /> 
    67         <field id="taum"         description="wind stress module"                                           unit="N/m2"     /> 
    68         <field id="wspd"         description="Wind speed module at 10 m"                                    unit="m/s"      /> 
    69          
    70         <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
    71         <field id="qns_oce"      description="Non solar Downward Heat Flux over open ocean"                 unit="W/m2"     /> 
    72         <field id="qlw_oce"      description="Longwave Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    73         <field id="qsb_oce"      description="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    74         <field id="qla_oce"      description="Latent Downward Heat Flux over open ocean"                    unit="W/m2"     /> 
    75         <field id="taum_oce"     description="wind stress module over open ocean"                           unit="N/m2"     /> 
    76  
    77         <field id="ice_cover"    description="Ice fraction"                                                 unit="1"        /> 
    78  
    79         <field id="ioceflxb"     description="Oceanic flux at the ice base"                                 unit="W/m2"     /> 
    80         <field id="qsr_ai_cea"   description="Air-Ice downward solar heat flux (cell average)"              unit="W/m2"     /> 
    81         <field id="qns_ai_cea"   description="Air-Ice downward non-solar heat flux (cell average)"          unit="W/m2"     /> 
    82         <field id="qla_ai_cea"   description="Air-Ice downward Latent heat flux (cell average)"             unit="W/m2"     /> 
    83          
    84         <field id="qsr_io_cea"   description="Ice-Oce downward solar heat flux (cell average)"              unit="W/m2"     /> 
    85         <field id="qns_io_cea"   description="Ice-Oce downward non-solar heat flux (cell average)"          unit="W/m2"     /> 
    86          
    87         <field id="snowthic_cea" description="Snow thickness (cell average)"                                unit="m"        /> 
    88         <field id="icethic_cea"  description="Ice thickness (cell average)"                                 unit="m"        /> 
    89         <field id="iceprod_cea"  description="Ice production (cell average)"                                unit="m/s"      /> 
    90         <field id="sim"          description="Sea Ice Plus Surface Snow Amount"                             unit="kg/m2"    /> 
    91          
    92         <field id="ice_pres"     description="Ice presence"                                                 unit="-"        /> 
    93         <field id="ist_cea"      description="Ice surface temperature (cell average)"                       unit="degC"     /> 
    94         <field id="ist_ipa"      description="Ice surface temperature (ice presence average)"               unit="degC"     />       
    95         <field id="uice_ipa"     description="Ice velocity along i-axis at I-point (ice presence average)"  unit="m/s"      />       
    96         <field id="vice_ipa"     description="Ice velocity along j-axis at I-point (ice presence average)"  unit="m/s"      />       
    97          
    98         <field id="utau_ice"     description="Wind stress along i-axis over the ice at i-point"             unit="N/m2"     /> 
    99         <field id="vtau_ice"     description="Wind stress along j-axis over the ice at i-point"             unit="N/m2"     /> 
    100          
    101         <field id="u_imasstr"    description="Sea-ice mass transport along i-axis"                          unit="kg/s"     /> 
    102         <field id="v_imasstr"    description="Sea-ice mass transport along j-axis"                          unit="kg/s"     /> 
    103  
    104         <!-- available key_coupled --> 
    105         <field id="snow_ao_cea"  description="Snow over ice-free ocean (cell average)"                      unit="kg/m2/s"  /> 
    106         <field id="snow_ai_cea"  description="Snow over sea-ice (cell average)"                             unit="kg/m2/s"  /> 
    107         <field id="subl_ai_cea"  description="Sublimation over sea-ice (cell average)"                      unit="kg/m2/s"  /> 
    108         <field id="icealb_cea"   description="Ice albedo (cell average)"                                    unit="1"        /> 
    109         <field id="calving"      description="Calving"                                                      unit="kg/m2/s"  /> 
    110         <!-- available if key_coupled + conservative method --> 
    111         <field id="rain"         description="Liquid precipitation"                                         unit="kg/m2/s"  /> 
    112         <field id="evap_ao_cea"  description="Evaporation over ice-free ocean (cell average)"               unit="kg/m2/s"  /> 
    113         <!-- variables available with key_diaar5 --> 
    114         <field id="isnwmlt_cea"   description="Snow over Ice melting (cell average)"                        unit="kg/m2/s"  /> 
    115         <field id="fsal_virt_cea" description="Virtual salt flux due to ice formation (cell average)"       unit="kg/m2/s"  /> 
    116         <field id="fsal_real_cea" description="Real salt flux due to ice formation (cell average)"          unit="kg/m2/s"  /> 
    117         <field id="hflx_rain_cea" description="heat flux due to rainfall"                                   unit="W/m2"     /> 
    118         <field id="hflx_evap_cea" description="heat flux due to evaporation"                                unit="W/m2"     /> 
    119         <field id="hflx_snow_cea" description="heat flux due to snow falling over ice-free ocean"           unit="W/m2"     /> 
    120         <field id="hflx_ice_cea"  description="heat flux due to ice thermodynamics"                         unit="W/m2"     /> 
    121         <field id="hflx_rnf_cea"  description="heat flux due to runoffs"                                    unit="W/m2"     /> 
    122         <field id="hflx_cal_cea"  description="heat flux due to calving"                                    unit="W/m2"     /> 
    123         <field id="bicemel_cea"  description="Rate of Melt at Sea Ice Base (cell average)"                  unit="kg/m2/s"  /> 
    124         <field id="licepro_cea"  description="Lateral Sea Ice Growth Rate (cell average)"                   unit="kg/m2/s"  /> 
    125         <field id="snowmel_cea"  description="Snow Melt Rate (cell average)"                                unit="kg/m2/s"  /> 
    126         <field id="sntoice_cea"  description="Snow-Ice Formation Rate (cell average)"                       unit="kg/m2/s"  /> 
    127         <field id="ticemel_cea"  description="Rate of Melt at Upper Surface of Sea Ice (cell average)"      unit="kg/m2/s"  /> 
    128  
    129       </group> 
    130  
    131       <!-- U grid --> 
    132        
    133       <group id="grid_U"  axis_ref="depthu" grid_ref="grid_U"> 
    134         <field id="utau"         description="Wind Stress along i-axis"                    unit="N/m2" axis_ref="none" /> 
    135         <field id="uoce"         description="ocean current along i-axis"                  unit="m/s"                  /> 
    136         <field id="uoce_eff"     description="Effective ocean current along i-axis"        unit="m/s"                  /> 
    137         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    138         <field id="uoce_eiv"     description="EIV ocean current along i-axis"              unit="m/s"                  /> 
    139         <!-- uoce_eiv: available with key_trabbl_adv --> 
    140         <field id="uoce_bbl"     description="BBL ocean current along i-axis"              unit="m/s"                  /> 
    141         <!-- variables available with key_diaar5 --> 
    142         <field id="u_masstr"     description="ocean eulerian mass transport along i-axis"  unit="kg/s"                 /> 
    143         <field id="u_heattr"     description="ocean eulerian heat transport along i-axis"  unit="W"    axis_ref="none" /> 
    144         <field id="ueiv_heattr"  description="ocean bolus heat transport along i-axis"     unit="W"    axis_ref="none" /> 
    145         <field id="udiff_heattr" description="ocean diffusion heat transport along i-axis" unit="W"    axis_ref="none" /> 
    146      </group> 
    147        
    148       <!-- V grid --> 
    149        
    150       <group id="grid_V"  axis_ref="depthv" grid_ref="grid_V"> 
    151         <field id="vtau"         description="Wind Stress along j-axis"                    unit="N/m2" axis_ref="none" /> 
    152         <field id="voce"         description="ocean current along j-axis"                  unit="m/s"                  /> 
    153         <field id="voce_eff"     description="Effective ocean current along j-axis"        unit="m/s"                  /> 
    154         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    155         <field id="voce_eiv"     description="EIV ocean current along j-axis"              unit="m/s"                  /> 
    156         <!-- voce_eiv: available with key_trabbl_adv --> 
    157         <field id="voce_bbl"     description="BBL ocean current along j-axis"              unit="m/s"                  /> 
    158         <!-- variables available with key_diaar5 --> 
    159         <field id="v_masstr"     description="ocean eulerian mass transport along j-axis"  unit="kg/s"                 /> 
    160         <field id="v_heattr"     description="ocean eulerian heat transport along j-axis"  unit="W"    axis_ref="none" /> 
    161         <field id="veiv_heattr"  description="ocean bolus heat transport along j-axis"     unit="W"    axis_ref="none" /> 
    162         <field id="vdiff_heattr" description="ocean diffusion heat transport along j-axis" unit="W"    axis_ref="none" /> 
    163       </group> 
    164        
    165       <!-- W grid --> 
    166        
    167       <group id="grid_W"  axis_ref="depthw" grid_ref="grid_W"> 
    168         <field id="woce"         description="ocean vertical velocity"                     unit="m/s"                  /> 
    169         <field id="woce_eff"     description="effective ocean vertical velocity"           unit="m/s"                  /> 
    170         <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    171         <field id="woce_eiv"     description="EIV ocean vertical velocity"                 unit="m/s"                  /> 
    172         <!-- woce_eiv: available with key_trabbl_adv --> 
    173         <field id="woce_bbl"     description="BBL ocean vertical velocity"                 unit="m/s"                  /> 
    174         <field id="avt"          description="vertical eddy diffusivity"                   unit="m2/s"                 /> 
    175         <field id="avm"          description="vertical eddy viscosity"                     unit="m2/s"                 /> 
    176         <!-- avs: available with key_zdfddm --> 
    177         <field id="avs"          description="salt vertical eddy diffusivity"              unit="m2/s"                 /> 
    178         <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    179         <field id="avt_evd"      description="enhanced vertical diffusivity"               unit="m2/s"                 /> 
    180         <field id="avm_evd"      description="enhanced vertical viscosity"                 unit="m2/s"                 /> 
    181         <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d --> 
    182         <field id="aht2d"        description="lateral eddy diffusivity"                    unit="m2/s" axis_ref="none" /> 
    183         <field id="aht2d_eiv"    description="EIV lateral eddy diffusivity"                unit="m2/s" axis_ref="none" /> 
    184         <!-- avt_tide: available with key_zdftmx --> 
    185         <field id="av_tide"      description="tidal vertical diffusivity"                  unit="m2/s"                 /> 
    186         <!-- variables available with key_diaar5 -->    
    187         <field id="w_masstr"     description="vertical mass trasport"                      unit="kg/s"                 /> 
    188         <field id="w_masstr2"    description="square of vertical mass trasport"            unit="kg2/s2"               /> 
    189       </group> 
    190            
    191       <!-- scalar --> 
    192        
    193       <!-- variables available with key_diaar5 --> 
    194       <group id="scalar" axis_ref="none" grid_ref="scalarpoint" zoom_ref="1point" > 
    195         <field id="voltot"     description="global mean volume"                         unit="m3"   /> 
    196         <field id="sshtot"     description="global mean ssh"                            unit="m"    /> 
    197         <field id="sshsteric"  description="global mean ssh steric"                     unit="m"    /> 
    198         <field id="sshthster"  description="global mean ssh thermosteric"               unit="m"    /> 
    199         <field id="masstot"    description="global mean mass"                           unit="kg"   /> 
    200         <field id="temptot"    description="global mean temperature"                    unit="degC" /> 
    201         <field id="saltot"     description="global mean salinity"                       unit="psu"  /> 
    202         <field id="fram_trans" description="Sea Ice Mass Transport Through Fram Strait" unit="kg/s" /> 
    203       </group> 
    204  
    205       <!-- ptrc on T grid --> 
    206            
    207      <group id="ptrc_T" axis_ref="deptht" grid_ref="grid_T">   
    208        <field id="DIC"      description="Dissolved inorganic Concentration"        unit="molC/L" /> 
    209        <field id="Alkalini" description="Total Alkalinity Concentration"           unit="eq/L"    /> 
    210        <field id="O2"       description="Oxygen Concentration"                     unit="molO2/L" /> 
    211        <field id="CaCO3"    description="Calcite Concentration"                    unit="molC/L" /> 
    212        <field id="PO4"      description="Phosphate Concentration"                  unit="molC/L" /> 
    213        <field id="POC"      description="Small organic carbon Concentration"       unit="molC/L" /> 
    214        <field id="Si"       description="Silicate Concentration"                   unit="molSi/L" /> 
    215        <field id="PHY"      description="Nanophytoplankton Concentration"          unit="molC/L" /> 
    216        <field id="ZOO"      description="Microzooplankton Concentration"           unit="molC/L" /> 
    217        <field id="DOC"      description="Dissolved organic Concentration"          unit="molC/L" /> 
    218        <field id="PHY2"     description="Diatoms Concentration"                    unit="molC/L" /> 
    219        <field id="ZOO2"     description="Mesozooplankton Concentration"            unit="molC/L" /> 
    220        <field id="BSi"      description="Diatoms Silicate Concentration"           unit="molC/L" /> 
    221        <field id="Fer"      description="Dissolved Iron Concentration"             unit="molFe/L" /> 
    222        <field id="BFe"      description="Big iron particles Concentration"         unit="molFe/L" /> 
    223        <field id="GOC"      description="Big organic carbon Concentration"         unit="molC/L" /> 
    224        <field id="SFe"      description="Small iron particles Concentration"       unit="molFe/L" /> 
    225        <field id="DFe"      description="Diatoms iron  Concentration"              unit="molFe/L" /> 
    226        <field id="DSi"      description="Sinking biogenic Silicate Concentration"  unit="molC/L" /> 
    227        <field id="NFe"      description="Nano iron Concentration"                  unit="molC/L" /> 
    228        <field id="NCHL"     description="Nano chlorophyl Concentration"            unit="gChl/L" /> 
    229        <field id="DCHL"     description="Diatoms chlorophyl Concentration"         unit="gChl/L" /> 
    230        <field id="NO3"      description="Nitrates Concentration"                   unit="molC/L" /> 
    231        <field id="NH4"      description="Ammonium Concentration"                   unit="molC/L" /> 
    232      </group> 
    233  
    234       <!-- diad on T grid : : variables available with key_trc_diaadd && key_trc_dia3d --> 
    235  
    236      <group id="diad_T" axis_ref="none" grid_ref="grid_T"> 
    237        <field id="PH"          description="PH"                                       unit="-"           axis_ref="deptht" /> 
    238        <field id="CO3"         description="Bicarbonates"                             unit="mol/L"       axis_ref="deptht" /> 
    239        <field id="CO3sat"      description="CO3 saturation"                           unit="mol/L"       axis_ref="deptht" /> 
    240        <field id="PAR"         description="Photosynthetically Available Radiation"   unit="W/m2"        axis_ref="deptht" /> 
    241        <field id="PPPHY"       description="Primary production of nanophyto"          unit="molC/m3/s"   axis_ref="deptht" /> 
    242        <field id="PPPHY2"      description="Primary production of diatoms"            unit="molC/m3/s"   axis_ref="deptht" /> 
    243        <field id="PPNEWN"      description="New Primary production of nanophyto"      unit="molC/m3/s"   axis_ref="deptht" /> 
    244        <field id="PPNEWD"      description="New Primary production of diatoms"        unit="molC/m3/s"   axis_ref="deptht" /> 
    245        <field id="PBSi"        description="Primary production of Si diatoms"         unit="molSi/m3/s"  axis_ref="deptht" /> 
    246        <field id="PFeN"        description="Primary production of nano iron"          unit="molFe/m3/s"  axis_ref="deptht" /> 
    247        <field id="PFeD"        description="Primary production of diatoms iron"       unit="molFe/m3/s"  axis_ref="deptht" /> 
    248        <field id="PCAL"        description="Calcite production"                       unit="molC/m3/s"   axis_ref="deptht" /> 
    249        <field id="DCAL"        description="Calcite dissolution"                      unit="molC/m3/s"   axis_ref="deptht" /> 
    250        <field id="GRAZ"        description="Grazing by zooplankton"                   unit="molC/m3/s"   axis_ref="deptht" /> 
    251        <field id="Nfix"        description="Nitrogen fixation at surface"             unit="molN/m2/s"      /> 
    252        <field id="EPC100"      description="Export of carbon particles at 100 m"      unit="molC/m2/s"      /> 
    253        <field id="EPFE100"     description="Export of biogenic iron at 100 m"         unit="molFe/m2/s"     /> 
    254        <field id="EPSI100"     description="Export of Silicate at 100 m"              unit="molSi/m2/s"     /> 
    255        <field id="EPCAL100"    description="Export of Calcite at 100 m"               unit="molC/m2/s"      /> 
    256        <field id="Cflx"        description="DIC flux"                                 unit="molC/m2/s"      /> 
    257        <field id="Oflx"        description="Oxygen flux"                              unit="molC/m2/s"      /> 
    258        <field id="Kg"          description="Gas transfer"                             unit="molC/m2/s/uatm" /> 
    259        <field id="Dpco2"       description="Delta CO2"                                unit="uatm"            /> 
    260        <field id="Dpo2"        description="Delta O2"                                 unit="uatm"            /> 
    261        <field id="Heup"        description="Euphotic layer depth"                     unit="m"               /> 
    262        <field id="Irondep"     description="Iron deposition"                          unit="molFe/m2/s"     /> 
    263      </group> 
    264  
    265       <!-- dbio_T on T grid : variables available with key_diaar5 --> 
    266      <group id="dbio_T" axis_ref="none" grid_ref="grid_T"> 
    267        <field id="PHYT"        description="Total phytoplankton Concentration"                   unit="molC/L"     axis_ref="deptht" /> 
    268        <field id="ZOOT"        description="Total zooplankton Concentration"                     unit="molC/L"     axis_ref="deptht" /> 
    269        <field id="CHLT"        description="Total chlorophyll Concentration"                     unit="gChl/L"     axis_ref="deptht" /> 
    270        <field id="POCT"        description="Total carbon particles Concentration"                unit="molC/L"     axis_ref="deptht" /> 
    271        <field id="PFET"        description="Total biogenic Iron Concentration"                   unit="molFe/L"     axis_ref="deptht" /> 
    272        <field id="TPP"         description="Total Primary production of phyto"                   unit="molC/m3/s"  axis_ref="deptht" /> 
    273        <field id="TPNEW"       description="New Primary production of phyto"                     unit="molC/m3/s"  axis_ref="deptht" /> 
    274        <field id="TPBFE"       description="Total biogenic iron production"                      unit="molFe/m3/s" axis_ref="deptht" /> 
    275        <field id="EXPC"        description="Export of carbon"                                    unit="molC/m2/s"  axis_ref="deptht" /> 
    276        <field id="EXPFE"       description="Export of biogenic iron"                             unit="molFe/m2/s" axis_ref="deptht" /> 
    277        <field id="EXPSI"       description="Export of Silicate"                                  unit="molSi/m2/s" axis_ref="deptht" /> 
    278        <field id="EXPCAL"      description="Export of Calcite"                                   unit="molC/m2/s"  axis_ref="deptht" /> 
    279        <field id="PHSFC"       description="PH at surface"                                       unit="-"       /> 
    280        <field id="INTDIC"      description="DIC content"                                         unit="kg/m2"   /> 
    281        <field id="O2MIN"       description="Oxygen minimum concentration"                        unit="molC/L" /> 
    282        <field id="ZO2MIN"      description="Depth of oxygen minimum concentration"               unit="m"       /> 
    283        <field id="INTNFIX"     description="Nitrogen fixation rate : vert. integrated"           unit="molN/m2/s"  /> 
    284        <field id="INTPPPHY"    description="Vertically integrated primary production by nanophy" unit="molC/m2/s"  /> 
    285        <field id="INTPPPHY2"   description="Vertically integrated primary production by diatom"  unit="molC/m2/s"  /> 
    286        <field id="INTPP"       description="Vertically integrated primary production by phyto"   unit="molC/m2/s"  /> 
    287        <field id="INTPNEW"     description="Vertically integrated new primary production"        unit="molC/m2/s"  /> 
    288        <field id="INTPBFE"     description="Vertically integrated of biogenic iron production"   unit="molFe/m2/s" /> 
    289        <field id="INTPBSI"     description="Vertically integrated of biogenic Si production"     unit="molSi/m2/s" /> 
    290        <field id="INTPCAL"     description="Vertically integrated of calcite production"         unit="molSi/m2/s" /> 
    291      </group> 
    292  
    293     </field_definition> 
    294      
     14    <field_definition src="./field_def.xml"/> 
    29515    <!--  
    29616============================================================================================================ 
    297 =                                                                                                          = 
    29817=                                           output files definition                                        = 
    299 =                                                                                                          = 
    300 =                  http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Oclim                 = 
    301 =                  http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Omon                  = 
    302 =                  http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Oimon                 =               
    303 =          http://www.clivar.org/organization/wgomd/references/WGOMD_CMIP5_ocean_fields.June2009.pdf       = 
    304 =                                                                                                          = 
     18=                                            Define your own files                                         = 
     19=                                         put the variables you want...                                    = 
    30520============================================================================================================ 
    30621    --> 
    30722     
    308     <file_definition output_level="<OUTPUT_LEVEL>" > 
    309       <!-- 
    310 +++++++++++++++++++++++++++++++++++++++++++++++   daily   ++++++++++++++++++++++++++++++++++++++++++++++++++ 
    311         --> 
    312       <group id="1d" output_freq="86400"  enabled="<1D_ENABLE>">                             <!-- 1d files --> 
    313         <!-- 
    314 ..............................................    grid T   ................................................. 
    315           --> 
    316         <file id="ORCA2_1d_grid_T" description="table 2.2" >                                 <!-- grid T  Table 2.2 --> 
    317           <field ref="sst"          name="sosstsst" description="sea_surface_temperature"                              level="2"                      /> 
    318           <field ref="sst2"         name="tossq"    description="square_of_sea_surface_temperature"                    level="2"                      /> 
    319           <field ref="mldkz5"       name="omldamax" description="ocean_mixed_layer_thickness_defined_by_mixing_scheme" level="2" operation="t_max(X)" /> 
    320          <!--  IPSL variables --> 
    321           <field ref="sss"          name="sosaline" description="sea_surface_salinity"                                 level="3"                      /> 
    322           <field ref="ssh"          name="sossheig" description="sea_surface_height_above_geoid"                       level="3"                      /> 
    323           <field ref="toce"         name="votemper" description="sea_water_potential_temperature"                      level="3"                      /> 
    324           <field ref="soce"         name="vosaline" description="sea_water_salinity"                                   level="3"                      /> 
    325         </file>  
    326         <!-- 
    327 ..............................................    grid U   ................................................. 
    328           --> 
    329         <file id="ORCA2_1d_grid_U" description="ocean U grid variables : table 2.3" >                       <!-- grid U --> 
    330  
    331           <!--  table 2.3 --> 
    332             <field ref="uoce"         name="vozocrtx" description="sea_water_x_velocity"                                level="3"       /> 
    333  
    334         </file>  
    335  
    336         <!-- 
    337 ..............................................    grid V   ................................................. 
    338           --> 
    339         <file id="ORCA2_1d_grid_V" description="ocean V grid variables : table 2.3" >                      <!-- grid V --> 
    340  
    341           <!-- table 2.3 --> 
    342             <field ref="voce"         name="vomecrty" description="sea_water_y_velocity"                                level="3"       /> 
    343  
    344         </file>    
     23    <file_definition type="multiple_file" sync_freq="1d" min_digits="4"> 
     24     
     25      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> <!-- 1h files --> 
     26      <file_group id="2h" output_freq="2h"  output_level="10" enabled=".TRUE."/> <!-- 2h files --> 
     27      <file_group id="3h" output_freq="3h"  output_level="10" enabled=".TRUE."/> <!-- 3h files -->      
     28      <file_group id="4h" output_freq="4h"  output_level="10" enabled=".TRUE."/> <!-- 4h files --> 
     29      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files --> 
     30      
     31      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."> <!-- 1d files --> 
     32 
     33        <file id="1d_grid_T" name="auto" description="ocean T grid variables" > 
     34          <field field_ref="sst"          name="sosstsst"  /> 
     35          <field field_ref="sss"          name="sosaline"  /> 
     36          <field field_ref="ssh"          name="sossheig"  /> 
     37        </file> 
     38 
     39        <file id="1d_grid_U" name="auto" description="ocean U grid variables" > 
     40          <field field_ref="suoce"         name="vozocrtx"  /> 
     41        </file> 
     42         
     43        <file id="1d_grid_V" name="auto" description="ocean V grid variables" > 
     44          <field field_ref="svoce"         name="vomecrty"  /> 
     45        </file> 
     46         
     47      </file_group> 
     48      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->     
     49      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE.">  <!-- 5d files -->    
     50 
     51        <file id="5d_grid_T" name="auto" description="ocean T grid variables" > 
     52          <field field_ref="toce"         name="votemper"  /> 
     53          <field field_ref="soce"         name="vosaline"  /> 
     54          <field field_ref="sst"          name="sosstsst"  /> 
     55          <field field_ref="sss"          name="sosaline"  /> 
     56          <field field_ref="ssh"          name="sossheig"  /> 
     57          <field field_ref="empmr"        name="sowaflup"  /> 
     58          <field field_ref="qsr"          name="soshfldo"  /> 
     59          <field field_ref="saltflx"      name="sosfldow"  /> 
     60          <field field_ref="qt"           name="sohefldo"  /> 
     61          <field field_ref="mldr10_1"     name="somxl010"  /> 
     62          <field field_ref="mldkz5"       name="somixhgt"  /> 
     63          <field field_ref="ice_cover"    name="soicecov"  /> 
     64          <field field_ref="wspd"         name="sowindsp"  /> 
     65          <field field_ref="erp"          name="sowafldp"  /> 
     66          <field field_ref="ist_ipa"      name="soicetem"  /> 
     67          <field field_ref="icealb_cea"   name="soicealb"  />    
     68        </file> 
     69         
     70        <file id="5d_grid_U" name="auto" description="ocean U grid variables" > 
     71          <field field_ref="uoce"         name="vozocrtx"  /> 
     72          <field field_ref="uoce_eiv"     name="vozoeivu"  /> 
     73          <field field_ref="utau"         name="sozotaux"  /> 
     74        </file> 
     75         
     76        <file id="5d_grid_V" name="auto" description="ocean V grid variables" > 
     77          <field field_ref="voce"         name="vomecrty"  /> 
     78          <field field_ref="voce_eiv"     name="vomeeivv"  /> 
     79          <field field_ref="vtau"         name="sometauy"  /> 
     80        </file> 
     81         
     82        <file id="5d_grid_W" name="auto" description="ocean W grid variables" > 
     83          <field field_ref="woce"         name="vovecrtz" /> 
     84          <field field_ref="avt"          name="votkeavt" /> 
     85          <field field_ref="avs"          name="voddmavs" /> 
     86          <field field_ref="aht2d_eiv"    name="soleaeiw" /> 
     87        </file> 
     88         
     89        <file id="5d_icemod" name="auto" description="ice variables" > 
     90          <field field_ref="ice_pres"                     /> 
     91          <field field_ref="snowthic_cea" name="isnowthi" /> 
     92          <field field_ref="icethic_cea"  name="iicethic" /> 
     93          <field field_ref="iceprod_cea"  name="iiceprod" /> 
     94          <field field_ref="ist_ipa"      name="iicetemp" /> 
     95          <field field_ref="ioceflxb"     name="ioceflxb" /> 
     96          <field field_ref="uice_ipa"     name="iicevelu" /> 
     97          <field field_ref="vice_ipa"     name="iicevelv" /> 
     98          <field field_ref="utau_ice"     name="iicestru" /> 
     99          <field field_ref="vtau_ice"     name="iicestrv" /> 
     100          <field field_ref="qsr_io_cea"   name="iicesflx" /> 
     101          <field field_ref="qns_io_cea"   name="iicenflx" /> 
     102          <field field_ref="snowpre"      name="isnowpre" /> 
     103        </file> 
     104         
     105      </file_group> 
     106 
     107      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files --> 
     108 
     109        <file id="1m_ptrc_T" name="auto" description="pisces sms variables" > 
     110          <field field_ref="DIC"      /> 
     111          <field field_ref="Alkalini" /> 
     112          <field field_ref="O2"       /> 
     113          <field field_ref="PO4"      /> 
     114          <field field_ref="Si"       /> 
     115          <field field_ref="Fer"      /> 
     116          <field field_ref="NCHL"     /> 
     117          <field field_ref="DCHL"     /> 
     118          <field field_ref="NO3"      /> 
     119        </file> 
     120         
     121        <file id="1m_diad_T" name="auto" description="additional pisces diagnostics" > 
     122          <field field_ref="Cflx"     /> 
     123          <field field_ref="Dpco2"    /> 
     124        </file> 
     125 
     126      </file_group> 
     127      <file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files --> 
     128      <file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files --> 
     129      <file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files --> 
     130      <file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files --> 
     131 
     132      <file_group id="1y"  output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files --> 
     133 
     134        <file id="1y_ptrc_T" name="auto" description="pisces sms variables" > 
     135          <field field_ref="DIC"      /> 
     136          <field field_ref="Alkalini" /> 
     137          <field field_ref="O2"       /> 
     138          <field field_ref="CaCO3"    /> 
     139          <field field_ref="PO4"      /> 
     140          <field field_ref="POC"      /> 
     141          <field field_ref="Si"       /> 
     142          <field field_ref="PHY"      /> 
     143          <field field_ref="ZOO"      /> 
     144          <field field_ref="DOC"      /> 
     145          <field field_ref="PHY2"     /> 
     146          <field field_ref="ZOO2"     /> 
     147          <field field_ref="DSi"      /> 
     148          <field field_ref="Fer"      /> 
     149          <field field_ref="BFe"      /> 
     150          <field field_ref="GOC"      /> 
     151          <field field_ref="SFe"      /> 
     152          <field field_ref="DFe"      /> 
     153          <field field_ref="GSi"      /> 
     154          <field field_ref="NFe"      /> 
     155          <field field_ref="NCHL"     /> 
     156          <field field_ref="DCHL"     /> 
     157          <field field_ref="NO3"      /> 
     158          <field field_ref="NH4"      /> 
     159        </file> 
     160 
     161        <file id="1y_diad_T" name="auto" description="additional pisces diagnostics" > 
     162          <field field_ref="PH"       /> 
     163          <field field_ref="CO3"      /> 
     164          <field field_ref="CO3sat"   /> 
     165          <field field_ref="PAR"      /> 
     166          <field field_ref="PPPHY"    /> 
     167          <field field_ref="PPPHY2"   /> 
     168          <field field_ref="PPNEWN"   /> 
     169          <field field_ref="PPNEWD"   /> 
     170          <field field_ref="PBSi"     /> 
     171          <field field_ref="PFeN"     /> 
     172          <field field_ref="PFeD"     /> 
     173          <field field_ref="xfracal"  /> 
     174          <field field_ref="PCAL"     /> 
     175          <field field_ref="DCAL"     /> 
     176          <field field_ref="GRAZ1"    /> 
     177          <field field_ref="GRAZ2"    /> 
     178          <field field_ref="EPC100"   /> 
     179          <field field_ref="EPFE100"  /> 
     180          <field field_ref="EPSI100"  /> 
     181          <field field_ref="EPCAL100" /> 
     182          <field field_ref="Cflx"     /> 
     183          <field field_ref="Oflx"     /> 
     184          <field field_ref="Kg"       /> 
     185          <field field_ref="Dpco2"    /> 
     186          <field field_ref="Dpo2"     /> 
     187          <field field_ref="Heup"     /> 
     188          <field field_ref="Irondep"  /> 
     189          <field field_ref="Ironsed"  /> 
     190          <field field_ref="Ironice"  /> 
     191          <field field_ref="Nfix"     /> 
     192          <field field_ref="MuN"      /> 
     193          <field field_ref="MuD"      /> 
     194          <field field_ref="LNnut"    /> 
     195          <field field_ref="LDnut"    /> 
     196          <field field_ref="LNFe"     /> 
     197          <field field_ref="LDFe"     /> 
     198          <field field_ref="LNlight"  /> 
     199          <field field_ref="LDlight"  /> 
     200          <field field_ref="pdust"    /> 
     201          <field field_ref="Fe2"      /> 
     202          <field field_ref="Fe3"      /> 
     203          <field field_ref="FeL1"     /> 
     204          <field field_ref="FeL2"     /> 
     205          <field field_ref="FeP"      /> 
     206          <field field_ref="TL1"      /> 
     207          <field field_ref="TL2"      /> 
     208          <field field_ref="Sdenit"   /> 
     209          <field field_ref="Totlig"   /> 
     210        </file> 
     211 
     212      </file_group> 
     213      <file_group id="2y"  output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files --> 
     214      <file_group id="5y"  output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files --> 
     215      <file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files --> 
     216 
     217   </file_definition> 
     218     
     219    <!--  
     220============================================================================================================ 
     221= grid definition = = DO NOT CHANGE = 
     222============================================================================================================ 
     223    --> 
     224     
     225   <axis_definition>   
     226      <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 
     227      <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 
     228      <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 
     229      <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 
     230      <axis id="nfloat" long_name="Float number"      unit="-"  /> 
     231   </axis_definition>  
     232     
     233   <domain_definition src="./domain_def.xml"/> 
     234    
     235   <grid_definition>     
     236     <grid id="grid_T_2D" domain_ref="grid_T"/> 
     237     <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/> 
     238     <grid id="grid_U_2D" domain_ref="grid_U"/> 
     239     <grid id="grid_U_3D" domain_ref="grid_U" axis_ref="depthu"/> 
     240     <grid id="grid_V_2D" domain_ref="grid_V"/> 
     241     <grid id="grid_V_3D" domain_ref="grid_V" axis_ref="depthv"/> 
     242     <grid id="grid_W_2D" domain_ref="grid_W"/> 
     243     <grid id="grid_W_3D" domain_ref="grid_W" axis_ref="depthw"/> 
     244    </grid_definition>   
    345245   
    346         <!-- 
    347 ..............................................    icemod   ................................................. 
    348           --> 
    349         <file id="ORCA2_1d_icemod" description=" ice variables" >                                           <!-- icemod --> 
    350  
    351          <!-- table 2.2 --> 
    352           <field ref="uice_ipa"     name="iicevelu"   description="Ice velocity along i-axis at I-point (ice presence average)"  level="2"    /> 
    353           <field ref="vice_ipa"     name="iicevelv"   description="Ice velocity along j-axis at I-point (ice presence average)"  level="2"    /> 
    354           <field ref="ice_cover"    name="ileadfra"   description="sea_ice_area_fraction"                                        level="2"    /> 
    355           <field ref="icethic_cea"  name="iicethic"   description="sea_ice_thickness"                                            level="2"    /> 
    356         </file>    
    357  
    358       </group> 
    359         <!-- 
    360 ++++++++++++++++++++++++++++++++++++++++++++++   monthly   +++++++++++++++++++++++++++++++++++++++++++++++++ 
    361           --> 
    362       <group id="1m" output_freq="-1"     enabled="<1M_ENABLE>">                             <!-- real monthly files --> 
    363         <!-- 
    364 ..............................................    grid T   ................................................. 
    365           --> 
    366         <file id="ORCA2_1m_grid_T" description="ocean T grid variables table 2.2 2.5 2.6 2.7"> <!-- grid T --> 
    367  
    368           <!-- table 2.2 --> 
    369             <field ref="botpres"      name="pbo"      description="sea_water_pressure_at_sea_floor"                               /> 
    370             <!-- pso : sea_water_pressure_at_sea_water_surface = 0 --> 
    371             <field ref="ssh"          name="sossheig" description="sea_surface_height_above_geoid"                                 /> 
    372             <field ref="ssh2"         name="zossq"    description="square_of_sea_surface_height_above_geoid"             level="2" /> 
    373             <!-- masscello : sea_water_mass_per_unit_area = cellthc*rau0                                          no time changes --> 
    374             <field ref="cellthc"      name="thkcello" description="cell_thickness"                                                 /> <!-- no time changes --> 
    375             <field ref="toce"         name="votemper" description="sea_water_potential_temperature"                                /> 
    376             <field ref="sst"          name="sosstsst" description="sea_surface_temperature"                              level="1" /> 
    377             <field ref="sst2"         name="tossq"    description="square_of_sea_surface_temperature"                    level="2" /> 
    378             <field ref="soce"         name="vosaline" description="sea_water_salinity"                                             /> 
    379             <field ref="sss"          name="sosaline" description="sea_surface_salinity"                                 level="1" /> 
    380             <field ref="rhop"         name="rhopoto"  description="sea_water_potential_density"                          level="2" /> 
    381             <!-- no agessc : sea_water_age_since_surface_contact --> 
    382             <!-- no cfc11  : moles_per_unit_mass_of_cfc11_in_sea_water --> 
    383             <!-- msftbarot : ocean_barotropic_mass_streamfunction : offline --> 
    384             <!-- mlotst    :           ocean_mixed_layer_thickness_defined_by_sigma_t : must be done offline --> 
    385             <!-- mlotstsq  : square_of_ocean_mixed_layer_thickness_defined_by_sigma_t : must be done offline --> 
    386             <field ref="mldkz5"       name='omlmax'   description="ocean_mixed_layer_thickness_defined_by_mixing_scheme" level="2" operation="t_max(X)" /> 
    387  
    388           <!-- table 2.5 --> 
    389             <field ref="rain"         name="pr"       description="rainfall_flux"                                           level="1" /> 
    390             <field ref="snow_ao_cea"  name="prsn"     description="snowfall_flux"                                           level="1" /> 
    391             <field ref="evap_ao_cea"  name="evs"      description="water_evaporation_flux"                                  level="1" /> 
    392             <field ref="runoffs"      name="friver"   description="water_flux_into_sea_water_from_rivers"                   level="1" /> 
    393             <field ref="calving"      name="ficeberg" description="water_flux_into_sea_water_from_icebergs"                 level="1" /> 
    394             <field ref="isnwmlt_cea"  name="fsitherm" description="water_flux_into_sea_water_due_to_sea_ice_thermodynamics" level="1" /> 
    395             <field ref="emp"          name="sowaflup" description="water_flux_into_sea_water"                               level="1" /> 
    396             <!-- wfonocorr : water_flux_into_sea_water_without_flux_correction : emp - erp -->   
    397             <field ref="erp"          name="wfcorr"   description="water_flux_correction"                                   level="1" /> <!-- usually = 0 --> 
    398  
    399           <!-- table2.6 --> 
    400             <!-- vsfpr    : virtual_salt_flux_into_sea_water_due_to_rainfall    = 0 -->  
    401             <!-- vsfevap  : virtual_salt_flux_into_sea_water_due_to_evaporation = 0 -->  
    402             <!-- vsfriver : virtual_salt_flux_into_sea_water_from_rivers        = 0 -->            
    403             <field ref="fsal_virt_cea" name="vsfsit" description="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics" level="1" /> 
    404             <!-- vsf      : virtual_salt_flux_into_sea_water = fsal_virtual + fsal_real --> 
    405             <!-- wfcorr   : virtual_salt_flux_correction                        = 0 -->  
    406             <field ref="fsal_virt_cea" name="sfdsi"  description="downward_sea_ice_basal_salt_flux"                               level="1" /> 
    407             <!-- sfriver  : salt_flux_into_sea_water_from_rivers                = 0 -->  
    408  
    409           <!--  table 2.7 --> 
    410             <!-- hfgeou : upward_geothermal_heat_flux_at_sea_floor : cte, see nambbc and trabbc.F90 --> 
    411             <field ref="hflx_rain_cea" name="hfrainds"    description="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"      level="1" /> 
    412             <field ref="hflx_evap_cea" name="hfevapds"    description="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" level="1" /> 
    413             <field ref="hflx_rnf_cea"  name="hfrunoffds"  description="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water"        level="1" /> 
    414             <field ref="hflx_snow_cea" name="hfsnthermds" description="heat_flux_into_sea_water_due_to_snow_thermodynamics"                         level="1" /> 
    415             <field ref="hflx_ice_cea"  name="hfsithermds" description="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics"                      level="1" /> 
    416             <field ref="hflx_cal_cea"  name="hfibthermds" description="heat_flux_into_sea_water_due_to_iceberg_thermodynamics"                      level="1" /> 
    417             <!-- rlds   : surface_net_downward_longwave_flux : not available --> 
    418             <!-- hfls   : surface_downward_latent_heat_flux  : not available --> 
    419             <!-- hfss   : surface_downward_sensible_heat_flux: not available --> 
    420             <field ref="qns"           name="nshfls"      description="surface_net_downward_non_solar_flux"                                         level="1" /> 
    421             <field ref="qsr"           name="soshfldo"    description="surface_net_downward_shortwave_flux"                                         level="1" /> 
    422             <field ref="qsr3d"         name="rsds"        description="downwelling_shortwave_flux_in_sea_water"                                     level="1" /> 
    423             <field ref="qrp"           name="hfcorr"      description="heat_flux_correction"                                                        level="1" /> 
    424  
    425          <!--  IPSL variables --> 
    426         <field ref="qns+qsr"           name="sohefldo"    description="Net Downward Heat Flux"                                   /> 
    427         <field ref="20d"               name="so20chgt"    description="Depth of 20C isotherm"                                    /> 
    428         <field ref="28d"               name="so28chgt"    description="Depth of 28C isotherm"                                    /> 
    429         <field ref="hc300"             name="sohtc300"    description="Heat content 300 m"                                       /> 
    430         <field ref="mldr10_1"          name="somxl010"    description="Mixed Layer Depth 0.01 ref.10m"                           /> 
    431         <field ref="mldr10_3"          name="somx3010"    description="Mixed Layer Depth dr=0.03 (ref.10m)"                      /> 
    432         <field ref="mldkz5"            name="somixhgt"    description="mixing layer depth (Turbocline)"                          /> 
    433         <field ref="emps"              name="sowaflcd"    description="concentration/dilution water flux"                        /> 
    434         <field ref="taum"              name="sotaumod"    description="Wind stress module"                                       /> 
    435         <field ref="wspd"              name="sowindsp"    description="Wind speed module at 10 m"                                /> 
    436         <field ref="ice_cover"         name="soicecov"    description="sea_ice_area_fraction"                                    /> 
    437         <field ref="BLT"               name="blt"         description="Barrier Layer Thickness"                    level="3"     /> 
    438         <field ref="mld|dt|"           name="mld_dt02"    description="Mixed Layer Depth |dt|=0.2 (ref.10m)"       level="3"     /> 
    439  
    440           </file>    
    441         <!-- 
    442 ..............................................    grid U   ................................................. 
    443           --> 
    444         <file id="ORCA2_1m_grid_U" description="ocean U grid variables : table 2.3 2.8" >                       <!-- grid U --> 
    445  
    446           <!--  table 2.3 --> 
    447             <field ref="uoce"         name="vozocrtx" description="sea_water_x_velocity"                                       /> 
    448             <field ref="u_masstr"     name="umo"      description="ocean_mass_x_transport"                           level="1" /> 
    449             <field ref="u_heattr"     name="hfx"      description="ocean_heat_x_transport"                           level="1" /> 
    450             <field ref="ueiv_heattr"  name="hfxba"    description="ocean_heat_x_transport_due_to_bolus_advection"    level="2" /> 
    451             <field ref="udiff_heattr" name="hfxdiff"  description="ocean_heat_x_transport_due_to_diffusion"          level="2" /> 
    452  
    453           <!-- table 2.8 --> 
    454             <field ref="utau"         name="sozotaux" description="surface_downward_x_stress"                        level="1" /> 
    455             <!-- tauucorr : surface_downward_x_stress_correction = 0 --> 
    456  
    457          <!--  IPSL variables --> 
    458         <field ref="uoce_eiv"         name="vozoeivu" description="EIV ocean current along i-axis"                             /> 
    459  
    460         </file>    
    461  
    462         <!-- 
    463 ..............................................    grid V   ................................................. 
    464           --> 
    465         <file id="ORCA2_1m_grid_V" description="ocean V grid variables : table 2.3 2.8" >                      <!-- grid V --> 
    466  
    467           <!-- table 2.3" --> 
    468             <field ref="voce"         name="vomecrty" description="sea_water_y_velocity"                                       /> 
    469             <field ref="v_masstr"     name="vmo"      description="ocean_mass_y_transport"                           level="1" /> 
    470             <field ref="v_heattr"     name="hfy"      description="ocean_heat_y_transport"                           level="1" /> 
    471             <field ref="veiv_heattr"  name="hfyba"    description="ocean_heat_y_transport_due_to_bolus_advection"    level="2" /> 
    472             <field ref="vdiff_heattr" name="hfydiff"  description="ocean_heat_y_transport_due_to_diffusion"          level="2" /> 
    473  
    474           <!-- table 2.8 --> 
    475             <field ref="vtau"         name="sometauy" description="surface_downward_y_stress"                        level="1" /> 
    476             <!-- tauvcorr : surface_downward_y_stress_correction = 0 --> 
    477  
    478          <!--  IPSL variables --> 
    479         <field ref="voce_eiv"         name="vomeeivv" description="EIV ocean current along j-axis"                             /> 
    480  
    481         </file>    
    482  
    483         <!-- 
    484 ..............................................    grid W   ................................................. 
    485           --> 
    486         <file id="ORCA2_1m_grid_W" description="ocean W grid variables table 2.3 2.9 2.10" >                   <!-- grid W --> 
    487            
    488           <!-- table 2.3 --> 
    489             <field ref="w_masstr"   name="wmo"   description="upward_ocean_mass_transport"                                    /> 
    490             <field ref="w_masstr2"  name="wmosq" description="square_pf_upward_ocean_mass_transport"                          /> 
    491            
    492           <!-- table 2.9 --> 
    493             <field ref="avt"     name="votkeavt" description="ocean_vertical_heat_diffusivity"                  level="2" /> 
    494             <field ref="avs"     name="voddmavs" description="ocean_vertical_salt_diffusivity"                  level="2" /> 
    495             <!-- difvtrbo : ocean_vertical_tracer_diffusivity_due_to_background : cte with time, see namelist parameters nn_avb and  nn_havtb --> 
    496             <field ref="av_tide" name="difvtrto" description="ocean_vertical_tracer_diffusivity_due_to_tides"   level="2" /> 
    497             <!-- tnpeo     : tendency_of_ocean_potential_energy_content                   : not available --> 
    498             <!-- tnpeot    : tendency_of_ocean_potential_energy_content_due_to_tides      : not available --> 
    499             <!-- tnpeotb   : tendency_of_ocean_potential_energy_content_due_to_background : not available --> 
    500             <field ref="avm"     name="difvmo"   description="ocean_vertical_momentum_diffusivity"              level="2" /> 
    501             <!-- difvmbo : ocean_vertical_momentum_diffusivity_due_to_background : cte with time, see namelist parameters nn_avb --> 
    502             <field ref="av_tide" name="difvmto"  description="ocean_vertical_momentum_diffusivity_due_to_tides" level="2" /> <!-- same as tracer --> 
    503             <!-- difvmfdo : ocean_vertical_momentum_diffusivity_due_to_form_drag : ??? --> 
    504             <!-- dispkevfo : ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction : not available --> 
    505  
    506           <!-- table 2.10 --> 
    507             <!-- if ln_traldf_lap    =  .true. --> 
    508             <field ref="aht2d_eiv" name="soleaeiw" description="ocean_tracer_bolus_laplacian_diffusivity"                     level="2" /> 
    509             <!-- diftrelo : ocean_tracer_epineutral_laplacian_diffusivity : cte with time, see ln_traldf_iso --> 
    510             <!-- diftrxylo : ocean_tracer_xy_laplacian_diffusivity : cte with time --> 
    511             <field ref="aht2d"     name="soleahtw" description="ocean_tracer_xy_laplacian_diffusivity"                        level="2" /> 
    512             <!-- if ln_traldf_bilap  =  .true. --> 
    513             <!-- field ref="diftrbbo" name="aht2d_eiv" description="ocean_tracer_bolus_biharmonic_diffusivity"                level="2" /--> 
    514             <!-- diftrebo : ocean_tracer_epineutral_biharmonic_diffusivity : cte with time, see ln_traldf_iso --> 
    515             <!-- diftrxybo : ocean_tracer_xy_biharmonic_diffusivity : cte with time --> 
    516             <!-- tnkebto : tendency_of_ocean_eddy_kinetic_energy_content_due_to_bolus_transport : not available --> 
    517             <!-- difmxylo : ocean_momentum_xy_laplacian_diffusivity : cte with time, see ln_dynldf_lap --> 
    518             <!-- difmxybo : ocean_momentum_xy_biharmonic_diffusivity : cte with time, see ln_dynldf_bilap --> 
    519             <!-- dispkexyfo : ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction : not available --> 
    520  
    521          <!--  IPSL variables --> 
    522         <field ref="woce"             name="vovecrtz" description="ocean vertical velocity"                          /> 
    523         <field ref="woce_eiv"         name="voveeivw" description="EIV ocean vertical velocity"                      /> 
    524  
    525  
    526           </file>    
    527  
    528         <!-- 
    529 ..............................................    icemod   ................................................. 
    530           --> 
    531         <file id="ORCA2_1m_icemod" description=" ice variables" >                                           <!-- icemod --> 
    532          <!-- table 2.2 --> 
    533           <field ref="ice_pres"                     /> 
    534           <field ref="ice_cover"    name="soicecov"   description="sea_ice_area_fraction"                            /> 
    535           <field ref="icethic_cea"  name="iicethic"   description="sea_ice_thickness"                                /> 
    536           <field ref="subl_ai_cea"  name="evap"       description="water_evaporation_flux"                           /> 
    537           <field ref="snowthic_cea" name="isnowthi"   description="surface_snow_thickness"                           /> 
    538           <field ref="sim"          name="sim"        description="sea_ice_and_surface_snow_amount"                  /> 
    539           <!-- snc : surface_snow_area_fraction : same as sic --> 
    540           <field ref="icealb_cea"   name="ialb"       description="Bare Sea Ice Albedo"                              /> 
    541           <!-- ssi : Sea Ice Salinity : cte = sice = 6, expect in baltic sea = 2 --> 
    542           <field ref="ist_cea"      name="tsice"      description="Surface Temperature of Sea Ice"                   /> 
    543           <!-- tsnint : Temperature at Interface Between Sea Ice and Snow : not available --> 
    544           <!-- pr : Surface Rainfall Rate into the Sea Ice Portion of the Grid Cell = 0 --> 
    545           <field ref="snow_ai_cea"  name="prsn"       description="Surface Snowfall Rate into the Sea Ice Portion of the Grid Cell" /> 
    546           <!-- ageice   : Age of Sea Ice : not available ??? --> 
    547           <!-- grFrazil : Frazil Sea Ice Growth (Leads) Rate ??? --> 
    548           <!-- grCongel : Congelation Sea Ice Growth Rate ??? --> 
    549           <field ref="licepro_cea"  name="grLateral"  description="Lateral Sea Ice Growth Rate"                      /> 
    550           <field ref="sntoice_cea"  name="snoToIce"   description="Snow-Ice Formation Rate"                          /> 
    551           <field ref="snowmel_cea"  name="snomelt"    description="Snow Melt Rate"                                   /> 
    552           <field ref="ticemel_cea"  name="tmelt"      description="Rate of Melt at Upper Surface of Sea Ice"         /> 
    553           <field ref="bicemel_cea"  name="bmelt"      description="Rate of Melt at Sea Ice Base"                     /> 
    554           <!-- hcice  : Sea Ice Total Heat Content ??? --> 
    555           <field ref="qsr_ai_cea"   name="rsdssi"     description="surface_downwelling_shortwave_flux_in_air"        /> 
    556           <!-- rsussi : surface_upwelling_shortwave_flux_in_air ??? --> 
    557           <!-- rldssi : surface_downwelling_longwave_flux_in_air ??? --> 
    558           <!-- rlussi : surface_upwelling_longwave_flux_in_air ??? --> 
    559           <!-- hfssi  : surface_upward_sensible_heat_flux ??? --> 
    560           <!-- hflssi : surface_upward_latent_heat_flux ??? --> 
    561           <field ref="subl_ai_cea"  name="sblsi"      description="surface_snow_and_ice_sublimation_flux"             /> 
    562           <field ref="u_imasstr"    name="transix"    description="Eastward Sea Ice Transport"                        /> 
    563           <field ref="v_imasstr"    name="transiy"    description="Northward Sea Ice Transport"                       /> 
    564           <field ref="fram_trans"   name="transifs"   description="Sea Ice Mass Transport Through Fram Strait"        /> 
    565           <field ref="utau_ice"     name="strairx"    description="Eastward Atmospheric Stress On Sea Ice"            /> 
    566           <field ref="vtau_ice"     name="strairy"    description="Northward Atmospheric Stress On Sea Ice"           /> 
    567           <!-- strocnx : Eastward Ocean Stress On Sea Ice ??? --> 
    568           <!-- strocny : Northward Ocean Stress On Sea Ice ??? --> 
    569           <!-- streng  : Compressive Sea Ice Strength ??? --> 
    570           <!-- divice  : Strain Rate Divergence of Sea Ice ??? --> 
    571           <!-- shrice  : Strain Rate Shear of Sea Ice ??? --> 
    572           <!-- ridgice : Sea Ice Ridging Rate ??? --> 
    573          <!-- IPSL variables --> 
    574           <field ref="ist_ipa"    name="iicetemp"   description="Ice surface temperature (ice presence average)"               /> 
    575           <field ref="uice_ipa"   name="iicevelu"   description="Ice velocity along i-axis at I-point (ice presence average)"  /> 
    576           <field ref="vice_ipa"   name="iicevelv"   description="Ice velocity along j-axis at I-point (ice presence average)"  /> 
    577           <field ref="ice_cover"  name="ileadfra"   description="sea_ice_area_fraction"                                        /> 
    578  
    579         </file>    
    580  
    581         <!-- 
    582 ..............................................    scalar   ................................................. 
    583           --> 
    584         <file id="ORCA2_1m_scalar" description="Ocean scalar" >                                <!-- scalar --> 
    585          <!-- table 2.2 --> 
    586           <field ref="masstot"    name="masso"    description="sea_water_mass"                                /> 
    587           <field ref="voltot"     name="volo"     description="sea_water_volume"                              /> 
    588           <field ref="sshtot"     name="zosga"    description="global_average_sea_level_change"               /> 
    589           <field ref="sshsteric"  name="zossga"   description="global_average_steric_sea_level_change"        /> 
    590           <field ref="sshthster"  name="zostoga"  description="global_average_thermosteric_sea_level_change"  /> 
    591           <field ref="temptot"    name="thetaoga" description="sea_water_potential_temperature"               /> 
    592           <field ref="saltot"     name="soga"     description="sea_water_salinity"                            /> 
    593         </file>    
    594  
    595         <!-- 
    596 ..............................................   ptrc_T   ................................................. 
    597           --> 
    598         <file id="ORCA2_1m_ptrc_T" description="Marine Biogeochemical 3-D Fields" >            <!-- ptrc_T --> 
    599           <field ref="DIC"      /> 
    600           <field ref="Alkalini" /> 
    601           <field ref="O2"       /> 
    602           <field ref="CaCO3"    /> 
    603           <field ref="PO4"      /> 
    604           <field ref="POC"      /> 
    605           <field ref="Si"       /> 
    606           <field ref="PHY"      /> 
    607           <field ref="ZOO"      /> 
    608           <field ref="DOC"      /> 
    609           <field ref="PHY2"     /> 
    610           <field ref="ZOO2"     /> 
    611           <field ref="BSi"      /> 
    612           <field ref="Fer"      /> 
    613           <field ref="BFe"      /> 
    614           <field ref="GOC"      /> 
    615           <field ref="SFe"      /> 
    616           <field ref="DFe"      /> 
    617           <field ref="DSi"      /> 
    618           <field ref="NFe"      /> 
    619           <field ref="NCHL"     /> 
    620           <field ref="DCHL"     /> 
    621           <field ref="NO3"      /> 
    622           <field ref="NH4"      /> 
    623        </file> 
    624  
    625         <!-- 
    626 ..............................................  diad_T   ................................................. 
    627           --> 
    628        <file id="ORCA2_1m_diad_T" description="Marine Biogeochemical Fields : Additional diagnostics" > 
    629           <field ref="PH"       /> 
    630           <field ref="CO3"      /> 
    631           <field ref="CO3sat"   /> 
    632           <field ref="PAR"      /> 
    633           <field ref="PPPHY"    /> 
    634           <field ref="PPPHY2"   /> 
    635           <field ref="PPNEWN"   /> 
    636           <field ref="PPNEWD"   /> 
    637           <field ref="PBSi"     /> 
    638           <field ref="PFeN"     /> 
    639           <field ref="PFeD"     /> 
    640           <field ref="PCAL"     /> 
    641           <field ref="DCAL"     /> 
    642           <field ref="GRAZ"     /> 
    643           <field ref="EPC100"   /> 
    644           <field ref="EPFE100"  /> 
    645           <field ref="EPSI100"  /> 
    646           <field ref="EPCAL100" /> 
    647           <field ref="Cflx"     /> 
    648           <field ref="Oflx"     /> 
    649           <field ref="Kg"       /> 
    650           <field ref="Dpco2"    /> 
    651           <field ref="Dpo2"     /> 
    652           <field ref="Heup"     /> 
    653           <field ref="Irondep"  /> 
    654           <field ref="Nfix"     /> 
    655        </file> 
    656  
    657         <!-- 
    658 ..............................................  dbio_T   ................................................. 
    659           --> 
    660        <file id="ORCA2_1m_dbio_T" description="Marine Biogeochemical Fields : Additional AR5 variables" enabled="<DBIO_ENABLE>" > 
    661           <field ref="PHYT"        /> 
    662           <field ref="ZOOT"        /> 
    663           <field ref="CHLT"        /> 
    664           <field ref="POCT"        /> 
    665           <field ref="PFET"        /> 
    666           <field ref="TPP"         /> 
    667           <field ref="TPNEW"       /> 
    668           <field ref="TPBFE"       /> 
    669           <field ref="EXPC"        /> 
    670           <field ref="EXPFE"       /> 
    671           <field ref="EXPSI"       /> 
    672           <field ref="EXPCAL"      /> 
    673           <field ref="INTDIC"      /> 
    674           <field ref="O2MIN"       /> 
    675           <field ref="ZO2MIN"      /> 
    676           <field ref="INTNFIX"     /> 
    677           <field ref="INTPPPHY"    /> 
    678           <field ref="INTPPPHY2"   /> 
    679           <field ref="INTPP"       /> 
    680           <field ref="INTPNEW"     /> 
    681           <field ref="INTPBFE"     /> 
    682           <field ref="INTPBSI"     /> 
    683           <field ref="INTPCAL"     /> 
    684        </file> 
    685  
    686       </group>    
    687  
    688       <!--  
    689 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
    690         --> 
    691     </file_definition> 
    692      
    693     <!--  
    694 ============================================================================================================ 
    695 =                                           grid definition                                                = 
    696 =                                            DO NOT CHANGE                                                 = 
    697 ============================================================================================================ 
    698     --> 
    699      
    700     <axis_definition>   
    701       <axis id="deptht" description="Vertical T levels" unit="m" positive=".false." /> 
    702       <axis id="depthu" description="Vertical U levels" unit="m" positive=".false." /> 
    703       <axis id="depthv" description="Vertical V levels" unit="m" positive=".false." /> 
    704       <axis id="depthw" description="Vertical W levels" unit="m" positive=".false." /> 
    705       <axis id="none" description="axe non defini" unit="none" size="1" /> 
    706     </axis_definition>  
    707      
    708     <grid_definition> 
    709       <grid id="grid_T" description="grid T" /> 
    710       <grid id="grid_U" description="grid U" /> 
    711       <grid id="grid_V" description="grid V" /> 
    712       <grid id="grid_W" description="grid W" /> 
    713       <grid id="scalarpoint" description="scalar" > 
    714         <zoom id="1point" ibegin="1" jbegin="1" ni="1" nj="1" /> 
    715       </grid> 
    716     </grid_definition>     
    717      
     246  </context> 
     247   
     248 
     249  <context id="xios"> 
     250 
     251      <variable_definition> 
     252         
     253     <!--  
     254        We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 
     255--> 
     256          <variable id="buffer_size"               type="integer">50000000</variable> 
     257          <variable id="buffer_server_factor_size" type="integer">2</variable> 
     258          <variable id="info_level"                type="integer">0</variable> 
     259          <variable id="using_server"              type="boolean">false</variable> 
     260          <variable id="using_oasis"               type="boolean">false</variable> 
     261          <variable id="oasis_codes_id"            type="string" >oceanx</variable> 
     262         
     263      </variable_definition> 
     264                
    718265  </context> 
    719266   
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2

    r970 r2061  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    22!! NEMO/OPA  :  1 - run manager      (namrun) 
    3 !! namelists    2 - Domain           (namzgr, namzgr_sco, namdom) 
     3!! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
    5 !!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 
     5!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
     6!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    67!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
    78!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
     
    910!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    1011!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
    11 !!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
    12 !!              9 - miscellaneous    (namsol, nammpp, nammpp_dyndist, namctl) 
     12!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namptr, namhsb) 
     13!!             10 - miscellaneous    (namsol, nammpp, nammpp_dyndist, namctl) 
     14!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1315!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    14 !  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. 
    1516 
    1617!!====================================================================== 
    1718!!                   ***  Run management namelists  *** 
    1819!!====================================================================== 
    19 !!   namrun            parameters of the run 
    20 !!====================================================================== 
    21  
     20!!   namrun        parameters of the run 
     21!!====================================================================== 
     22! 
    2223!----------------------------------------------------------------------- 
    2324&namrun        !   parameters of the run 
    2425!----------------------------------------------------------------------- 
    25    nn_no       =       0   !  job number 
    26    cn_exp      =  "ORCA2"  !  AUTO - experience name  
     26   nn_no       =       0   !  job number (no more used...) 
     27   cn_exp      =  "ORCA2"  !  AUTO - experience name 
    2728   nn_it000    =       1   !  AUTO - first time step 
    28    nn_itend    =    5475   !  AUTO - last  time step(std 5475) 
    29    nn_date0    =  010101   !  AUTO - initial calendar date yymmdd (used if  nn_rstctl=1) 
     29   nn_itend    =    5475   !  AUTO - last  time step (std 5475) 
     30   nn_date0    =  010101   !  AUTO - date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    3031   nn_leapy    =       0   !  AUTO - Leap year calendar (1) or not (0) 
     32   ln_rstart   = .false.   !  AUTO - start from rest (F) or from a restart file (T) 
     33   nn_rstctl   =       0   !  AUTO - restart control => activated only if ln_rstart = T 
     34                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
     35                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     36                           !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
     37   cn_ocerst_in  = "restartopa"   !  suffix of ocean restart name (input) 
     38   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    3139   nn_istate   =       0   !  output the initial state (1) or not (0) 
    3240   nn_stock    =    5475   !  AUTO - frequency of creation of a restart file (modulo referenced to 1) 
    3341   nn_write    =    5475   !  AUTO - frequency of write in the output file   (modulo referenced to nn_it000) 
    3442   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    35    ln_mskland  = .true.    !  mask land points in NetCDF outputs (costly: + ~15%) 
     43   ln_mskland  = .true.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    3644   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
    37    nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (working only with iom_nf90 routines) 
    38    ln_rstart   = .false.   !  AUTO - start from rest (F) or from a restart file (T) 
    39    nn_rstctl   =       0   !  AUTO - restart control = 0 nn_it000 is not compared to the restart file value 
    40                            !                  = 1 use nn_date0 in namelist (not the value in the restart file) 
    41                            !                  = 2 calendar parameters read in the restart file 
    42    cn_ocerst_in  = "restartopa"!  suffix of ocean restart name (input) 
    43    cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    44 / 
     45   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
     46/ 
     47 
    4548!!====================================================================== 
    4649!!                      ***  Domain namelists  *** 
     
    4952!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    5053!!   namdom       space and time domain (bathymetry, mesh, timestep) 
    51 !!====================================================================== 
    52  
     54!!   namtsd       data: temperature & salinity 
     55!!====================================================================== 
     56! 
    5357!----------------------------------------------------------------------- 
    5458&namzgr        !   vertical coordinate 
     
    6165&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
    6266!----------------------------------------------------------------------- 
    63    rn_sbot_min =  300.     !  minimum depth of s-bottom surface (>0) (m) 
    64    rn_sbot_max = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     67   ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
     68   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
     69   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
     70                           !  stretching coefficients for all functions 
     71   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
     72   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     73   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
     74                        !!!!!!!  Envelop bathymetry 
     75   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
     76                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    6577   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    66    rn_thetb    =    0.75   !  bottom control parameter  (0<=thetb<= 1) 
    67    rn_rmax     =    0.15   !  maximum cut-off r-value allowed (0<r_max<1) 
    68    ln_s_sigma  = .false.   !  hybrid s-sigma coordinates 
    69    rn_bb       =    0.8    !  stretching with s-sigma 
    70    rn_hc       =  150.0    !  critical depth with s-sigma  
     78   rn_bb       =    0.8    !  stretching with SH94 s-sigma    
     79                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
     80   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     81   rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
     82   rn_zs       =    1.0    !  depth of surface grid box 
     83                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
     84   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
     85   rn_zb_b     =   -0.2    !  offset for calculating Zb 
     86                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
     87   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
    7188/ 
    7289!----------------------------------------------------------------------- 
    7390&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    7491!----------------------------------------------------------------------- 
    75    nn_bathy    =    1      !  compute (=0) or read(=1) the bathymetry file 
    76    nn_closea   =    1      !  closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 
    77    nn_msh      =    0      !  AUTO create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 
    78    rn_e3zps_min=   20.     !  the thickness of the partial step is set larger than the minimum 
    79    rn_e3zps_rat=    0.1    !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1) 
     92   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
     93   nn_closea    =   0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
     94   nn_msh      =    0      !  AUTO - create (=1) a mesh file or not (=0) 
     95   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
     96   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
     97   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    8098                           ! 
    81    rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nacc=0)   ==> 5760 
    82    nn_baro     =   64      !  number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts") 
     99   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
     100   nn_baro     =   64      !  number of barotropic time step           ("key_dynspg_ts") 
    83101   rn_atfp     =    0.1    !  asselin time filter parameter 
    84102   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    85103                                 !                          =0, not used, rdt = rdttra 
    86    rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nacc=1) 
    87    rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nacc=1) 
    88    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nacc=1) 
     104   rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1) 
     105   rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
     106   rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
     107/ 
     108!----------------------------------------------------------------------- 
     109&namtsd    !   data : Temperature  & Salinity 
     110!----------------------------------------------------------------------- 
     111!          ! file name ! frequency (hours)    ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     112!          !           !  (if <0  months)     !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     113   sn_tem  = 'data_1m_potential_temperature_nomask', -1,'votemper',  .true.  , .true., 'yearly'   , ' '      , ' ' 
     114   sn_sal  = 'data_1m_salinity_nomask'             , -1,'vosaline',  .true.  , .true., 'yearly'   , ''       , ' ' 
     115   ! 
     116   cn_dir        = './'     !  root directory for the location of the runoff files 
     117   ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F) 
     118   ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
    89119/ 
    90120!!====================================================================== 
    91121!!            ***  Surface Boundary Condition namelists  *** 
    92122!!====================================================================== 
    93 !!   namsbc        surface boundary condition 
    94 !!   namsbc_ana    analytical         formulation 
    95 !!   namsbc_flx    flux               formulation 
    96 !!   namsbc_clio   CLIO bulk formulea formulation 
    97 !!   namsbc_core   CORE bulk formulea formulation 
    98 !!   namsbc_cpl    CouPLed            formulation                       ("key_coupled") 
    99 !!   namtra_qsr    penetrative solar radiation 
    100 !!   namsbc_rnf    river runoffs 
    101 !!   namsbc_ssr    sea surface restoring term (for T and/or S) 
    102 !!   namsbc_alb    albedo parameters 
    103 !!====================================================================== 
    104  
     123!!   namsbc          surface boundary condition 
     124!!   namsbc_ana      analytical         formulation 
     125!!   namsbc_flx      flux               formulation 
     126!!   namsbc_clio     CLIO bulk formulae formulation 
     127!!   namsbc_core     CORE bulk formulae formulation 
     128!!   namsbc_mfs      MFS  bulk formulae formulation 
     129!!   namsbc_cpl      CouPLed            formulation                     ("key_coupled") 
     130!!   namtra_qsr      penetrative solar radiation 
     131!!   namsbc_rnf      river runoffs 
     132!!   namsbc_apr      Atmospheric Pressure 
     133!!   namsbc_ssr      sea surface restoring term (for T and/or S) 
     134!!   namsbc_alb      albedo parameters 
     135!!====================================================================== 
     136! 
    105137!----------------------------------------------------------------------- 
    106138&namsbc        !   Surface Boundary Condition (surface module) 
    107139!----------------------------------------------------------------------- 
    108    nn_fsbc     = 5         !  frequency of surface boundary condition computation  
    109                            !               (= the frequency of sea-ice model call) 
    110    ln_ana      = .false.   !  analytical formulation (T => fill namsbc_ana )  
    111    ln_flx      = .false.   !  flux formulation       (T => fill namsbc_flx ) 
    112    ln_blk_clio = .false.   !  CLIO bulk formulation  (T => fill namsbc_clio)  
    113    ln_blk_core = .false.   !  CORE bulk formulation  (T => fill namsbc_core)  
    114    ln_cpl      = .true.    !  Coupled formulation    (T => fill namsbc_cpl ) 
     140   nn_fsbc     = 5         !  frequency of surface boundary condition computation 
     141                           !     (also = the frequency of sea-ice model call) 
     142   ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana ) 
     143   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
     144   ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio) 
     145   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     146   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
     147   ln_cpl      = .true.    !  Coupled formulation                       (T => fill namsbc_cpl ) 
     148   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    115149   nn_ice      = 2         !  =0 no ice boundary condition   , 
    116150                           !  =1 use observed ice-cover      , 
    117                            !  =2 ice-model used                             ("key_lim3" or "key_lim2) 
    118    nn_ico_cpl  = 0         !  ice-ocean coupling : =0 each nn_fsbc  
    119                            !                       =1 stresses recomputed each ocean time step ("key_lim3" only) 
    120                            !                       =2 combination of 0 and 1 cases             ("key_lim3" only) 
    121    ln_dm2dc    = .false.   !  daily mean to diurnal cycle short wave (qsr) 
    122    ln_rnf      = .true.    !  runoffs (T => fill namsbc_rnf) 
    123    ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 
    124    nn_fwb      = 0         !  FreshWater Budget: =0 unchecked  
    125                            !                     =1 global mean of e-p-r set to zero at each time step  
    126                            !                     =2 annual global mean of e-p-r set to zero 
    127                            !                     =3 global emp set to zero and spread out over erp area 
     151                           !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     152   nn_ice_embd = 0         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
     153                           !  =1 levitating ice with mass and salt exchange but no presure effect 
     154                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
     155   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
     156   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     157   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
     158   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
     159                           !     =1 global mean of e-p-r set to zero at each time step 
     160                           !     =2 annual global mean of e-p-r set to zero 
     161                           !     =3 global emp set to zero and spread out over erp area 
     162   ln_wave = .false.       !  Activate coupling with wave (either Stokes Drift or Drag coefficient, or both)  (T => fill namsbc_wave) 
     163   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
     164   ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave) 
    128165/ 
    129166!----------------------------------------------------------------------- 
     
    140177&namsbc_flx    !   surface boundary condition : flux formulation 
    141178!----------------------------------------------------------------------- 
    142 !              !   file name  ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    143 !              !              !  (if <0  months)  !    name    !    (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    144    sn_utau     = 'utau'       ,        24         ,    'utau'  ,     .false.    , .false. , 'yearly'  , ''       , '' 
    145    sn_vtau     = 'vtau'       ,        24         ,    'vtau'  ,     .false.    , .false. , 'yearly'  , ''       , '' 
    146    sn_qtot     = 'qtot'       ,        24         ,    'qtot'  ,     .false.    , .false. , 'yearly'  , ''       , '' 
    147    sn_qsr      = 'qsr'        ,        24         ,    'qsr'   ,     .false.    , .false. , 'yearly'  , ''       , '' 
    148    sn_emp      = 'emp'        ,        24         ,    'emp'   ,     .false.    , .false. , 'yearly'  , ''       , '' 
    149 ! 
     179!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     180!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     181   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , '' 
     182   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , '' 
     183   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , '' 
     184   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , '' 
     185   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , '' 
     186 
    150187   cn_dir      = './'      !  root directory for the location of the flux files 
    151 /       
    152 !----------------------------------------------------------------------- 
    153 &namsbc_clio   !   namsbc_clio  CLIO bulk formulea 
    154 !----------------------------------------------------------------------- 
    155 !              !   file name  ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    156 !              !              !  (if <0  months)  !    name    !    (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    157    sn_utau     = 'taux_1m'    ,       -1          , 'sozotaux' ,    .true.      , .true. , 'yearly'  , ''       , '' 
    158    sn_vtau     = 'tauy_1m'    ,       -1          , 'sometauy' ,    .true.      , .true. , 'yearly'  , ''       , '' 
    159    sn_wndm     = 'flx'        ,       -1          , 'socliowi' ,    .true.      , .true. , 'yearly'  , ''       , '' 
    160    sn_tair     = 'flx'        ,       -1          , 'socliot2' ,    .true.      , .true. , 'yearly'  , ''       , '' 
    161    sn_humi     = 'flx'        ,       -1          , 'socliohu' ,    .true.      , .true. , 'yearly'  , ''       , '' 
    162    sn_ccov     = 'flx'        ,       -1          , 'socliocl' ,    .false.     , .true. , 'yearly'  , ''       , '' 
    163    sn_prec     = 'flx'        ,       -1          , 'socliopl' ,    .false.     , .true. , 'yearly'  , ''       , '' 
    164 ! 
     188/ 
     189!----------------------------------------------------------------------- 
     190&namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
     191!----------------------------------------------------------------------- 
     192!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     193!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     194   sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , '' 
     195   sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , '' 
     196   sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , '' 
     197   sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , '' 
     198   sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , '' 
     199   sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , '' 
     200   sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , '' 
     201 
    165202   cn_dir      = './'      !  root directory for the location of the bulk files are 
    166203/ 
    167204!----------------------------------------------------------------------- 
    168 &namsbc_core   !   namsbc_core  CORE bulk formulea 
    169 !----------------------------------------------------------------------- 
    170 !              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights              ! rotation ! 
    171 !              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename            ! pairing  ! 
    172    sn_wndi     = 'u10_core'   ,       -1          , 'u10'      ,    .true.      , .true.  , 'yearly'  ,'bicubic_weights_orca2.nc' , 'U1' 
    173    sn_wndj     = 'v10_core'   ,       -1          , 'v10'      ,    .true.      , .true.  , 'yearly'  ,'bicubic_weights_orca2.nc' , 'V1' 
    174    sn_qsr      = 'qsw_core'   ,       -1          , 'swdn'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    175    sn_qlw      = 'qlw_core'   ,       -1          , 'lwdn'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    176    sn_tair     = 't2_core'    ,       -1          , 't2'       ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    177    sn_humi     = 'q2_core'    ,       -1          , 'q2'       ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    178    sn_prec     = 'precip_core',       -1          , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    179    sn_snow     = 'snow_core'  ,       -1          , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    180    sn_tdif     = 'taudif_core',       24          , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', '' 
    181 ! 
     205&namsbc_core   !   namsbc_core  CORE bulk formulae 
     206!----------------------------------------------------------------------- 
     207!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     208!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     209   sn_wndi     = 'u_10.15JUNE2009_orca2'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , ''       , 'Uwnd' 
     210   sn_wndj     = 'v_10.15JUNE2009_orca2'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , ''       , 'Vwnd' 
     211   sn_qsr      = 'ncar_rad.15JUNE2009_orca2'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
     212   sn_qlw      = 'ncar_rad.15JUNE2009_orca2'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
     213   sn_tair     = 't_10.15JUNE2009_orca2'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
     214   sn_humi     = 'q_10.15JUNE2009_orca2'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
     215   sn_prec     = 'ncar_precip.15JUNE2009_orca2',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , ''       , '' 
     216   sn_snow     = 'ncar_precip.15JUNE2009_orca2',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , ''       , '' 
     217   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , ''       , '' 
     218 
    182219   cn_dir      = './'      !  root directory for the location of the bulk files 
    183    ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    184    ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ? 
     220   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     221   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    185222   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    186223/ 
    187224!----------------------------------------------------------------------- 
    188 &namsbc_cpl    !   coupled ocean/atmosphere model                        ("key_coupled") 
    189 !----------------------------------------------------------------------- 
    190                                        ! send 
    191 cn_snd_temperature= 'weighted oce and ice'  ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 
    192 cn_snd_albedo     = 'weighted ice'          ! 'none' 'weighted ice' 'mixed oce-ice' 
    193 cn_snd_thickness  = 'none'                  ! 'none' 'weighted ice and snow' 
    194 cn_snd_crt_nature = 'mixed oce-ice'         ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 
    195 cn_snd_crt_refere = 'cartesian'             ! 'spherical' 'cartesian' 
    196 cn_snd_crt_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
    197 cn_snd_crt_grid   = 'T'                     ! 'T' 
    198                                        ! receive 
    199 cn_rcv_w10m       = 'coupled'               ! 'none' 'coupled' 
    200 cn_rcv_taumod     = 'none'                  ! 'none' 'coupled' 
    201 cn_rcv_tau_nature = 'mixed oce-ice'         ! 'oce only' 'oce and ice' 'mixed oce-ice' 
    202 cn_rcv_tau_refere = 'cartesian'             ! 'spherical' 'cartesian' 
    203 cn_rcv_tau_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid' 
    204 cn_rcv_tau_grid   = 'U,V'                   ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 
    205 cn_rcv_dqnsdt     = 'coupled'               ! 'none' 'coupled' 
    206 cn_rcv_qsr        = 'conservative'          ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    207 cn_rcv_qns        = 'conservative'          ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    208 cn_rcv_emp        = 'conservative'          ! 'conservative' 'oce and ice' 'mixed oce-ice' 
    209 cn_rcv_rnf        = 'coupled'               ! 'coupled' 'climato' 'mixed' 
    210 cn_rcv_cal        = 'coupled'               ! 'none' 'coupled' 
    211 / 
    212 !----------------------------------------------------------------------- 
    213 &namsbc_cpl_co2   !   coupled ocean/biogeo/atmosphere model            ("key_cpl_carbon_cycle") 
    214 !----------------------------------------------------------------------- 
    215 cn_snd_co2        = 'coupled'         ! send    :  'none' 'coupled' 
    216 cn_rcv_co2        = 'coupled'         ! receive : 'none' 'coupled' 
     225&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
     226!----------------------------------------------------------------------- 
     227!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     228!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     229   sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
     230   sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
     231   sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', '' 
     232   sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
     233   sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
     234   sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', '' 
     235   sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , '' 
     236 
     237   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files 
     238/ 
     239!----------------------------------------------------------------------- 
     240&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_coupled") 
     241!----------------------------------------------------------------------- 
     242!                    !     description       !  multiple  !    vector   !      vector          ! vector ! 
     243!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
     244! send 
     245sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     246sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     247sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     248sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     249sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     250! receive 
     251sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     252sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     253sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     254sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     255sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     256sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     257sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     258sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     259sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     260sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    217261/ 
    218262!----------------------------------------------------------------------- 
    219263&namtra_qsr    !   penetrative solar radiation 
    220264!----------------------------------------------------------------------- 
    221 !              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    222 !              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    223    sn_chl      = 'chlorophyll',        -1         , 'CHLA'     ,    .true.      , .true. , 'yearly'  , ''       , '' 
    224   
     265!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     266!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     267   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , '' 
     268 
    225269   cn_dir      = './'      !  root directory for the location of the runoff files 
    226270   ln_traqsr   = .true.    !  Light penetration (T) or not (F) 
     
    228272   ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    229273   ln_qsr_bio  = .true.    !  bio-model light penetration 
    230    nn_chldta   =      0    !  RGB : Chl 2D data (=1), 3D data (=2) or cst value (=0) 
     274   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0) 
    231275   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1) 
    232276   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction 
    233277   rn_si1      =   23.0    !  2 bands: longest depth of extinction 
    234    rn_si2      =   62.0    !  3 bands: longest depth of extinction (for blue waveband & 0.01 mg/m2 Chl) 
    235278/ 
    236279!----------------------------------------------------------------------- 
    237280&namsbc_rnf    !   runoffs namelist surface boundary condition 
    238281!----------------------------------------------------------------------- 
    239 !              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    240 !              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
    241    sn_rnf      = 'runoff_1m_nomask' ,  -1         , 'sorunoff' ,    .true.      , .true.  , 'yearly'  , ''       , '' 
    242    sn_cnf      = 'runoff_1m_nomask' ,   0         , 'socoefr'  ,    .false.     , .true.  , 'yearly'  , ''       , '' 
    243   
     282!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     283!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     284   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , '' 
     285   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , '' 
     286   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , '' 
     287   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , '' 
     288   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , '' 
     289 
    244290   cn_dir       = './'      !  root directory for the location of the runoff files 
    245    ln_rnf_emp   =   .true.  !  runoffs included into precipitation field (T) or into a file (F) 
    246    ln_rnf_mouth =   .true.  !  specific treatment at rivers mouths 
    247    rn_hrnf      =   15.e0   !  depth over which enhanced vertical mixing is used 
     291   ln_rnf_emp   = .false.   !  runoffs included into precipitation field (T) or into a file (F) 
     292   ln_rnf_mouth = .true.    !  specific treatment at rivers mouths 
     293   rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used 
    248294   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] 
    249295   rn_rfact     =   1.e0    !  multiplicative factor for runoff 
     296   ln_rnf_depth = .false.   !  read in depth information for runoff 
     297   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
     298   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
     299/ 
     300!----------------------------------------------------------------------- 
     301&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
     302!----------------------------------------------------------------------- 
     303!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     304!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     305   sn_apr      = 'patm'    ,         -1        ,'somslpre',    .true.      , .true.  , 'yearly'  ,  ''      ,   '' 
     306 
     307   cn_dir      = './'       !  root directory for the location of the bulk files 
     308   rn_pref     = 101000._wp !  reference atmospheric pressure   [N/m2]/ 
     309   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F) 
     310   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data 
    250311/ 
    251312!----------------------------------------------------------------------- 
    252313&namsbc_ssr    !   surface boundary condition : sea surface restoring 
    253314!----------------------------------------------------------------------- 
    254 !              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    255 !              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    256    sn_sst      = 'sst_data'   ,        24         ,  'sst'     ,     .false.    , .false. , 'yearly'  , ''       , '' 
    257    sn_sss      = 'sss_data'   ,        -1         ,  'sss'     ,     .true.     , .true. , 'yearly'  , ''       , '' 
    258      
     315!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     316!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     317   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , '' 
     318   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , '' 
     319 
    259320   cn_dir      = './'      !  root directory for the location of the runoff files 
    260321   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0) 
    261    nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)  
     322   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2) 
    262323                           !  or to SSS only (=1) or no damping term (=0) 
    263324   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K] 
    264    rn_deds     =   -27.7   !  magnitude of the damping on salinity   [mm/day/psu] 
     325   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    265326   ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2) 
    266327   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    267 /       
     328/ 
    268329!----------------------------------------------------------------------- 
    269330&namsbc_alb    !   albedo parameters 
    270331!----------------------------------------------------------------------- 
    271    rn_cloud    =    0.0    !  cloud correction to snow and ice albedo  
    272    rn_albice   =    0.5    !  albedo of melting ice in the arctic and antarctic 
     332   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo 
     333   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic 
    273334   rn_alphd    =    0.80   !  coefficients for linear interpolation used to 
    274    rn_alphc    =    0.65   !  compute albedo between two extremes values  
     335   rn_alphc    =    0.65   !  compute albedo between two extremes values 
    275336   rn_alphdi   =    0.72   !  (Pyane, 1972) 
     337/ 
     338!----------------------------------------------------------------------- 
     339&namberg       !   iceberg parameters 
     340!----------------------------------------------------------------------- 
     341      ln_icebergs              = .false. 
     342      ln_bergdia               = .true.               ! Calculate budgets 
     343      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     344      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     345      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     346                                                      ! Initial mass required for an iceberg of each class 
     347      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     348                                                      ! Proportion of calving mass to apportion to each class   
     349      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     350                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     351                                                      ! i.e. number of icebergs represented at a point          
     352      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     353                                                      ! thickness of newly calved bergs (m) 
     354      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     355      rn_rho_bergs             = 850.                 ! Density of icebergs 
     356      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     357      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     358      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     359      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     360      ln_passive_mode          = .false.              ! iceberg - ocean decoupling    
     361      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     362                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     363      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     364      rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
     365 
     366               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     367               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     368      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     369    
     370      cn_dir = './'  
    276371/ 
    277372 
     
    282377!!   namcla        cross land advection 
    283378!!   namobc        open boundaries parameters                           ("key_obc") 
    284 !!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
     379!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    285380!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    286381!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides") 
    287382!!====================================================================== 
    288  
     383! 
    289384!----------------------------------------------------------------------- 
    290385&namlbc        !   lateral momentum boundary condition 
     
    292387   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    293388                           !  free slip  !   partial slip  !   no slip   ! strong slip 
     389   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs. 
    294390/ 
    295391!----------------------------------------------------------------------- 
     
    301397&namobc        !   open boundaries parameters                           ("key_obc") 
    302398!----------------------------------------------------------------------- 
    303     ln_obc_clim= .false.   !  climatological obc data files (T) or not (F) 
    304     ln_vol_cst = .true.    !  impose the total volume conservation (T) or not (F) 
    305     ln_obc_fla = .false.   !  Flather open boundary condition  
    306     nn_obcdta  =    1      !  = 0 the obc data are equal to the initial state 
     399   ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
     400   ln_vol_cst = .true.    !  impose the total volume conservation (T) or not (F) 
     401   ln_obc_fla  = .false.   !  Flather open boundary condition 
     402   nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    307403                           !  = 1 the obc data are read in 'obc.dta' files 
    308     cn_obcdta  = 'annual'  !  set to annual if obc datafile hold 1 year of data 
     404   cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    309405                           !  set to monthly if obc datafile hold 1 month of data 
    310     rn_dpein   =    1.     !  damping time scale for inflow at east  open boundary 
    311     rn_dpwin   =    1.     !     -           -         -       west    -      - 
    312     rn_dpnin   =    1.     !     -           -         -       north   -      - 
    313     rn_dpsin   =    1.     !     -           -         -       south   -      - 
    314     rn_dpeob   = 3000.     !  time relaxation (days) for the east  open boundary 
    315     rn_dpwob   =   15.     !     -           -         -     west    -      - 
    316     rn_dpnob   = 3000.     !     -           -         -     north   -      - 
    317     rn_dpsob   =   15.     !     -           -         -     south   -      - 
    318     rn_volemp  =  1.       !  = 0 the total volume change with the surface flux (E-P-R) 
     406   rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
     407   rn_dpwin    =    1.     !     -           -         -       west    -      - 
     408   rn_dpnin    =    1.     !     -           -         -       north   -      - 
     409   rn_dpsin    =    1.     !     -           -         -       south   -      - 
     410   rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
     411   rn_dpwob    =   15.     !     -           -         -     west    -      - 
     412   rn_dpnob    = 3000.     !     -           -         -     north   -      - 
     413   rn_dpsob    =   15.     !     -           -         -     south   -      - 
     414   rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    319415                           !  = 1 the total volume remains constant 
    320416/ 
     
    322418&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    323419!----------------------------------------------------------------------- 
    324     nn_cln_update = 3      !  baroclinic update frequency 
    325     ln_spc_dyn    = .true. !  use 0 as special value for dynamics 
    326     rn_sponge_tra = 2880.  !  coefficient for tracer   sponge layer [s] 
    327     rn_sponge_dyn = 2880.  !  coefficient for dynamics sponge layer [s] 
     420   nn_cln_update =    3    !  baroclinic update frequency 
     421   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
     422   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
     423   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
     424/ 
     425!----------------------------------------------------------------------- 
     426&nam_tide      !   tide parameters (#ifdef key_tide) 
     427!----------------------------------------------------------------------- 
     428   ln_tide_pot   = .true.   !  use tidal potential forcing 
     429   nb_harmo      =    11    !  number of constituents used 
     430   clname(1)     =   'M2'   !  name of constituent 
     431   clname(2)     =   'S2' 
     432   clname(3)     =   'N2' 
     433   clname(4)     =   'K1' 
     434   clname(5)     =   'O1' 
     435   clname(6)     =   'Q1' 
     436   clname(7)     =   'M4' 
     437   clname(8)     =   'K2' 
     438   clname(9)     =   'P1' 
     439   clname(10)    =   'Mf' 
     440   clname(11)    =   'Mm' 
    328441/ 
    329442!----------------------------------------------------------------------- 
    330443&nambdy        !  unstructured open boundaries                          ("key_bdy") 
    331444!----------------------------------------------------------------------- 
    332     filbdy_mask    =  ''                  !  name of mask file (if ln_bdy_mask=.TRUE.) 
    333     filbdy_data_T  = 'bdydata_grid_T.nc'  !  name of data file (T-points) 
    334     filbdy_data_U  = 'bdydata_grid_U.nc'  !  name of data file (U-points) 
    335     filbdy_data_V  = 'bdydata_grid_V.nc'  !  name of data file (V-points) 
    336     ln_bdy_clim    = .false.              !  contain 1 (T) or 12 (F) time dumps and be cyclic 
    337     ln_bdy_vol     = .true.               !  total volume correction (see volbdy parameter) 
    338     ln_bdy_mask    = .false.              !  boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) 
    339     ln_bdy_tides   = .true.               !  Apply tidal harmonic forcing with Flather condition 
    340     ln_bdy_dyn_fla = .true.               !  Apply Flather condition to velocities 
    341     ln_bdy_tra_frs = .false.              !  Apply FRS condition to temperature and salinity  
    342     ln_bdy_dyn_frs = .false.              !  Apply FRS condition to velocities 
    343     nbdy_dta       =  1                   !  = 0, bdy data are equal to the initial state 
     445    nb_bdy = 1                            !  number of open boundary sets 
     446    ln_coords_file = .true.               !  =T : read bdy coordinates from file 
     447    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
     448    ln_mask_file = .false.                !  =T : read mask from file 
     449    cn_mask_file = ''                     !  name of mask file (if ln_mask_file=.TRUE.) 
     450    nn_dyn2d      =  2                    !  boundary conditions for barotropic fields 
     451    nn_dyn2d_dta  =  3                    !  = 0, bdy data are equal to the initial state 
    344452                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    345     nb_rimwidth    = 9                    !  width of the relaxation zone 
    346     volbdy         = 0                    !  = 0, the total water flux across open boundaries is zero 
    347                                           !  = 1, the total volume of the system is conserved 
    348 / 
    349 !----------------------------------------------------------------------- 
    350 &nambdy_tide     ! tidal forcing at unstructured boundaries               
    351 !----------------------------------------------------------------------- 
    352     filtide      = 'bdytide_'           !  file name root of tidal forcing files 
    353     tide_cpt     = 'M2','S1'            !  names of tidal components used 
    354     tide_speed   = 28.984106, 15.000001 !  phase speeds of tidal components (deg/hour) 
    355     ln_tide_date = .false.              !  adjust tidal harmonics for start date of run 
    356 / 
    357  
     453                                          !  = 2, use tidal harmonic forcing data from files 
     454                                          !  = 3, use external data AND tidal harmonic forcing 
     455    nn_dyn3d      =  0                    !  boundary conditions for baroclinic velocities 
     456    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
     457                           !  = 1, bdy data are read in 'bdydata   .nc' files 
     458    nn_tra        =  1                    !  boundary conditions for T and S 
     459    nn_tra_dta    =  1                    !  = 0, bdy data are equal to the initial state 
     460                           !  = 1, bdy data are read in 'bdydata   .nc' files 
     461    nn_rimwidth  = 10                      !  width of the relaxation zone 
     462    ln_vol     = .false.                  !  total volume correction (see nn_volctl parameter) 
     463    nn_volctl  = 1                        !  = 0, the total water flux across open boundaries is zero 
     464/ 
     465!----------------------------------------------------------------------- 
     466&nambdy_dta      !  open boundaries - external data           ("key_bdy") 
     467!----------------------------------------------------------------------- 
     468!              !   file name    ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     469!              !                !  (if <0  months)  !    name    !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     470   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     471   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     472   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     473   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     474   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     475   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     476   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     477   cn_dir  =    'bdydta/' 
     478   ln_full_vel = .false. 
     479/ 
     480!----------------------------------------------------------------------- 
     481&nambdy_tide     ! tidal forcing at open boundaries 
     482!----------------------------------------------------------------------- 
     483   filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     484    tide_cpt(1)   ='Q1'  !  names of tidal components used 
     485    tide_cpt(2)   ='O1'  !  names of tidal components used 
     486    tide_cpt(3)   ='P1'  !  names of tidal components used 
     487    tide_cpt(4)   ='S1'  !  names of tidal components used 
     488    tide_cpt(5)   ='K1'  !  names of tidal components used 
     489    tide_cpt(6)   ='2N2' !  names of tidal components used 
     490    tide_cpt(7)   ='MU2' !  names of tidal components used 
     491    tide_cpt(8)   ='N2'  !  names of tidal components used 
     492    tide_cpt(9)   ='NU2' !  names of tidal components used 
     493    tide_cpt(10)   ='M2'  !  names of tidal components used 
     494    tide_cpt(11)   ='L2'  !  names of tidal components used 
     495    tide_cpt(12)   ='T2'  !  names of tidal components used 
     496    tide_cpt(13)   ='S2'  !  names of tidal components used 
     497    tide_cpt(14)   ='K2'  !  names of tidal components used 
     498    tide_cpt(15)   ='M4'  !  names of tidal components used 
     499    tide_speed(1)   = 13.398661 !  phase speeds of tidal components (deg/hour) 
     500    tide_speed(2)   = 13.943036 !  phase speeds of tidal components (deg/hour) 
     501    tide_speed(3)   = 14.958932 !  phase speeds of tidal components (deg/hour) 
     502    tide_speed(4)   = 15.000001 !  phase speeds of tidal components (deg/hour) 
     503    tide_speed(5)   = 15.041069 !  phase speeds of tidal components (deg/hour) 
     504    tide_speed(6)   = 27.895355 !  phase speeds of tidal components (deg/hour) 
     505    tide_speed(7)   = 27.968210 !  phase speeds of tidal components (deg/hour) 
     506    tide_speed(8)   = 28.439730 !  phase speeds of tidal components (deg/hour) 
     507    tide_speed(9)   = 28.512585 !  phase speeds of tidal components (deg/hour) 
     508    tide_speed(10)   = 28.984106 !  phase speeds of tidal components (deg/hour) 
     509    tide_speed(11)   = 29.528479 !  phase speeds of tidal components (deg/hour) 
     510    tide_speed(12)   = 29.958935 !  phase speeds of tidal components (deg/hour) 
     511    tide_speed(13)   = 30.000002 !  phase speeds of tidal components (deg/hour) 
     512    tide_speed(14)   = 30.082138 !  phase speeds of tidal components (deg/hour) 
     513    tide_speed(15)   = 57.968212 !  phase speeds of tidal components (deg/hour) 
     514    ln_tide_date = .true.               !  adjust tidal harmonics for start date of run 
     515/ 
    358516!!====================================================================== 
    359517!!                 ***  Bottom boundary condition  *** 
    360518!!====================================================================== 
    361519!!   nambfr        bottom friction 
    362 !!   nambbc        bottom temperature boundary condition                ("key_trabbc") 
    363 !!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
    364 !!====================================================================== 
    365  
     520!!   nambbc        bottom temperature boundary condition 
     521!!   nambbl        bottom boundary layer scheme                         ("key_trabbl") 
     522!!====================================================================== 
     523! 
    366524!----------------------------------------------------------------------- 
    367525&nambfr        !   bottom friction 
     
    371529   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case) 
    372530   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case) 
    373    rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m^2/s^2) 
    374    ln_bfr2d    =   .false. !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    375    rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d = .true.) 
     531   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
     532   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
     533   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     534   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
    376535/ 
    377536!----------------------------------------------------------------------- 
    378537&nambbc        !   bottom temperature boundary condition 
    379538!----------------------------------------------------------------------- 
    380    nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux  
     539   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
     540   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
    381541                           !     = 1 constant flux 
    382                            !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)  
     542                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    383543   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    384544/ 
     
    386546&nambbl        !   bottom boundary layer scheme 
    387547!----------------------------------------------------------------------- 
    388 !                          !  diffusive bbl                             ("key_trabbl") 
    389 !                          !  advective bbl                             ("key_trabbl_adv") 
    390    rn_ahtbbl   =  10000.   !  lateral mixing coefficient in the bbl  [m2/s] 
    391 / 
     548   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0) 
     549   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0) 
     550   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s] 
     551   rn_gambbl   =  10.    !  advective bbl coefficient                 [s] 
     552/ 
     553 
    392554!!====================================================================== 
    393555!!                        Tracer (T & S ) namelists 
     
    396558!!   namtra_adv    advection scheme 
    397559!!   namtra_ldf    lateral diffusion scheme 
    398 !!   namtra_dmp    T & S newtonian damping                              ("key_tradmp") 
    399 !!====================================================================== 
    400  
     560!!   namtra_dmp    T & S newtonian damping 
     561!!====================================================================== 
     562! 
    401563!----------------------------------------------------------------------- 
    402564&nameos        !   ocean physical parameters 
    403565!----------------------------------------------------------------------- 
    404    nn_eos      =    0      !  type of equation of state and Brunt-Vaisala frequency 
     566   nn_eos      =   0       !  type of equation of state and Brunt-Vaisala frequency 
    405567                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
    406568                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    407569                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    408    rn_alpha    =    2.e-4  !  thermal expension coefficient (neos= 1 or 2) 
    409    rn_beta     =    0.001  !  saline  expension coefficient (neos= 2) 
    410 / 
    411 !----------------------------------------------------------------------- 
    412 &namtra_adv    !   advection scheme for tracer  
    413 !----------------------------------------------------------------------- 
    414    ln_traadv_cen2   =  .false.  !  2nd order centered scheme    
    415    ln_traadv_tvd    =  .true.   !  TVD scheme                 
    416    ln_traadv_muscl  =  .false.  !  MUSCL scheme              
    417    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
    418    ln_traadv_ubs    =  .false.  !  UBS scheme                  
    419 / 
    420 !----------------------------------------------------------------------- 
    421 &namtra_ldf    !   lateral diffusion scheme for tracer  
    422 !----------------------------------------------------------------------- 
    423                            !  Type of the operator :  
    424    ln_traldf_lap    =  .true.   !     laplacian operator        
    425    ln_traldf_bilap  =  .false.  !     bilaplacian operator      
    426                            !  Direction of action  : 
    427    ln_traldf_level  =  .false.  !     iso-level                 
    428    ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    429    ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
    430                            !  Coefficient 
    431    rn_aht_0         =  2000.    !     horizontal eddy diffusivity for tracers [m2/s] 
    432    rn_ahtb_0        =     0.    !     background eddy diffusivity for ldf_iso [m2/s] 
    433    rn_aeiv_0        =  2000.    !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
    434 / 
    435 !----------------------------------------------------------------------- 
    436 &namtra_dmp    !   tracer: T & S newtonian damping                      ('key_tradmp') 
    437 !----------------------------------------------------------------------- 
     570   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
     571   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     572/ 
     573!----------------------------------------------------------------------- 
     574&namtra_adv    !   advection scheme for tracer 
     575!----------------------------------------------------------------------- 
     576   ln_traadv_cen2   =  .false.  !  2nd order centered scheme 
     577   ln_traadv_tvd    =  .true.   !  TVD scheme 
     578   ln_traadv_muscl  =  .false.  !  MUSCL scheme 
     579   ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
     580   ln_traadv_ubs    =  .false.  !  UBS scheme 
     581   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
     582   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     583/ 
     584!---------------------------------------------------------------------------------- 
     585&namtra_ldf    !   lateral diffusion scheme for tracers 
     586!---------------------------------------------------------------------------------- 
     587   !                       !  Operator type: 
     588   ln_traldf_lap    =  .true.   !  laplacian operator 
     589   ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     590   !                       !  Direction of action: 
     591   ln_traldf_level  =  .false.  !  iso-level 
     592   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
     593   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
     594   !                       !  Griffies parameters              (all need "key_ldfslp") 
     595   ln_traldf_grif   =  .false.  !  use griffies triads 
     596   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
     597   ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
     598   ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
     599   !                       !  Coefficients 
     600   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
     601   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
     602   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
     603   rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s] 
     604   rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     605   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
     606   !                                           (normally=0; not used with Griffies) 
     607/ 
     608!----------------------------------------------------------------------- 
     609&namtra_dmp    !   tracer: T & S newtonian damping 
     610!----------------------------------------------------------------------- 
     611   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    438612   nn_hdmp     =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
    439613                           !                   =XX, damping poleward of XX degrees (XX>0) 
    440614                           !                      + F(distance-to-coast) + Red and Med Seas 
    441    nn_zdmp     =    1      !  vertical   shape =0    damping throughout the water column 
     615   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column 
    442616                           !                   =1 no damping in the mixing layer (kz  criteria) 
    443617                           !                   =2 no damping in the mixed  layer (rho crieria) 
     
    445619   rn_bot      =  360.     !  bottom  time scale of damping   [days] 
    446620   rn_dep      =  800.     !  depth of transition between rn_surf and rn_bot [meters] 
    447    nn_file     =    1      !  create a damping.coeff NetCDF file (=1) or not (=0) 
    448 / 
     621   nn_file     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     622/ 
     623 
    449624!!====================================================================== 
    450625!!                      ***  Dynamics namelists  *** 
     
    456631!!   namdyn_ldf    lateral diffusion scheme 
    457632!!====================================================================== 
    458  
     633! 
    459634!----------------------------------------------------------------------- 
    460635&namdyn_adv    !   formulation of the momentum advection 
    461636!----------------------------------------------------------------------- 
    462    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)   
     637   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    463638   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    464    ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme  
    465 /   
     639   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     640/ 
    466641!----------------------------------------------------------------------- 
    467642&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    468643!----------------------------------------------------------------------- 
    469    ln_dynvor_ene = .false. !  enstrophy conserving scheme   
    470    ln_dynvor_ens = .false. !  energy conserving scheme     
    471    ln_dynvor_mix = .false. !  mixed scheme                
    472    ln_dynvor_een = .true.  !  energy & enstrophy scheme   
     644   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     645   ln_dynvor_ens = .false. !  energy conserving scheme 
     646   ln_dynvor_mix = .false. !  mixed scheme 
     647   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    473648/ 
    474649!----------------------------------------------------------------------- 
    475650&namdyn_hpg    !   Hydrostatic pressure gradient option 
    476651!----------------------------------------------------------------------- 
    477    ln_hpg_zco  = .false.   !  z-coordinate - full steps                    
     652   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    478653   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    479654   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    480    ln_hpg_hel  = .false.   !  s-coordinate (helsinki modification) 
    481    ln_hpg_wdj  = .false.   !  s-coordinate (weighted density jacobian) 
    482655   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    483    ln_hpg_rot  = .false.   !  s-coordinate (ROTated axes scheme) 
    484    rn_gamma    = 0.e0      !  weighting coefficient (wdj scheme) 
     656   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    485657   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    486658                                 !           centered      time scheme  (F) 
    487    nn_dynhpg_rst =  0            ! =1 dynhpg restartable restart or not (=0) 
    488659/ 
    489660!----------------------------------------------------------------------- 
     
    497668&namdyn_ldf    !   lateral diffusion on momentum 
    498669!----------------------------------------------------------------------- 
    499                            !  Type of the operator :  
    500    ln_dynldf_lap    =  .true.   !  laplacian operator          
    501    ln_dynldf_bilap  =  .false.  !  bilaplacian operator     
    502                            !  Direction of action  :  
    503    ln_dynldf_level  =  .false.  !  iso-level                
     670   !                       !  Type of the operator : 
     671   ln_dynldf_lap    =  .true.   !  laplacian operator 
     672   ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     673   !                       !  Direction of action  : 
     674   ln_dynldf_level  =  .false.  !  iso-level 
    504675   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    505676   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
    506                            !  Coefficient 
    507    rn_ahm_0    = 40000.         !  horizontal eddy viscosity   [m2/s] 
    508    rn_ahmb_0   =     0.         !  background eddy viscosity for ldf_iso [m2/s] 
    509 / 
     677   !                       !  Coefficient 
     678   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
     679   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
     680   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
     681/ 
     682 
    510683!!====================================================================== 
    511684!!             Tracers & Dynamics vertical physics namelists 
    512685!!====================================================================== 
    513 !!       namzdf        vertical physics 
    514 !!       namzdf_ric    richardson number dependent vertical mixing      ("key_zdfric"      ) 
    515 !!       namzdf_tke    TKE dependent vertical mixing                    ("key_zdftke"      ) 
    516 !!       namzdf_kpp    KPP dependent vertical mixing                    ("key_zdfkpp"      ) 
    517 !!       namzdf_ddm    double diffusive mixing parameterization         ("key_zdfddm"      ) 
    518 !!       namzdf_tmx    tidal mixing parameterization                    ("key_zdftmx"      ) 
    519 !!====================================================================== 
    520  
     686!!    namzdf        vertical physics 
     687!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
     688!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
     689!!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp") 
     690!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
     691!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
     692!!====================================================================== 
     693! 
    521694!----------------------------------------------------------------------- 
    522695&namzdf        !   vertical physics 
     
    529702   nn_evdm     =    0      !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    530703   rn_avevd    =  100.     !  evd mixing coefficient [m2/s] 
    531    ln_zdfnpc   = .false.   !  Non-Penetrative algorithm (T) or not (F) 
     704   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F) 
    532705   nn_npc      =    1            !  frequency of application of npc 
    533706   nn_npcp     =  365            !  npc control print frequency 
     
    541714   rn_alp      =   5.      !  coefficient of the parameterization 
    542715   nn_ric      =   2       !  coefficient of the parameterization 
     716   rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation 
     717   rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m) 
     718   rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m) 
     719   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
     720   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
     721   ln_mldw     = .true.    !  Flag to use or not the mized layer depth param. 
    543722/ 
    544723!----------------------------------------------------------------------- 
     
    547726   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
    548727   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation 
    549    rn_ebb      =  60.      !  coef. of the surface input of tke 
     728   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 
    550729   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2] 
    551730   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2] 
    552    rn_bshear   =   1.e-20  !  background shear (>0) 
    553731   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
    554732                           !                 = 1 bounded by the local vertical scale factor 
    555733                           !                 = 2 first vertical derivative of mixing length bounded by 1 
    556                            !                 = 3 same criteria as case 2 but applied in a different way 
     734                           !                 = 3 as =2 with distinct disspipative an mixing length scale 
    557735   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    558    ln_mxl0     = .true.    !  mixing length scale surface value as function of wind stress (T) or not (F) 
    559    rn_lmin     =   0.001   !  interior buoyancy lenght scale minimum value 
    560    rn_lmin0    =   0.01    !  surface  buoyancy lenght scale minimum value 
    561    nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves 
    562                            !        = 0 no penetration ( O(2 km) resolution) 
    563                            !        = 1 additional tke source (rn_efr * en) 
    564                            !        = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 
    565                            !        = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 
    566    nn_htau     =   1       !  type of exponential decrease of tke penetration 
     736   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     737   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
     738   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
     739   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
     740   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
     741                           !        = 0 no penetration 
     742                           !        = 1 add a tke source below the ML 
     743                           !        = 2 add a tke source just at the base of the ML 
     744                           !        = 3 as = 1 applied on HF part of the stress    ("key_coupled") 
     745   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
     746   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
    567747                           !        = 0  constant 10 m length scale 
    568                            !        = 1  0.5m at the equator to 30m at high latitudes 
    569                            !        = 2  30 meters constant depth penetration 
    570                            !  otion used only if nn_etau = 1 or 2:  
    571    rn_efr      =   0.05    !     fraction of surface tke value which penetrates inside the ocean 
    572                            !  otion used only if nn_etau = 3: 
    573    rn_addhft   =  -1.e-3   !     add offset   applied to the "mean of stress module - module of mean stress" (always kept > 0) 
    574    rn_sclhft   =   1.      !     scale factor applied to the "mean of stress module - module of mean stress" 
    575    ln_lc       = .true.    !  Langmuir cell parameterisation 
    576    rn_lc       =   0.15    !  coef. associated to Langmuir cells 
     748                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    577749/ 
    578750!------------------------------------------------------------------------ 
    579 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionnally: 
     751&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    580752!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    581    ln_kpprimix = .true.    !  shear instability mixing  
     753   ln_kpprimix = .true.    !  shear instability mixing 
    582754   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    583755   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    584756   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    585757   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
    586    rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]  
    587    rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]  
     758   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2] 
     759   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s] 
    588760   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    589761   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
     762/ 
     763!----------------------------------------------------------------------- 
     764&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     765!----------------------------------------------------------------------- 
     766   rn_emin       = 1.e-6   !  minimum value of e   [m2/s2] 
     767   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3] 
     768   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988) 
     769   rn_clim_galp  = 0.53    !  galperin limit 
     770   ln_crban      = .true.  !  Use Craig & Banner (1994) surface wave mixing parametrisation 
     771   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case 
     772   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux 
     773   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
     774   nn_tkebc_surf =   1     !  surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
     775   nn_tkebc_bot  =   1     !  bottom tke condition (0/1=Dir/Neum) 
     776   nn_psibc_surf =   1     !  surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
     777   nn_psibc_bot  =   1     !  bottom psi condition (0/1=Dir/Neum) 
     778   nn_stab_func  =   2     !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
     779   nn_clos       =   1     !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
    590780/ 
    591781!----------------------------------------------------------------------- 
     
    601791   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1) 
    602792   rn_tfe      = 0.333     !  tidal dissipation efficiency 
    603    rn_me       = 0.2       !  mixing efficiency  
    604    ln_tmx_itf  = .TRUE.    !  ITF specific parameterisation 
     793   rn_me       = 0.2       !  mixing efficiency 
     794   ln_tmx_itf  = .true.    !  ITF specific parameterisation 
    605795   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
    606796/ 
    607 !!====================================================================== 
    608 !!                  ***  Miscelaneous namelists  *** 
     797 
     798!!====================================================================== 
     799!!                  ***  Miscellaneous namelists  *** 
    609800!!====================================================================== 
    610801!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    611 !!   nammpp_dyndist    Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist") 
    612802!!   namctl            Control prints & Benchmark 
    613 !!   namsol            elliptic solver / island / free surface  
    614 !!====================================================================== 
    615  
    616 !----------------------------------------------------------------------- 
    617 &namsol        !   elliptic solver / island / free surface  
     803!!   namsol            elliptic solver / island / free surface 
     804!!====================================================================== 
     805! 
     806!----------------------------------------------------------------------- 
     807&namsol        !   elliptic solver / island / free surface 
    618808!----------------------------------------------------------------------- 
    619809   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
     
    633823                           !  buffer blocking send or immediate non-blocking sends, resp. 
    634824   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    635 / 
    636 !----------------------------------------------------------------------- 
    637 &nammpp_dyndist !   Massively Parallel Distribution for AGRIF zoom      ("key_agrif" && "key_mpp_dyndist") 
    638 !----------------------------------------------------------------------- 
    639    jpni        =    1      !  jpni   number of processors following i 
    640    jpnj        =    1      !  jpnj   number of processors following j 
    641    jpnij       =    1      !  jpnij  number of local domains 
     825   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold 
     826   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1) 
     827   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1) 
     828   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1) 
    642829/ 
    643830!----------------------------------------------------------------------- 
     
    654841   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    655842                           !     (no physical validity of the results) 
    656    nn_bit_cmp  =    0      !  bit comparison mode (1/0): CAUTION use zero except for test 
    657                            !     of comparison between single and multiple processor runs 
    658 / 
    659  
    660 !!====================================================================== 
    661 !!                       ***  Diagnostics namelists  *** 
    662 !!====================================================================== 
     843   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
     844/ 
     845 
     846!!====================================================================== 
     847!!                  ***  Diagnostics namelists  *** 
     848!!====================================================================== 
     849!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    663850!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
    664 !!   namgap       level mean model-data gap                             ("key_diagap") 
    665851!!   namflo       float parameters                                      ("key_float") 
    666852!!   namptr       Poleward Transport Diagnostics 
    667 !!====================================================================== 
    668  
     853!!   namhsb       Heat and salt budgets 
     854!!====================================================================== 
     855! 
     856!----------------------------------------------------------------------- 
     857&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     858!----------------------------------------------------------------------- 
     859   nn_nchunks_i=   4       !  number of chunks in i-dimension 
     860   nn_nchunks_j=   4       !  number of chunks in j-dimension 
     861   nn_nchunks_k=   31      !  number of chunks in k-dimension 
     862                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
     863                           !  is optimal for postprocessing which works exclusively with horizontal slabs 
     864   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
     865                           !  (F) ignore chunking information and produce netcdf3-compatible files 
     866/ 
    669867!----------------------------------------------------------------------- 
    670868&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra") 
    671 !              !       or mixed-layer trends or barotropic vorticity    ('key_trdmld' or  "key_trdvor") 
     869!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor") 
    672870!----------------------------------------------------------------------- 
    673871   nn_trd      = 365       !  time step frequency dynamics and tracers trends 
     
    680878/ 
    681879!----------------------------------------------------------------------- 
    682 &namgap       !   level mean model-data gap                             ('key_diagap') 
    683 !----------------------------------------------------------------------- 
    684    nn_gap     =  15        !  time-step frequency of model-data gap computation 
    685    nn_prg     =  10        !  AUTO - time-step frequency of gap print in model output 
    686 / 
    687 !----------------------------------------------------------------------- 
    688880&namflo       !   float parameters                                      ("key_float") 
    689881!----------------------------------------------------------------------- 
    690     ln_rstflo = .false.    !  float restart (T) or not (F) 
    691     nn_writefl=      75    !  frequency of writing in float output file  
    692     nn_stockfl=    5475    !  frequency of creation of the float restart file  
    693     ln_argo   = .false.    !  Argo type floats (stay at the surface each 10 days) 
    694     ln_flork4 = .false.    !  trajectories computed with a 4th order Runge-Kutta (T) 
    695                            !  or computed with Blanke' scheme (F) 
     882   jpnfl         = 1          !  total number of floats during the run 
     883   jpnnewflo     = 0          !  number of floats for the restart 
     884   ln_rstflo     = .false.    !  float restart (T) or not (F) 
     885   nn_writefl    =      75    !  frequency of writing in float output file 
     886   nn_stockfl    =    5475    !  frequency of creation of the float restart file 
     887   ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days) 
     888   ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T) 
     889                              !  or computed with Blanke' scheme (F) 
     890   ln_ariane     = .true.     !  Input with Ariane tool convention(T) 
     891   ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
    696892/ 
    697893!----------------------------------------------------------------------- 
    698894&namptr       !   Poleward Transport Diagnostic 
    699895!----------------------------------------------------------------------- 
    700    ln_diaptr  = .true.     !  Poleward heat and salt transport (T) or not (F) 
     896   ln_diaptr  = .true.    !  Poleward heat and salt transport (T) or not (F) 
    701897   ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    702    ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not  
     898   ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    703899                           !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    704900   ln_ptrcomp = .true.     !  Add decomposition : overturning 
    705    nf_ptr     =  1         !  Frequency of ptr computation [time step] 
    706    nf_ptr_wri =  15        !  AUTO - Frequency of ptr outputs 
    707 / 
     901   nn_fptr    =  1         !  Frequency of ptr computation [time step] 
     902   nn_fwri    =  15        !  Frequency of ptr outputs [time step] 
     903/ 
     904!----------------------------------------------------------------------- 
     905&namhsb       !  Heat and salt budgets 
     906!----------------------------------------------------------------------- 
     907   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
     908/ 
     909!----------------------------------------------------------------------- 
     910&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm') 
     911!----------------------------------------------------------------------- 
     912    nit000_han = 1         ! First time step used for harmonic analysis 
     913    nitend_han = 75        ! Last time step used for harmonic analysis 
     914    nstep_han  = 15        ! Time step frequency for harmonic analysis 
     915    tname(1)   = 'M2'      ! Name of tidal constituents 
     916    tname(2)   = 'K1' 
     917/ 
     918!----------------------------------------------------------------------- 
     919&namdct        ! transports through sections 
     920!----------------------------------------------------------------------- 
     921    nn_dct      = 15       !  time step frequency for transports computing 
     922    nn_dctwri   = 15       !  time step frequency for transports writing 
     923    nn_secdebug = 112      !      0 : no section to debug 
     924                           !     -1 : debug all section 
     925                           !  0 < n : debug section number n 
     926/ 
     927 
     928!!====================================================================== 
     929!!            ***  Observation & Assimilation namelists *** 
     930!!====================================================================== 
     931!!   namobs       observation and model comparison                      ('key_diaobs') 
     932!!   nam_asminc   assimilation increments                               ('key_asminc') 
     933!!====================================================================== 
     934! 
     935!----------------------------------------------------------------------- 
     936&namobs       !  observation usage switch                               ('key_diaobs') 
     937!----------------------------------------------------------------------- 
     938   ln_t3d     = .false.    ! Logical switch for T profile observations 
     939   ln_s3d     = .false.    ! Logical switch for S profile observations 
     940   ln_ena     = .false.    ! Logical switch for ENACT insitu data set 
     941   !                       !     ln_cor                  Logical switch for Coriolis insitu data set 
     942   ln_profb   = .false.    ! Logical switch for feedback insitu data set 
     943   ln_sla     = .false.    ! Logical switch for SLA observations 
     944 
     945   ln_sladt   = .false.    ! Logical switch for AVISO SLA data 
     946 
     947   ln_slafb   = .false.    ! Logical switch for feedback SLA data 
     948                           !     ln_ssh                  Logical switch for SSH observations 
     949 
     950   ln_sst     = .true.     ! Logical switch for SST observations 
     951   ln_reysst  = .true.     !     ln_reysst               Logical switch for Reynolds observations 
     952   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
     953 
     954   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
     955                           !     ln_sss                  Logical switch for SSS observations 
     956                           !     ln_seaice               Logical switch for Sea Ice observations 
     957                           !     ln_vel3d                Logical switch for velocity observations 
     958                           !     ln_velavcur             Logical switch for velocity daily av. cur. 
     959                           !     ln_velhrcur             Logical switch for velocity high freq. cur. 
     960                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP 
     961                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP 
     962                           !     ln_velfb                Logical switch for feedback velocity data 
     963                           !     ln_grid_global          Global distribtion of observations 
     964                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table 
     965                           !     grid_search_file        Grid search lookup file header 
     966                           !     enactfiles              ENACT input observation file names 
     967                           !     coriofiles              Coriolis input observation file name 
     968   !                       ! profbfiles: Profile feedback input observation file name 
     969   profbfiles = 'profiles_01.nc' 
     970                           !     ln_profb_enatim         Enact feedback input time setting switch 
     971                           !     slafilesact             Active SLA input observation file name 
     972                           !     slafilespas             Passive SLA input observation file name 
     973   !                       ! slafbfiles: Feedback SLA input observation file name 
     974   slafbfiles = 'sla_01.nc' 
     975                           !     sstfiles                GHRSST input observation file name 
     976   !                       ! sstfbfiles: Feedback SST input observation file name 
     977   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
     978                           !     seaicefiles             Sea Ice input observation file name 
     979                           !     velavcurfiles           Vel. cur. daily av. input file name 
     980                           !     velhvcurfiles           Vel. cur. high freq. input file name 
     981                           !     velavadcpfiles          Vel. ADCP daily av. input file name 
     982                           !     velhvadcpfiles          Vel. ADCP high freq. input file name 
     983                           !     velfbfiles              Vel. feedback input observation file name 
     984                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS 
     985                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS 
     986                           !     n1dint                  Type of vertical interpolation method 
     987                           !     n2dint                  Type of horizontal interpolation method 
     988                           !     ln_nea                  Rejection of observations near land switch 
     989   nmsshc     = 0          ! MSSH correction scheme 
     990                           !     mdtcorr                 MDT  correction 
     991                           !     mdtcutoff               MDT cutoff for computed correction 
     992   ln_altbias = .false.    ! Logical switch for alt bias 
     993   ln_ignmis  = .true.     ! Logical switch for ignoring missing files 
     994                           !     endailyavtypes   ENACT daily average types 
     995   ln_grid_global = .true. 
     996   ln_grid_search_lookup = .false. 
     997/ 
     998!----------------------------------------------------------------------- 
     999&nam_asminc   !   assimilation increments                               ('key_asminc') 
     1000!----------------------------------------------------------------------- 
     1001    ln_bkgwri = .false.    !  Logical switch for writing out background state 
     1002    ln_trainc = .false.    !  Logical switch for applying tracer increments 
     1003    ln_dyninc = .false.    !  Logical switch for applying velocity increments 
     1004    ln_sshinc = .false.    !  Logical switch for applying SSH increments 
     1005    ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI) 
     1006    ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU) 
     1007    nitbkg    = 0          !  Timestep of background in [0,nitend-nit000-1] 
     1008    nitdin    = 0          !  Timestep of background for DI in [0,nitend-nit000-1] 
     1009    nitiaustr = 1          !  Timestep of start of IAU interval in [0,nitend-nit000-1] 
     1010    nitiaufin = 15         !  Timestep of end of IAU interval in [0,nitend-nit000-1] 
     1011    niaufn    = 0          !  Type of IAU weighting function 
     1012    ln_salfix = .false.    !  Logical switch for ensuring that the sa > salfixmin 
     1013    salfixmin = -9999      !  Minimum salinity after applying the increments 
     1014    nn_divdmp = 0          !  Number of iterations of divergence damping operator 
     1015/ 
     1016!----------------------------------------------------------------------- 
     1017&namsbc_wave   ! External fields from wave model 
     1018!----------------------------------------------------------------------- 
     1019!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     1020!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     1021   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' , .true.   , .false. , 'daily'  ,''         , '' 
     1022   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
     1023   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
     1024   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   , .true.   , .false. , 'daily'  ,''         , '' 
     1025! 
     1026   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
     1027/ 
     1028!----------------------------------------------------------------------- 
     1029&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
     1030!----------------------------------------------------------------------- 
     1031   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 
     1032   ln_neptsimp       = .false.  ! yes/no use simplified neptune 
     1033 
     1034   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep 
     1035   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator 
     1036   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole 
     1037   ! Specify whether to ramp down the Neptune velocity in shallow 
     1038   ! water, and if so the depth range controlling such ramping down 
     1039   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water 
     1040   rn_htrmin         =  100.0   ! min. depth of transition range 
     1041   rn_htrmax         =  200.0   ! max. depth of transition range 
     1042/ 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ice_ORCA2

    r705 r2061  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/LIM2 : 1 - dynamics/advection/thermo          (namicerun) 
     2!! NEMO/LIM-2 : 1 - dynamics/advection/thermo          (namicerun) 
    33!! namelists    2 - ice intialisation                  (namiceini) 
    44!!              3 - ice dynamic                        (namicedyn) 
    5 !!              5 - ice advection                      (namicetrp) 
    6 !!              6 - thermodynamic                      (namicethd) 
    7 !!              7 - ice salinity                       (namicesal) 
    8 !!              8 - mechanical redistribution of ice   (namiceitdme) 
    9 !!              3 - ice diagnostics                    (namicedia) 
    10 !!              9 - ice outputs                        (namiceout) 
     5!!              4 - ice advection                      (namicetrp) 
     6!!              5 - thermodynamic                      (namicethd) 
     7!!              6 - ice damping                        (namice_dmp) 
     8!!              7 - ice diagnostics                    (namicedia) 
     9!!              8 - ice outputs                        (namiceout) 
    1110!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1211 
     
    1615   cn_icerst_in  = "restart_ice_in"   !  suffix of ice restart name (input) 
    1716   cn_icerst_out = "restart_ice"      !  suffix of ice restart name (output) 
    18    ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F) 
    19    ln_limdmp   = .false.   !  restoring ice thickness and fraction leads (T) or not (F) 
    20    acrit       = 1.0e-06   , 1.0e-06   ! minimum fraction for leads in the Northern (Southern) Hemisphere 
    21    hsndif      =  0.0      !  computation of temperature in snow (=0.0) or not 
    22    hicdif      =  0.0      !  computation of temperature in ice  (=0.0) or not (=9999.0) 
     17   ln_limdyn     = .true.             !  ice dynamics (T) or thermodynamics only (F) 
     18   ln_limdmp     = .false.            !  restoring ice thickness and fraction leads   (T => fill  namice_dmp) 
     19   acrit         = 1.0e-06 , 1.0e-06  !  minimum lead fraction in the Northern & Southern Hemispheres 
     20   hsndif        =  0.e0              !  computation of temperature in snow (=0.0) or not 
     21   hicdif        =  0.e0              !  computation of temperature in ice  (=0.0) or not (=9999.0) 
    2322/ 
    2423!----------------------------------------------------------------------- 
    2524&namiceini     !   ice initialisation 
    2625!----------------------------------------------------------------------- 
    27    ln_limini   = .false.   !  read the ice initial state in the file 'Ice_initialization.nc' (T) or not (F) 
     26   ln_limini   = .false.   !  read the initial state in 'Ice_initialization.nc' (T) or not (F) 
    2827   ttest       =  2.0      !  threshold water temperature for initial sea ice 
    2928   hninn       =  0.5      !  initial snow thickness in the north 
     
    4948   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    5049   etamn       =   0.0e+07 !  minimun value for viscosity 
    51    creepl      =   2.0e-08 !  creep limit 
     50   creepl      =   1.0e-08 !  creep limit 
    5251   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    5352   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
     53   nevp        =   120     !  number of EVP subcycling iterations 
     54   telast      =   9600    !  timescale for EVP elastic waves 
     55   alphaevp    =   1.0     !  coefficient for the solution of EVP int. stresses 
    5456/ 
    5557!----------------------------------------------------------------------- 
    5658&namicetrp     !   ice transport 
    5759!----------------------------------------------------------------------- 
    58    bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip) 
     60   bound       =   0.0     !  boundary conditions (=0. no-slip, =1. free-slip) 
    5961/ 
    6062!----------------------------------------------------------------------- 
    6163&namicethd     !   ice thermodynamic 
    6264!----------------------------------------------------------------------- 
    63 !  hmelt       : maximum melting at the bottom 
    64 !  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
    65 !                caution 1.0, 1.0 best value to be used!!! (gilles G.) 
    66 !  hicmin      : ice thickness corr. to max. energy stored in brine pocket 
    67 !  hiclim      : minimum ice thickness 
    68 !  amax        : maximum lead fraction 
    69 !  swiqst      : energy stored in brine pocket (=1) or not (=0)  
    70 !  sbeta       : numerical caracteritic of the scheme for diffusion in ice 
    71 !               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) 
    72 !  parlat      : percentage of energy used for lateral ablation 
    73 !  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting 
    74 !  hibspl      : slope of distribution for Hibler's lateral melting 
    75 !  exld        : exponent for leads-closure rate 
    76 !  hakdif      : coefficient for diffusions of ice and snow 
    77 !  thth         : threshold thickness for comp. of eq. thermal conductivity 
    78 !  hnzst       : thickness of the surf. layer in temp. computation 
    79 !  parsub      : switch for snow sublimation or not 
    80 !  alphs       : coefficient for snow density when snow ice formation 
    81 !  
    8265   hmelt       = -0.15     !  maximum melting at the bottom 
    8366   hiccrit     = 0.3 , 0.3 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     67   !                       !  (caution 1.0, 1.0 best value to be used!!! (gilles G.)) 
    8468   hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket 
    8569   hiclim      = 0.05      !  minimum ice thickness 
     
    9579   thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity 
    9680   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
    97    parsub      = 0.0       !  switch for snow sublimation or not 
     81   parsub      = 1.0       !  switch for snow sublimation or not 
    9882   alphs       = 1.0       !  coefficient for snow density when snow ice formation 
     83/ 
     84!----------------------------------------------------------------------- 
     85&namice_dmp    !   damping of sea ice alone open boundaries 
     86!              !   (hard coded damping area: check if it fit your config) 
     87!----------------------------------------------------------------------- 
     88!              !  file name   ! frequency (hours) ! variable ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     89!              !              !  (if <0  months)  !   name   !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     90   sn_hicif    = 'ice_damping',  -1.              , 'hicif'  ,   .true.     , .true. , 'yearly'  , ''       , '' 
     91   sn_cnf      = 'ice_damping',  -1.              , 'frld'   ,   .true.     , .true. , 'yearly'  , ''       , '' 
     92! 
     93   cn_dir      = './'      !  root directory for the location of the runoff files 
    9994/ 
    10095!----------------------------------------------------------------------- 
    10196&namicedia     !   ice diagnostics 
    10297!----------------------------------------------------------------------- 
    103    fmtinf      ='1PE13.5 ' !  format of the output values 
    104    nfrinf      = 4         !  number of variables written in one line 
    105    ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
    106    ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     98   fmtinf      = '1PE13.5' !  format of the output values 
     99   nfrinf      =  4        !  number of variables written in one line 
     100   ntmoy       =  1        !  instantaneous values of ice evolution or averaging 
     101   ninfo       =  1        !  frequency of ouputs on file ice_evolu in case of averaging 
    107102/ 
    108103!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    109104&namiceout     !   parameters for outputs 
    110105!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     106!              !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
     107!              !                                      !           !           ! or not !  factor  !  factor  ! 
     108   field_1     = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
     109   field_2     = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
     110   field_3     = 'Ice produced                       ', 'iiceprod', 'm/kt    ',    1   ,  1.0     ,    0.0 
     111   field_4     = 'Ice concentration                  ', 'ileadfra', '%       ',    1   , -1.0     ,    1.0 
     112   field_5     = 'Ice temperature                    ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
     113   field_6     = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
     114   field_7     = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
     115   field_8     = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
     116   field_9     = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
     117   field_10    = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
     118   field_11    = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     119   field_12    = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     120   field_13    = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
     121   field_14    = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt',    1   ,  1.0     ,    0.0 
     122   field_15    = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
     123   field_16    = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
     124   field_17    = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     125   field_18    = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     126   field_19    = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ',    1   ,  1.0     ,    0.0 
     127! 
    111128   noumef      =   19      !  number of fields 
    112 ! 
    113 !           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
    114 !           !                                      !           !           ! or not !  factor  !  factor  ! 
    115    field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
    116    field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
    117    field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ',    1   ,  1.0     ,    0.0 
    118    field_4  = 'Ice concentration                  ', 'ileadfra', '%       ',    1   , -1.0     ,    1.0 
    119    field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
    120    field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
    121    field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
    122    field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
    123    field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
    124    field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
    125    field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    126    field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    127    field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
    128    field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt',    1   ,  1.0     ,    0.0 
    129    field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
    130    field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
    131    field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    132    field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    133    field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ',    1   ,  1.0     ,    0.0 
    134129/       
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_pisces_ORCA2

    r1775 r2061  
    1515&nampisext     !   air-sea exchange 
    1616!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     17   ln_co2int  =  .false. ! read atm pco2 from a file (T) or constant (F) 
    1718   atcco2     = 0.36886500E+03 ! AUTO - atmospheric pCO2 
     19   clname     =  'atcco2.txt'  ! Name of atm pCO2 file - ln_co2int = T 
     20   nn_offset  =  0       ! Offset model-data start year - ln_co2int = T 
     21!                        ! If your model year is iyy, nn_offset=(years(1)-iyy)  
     22!                        ! then the first atmospheric CO2 record read is at years(1) 
     23/ 
     24!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     25&nampisatm     !  Atmospheric prrssure  
     26!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     27!              !  file name   ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     28!              !              !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     29   sn_patm     = 'presatm'    ,     -1            , 'patm'     ,  .true.      , .true. ,   'yearly'  , ''       , '' 
     30   cn_dir      = './'     !  root directory for the location of the dynamical files 
     31! 
     32   ln_presatm  = .true.   ! constant atmopsheric pressure (F) or from a file (T) 
    1833/ 
    1934!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    2035&nampisbio     !   biological parameters 
    2136!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    22    part       =  0.85    ! part of calcite not dissolved in guts 
    23    nrdttrc    =  1       ! time step frequency for biology 
    24    wsbio      =  2.      ! POC sinking speed 
    25    xkmort     =  1.E-7   ! half saturation constant for mortality 
    26    ferat3     =  3.E-6   ! Fe/C in zooplankton  
    27    wsbio2     =  50.     ! Big particles sinking speed 
     37   nrdttrc    =  1        ! time step frequency for biology 
     38   wsbio      =  2.       ! POC sinking speed 
     39   xkmort     =  2.E-7    ! half saturation constant for mortality 
     40   ferat3     =  10.E-6   ! Fe/C in zooplankton  
     41   wsbio2     =  30.      ! Big particles sinking speed 
     42   niter1max  =  1        ! Maximum number of iterations for POC 
     43   niter2max  =  1        ! Maximum number of iterations for GOC 
    2844/ 
    2945!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    3046&nampislim     !   parameters for nutrient limitations 
    3147!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    32    conc0      =  2.e-6    ! Phosphate half saturation 
    33    conc1      =  10E-6    ! Phosphate half saturation for diatoms 
    34    conc2      =  0.02E-9  ! Iron half saturation for phyto 
    35    conc2m     =  0.08E-9  ! Max iron half saturation for phyto 
    36    conc3      =  0.1E-9   ! Iron half saturation for diatoms 
    37    conc3m     =  0.4E-9   ! Maxi iron half saturation for diatoms 
     48   concnno3   =  1.e-6    ! Nitrate half saturation of nanophytoplankton 
     49   concdno3   =  3.E-6   ! Phosphate half saturation for diatoms 
    3850   concnnh4   =  1.E-7    ! NH4 half saturation for phyto 
    39    concdnh4   =  5.E-7    ! NH4 half saturation for diatoms 
     51   concdnh4   =  3.E-7   ! NH4 half saturation for diatoms 
     52   concnfer   =  1.E-9     ! Iron half saturation for phyto 
     53   concdfer   =  3.E-9   ! Iron half saturation for diatoms 
     54   concbfe    =  1.E-11   ! Half-saturation for Fe limitation of Bacteria 
     55   concbnh4   =  2.5E-8   ! NH4 half saturation for phyto 
     56   concbno3   =  2.5E-7   ! Phosphate half saturation for diatoms 
     57   xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
     58   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     59   xsizern    =  3.0      ! Size ratio for nanophytoplankton 
     60   xsizerd    =  3.0      ! Size ratio for diatoms 
    4061   xksi1      =  2.E-6    ! half saturation constant for Si uptake 
    41    xksi2      =  3.33E-6  ! half saturation constant for Si/C 
     62   xksi2      =  20E-6  ! half saturation constant for Si/C 
    4263   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
    43    caco3r     =  0.15      ! mean rain ratio 
    44 / 
     64   qnfelim    =  7.E-6    ! Optimal quota of phyto 
     65   qdfelim    =  7.E-6    ! Optimal quota of diatoms 
     66   caco3r     =  0.3      ! mean rain ratio 
     67/ 
     68!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     69&nampisopt     !   parameters for optics 
     70!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     71!              !  file name        ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     72!              !                   !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     73   sn_par     = 'par.orca'       ,     24            , 'fr_par'     ,  .true.      , .true. ,   'yearly'  , ''       , '' 
     74   cn_dir      = './'      !  root directory for the location of the dynamical files 
     75   ln_varpar   =  .true.   ! boolean for PAR variable 
     76   parlux      =  0.43      ! Fraction of shortwave as PAR 
     77/  
    4578!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    4679&nampisprod     !   parameters for phytoplankton growth 
    4780!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    48    pislope    =  4.       ! P-I slope   
    49    pislope2   =  4.       ! P-I slope  for diatoms 
     81   pislope    =  2.       ! P-I slope 
     82   pislope2   =  2.       ! P-I slope  for diatoms 
    5083   excret     =  0.05     ! excretion ratio of phytoplankton 
    5184   excret2    =  0.05     ! excretion ratio of diatoms 
     85   ln_newprod =  .true.   ! Enable new parame. of production (T/F)  
     86   bresp      =  0.00333  ! Basal respiration rate 
    5287   chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton 
    5388   chlcdm     =  0.05     ! Minimum Chl/C in diatoms 
    54    fecnm      =  10E-6    ! Maximum Fe/C in nanophytoplankton 
    55    fecdm      =  15E-6    ! Minimum Fe/C in diatoms 
    56    grosip     =  0.151    ! mean Si/C ratio 
     89   chlcmin    =  0.004    ! Maximum Chl/c in phytoplankton 
     90   fecnm      =  40E-6    ! Maximum Fe/C in nanophytoplankton 
     91   fecdm      =  40E-6    ! Minimum Fe/C in diatoms 
     92   grosip     =  0.159    ! mean Si/C ratio 
    5793/ 
    5894!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    5995&nampismort     !   parameters for phytoplankton sinks 
    6096!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    61    wchl       =  0.001    ! quadratic mortality of phytoplankton 
    62    wchld      =  0.02     ! maximum quadratic mortality of diatoms 
     97   wchl       =  0.01    ! quadratic mortality of phytoplankton 
     98   wchld      =  0.01     ! maximum quadratic mortality of diatoms 
     99   wchldm     =  0.03     ! maximum quadratic mortality of diatoms 
    63100   mprat      =  0.01     ! phytoplankton mortality rate 
    64101   mprat2     =  0.01     ! Diatoms mortality rate 
    65    mpratm     =  0.01     ! Phytoplankton minimum mortality rate 
    66102/ 
    67103!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    68104&nampismes     !   parameters for mesozooplankton 
    69105!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    70    grazrat2   =  0.7      ! maximal mesozoo grazing rate 
     106   part2      =  0.75     ! part of calcite not dissolved in mesozoo guts 
     107   grazrat2   =  0.75     ! maximal mesozoo grazing rate 
    71108   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    72109   mzrat2     =  0.03     ! mesozooplankton mortality rate 
    73110   xprefc     =  1.       ! zoo preference for phyto 
    74    xprefp     =  0.2      ! zoo preference for POC 
     111   xprefp     =  0.3      ! zoo preference for POC 
    75112   xprefz     =  1.       ! zoo preference for zoo 
    76    xprefpoc   =  0.2      ! zoo preference for poc 
     113   xprefpoc   =  0.3      ! zoo preference for poc 
     114   xthresh2zoo = 1E-8     ! zoo feeding threshold for mesozooplankton  
     115   xthresh2dia = 1E-8     ! diatoms feeding threshold for mesozooplankton  
     116   xthresh2phy = 1E-8     ! nanophyto feeding threshold for mesozooplankton  
     117   xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton  
     118   xthresh2   =  3E-7     ! Food threshold for grazing 
    77119   xkgraz2    =  20.E-6   ! half sturation constant for meso grazing 
    78    epsher2    =  0.33     ! Efficicency of Mesozoo growth  
     120   epsher2    =  0.3      ! Efficicency of Mesozoo growth 
    79121   sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
    80122   unass2     =  0.3      ! non assimilated fraction of P by mesozoo 
    81    grazflux   =  5.e3     ! flux-feeding rate 
     123   grazflux   =  2.e3     ! flux-feeding rate 
    82124/ 
    83125!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    84126&nampiszoo     !   parameters for microzooplankton 
    85127!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    86    grazrat    =  4.0      ! maximal zoo grazing rate    
     128   part       =  0.5      ! part of calcite not dissolved in microzoo gutsa 
     129   grazrat    =  3.0      ! maximal zoo grazing rate 
    87130   resrat     =  0.03     ! exsudation rate of zooplankton 
    88    mzrat      =  0.0      ! zooplankton mortality rate 
    89    xpref2c    =  0.0      ! Microzoo preference for POM  
    90    xpref2p    =  0.5      ! Microzoo preference for Nanophyto 
     131   mzrat      =  0.004    ! zooplankton mortality rate 
     132   xpref2c    =  0.1      ! Microzoo preference for POM 
     133   xpref2p    =  1.       ! Microzoo preference for Nanophyto 
    91134   xpref2d    =  0.5      ! Microzoo preference for Diatoms 
    92    xkgraz     =  20.E-6   ! half sturation constant for grazing  
    93    epsher     =  0.33     ! Efficiency of microzoo growth 
     135   xthreshdia =  1.E-8    ! Diatoms feeding threshold for microzooplankton  
     136   xthreshphy =  1.E-8    ! Nanophyto feeding threshold for microzooplankton  
     137   xthreshpoc =  1.E-8    ! POC feeding threshold for microzooplankton  
     138   xthresh    =  3.E-7    ! Food threshold for feeding 
     139   xkgraz     =  20.E-6   ! half sturation constant for grazing 
     140   epsher     =  0.3      ! Efficiency of microzoo growth 
    94141   sigma1     =  0.6      ! Fraction of microzoo excretion as DOM 
    95142   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
    96143/ 
    97144!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     145&nampisfer     !   parameters for iron chemistry 
     146!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     147   ln_fechem =  .false.   ! complex iron chemistry ( T/F ) 
     148   ln_ligvar =  .true.   ! variable ligand concentration 
     149   xlam1     =  0.005     ! scavenging rate of Iron 
     150   xlamdust  =  150.0     ! Scavenging rate of dust 
     151   ligand    =  0.6E-9    ! Ligands concentration  
     152 
     153!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    98154&nampisrem     !   parameters for remineralization 
    99155!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    100    xremik    =  0.3       ! remineralization rate of DOC 
     156   xremik    =  0.35      ! remineralization rate of DOC 
    101157   xremip    =  0.025     ! remineralisation rate of POC 
    102158   nitrif    =  0.05      ! NH4 nitrification rate 
    103    xsirem    =  0.015     ! remineralization rate of Si 
    104    xlam1     =  0.005     ! scavenging rate of Iron 
    105    oxymin    =  1.E-6     ! Half-saturation constant for anoxia  
     159   xsirem    =  0.003     ! remineralization rate of Si 
     160   xsiremlab =  0.03      ! fast remineralization rate of Si 
     161   xsilab    =  0.5       ! Fraction of labile biogenic silica 
     162   oxymin    =  1.E-6     ! Half-saturation constant for anoxia 
    106163/ 
    107164!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    108165&nampiscal     !   parameters for Calcite chemistry 
    109166!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    110    kdca       =  0.327e3  ! calcite dissolution rate constant (1/time) 
     167   kdca       =  6.       ! calcite dissolution rate constant (1/time) 
    111168   nca        =  1.       ! order of dissolution reaction (dimensionless) 
    112169/ 
    113170!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    114 &nampissed     !   parameters for inputs deposition 
    115 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    116    ln_dustfer  =  .true.   ! boolean for dust input from the atmosphere 
    117    ln_river    =  .true.   ! AUTO boolean for river input of nutrients 
     171&nampissbc     !   parameters for inputs deposition 
     172!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     173!              !  file name        ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     174!              !                   !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     175   sn_dust     = 'dust.orca'       ,     -1            , 'dust'     ,  .true.      , .true. ,   'yearly'  , ''       , '' 
     176   sn_solub    = 'solubility.orca' ,    -12            , 'solubility1' ,  .false.      , .true. ,   'yearly'  , ''       , '' 
     177   sn_riverdic = 'river.orca'      ,    120            , 'riverdic' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     178   sn_riverdoc = 'river.orca'      ,    120            , 'riverdoc' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     179   sn_riverdin = 'river.orca'      ,    120            , 'riverdin' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     180   sn_riverdon = 'river.orca'      ,    120            , 'riverdon' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     181   sn_riverdip = 'river.orca'      ,    120            , 'riverdip' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     182   sn_riverdop = 'river.orca'      ,    120            , 'riverdop' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     183   sn_riverdsi = 'river.orca'      ,    120            , 'riverdsi' ,  .true.     , .true. ,   'yearly'  , ''       , '' 
     184   sn_ndepo    = 'ndeposition.orca',    -12            , 'ndep'     ,  .false.     , .true. ,   'yearly'  , ''       , '' 
     185   sn_ironsed  = 'bathy.orca'      ,    -12            , 'bathy'    ,  .false.     , .true. ,   'yearly'  , ''       , '' 
     186   sn_hydrofe  = 'hydrofe'         ,    -12            , 'epsdb'    ,  .false.     , .true. ,   'yearly'  , ''       , '' 
     187! 
     188   cn_dir      = './'      !  root directory for the location of the dynamical files 
     189   ln_dust     =  .true.   ! boolean for dust input from the atmosphere 
     190   ln_solub    =  .true.   ! boolean for variable solubility of atm. Iron 
     191   ln_river    =  .true.   ! boolean for river input of nutrients 
    118192   ln_ndepo    =  .true.   ! boolean for atmospheric deposition of N 
    119    ln_sedinput =  .true.   ! boolean for Fe input from sediments 
    120    sedfeinput  =  1E-9     ! Coastal release of Iron 
    121    dustsolub   =  0.014    ! Solubility of the dust 
     193   ln_ironsed  =  .true.   ! boolean for Fe input from sediments 
     194   ln_ironice  =  .true.   ! boolean for Fe input from sea ice 
     195   ln_hydrofe  =  .false.  ! boolean for from hydrothermal vents 
     196   sedfeinput  =  1.e-9    ! Coastal release of Iron 
     197   dustsolub   =  0.02     ! Solubility of the dust 
     198   wdust       =  2.0      ! Dust sinking speed  
     199   icefeinput  =  15.e-9   ! Iron concentration in sea ice 
     200   nitrfix     =  1.e-7    ! Nitrogen fixation rate 
     201   diazolight  =  50.      ! Diazotrophs sensitivity to light (W/m2) 
     202   concfediaz  =  1.e-10   ! Diazotrophs half-saturation Cste for Iron 
     203   hratio      =  9.e+5    ! Fe to 3He ratio assumed for vent iron supply  
    122204/ 
    123205!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     
    126208   xkr_eta      = 1.17    ! Sinking  exponent 
    127209   xkr_zeta     = 2.28    ! N content exponent 
     210   xkr_ncontent = 5.7E-6  ! N content factor 
    128211   xkr_mass_min = 0.0002  ! Minimum mass for Aggregates 
    129212   xkr_mass_max = 1.      ! Maximum mass for Aggregates 
     
    140223/ 
    141224!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    142 &nampisdia     !   additional 2D/3D tracers diagnostics ("key_trc_diaadd") 
    143 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    144    nwritedia   =  5475     !  time step frequency for tracers diagnostics 
    145 ! 
     225&nampisdia     !   additional 2D/3D tracers diagnostics  
     226!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    146227!              !    name   !           title of the field          !     units      ! 
    147228!              !           !                                       !                !   
     
    175256!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    176257   ln_pisdmp    =  .true.     !  Relaxation fo some tracers to a mean value 
    177    ln_pisclo    =  .true.    !  Restoring of tracer to initial value on closed sea ("key_dtatrc") 
    178 / 
     258   nn_pisdmp    =  5475       !  Frequency of Relaxation  
     259/ 
     260!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     261&nampismass     !  Mass conservation 
     262!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     263   ln_check_mass =  .true.    !  Check mass conservation 
     264/ 
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_top_ORCA2

    r1775 r2061  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/TOP1 :  1 - tracer definition                     (namtoptrc) 
    3 !! namelists    2 - dynamical tracer trends               (namtoptrd) 
    4 !!              6 - tracer advection                      (namtopadv) 
    5 !!              7 - tracer bottom boundary                (namtopbbl) 
    6 !!              8 - tracer lateral diffusion              (namtopldf) 
    7 !!              3 - tracer vertical physics               (namtopzdf) 
    8 !!              9 - tracer newtonian damping              (namtopdmp) 
     2!! NEMO/TOP1 :  1 - tracer definition                     (namtrc    ) 
     3!!              2 - tracer data initialisation            (namtrc_dta) 
     4!!              3 - tracer advection                      (namtrc_adv) 
     5!!              4 - tracer lateral diffusion              (namtrc_ldf) 
     6!!              5 - tracer vertical physics               (namtrc_zdf) 
     7!!              6 - tracer newtonian damping              (namtrc_dmp) 
     8!!              7 - dynamical tracer trends               (namtrc_trd) 
     9!!              8 - tracer output diagonstics             (namtrc_dia) 
    910!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1011!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    11 &namtoptrc     !   tracers definition 
     12&namtrc     !   tracers definition 
    1213!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    13    ndttrc      =  1        !  time step frequency for passive tracers       
    14    nwritetrc   =  5475     !  time step frequency for tracer outputs 
    15    ln_rsttr    = .false.   !  AUTO start from a restart file (T) or not (F) 
    16    nrsttr      =   0       !  AUTO restart control = 0 initial time step is not compared to the restart file value 
    17                            !                  = 1 do not use the value in the restart file 
    18                            !                  = 2 calendar parameters read in the restart file 
    19    cn_trcrst_in  = "restart_trc"   !  suffix of pass. tracer restart name (input) 
    20    cn_trcrst_out = "restart_trc"   !  suffix of pass. tracer restart name (output) 
     14   nn_dttrc      =  1        !  time step frequency for passive sn_tracers       
     15   nn_writetrc   =  5475     !  time step frequency for sn_tracer outputs 
     16   ln_top_euler  = .true.    !  use Euler time-stepping for TOP 
     17   ln_rsttr      = .false.   !  AUTO start from a restart file (T) or not (F) 
     18   nn_rsttr      =   0       !  AUTO restart control = 0 initial time step is not compared to the restart file value 
     19                             !                  = 1 do not use the value in the restart file 
     20                             !                  = 2 calendar parameters read in the restart file 
     21   cn_trcrst_in  = "restart_trc"   !  suffix of pass. sn_tracer restart name (input) 
     22   cn_trcrst_out = "restart_trc"   !  suffix of pass. sn_tracer restart name (output) 
     23   ln_trcdta     =   .true.  !  Initialisation from data input file (T) or not (F) 
     24   ln_trcdmp     =  .false.  !  add a damping termn (T) or not (F) 
     25   ln_trcdmp_clo =  .true.   !  restoring on close sea + baltic sea (T) or not (F) 
    2126! 
    22 !              !    name   !           title of the field              !   units    ! initial data ! save   ! 
    23 !              !           !                                           !            ! from file    ! or not !  
    24 !              !           !                                           !            ! or not       !        ! 
    25    tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'molC/L' ,  .true.     ,  .true. 
    26    tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true. 
    27    tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'molO2/L' ,  .true.     ,  .true. 
    28    tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'molC/L' ,  .false.    ,  .true. 
    29    tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'molC/L' ,  .true.     ,  .true. 
    30    tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'molC/L' ,  .false.    ,  .true. 
    31    tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'molSi/L' ,  .true.     ,  .true. 
    32    tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'molC/L' ,  .false.    ,  .true. 
    33    tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'molC/L' ,  .false.    ,  .true. 
    34    tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'molC/L'  ,  .false.    ,  .true. 
    35    tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'molC/L' ,  .false.    ,  .true. 
    36    tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'molC/L' ,  .false.    ,  .true. 
    37    tracer(13)  = 'BSi     ' , 'Diatoms Silicate Concentration         ',  'molC/L' ,  .false.    ,  .true. 
    38    tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'molFe/L' ,  .true.     ,  .true. 
    39    tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'molFe/L' ,  .false.    ,  .true. 
    40    tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'molC/L' ,  .false.    ,  .true. 
    41    tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'molFe/L' ,  .false.    ,  .true. 
    42    tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'molFe/L' ,  .false.    ,  .true. 
    43    tracer(19)  = 'DSi     ' , 'Sinking biogenic Silicate Concentration',  'molC/L' ,  .false.    ,  .true. 
    44    tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'molFe/L' ,  .false.    ,  .true. 
    45    tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'gChl/L' ,  .false.    ,  .true. 
    46    tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'gChl/L' ,  .false.    ,  .true. 
    47    tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'molC/L' ,  .true.     ,  .true. 
    48    tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'molC/L' ,  .false.    ,  .true. 
     27!                !    name   !           title of the field              ! initial data ! initial data ! save   ! 
     28!                !           !                                           !  units       ! from file    ! or not !  
     29!                !           !                                           !              ! or not       !        ! 
     30   sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true.     ,  .true. 
     31   sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true. 
     32   sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true.     ,  .true. 
     33   sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
     34   sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true.     ,  .true. 
     35   sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
     36   sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
     37   sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false.    ,  .true. 
     38   sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
     39   sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .true.     ,  .true. 
     40   sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
     41   sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
     42   sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
     43   sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true.     ,  .true. 
     44   sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
     45   sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
     46   sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
     47   sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false.    ,  .true. 
     48   sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false.    ,  .true. 
     49   sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false.    ,  .true. 
     50   sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
     51   sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
     52   sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
     53   sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false.    ,  .true. 
    4954/ 
    5055!----------------------------------------------------------------------- 
    51 &namtopadv    !   advection scheme for passive tracer  
     56&namtrc_dta      !    Initialisation from data input file 
    5257!----------------------------------------------------------------------- 
    53    ln_trcadv_cen2   =  .false.  !  2nd order centered scheme    
    54    ln_trcadv_tvd    =  .false.  !  TVD scheme    
    55    ln_trcadv_muscl  =  .true.   !  MUSCL scheme 
    56    ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    57    ln_trcadv_smolar =  .false.  !  SMOLAR scheme 
    58    rsc              =  1.       !  tuning coefficient for Smol-Car. scheme  
    59    ncortrc          =  1        !  number of corrective phases for Smol-Car. scheme  
    60    crosster         =  .false.  !  computes Smol-Car crossterms (T) or not (F) 
     58! 
     59!                !  file name               ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     60!                !                          !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     61   sn_trcdta(1)  = 'data_DIC_nomask'        ,        -12        ,  'DIC'     ,    .false.   , .true. , 'yearly'  , ''       , '' 
     62   sn_trcdta(2)  = 'data_Alkalini_nomask'   ,        -12        ,  'Alkalini',    .false.   , .true. , 'yearly'  , ''       , '' 
     63   sn_trcdta(3)  = 'data_O2_nomask'         ,        -1         ,  'O2'      ,    .true.    , .true. , 'yearly'  , ''       , '' 
     64   sn_trcdta(5)  = 'data_PO4_nomask'        ,        -1         ,  'PO4'     ,    .true.    , .true. , 'yearly'  , ''       , '' 
     65   sn_trcdta(7)  = 'data_Si_nomask'         ,        -1         ,  'Si'      ,    .true.    , .true. , 'yearly'  , ''       , '' 
     66   sn_trcdta(10) = 'data_DOC_nomask'        ,        -12        ,  'DOC'     ,    .false.   , .true. , 'yearly'  , ''       , '' 
     67   sn_trcdta(14) = 'data_Fer_nomask'        ,        -12        ,  'Fer'     ,    .false.   , .true. , 'yearly'  , ''       , '' 
     68   sn_trcdta(23) = 'data_NO3_nomask'        ,        -1         ,  'NO3'     ,    .true.    , .true. , 'yearly'  , ''       , '' 
     69! 
     70   cn_dir        =  './'      !  root directory for the location of the data files 
     71   rn_trfac(1)   =   1.0e-06  !  multiplicative factor 
     72   rn_trfac(2)   =   1.0e-06  !  -      -      -     - 
     73   rn_trfac(3)   =  44.6e-06  !  -      -      -     - 
     74   rn_trfac(5)   = 122.0e-06  !  -      -      -     - 
     75   rn_trfac(7)   =   1.0e-06  !  -      -      -     - 
     76   rn_trfac(10)  =   1.0      !  -      -      -     - 
     77   rn_trfac(14)  =   1.0      !  -      -      -     - 
     78   rn_trfac(23)  =   7.6e-06  !  -      -      -     - 
    6179/ 
    6280!----------------------------------------------------------------------- 
    63 &namtopbbl    !   bottom boundary layer scheme for passive tracer  
     81&namtrc_adv    !   advection scheme for passive tracer  
    6482!----------------------------------------------------------------------- 
    65    atrcbbl          = 1000.     !  lateral mixing coefficient in the bbl  [m2/s] 
     83   ln_trcadv_cen2    =  .false.  !  2nd order centered scheme    
     84   ln_trcadv_tvd     =  .false.  !  TVD scheme 
     85   ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     86   ln_trcadv_muscl2  =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
     87   ln_trcadv_ubs     =  .false.  !  UBS scheme 
     88   ln_trcadv_qck     =  .false.  !  QUICKEST scheme 
     89   ln_trcadv_msc_ups =  .false.  !  use upstream scheme within muscl 
    6690/ 
    6791!----------------------------------------------------------------------- 
    68 &namtopldf    !   lateral diffusion scheme for passive tracer  
     92&namtrc_ldf    !   lateral diffusion scheme for passive tracer  
    6993!----------------------------------------------------------------------- 
    7094   ln_trcldf_diff   =  .true.   !  performs lateral diffusion (T) or not (F) 
     
    77101   ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
    78102!                               !  Coefficient 
    79    ahtrc0      =  2000.         !     horizontal eddy diffusivity for tracers [m2/s] 
    80    ahtrb0      =     0.         !     background eddy diffusivity for ldf_iso [m2/s] 
    81    aeivtr0     =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_trcldf_eiv") 
    82    trcrat      =  1.            !  ratio betweeen passive and active tracer diffusion coeff 
     103   rn_ahtrc_0       =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     104   rn_ahtrb_0       =     0.    !     background eddy diffusivity for ldf_iso [m2/s] 
    83105/ 
    84106!----------------------------------------------------------------------- 
    85 &namtopzdf        !   vertical physics 
     107&namtrc_zdf        !   vertical physics 
    86108!----------------------------------------------------------------------- 
    87109   ln_trczdf_exp   =  .false.  !  split explicit (T) or implicit (F) time stepping 
    88    n_trczdf_exp    =   3       !  number of sub-timestep for ln_trczdfexp=T 
     110   nn_trczdf_exp   =   3       !  number of sub-timestep for ln_trczdfexp=T 
    89111/ 
    90112!----------------------------------------------------------------------- 
    91 &namtoprad        !  treatment of negative concentrations  
     113&namtrc_rad        !  treatment of negative concentrations  
    92114!----------------------------------------------------------------------- 
    93115   ln_trcrad   =  .true.  !  artificially correct negative concentrations (T) or not (F) 
    94116/ 
    95117!----------------------------------------------------------------------- 
    96 &namtopdmp        !   passive tracer newtonian damping                 ('key_trcdmp') 
     118&namtrc_dmp    !   passive tracer newtonian damping    
    97119!----------------------------------------------------------------------- 
    98    ndmptr      =   20      !  type of damping in passive tracers 
    99                            !     ='latitude', damping poleward of 'ndmp' degrees and function  
    100                            !                  of the distance-to-coast. Red and Med Seas as ndmptr=-1 
    101                            !     =-1 damping only in Med and Red Seas 
    102    ndmpftr     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
    103    nmldmptr    =    1      !  type of damping: =0 damping throughout the water column 
    104                            !                   =1 no damping in the mixed layer defined by avt >5cm2/s ) 
    105                            !                   =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    106    sdmptr      =   50.     !  surface time scale for internal damping (days) 
    107    bdmptr      =  360.     !  bottom  time scale for internal damping (days) 
    108    hdmptr      =  800.     !  depth of transition between sdmptr and bdmptr (meters) 
     120   nn_hdmp_tr  =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
     121                           !                   =XX, damping poleward of XX degrees (XX>0) 
     122                           !                      + F(distance-to-coast) + Red and Med Seas 
     123   nn_zdmp_tr  =    1      !  vertical   shape =0    damping throughout the water column 
     124                           !                   =1 no damping in the mixing layer (kz  criteria) 
     125                           !                   =2 no damping in the mixed  layer (rho crieria) 
     126   rn_surf_tr  =   50.     !  surface time scale of damping   [days] 
     127   rn_bot_tr   =  360.     !  bottom  time scale of damping   [days] 
     128   rn_dep_tr   =  800.     !  depth of transition between rn_surf and rn_bot [meters] 
     129   nn_file_tr  =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
    109130/ 
    110131!----------------------------------------------------------------------- 
    111 &namtoptrd       !   diagnostics on tracer trends 
     132&namtrc_trd       !   diagnostics on tracer trends        ('key_trdtrc') 
    112133!                          or mixed-layer trends          ('key_trdmld_trc') 
    113134!---------------------------------------------------------------------- 
    114    ntrd_trc   =  5475      !  time step frequency and tracers trends 
    115    nctls_trc =   0        !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    116    ucf_trc    =   1        !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     135   nn_trd_trc  =  5475      !  time step frequency and tracers trends 
     136   nn_ctls_trc =   0        !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     137   rn_ucf_trc  =   1        !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    117138   ln_trdmld_trc_restart = .false.  !  restart for ML diagnostics 
    118139   ln_trdmld_trc_instant = .true.  !  flag to diagnose trends of instantantaneous or mean ML T/S 
    119    luttrd(1)  =   .true. 
    120    luttrd(2)  =   .true. 
    121    luttrd(3)  =   .false. 
    122    luttrd(4)  =   .false. 
    123    luttrd(5)  =   .false. 
    124    luttrd(6)  =   .false. 
    125    luttrd(7)  =   .false. 
    126    luttrd(8)  =   .false. 
    127    luttrd(9)  =   .false. 
    128    luttrd(10) =   .false. 
    129    luttrd(11) =   .false. 
    130    luttrd(12) =   .false. 
    131    luttrd(13) =   .false. 
    132    luttrd(14) =   .false. 
    133    luttrd(15) =   .false. 
    134    luttrd(16) =   .false. 
    135    luttrd(17) =   .false. 
    136    luttrd(18) =   .false. 
    137    luttrd(19) =   .false. 
    138    luttrd(20) =   .false. 
    139    luttrd(21) =   .false. 
    140    luttrd(22) =   .false. 
    141    luttrd(23) =   .true. 
    142    luttrd(24) =   .false. 
     140   ln_trdtrc(1)  =   .true. 
     141   ln_trdtrc(2)  =   .true. 
     142   ln_trdtrc(23) =   .true. 
    143143/ 
     144!----------------------------------------------------------------------- 
     145&namtrc_dia       !   parameters for passive tracer additional diagnostics 
     146!---------------------------------------------------------------------- 
     147   ln_diatrc     =  .true.   !  save additional diag. (T) or not (F) 
     148   nn_writedia   =  5475     !  time step frequency for diagnostics 
     149/ 
Note: See TracChangeset for help on using the changeset viewer.