Ignore:
Timestamp:
03/21/23 14:42:28 (16 months ago)
Author:
cetlod
Message:

IPSLCM7_WORK : 1st step of coupling DYNAMICO-LMDZ-ORCHIDEE-NEMO4.2 + Oasis3-mct.5.0

Location:
CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM
Files:
3 added
10 deleted
3 edited

Legend:

Unmodified
Added
Removed

  • CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/namcouple_ORCA1xICO40

    r6329 r6346  
    5252#  
    5353O_SSTSST  SISUTESW 1 <freq_coupling>  2  sstoc.nc  <output_mode>  
    54 362 332 1 16002 torc  tico  LAG=<lag_oce> 
     54360 331 1 16002 torc  tico  LAG=<lag_oce> 
    5555P  2 P  0 
    5656LOCTRANS MAPPING 
    … …  
    6060# Mozaic: 1) mapping filename 2) connected unit 3) dataset rank 4) Maximum 
    6161#         number of overlapped neighbors 
    62 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     62rmp_torc_to_tico_TempIceAlb.nc dst 
    6363# CHECKOUT: indicate computation of global, land and sea field integrals. 
    6464# 
    … …  
    6868# 
    6969OIceFrc SIICECOV 44 <freq_coupling>  2  sstoc.nc <output_mode> 
    70 362 332 1 16002 torc  tico   LAG=<lag_oce> 
    71 P  2 P  0 
    72 # 
    73 LOCTRANS MAPPING 
    74  AVERAGE 
    75 # CHECKIN: indicate computation of global, land and sea field integrals. 
    76 # 
    77 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     70360 331 1 16002 torc  tico   LAG=<lag_oce> 
     71P  2 P  0 
     72# 
     73LOCTRANS MAPPING 
     74 AVERAGE 
     75# CHECKIN: indicate computation of global, land and sea field integrals. 
     76# 
     77rmp_torc_to_tico_TempIceAlb.nc dst 
    7878# CHECKOUT: indicate computation of global, land and sea field integrals. 
    7979# 
    … …  
    8484# 
    8585O_TepIce  SIICTEMW 34 <freq_coupling>  2  sstoc.nc  <output_mode> 
    86 362 332 1 16002 torc  tico   LAG=<lag_oce> 
    87 P  2 P  0 
    88 LOCTRANS MAPPING 
    89  AVERAGE 
    90 # CHECKIN: indicate computation of global, land and sea field integrals. 
    91 #  
    92 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     86360 331 1 16002 torc  tico   LAG=<lag_oce> 
     87P  2 P  0 
     88LOCTRANS MAPPING 
     89 AVERAGE 
     90# CHECKIN: indicate computation of global, land and sea field integrals. 
     91#  
     92rmp_torc_to_tico_TempIceAlb.nc dst 
    9393# CHECKOUT: indicate computation of global, land and sea field integrals. 
    9494# 
    … …  
    9898# 
    9999O_AlbIce  SIICEALW 17 <freq_coupling>  2  sstoc.nc  <output_mode> 
    100 362 332 1 16002 torc  tico  LAG=<lag_oce>   
    101 P  2 P  0 
    102 # 
    103 LOCTRANS MAPPING 
    104  AVERAGE 
    105 # CHECKIN: indicate computation of global, land and sea field integrals. 
    106 # 
    107  rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     100360 331 1 16002 torc  tico  LAG=<lag_oce>   
     101P  2 P  0 
     102# 
     103LOCTRANS MAPPING 
     104 AVERAGE 
     105# CHECKIN: indicate computation of global, land and sea field integrals. 
     106# 
     107 rmp_torc_to_tico_TempIceAlb.nc dst 
    108108# CHECKOUT: indicate computation of global, land and sea field integrals. 
    109109# 
    … …  
    113113# Field 5 : Current surface (o->a 5) 
    114114O_OCurx1 CURRENTX 321 <freq_coupling>  2  sstoc.nc  <output_mode> 
    115 362 332 1 16002 torc    tico  LAG=<lag_oce> 
    116 P  2 P  0 
    117 LOCTRANS MAPPING 
    118  AVERAGE 
    119 # CHECKIN: indicate computation of global, land and sea field integrals. 
    120 # 
    121 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     115360 331 1 16002 torc    tico  LAG=<lag_oce> 
     116P  2 P  0 
     117LOCTRANS MAPPING 
     118 AVERAGE 
     119# CHECKIN: indicate computation of global, land and sea field integrals. 
     120# 
     121rmp_torc_to_tico_TempIceAlb.nc dst 
    122122# CHECKOUT: indicate computation of global, land and sea field integrals. 
    123123# 
    … …  
    127127# Field 6 : Current surface (o->a 6) 
    128128O_OCury1 CURRENTY 321 <freq_coupling>  2  sstoc.nc  <output_mode> 
    129 362 332 1 16002 torc    tico  LAG=<lag_oce> 
    130 P  2 P  0 
    131 LOCTRANS MAPPING 
    132  AVERAGE 
    133 # CHECKIN: indicate computation of global, land and sea field integrals. 
    134 # 
    135 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     129360 331 1 16002 torc    tico  LAG=<lag_oce> 
     130P  2 P  0 
     131LOCTRANS MAPPING 
     132 AVERAGE 
     133# CHECKIN: indicate computation of global, land and sea field integrals. 
     134# 
     135rmp_torc_to_tico_TempIceAlb.nc dst 
    136136# CHECKOUT: indicate computation of global, land and sea field integrals. 
    137137#   
    … …  
    140140# Field 7 : Current surface (o->a 7) 
    141141O_OCurz1 CURRENTZ 321 <freq_coupling>  2  sstoc.nc  <output_mode> 
    142 362 332 1 16002 torc    tico  LAG=<lag_oce> 
    143 P  2 P  0 
    144 LOCTRANS MAPPING 
    145  AVERAGE 
    146 # CHECKIN: indicate computation of global, land and sea field integrals. 
    147 # 
    148 rmp_teORCA1.2_to_tICO40_TempIceAlb_1stOrder_v3.nc dst 
     142360 331 1 16002 torc    tico  LAG=<lag_oce> 
     143P  2 P  0 
     144LOCTRANS MAPPING 
     145 AVERAGE 
     146# CHECKIN: indicate computation of global, land and sea field integrals. 
     147# 
     148rmp_torc_to_tico_TempIceAlb.nc dst 
    149149# CHECKOUT: indicate computation of global, land and sea field integrals. 
    150150# 
    … …  
    159159# 
    160160COTAUXXU O_OTaux1 23 <freq_coupling>  1  flxat.nc   <output_mode> 
    161 1 16002 362 332 tico    uorc  LAG=<lag_atm>  
    162 P  0 P  2 
    163 MAPPING 
    164 # CHECKIN: indicate computation of global, land and sea field integrals. 
    165 # Interpolation method ou parametres mozaic 
    166 rmp_tICO40_to_ueORCA1.2_WindStress_2ndOrder_v3.nc src 
     1611 16002 360 331 tico    uorc  LAG=<lag_atm>  
     162P  0 P  2 
     163MAPPING 
     164# CHECKIN: indicate computation of global, land and sea field integrals. 
     165# Interpolation method ou parametres mozaic 
     166rmp_tico_to_uorc_WindStress.nc src 
    167167# CHECKOUT: indicate computation of global, land and sea field integrals. 
    168168# 
    … …  
    172172# 
    173173COTAUYYU O_OTauy1 23 <freq_coupling>  1  flxat.nc  <output_mode> 
    174 1 16002 362 332 tico    uorc  LAG=<lag_atm>  
    175 P  0 P  2 
    176 MAPPING 
    177 # CHECKIN: indicate computation of global, land and sea field integrals. 
    178 # Interpolation method ou parametres mozaic 
    179 rmp_tICO40_to_ueORCA1.2_WindStress_2ndOrder_v3.nc src 
     1741 16002 360 331 tico    uorc  LAG=<lag_atm>  
     175P  0 P  2 
     176MAPPING 
     177# CHECKIN: indicate computation of global, land and sea field integrals. 
     178# Interpolation method ou parametres mozaic 
     179rmp_tico_to_uorc_WindStress.nc src 
    180180# CHECKOUT: indicate computation of global, land and sea field integrals. 
    181181# 
    … …  
    185185# 
    186186COTAUZZU O_OTauz1 23 <freq_coupling>  1  flxat.nc  <output_mode> 
    187 1 16002 362 332 tico    uorc  LAG=<lag_atm>  
    188 P  0 P  2 
    189 MAPPING 
    190 # CHECKIN: indicate computation of global, land and sea field integrals. 
    191 # Interpolation method ou parametres mozaic 
    192 rmp_tICO40_to_ueORCA1.2_WindStress_2ndOrder_v3.nc src 
     1871 16002 360 331 tico    uorc  LAG=<lag_atm>  
     188P  0 P  2 
     189MAPPING 
     190# CHECKIN: indicate computation of global, land and sea field integrals. 
     191# Interpolation method ou parametres mozaic 
     192rmp_tico_to_uorc_WindStress.nc src 
    193193# CHECKOUT: indicate computation of global, land and sea field integrals. 
    194194# 
    … …  
    198198# 
    199199COTAUXXV O_OTaux2  24 <freq_coupling>  1    flxat.nc  <output_mode> 
    200 1 16002 362 332 tico    vorc  LAG=<lag_atm>  
     2001 16002 360 331 tico    vorc  LAG=<lag_atm>  
    201201P  0 P  2 
    202202#SCRIPR 
    … …  
    205205MAPPING 
    206206# Interpolation method or mozaic parameters 
    207 rmp_tICO40_to_veORCA1.2_WindStress_2ndOrder_v3.nc src 
     207rmp_tico_to_vorc_WindStress.nc src 
    208208# CHECKOUT: indicate computation of global, land and sea field integrals. 
    209209# 
    … …  
    213213# 
    214214COTAUYYV O_OTauy2  24 <freq_coupling>  1    flxat.nc  <output_mode> 
    215 1 16002 362 332 tico    vorc  LAG=<lag_atm> 
     2151 16002 360 331 tico    vorc  LAG=<lag_atm> 
    216216P  0 P  2 
    217217#SCRIPR 
    … …  
    219219MAPPING 
    220220# Interpolation method or mozaic parameters 
    221 rmp_tICO40_to_veORCA1.2_WindStress_2ndOrder_v3.nc src 
     221rmp_tico_to_vorc_WindStress.nc src 
    222222# 
    223223# CHECKOUT: indicate computation of global, land and sea field integrals. 
    … …  
    228228# 
    229229COTAUZZV O_OTauz2  24 <freq_coupling>  1    flxat.nc  <output_mode> 
    230 1 16002 362 332 tico    vorc  LAG=<lag_atm> 
     2301 16002 360 331 tico    vorc  LAG=<lag_atm> 
    231231P  0 P  2 
    232232#SCRIPR 
    … …  
    234234MAPPING 
    235235# Interpolation method or mozaic parameters 
    236 rmp_tICO40_to_veORCA1.2_WindStress_2ndOrder_v3.nc src 
     236rmp_tico_to_vorc_WindStress.nc src 
    237237# CHECKOUT: indicate computation of global, land and sea field integrals. 
    238238# 
    … …  
    241241# 
    242242COWINDSP O_Wind10  56 <freq_coupling>  1    flxat.nc  <output_mode> 
    243 1 16002 362 332 tico    torc  LAG=<lag_atm> 
    244 P  0 P  2 
    245 MAPPING 
    246 # CHECKIN: indicate computation of global, land and sea field integrals. 
    247 # Interpolation method ou parametres mozaic 
    248 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     2431 16002 360 331 tico    torc  LAG=<lag_atm> 
     244P  0 P  2 
     245MAPPING 
     246# CHECKIN: indicate computation of global, land and sea field integrals. 
     247# Interpolation method ou parametres mozaic 
     248rmp_tico_to_torc_HeatWaterFluxes.nc src 
    249249# CHECKOUT: indicate computation of global, land and sea field integrals. 
    250250#  
    … …  
    254254# 
    255255COTOTRAI OTotRain 26 <freq_coupling>  1   flxat.nc   <output_mode> 
    256 1 16002 362 332 tico    torc  LAG=<lag_atm>  
    257 P  0 P  2 
    258 MAPPING 
    259 # CHECKIN: indicate computation of global, land and sea field integrals. 
    260 # Interpolation method ou parametres mozaic 
    261 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     2561 16002 360 331 tico    torc  LAG=<lag_atm>  
     257P  0 P  2 
     258MAPPING 
     259# CHECKIN: indicate computation of global, land and sea field integrals. 
     260# Interpolation method ou parametres mozaic 
     261rmp_tico_to_torc_HeatWaterFluxes.nc src 
    262262# CHECKOUT: indicate computation of global, land and sea field integrals. 
    263263# 
    … …  
    267267# 
    268268COTOTSNO  OTotSnow 28 <freq_coupling>  1   flxat.nc   <output_mode> 
    269 1 16002 362 332 tico    torc  LAG=<lag_atm>  
    270 P  0 P  2 
    271 MAPPING 
    272 # CHECKIN: indicate computation of global, land and sea field integrals. 
    273 # Interpolation method ou parametres mozaic 
    274 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     2691 16002 360 331 tico    torc  LAG=<lag_atm>  
     270P  0 P  2 
     271MAPPING 
     272# CHECKIN: indicate computation of global, land and sea field integrals. 
     273# Interpolation method ou parametres mozaic 
     274rmp_tico_to_torc_HeatWaterFluxes.nc src 
    275275# CHECKOUT: indicate computation of global, land and sea field integrals. 
    276276# 
    … …  
    280280# 
    281281COTOTEVA  OTotEvap 25 <freq_coupling>  1   flxat.nc   <output_mode> 
    282 1 16002 362 332 tico    torc  LAG=<lag_atm>  
    283 P  0 P  2 
    284 MAPPING 
    285 # CHECKIN: indicate computation of global, land and sea field integrals. 
    286 # Interpolation method ou parametres mozaic 
    287 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     2821 16002 360 331 tico    torc  LAG=<lag_atm>  
     283P  0 P  2 
     284MAPPING 
     285# CHECKIN: indicate computation of global, land and sea field integrals. 
     286# Interpolation method ou parametres mozaic 
     287rmp_tico_to_torc_HeatWaterFluxes.nc src 
    288288# CHECKOUT: indicate computation of global, land and sea field integrals. 
    289289# 
    … …  
    293293# 
    294294COICEVAP OIceEvap 41 <freq_coupling>  1   flxat.nc   <output_mode> 
    295 1 16002 362 332 tico    torc  LAG=<lag_atm>  
    296 P  0 P  2 
    297 MAPPING 
    298 # CHECKIN: indicate computation of global, land and sea field integrals. 
    299 # Interpolation method ou parametres mozaic 
    300 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     2951 16002 360 331 tico    torc  LAG=<lag_atm>  
     296P  0 P  2 
     297MAPPING 
     298# CHECKIN: indicate computation of global, land and sea field integrals. 
     299# Interpolation method ou parametres mozaic 
     300rmp_tico_to_torc_HeatWaterFluxes.nc src 
    301301# CHECKOUT: indicate computation of global, land and sea field integrals. 
    302302# 
    … …  
    306306# 
    307307COQSRMIX O_QsrMix  7 <freq_coupling> 1   flxat.nc  <output_mode> 
    308 1 16002 362 332 tico    torc  LAG=<lag_atm> 
     3081 16002 360 331 tico    torc  LAG=<lag_atm> 
    309309P  0 P  2 
    310310MAPPING 
    311311# CHECKIN: indicate computation of global, land and sea field integrals. 
    312312# Interpolation method or mozaic parameters 
    313 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     313rmp_tico_to_torc_HeatWaterFluxes.nc src 
    314314# CHECKOUT: indicate computation of global, land and sea field integrals. 
    315315# 
    … …  
    319319# 
    320320COQNSMIX O_QnsMix 6 <freq_coupling>  1   flxat.nc  <output_mode> 
    321 1 16002 362 332 tico    torc   LAG=<lag_atm> 
    322 P  0 P  2 
    323 MAPPING 
    324 # CHECKIN: indicate computation of global, land and sea field integrals. 
    325 # Interpolation method ou parametres mozaic 
    326 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     3211 16002 360 331 tico    torc   LAG=<lag_atm> 
     322P  0 P  2 
     323MAPPING 
     324# CHECKIN: indicate computation of global, land and sea field integrals. 
     325# Interpolation method ou parametres mozaic 
     326rmp_tico_to_torc_HeatWaterFluxes.nc src 
    327327# CHECKOUT: indicate computation of global, land and sea field integrals. 
    328328# 
    … …  
    332332# 
    333333COSHFICE O_QsrIce  7 <freq_coupling> 1   flxat.nc  <output_mode> 
    334 1 16002 362 332 tico    torc  LAG=<lag_atm> 
     3341 16002 360 331 tico    torc  LAG=<lag_atm> 
    335335P  0 P  2 
    336336MAPPING 
    337337# CHECKIN: indicate computation of global, land and sea field integrals. 
    338338# Interpolation method or mozaic parameters 
    339 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     339rmp_tico_to_torc_HeatWaterFluxes.nc src 
    340340# CHECKOUT: indicate computation of global, land and sea field integrals. 
    341341# 
    … …  
    345345# 
    346346CONSFICE O_QnsIce 6 <freq_coupling>  1  flxat.nc  <output_mode> 
    347 1 16002 362 332 tico    torc  LAG=<lag_atm> 
    348 P  0 P  2 
    349 MAPPING 
    350 # CHECKIN: indicate computation of global, land and sea field integrals. 
    351 # Interpolation method ou parametres mozaic 
    352 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     3471 16002 360 331 tico    torc  LAG=<lag_atm> 
     348P  0 P  2 
     349MAPPING 
     350# CHECKIN: indicate computation of global, land and sea field integrals. 
     351# Interpolation method ou parametres mozaic 
     352rmp_tico_to_torc_HeatWaterFluxes.nc src 
    353353# CHECKOUT: indicate computation of global, land and sea field integrals. 
    354354# 
    … …  
    358358# 
    359359CODFLXDT O_dQnsdT  35 <freq_coupling>  1   flxat.nc  <output_mode> 
    360 1 16002 362 332 tico    torc  LAG=<lag_atm> 
    361 P  0 P  2 
    362 MAPPING 
    363 # CHECKIN: indicate computation of global, land and sea field integrals. 
    364 # Interpolation method ou parametres mozaic 
    365 rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_2ndOrder_v3.nc src 
     3601 16002 360 331 tico    torc  LAG=<lag_atm> 
     361P  0 P  2 
     362MAPPING 
     363# CHECKIN: indicate computation of global, land and sea field integrals. 
     364# Interpolation method ou parametres mozaic 
     365rmp_tico_to_torc_HeatWaterFluxes.nc src 
    366366# CHECKOUT: indicate computation of global, land and sea field integrals. 
    367367# 
    … …  
    371371# 
    372372COCALVIN OCalving  36 <freq_coupling_roff_calv>  2  icbrg.nc  <output_mode> 
    373 1 16002 362 332 tico    torc  LAG=<lag_atm>  
     3731 16002 360 331 tico    torc  LAG=<lag_atm>  
    374374P  0 P  2 
    375375LOCTRANS MAPPING 
    … …  
    377377# CHECKIN: indicate computation of global, land and sea field integrals. 
    378378# Interpolation method ou parametres mozaic 
    379 rmp_tICO40_to_teORCA1.2_calving_nosouth_v3.nc dst 
     379rmp_tico_to_torc_calving_nosouth.nc dst 
    380380# CHECKOUT: indicate computation of global, land and sea field integrals. 
    381381# 
    … …  
    385385# 
    386386COCALVIN OIceberg  36 <freq_coupling_roff_calv>  3  icbrg.nc  <output_mode> 
    387 1 16002 362 332 tico    torc  LAG=<lag_atm>  
     3871 16002 360 331 tico    torc  LAG=<lag_atm>  
    388388P  0 P  2 
    389389LOCTRANS MAPPING BLASNEW 
    … …  
    391391# CHECKIN: indicate computation of global, land and sea field integrals. 
    392392# Interpolation method ou parametres mozaic 
    393 rmp_tICO40_to_teORCA1.2_calving_iceberg_v3.nc dst 
     393rmp_tico_to_torc_calving_iceberg.nc dst 
    3943940.5 0 
    395395# CHECKOUT: indicate computation of global, land and sea field integrals. 
    … …  
    399399# 
    400400COCALVIN OIcshelf  36 <freq_coupling_roff_calv>  3  icshf.nc  <output_mode> 
    401 1 16002 362 332 tico    torc  LAG=<lag_atm>  
     4011 16002 360 331 tico    torc  LAG=<lag_atm>  
    402402P  0 P  2 
    403403LOCTRANS MAPPING BLASNEW 
    … …  
    405405# CHECKIN: indicate computation of global, land and sea field integrals. 
    406406# Interpolation method ou parametres mozaic 
    407 rmp_tICO40_to_teORCA1.2_calving_iceshelf_v3.nc dst 
     407rmp_tico_to_torc_calving_iceshelf.nc dst 
    4084080.5 0 
    409409# CHECKOUT: indicate computation of global, land and sea field integrals. 
    … …  
    413413# 
    414414COLIQRUN  O_Runoff 32 <freq_coupling_roff_calv>  4  flxat.nc   <output_mode> 
    415 1 16002 362 332 oico  torc  LAG=<lag_atm>  
     4151 16002 360 331 oico  torc  LAG=<lag_atm>  
    416416P  0 P  2 
    417417LOCTRANS MAPPING CONSERV BLASNEW 
    … …  
    419419# Interpolation method ou parametres mozaic 
    420420# weights convert from kg/s to kg/m^2/s 
    421 rmp_tICO40_to_teORCA1.2_runoff_Quantity_to_Surfacic_v3.nc dst 
    422 #rmp_tICO40_to_teORCA1.2_Quantity_v3.nc src 
    423 #rmp_tICO40_to_teORCA1.2_HeatWaterFluxes_v3.nc src 
     421rmp_tico_to_torc_runoff.nc dst 
    424422# CONSERV 
    425423GLOBAL bfb 
    … …  
    432430# 
    433431##COTAUMOD  O_TauMod 466 <freq_coupling>  6   flxat.nc   <output_mode> 
    434 ##tlmd    torc  LAG=<lag_atm>  
     432##tico    torc  LAG=<lag_atm>  
    435433##P  0 P  2 
    436434##INVERT CHECKIN MASK EXTRAP INTERP CHECKOUT 

  • CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/namelist_ORCA1_cfg

    r5479 r6346  
    1  
    21!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    3 !! NEMO/OPA  Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref 
     2!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref 
    43!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     4!!             ORCA2 - ICE - PISCES configuration                     !! 
     5!!====================================================================== 
     6!!              ***  Domain & Run management namelists  ***           !! 
     7!!                                                                    !! 
     8!!   namrun       parameters of the run 
     9!!   namdom       space and time domain 
     10!!   namcfg       parameters of the configuration                       (default: user defined GYRE) 
     11!!   namwad       Wetting and drying                                    (default: OFF) 
     12!!   namtsd       data: temperature & salinity                          (default: OFF) 
     13!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
     14!!   namc1d       1D configuration options                              ("key_c1d") 
     15!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
     16!!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
     17!!====================================================================== 
    518! 
    619!----------------------------------------------------------------------- 
    … …  
    1326   nn_leapy    = _AUTOBLOCKER_    !  Leap year calendar (1) or not (0) 
    1427   ln_rstart   = _AUTOBLOCKER_    !  start from rest (F) or from a restart file (T) 
    15    ln_rstart_ts = _AUTOBLOCKER_   !  start from rest for current only (F) or from a restart file (T) 
    1628   nn_rstctl   = _AUTOBLOCKER_    !  Restart control => activated only if ln_rstart = T 
    1729                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    … …  
    2234   cn_ocerst_out = "restart"      !  Suffix of ocean restart name (output) 
    2335   cn_ocerst_outdir = "."      !  directory in which to write output ocean restarts 
    24    nn_istate   =       0   !  Output the initial state (1) or not (0) 
    2536   nn_stock    =  _AUTOBLOCKER_   !  Frequency of creation of a restart file (modulo referenced to 1) 
    26    nn_write    =    5475   !  Requency of write in the output file   (modulo referenced to nn_it000) 
    2737   ln_mskland  = .true.    !  Masks land points in NetCDF outputs 
    28    ln_mskutil  = .true.    !  Outputs without halos 
    2938   ln_cfmeta   = .true.    !  output additional data to netCDF files required for compliance with the CF metadata standard 
    3039/ 
    3140!----------------------------------------------------------------------- 
    32 &namcfg        !   parameters of the configuration 
    33 !----------------------------------------------------------------------- 
    34    cp_cfg      =  "orca"               !  name of the configuration 
    35    jp_cfg      =       1               !  resolution of the configuration 
    36    jpidta      =     362               !  1st lateral dimension ( >= jpi ) 
    37    jpjdta      =     332               !  2nd    "         "    ( >= jpj ) 
    38    jpkdta      =      75               !  number of levels      ( >= jpk ) 
    39    jpiglo      =     362               !  1st dimension of global domain --> i =jpidta 
    40    jpjglo      =     332               !  2nd    -                  -    --> j  =jpjdta 
    41    jperio      =       6               !  lateral cond. type (between 0 and 6) 
    42 / 
    43 !----------------------------------------------------------------------- 
    44 &namzgr        !   vertical coordinate 
    45 !----------------------------------------------------------------------- 
    46 / 
    47 !----------------------------------------------------------------------- 
    48 &namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
    49 !----------------------------------------------------------------------- 
    50 / 
    51 !----------------------------------------------------------------------- 
    52 &namdom        !   space and time domain (bathymetry, mesh, timestep) 
    53 !----------------------------------------------------------------------- 
    54    nn_closea    =   1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    55    ln_clodyn    =  _AUTOBLOCKER_      !   
    56    ! 
    57    jphgr_msh   =       0               !  type of horizontal mesh 
    58    ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    59    ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
    60    ppe1_deg    =  999999.0             !  zonal      grid-spacing (degrees) 
    61    ppe2_deg    =  999999.0             !  meridional grid-spacing (degrees) 
    62    ppe1_m      =  999999.0             !  zonal      grid-spacing (degrees) 
    63    ppe2_m      =  999999.0             !  meridional grid-spacing (degrees) 
    64    ppsur       =   -3958.951371276829  !  ORCA r4, r2 and r05 coefficients 
    65    ppa0        =     103.9530096000000 ! (default coefficients) 
    66    ppa1        =       2.415951269000000   ! 
    67    ppkth       =      15.35101370000000    ! 
    68    ppacr       =       7.0             ! 
    69    ppdzmin     =  999999.0             !  Minimum vertical spacing 
    70    pphmax      =  999999.0             !  Maximum depth 
    71    ppa2        =     100.7609285000000 !  Double tanh function parameters 
    72    ppkth2      =      48.02989372000000    ! 
    73    ppacr2      =      13.00000000000   ! 
    74    rn_rdt      = 2700.     !  time step for the dynamics (and tracer if nn_acc=0) 
    75    rn_hmin     =   20. 
    76    nn_msh      = _AUTO_    !  AUTO - Create (=1) a mesh file or not (=0) 
    77 / 
    78 !----------------------------------------------------------------------- 
    79 &namsplit        
    80 !-----------------------------------------------------------------------  
    81    ln_bt_fw      =    .FALSE.          !  leap-frog integration of barotropic equations 
    82    ln_bt_av      =    .TRUE.           !  Time filtering of barotropic variables 
    83    ln_bt_nn_auto =    .TRUE.           !  Set nn_baro automatically to be just below 
    84                                        !  a user defined maximum courant number (rn_bt_cmax) 
    85    nn_baro       =    30               !  Number of iterations of barotropic mode 
    86                                        !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
    87    rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T  
    88    nn_bt_flt     =    1                !  Time filter choice 
    89                                        !  = 0 None 
    90                                        !  = 1 Boxcar over   nn_baro barotropic steps 
    91                                        !  = 2 Boxcar over 2*nn_baro     "         
    92 / 
    93 !----------------------------------------------------------------------- 
    94 &namcrs        !   Grid coarsening for dynamics output and/or 
    95                !   passive tracer coarsened online simulations 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    99 &namtsd    !   data : Temperature  & Salinity 
    100 !----------------------------------------------------------------------- 
    101    ln_tsd_tradmp = .false.  !  damping of ocean T & S toward T &S input data (T) or not (F) 
    102    sn_tem  = 'conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'votemper' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc'  ,   ''    ,    '' 
    103    sn_sal  = 'absolute_salinity_WOA13_decav_Reg1L75_clim'       , -1 ,'vosaline' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc'  ,   ''    ,    '' 
    104 / 
    105 !----------------------------------------------------------------------- 
    106 &namsbc        !   Surface Boundary Condition (surface module) 
    107 !----------------------------------------------------------------------- 
    108    nn_fsbc     =  2        !  frequency of surface boundary condition computation 
    109                            !     (also = the frequency of sea-ice model call) 
    110    ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     41&namdom        !   time and space domain 
     42!----------------------------------------------------------------------- 
     43   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
     44   ! 
     45   rn_Dt      = 2700.     !  time step for the dynamics and tracer 
     46   ln_meshmask =  _AUTO_ 
     47/ 
     48!----------------------------------------------------------------------- 
     49&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
     50!----------------------------------------------------------------------- 
     51   ln_read_cfg = .true.    !  (=T) read the domain configuration file 
     52      cn_domcfg = "domain_cfg.nc"    ! domain configuration filename 
     53      ! 
     54      ln_closea    = .true.    !  F => suppress closed seas (defined by closea_mask field)  
     55      !                         !       from the bathymetry at runtime. 
     56/ 
     57!----------------------------------------------------------------------- 
     58&namclo        !   parameters of the closed sea (cs) behavior                (default: OFF) 
     59!----------------------------------------------------------------------- 
     60   ln_maskcs = .false.        ! (=T) cs are masked ; So, in this case ln_mask_csundef and ln_clo_rnf have no effect. 
     61      !                       ! (=F => set ln_mask_csundef and ln_clo_rnf) 
     62      !                       ! cs masks are read and net evap/precip over closed sea spread out depending on domain_cfg.nc masks. 
     63      !                       ! See ln_mask_csundef and ln_clo_rnf for specific option related to this case 
     64      ! 
     65      ln_mask_csundef = .true.   ! (=T) undefined closed seas are masked ; 
     66      !                          ! (=F) undefined closed seas are kept and no specific treatment is done for these closed seas 
     67      ! 
     68      ln_clo_rnf = .true.        ! (=T) river mouth specified in domain_cfg.nc masks (rnf and emp case) are added to the runoff mask. 
     69      !                          !      allow the treatment of closed sea outflow grid-points to be the same as river mouth grid-points 
     70/ 
     71!----------------------------------------------------------------------- 
     72&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF) 
     73!----------------------------------------------------------------------- 
     74   !                       ! =T  read T-S fields for: 
     75   ln_tsd_init = .true.          !  ocean initialisation 
     76   ln_tsd_dmp  = .false.          !  T-S restoring   (see namtra_dmp) 
     77 
     78   cn_dir = './'     !  root directory for the T-S data location 
     79   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     80   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     81   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     82   sn_tem  = 'conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'votemper' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_bilinear.nc'  ,   ''    ,    '' 
     83   sn_sal  = 'absolute_salinity_WOA13_decav_Reg1L75_clim'       , -1 ,'vosaline' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_bilinear.nc'  ,   '' 
     84/ 
     85!!====================================================================== 
     86!!            ***  Surface Boundary Condition namelists  ***          !! 
     87!!                                                                    !! 
     88!!   namsbc          surface boundary condition manager                 (default: NO selection) 
     89!!   namsbc_flx      flux               formulation                     (ln_flx     =T) 
     90!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T) 
     91!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
     92!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only) 
     93!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   ) 
     94!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
     95!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
     96!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
     97!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
     98!!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
     99!!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
     100!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
     101!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     102!!====================================================================== 
     103! 
     104!----------------------------------------------------------------------- 
     105&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
     106!----------------------------------------------------------------------- 
     107   nn_fsbc     = 2         !  frequency of SBC module call 
     108                           !     (also = the frequency of sea-ice & iceberg model call) 
     109                     ! Type of air-sea fluxes  
    111110   ln_cpl      = .true.    !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    112    nn_limflx =   2         !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
    113                            !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
    114                            !  = 0  Average per-category fluxes (forced and coupled mode) 
    115                            !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
    116                            !  = 2  Redistribute a single flux over categories (coupled mode only) 
    117    nn_ice_embd = 1         !  AUTO -  
    118                            !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    119                            !  =1 levitating ice with mass and salt exchange but no presure effect 
    120                            !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    121    ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    122    ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    123    nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
    124    nn_isf      = 3         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
    125                            !  3 = rnf file for isf 
    126 !----------------------------------------------------------------------- 
    127 &namsbc_core   !   namsbc_core  CORE bulk formulae 
    128 !----------------------------------------------------------------------- 
    129 / 
    130 !----------------------------------------------------------------------- 
    131 &namtra_qsr    !   penetrative solar radiation 
    132 !----------------------------------------------------------------------- 
    133    sn_chl      ='merged_ESACCI_BIOMER4V1R1_CHL_REG05',  -1  , 'CHLA' , .true. , .true. , 'yearly' , 'weights_reg05_2_eorca1_bilinear.nc' , '' , '' 
    134    ln_traqsr   = .true.    !  Light penetration (T) or not (F) 
    135    ln_qsr_rgb  = .false.    !  RGB (Red-Green-Blue) light penetration 
    136    ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    137    ln_qsr_bio  = .true.   !  bio-model light penetration 
    138 / 
    139 !----------------------------------------------------------------------- 
    140 &namsbc_rnf    !   runoffs namelist surface boundary condition 
    141 !----------------------------------------------------------------------- 
    142 !              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    143 !              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    144    sn_rnf      = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    145    sn_i_rnf    = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc',        -1         , 'Icb_flux',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    146    sn_cnf      = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc',         0         , 'socoeff' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    147    sn_s_rnf    = 'runoffs'                         ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    148    sn_t_rnf    = 'runoffs'                         ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    149    sn_dep_rnf  = 'runoffs_eORCA1.0_depths.nc'      ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
     111                     ! Sea-ice : 
     112   nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
     113                     ! Misc. options of sbc :  
     114   ln_traqsr   = .true.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
     115/ 
     116!----------------------------------------------------------------------- 
     117&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
     118!----------------------------------------------------------------------- 
     119   nn_cplmodel       =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
     120   ln_usecplmask     = .false. !  use a coupling mask file to merge data received from several models 
     121   !                           !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
     122   ln_scale_ice_flux = .false. !  use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 
     123   nn_cats_cpl       =     1   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
     124   !_____________!__________________________!____________!_____________!______________________!________! 
     125   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
     126   !             !                          ! categories !  reference  !     orientation      ! grids  ! 
     127!***   send    *** 
     128   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     129   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     130   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     131   sn_snd_crt    =   'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward' ,  'T' 
     132   sn_snd_co2    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     133   sn_snd_crtw   =   'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V' 
     134   sn_snd_ifrac  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     135   sn_snd_wlev   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     136   sn_snd_cond   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     137   sn_snd_thick1 =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     138   sn_snd_mpnd   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     139   sn_snd_sstfrz =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     140   sn_snd_ttilyr =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     141!***  receive  *** 
     142   sn_rcv_w10m   =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     143   sn_rcv_taumod =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     144   sn_rcv_tau    =   'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward' ,  'U,V' 
     145   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     146   sn_rcv_qsr    =   'conservative'         ,    'no'    ,     ''      ,         ''           ,   '' 
     147   sn_rcv_qns    =   'conservative'         ,    'no'    ,     ''      ,         ''           ,   '' 
     148   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''           ,   '' 
     149   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     150   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     151   sn_rcv_co2    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     152   sn_rcv_iceflx =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     153   sn_rcv_mslp   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     154   sn_rcv_ts_ice =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     155   sn_rcv_isf    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     156   sn_rcv_icb    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     157   sn_rcv_hsig   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     158   sn_rcv_phioc  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     159   sn_rcv_sdrfx  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     160   sn_rcv_sdrfy  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     161   sn_rcv_wper   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     162   sn_rcv_wnum   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     163   sn_rcv_wstrf  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     164   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     165   sn_rcv_charn  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     166   sn_rcv_taw    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'U,V' 
     167   sn_rcv_bhd    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     168   sn_rcv_tusd   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     169   sn_rcv_tvsd   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   'T' 
     170/ 
     171!----------------------------------------------------------------------- 
     172&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T) 
     173!----------------------------------------------------------------------- 
     174   ln_qsr_bio  = .false.   !  bio-model light penetration 
     175   !                       ! type of penetration                        (default: NO selection) 
     176   ln_qsr_rgb  = .true.       !  RGB light penetration (Red-Green-Blue) 
     177   ! 
     178   nn_chldta   =      1       !  RGB : Chl data (=1) or cst value (=0) 
    150179 
    151    ln_rnf_icb   = .false.    !  read in iceberg flux 
    152    ln_rnf_mouth = .false.    !  specific treatment at rivers mouths 
    153    ln_rnf_depth = .true.     !  read in depth information for runoff 
    154    ln_rnf_tem   = .false.    !  read in temperature information for runoff 
    155    ln_rnf_sal   = .false.    !  read in salinity information for runoff 
    156    ln_rnf_depth_ini = .false.!  compute depth at initialisation from runoff file 
    157    rn_rnf_max   = 0.05       !  max value of the runoff climatology over global domain ( if ln_rnf_depth_ini = .true ) 
    158    rn_dep_max = 150.         !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
    159    nn_rnf_depth_file =  _AUTO_   ! create (=1) a runoff depth file or not (=0) 
    160 / 
    161 !----------------------------------------------------------------------- 
    162 &namsbc_isf    !  Top boundary layer (ISF) 
    163 !----------------------------------------------------------------------- 
    164 !              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    165 !              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
    166 !              ! 
    167    sn_rnfisf     = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
    168    sn_depmax_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12      ,'sodepmax_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
    169    sn_depmin_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12      ,'sodepmin_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
    170 / 
    171 !----------------------------------------------------------------------- 
    172 &namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
    173 !----------------------------------------------------------------------- 
    174 / 
    175 !----------------------------------------------------------------------- 
    176 &namsbc_ssr    !   surface boundary condition : sea surface restoring 
    177 !----------------------------------------------------------------------- 
    178 !              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    179 !              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
    180    sn_sss      = 'sss_absolute_salinity_WOA13_decav_Reg1L75_clim', -1. , 'sosaline', .true. , .true. , 'yearly' , 'weights_WOA13d1_2_eorca1_bilinear.nc' , '' 
    181 / 
    182 !----------------------------------------------------------------------- 
    183 &namsbc_alb    !   albedo parameters 
    184 !----------------------------------------------------------------------- 
    185    nn_ice_alb   =    1   !  parameterization of ice/snow albedo 
    186                          !     0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo 
    187                          !     1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004), 
    188                          !        giving cloud-sky albedo 
    189    rn_alb_sdry  =  0.87  !  dry snow albedo         : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) 
    190    rn_alb_smlt  =  0.82  !  melting snow albedo     : 0.65 ( '' )          ; 0.75 ( '' )          ; obs 0.72-0.82 ( '' ) 
    191    rn_alb_idry  =  0.65  !  dry ice albedo          : 0.72 ( '' )          ; 0.60 ( '' )          ; obs 0.54-0.65 ( '' ) 
    192    rn_alb_imlt  =  0.58  !  bare puddled ice albedo : 0.53 ( '' )          ; 0.50 ( '' )          ; obs 0.49-0.58 ( '' ) 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 &namsbc_cpl    !   coupling parameters 
    196 !----------------------------------------------------------------------- 
    197 !                    !     description       !  multiple  !    vector   !      vector          ! vector ! 
    198 !                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    199 ! send 
    200 sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    201 sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    202 sn_snd_thick  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    203 sn_snd_crt    =       'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward' ,  'T' 
    204 sn_snd_co2    =       'none'              ,    'no'    ,     ''      ,         ''           ,   '' 
    205 ! receive 
    206 sn_rcv_w10m   =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    207 sn_rcv_taumod =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    208 sn_rcv_tau    =       'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    209 sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    210 sn_rcv_qsr    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    211 sn_rcv_qns    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    212 sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    213 sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    214 sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    215 sn_rcv_co2    =       'none'              ,    'no'    ,     ''      ,         ''          ,   '' 
    216 sn_rcv_icb    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    217 sn_rcv_isf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     180   cn_dir = './'  !  root directory for the chlorophyl data location 
     181   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     182   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     183   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     184  sn_chl      ='merged_ESACCI_BIOMER4V1R1_CHL_REG05',  -1  , 'CHLA' , .true. , .true. , 'yearly' , 'weights_reg05_2_bilinear.nc' , '' , '' 
     185/ 
     186!----------------------------------------------------------------------- 
     187&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
     188!----------------------------------------------------------------------- 
     189/ 
     190!----------------------------------------------------------------------- 
     191&namsbc_rnf    !   runoffs                                              (ln_rnf =T) 
     192!----------------------------------------------------------------------- 
     193   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths 
     194      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
     195      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
     196   rn_rfact    =   1.e0    !  multiplicative factor for runoff 
     197   ln_rnf_depth = .false.   !  read in depth information for runoff 
     198   ln_rnf_tem  = .false.   !  read in temperature information for runoff 
     199   ln_rnf_sal  = .false.   !  read in salinity information for runoff 
     200   ln_rnf_icb  = .false.   !  read iceberg flux 
     201   ln_rnf_depth_ini = .true. ! compute depth at initialisation from runoff file 
     202      rn_rnf_max  = 0.05  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
     203      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
     204      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0) 
    218205 
    219 / 
    220 !----------------------------------------------------------------------- 
    221 &namberg       !   iceberg parameters 
    222 !----------------------------------------------------------------------- 
    223       ln_icebergs              = .false. 
    224       ln_bergdia               = .false.              ! Calculate budgets 
    225       nn_verbose_level         = 0                    ! Turn on more verbose output if level > 0 
    226       nn_verbose_write         = 120                  ! Timesteps between verbose messages 
    227       nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
    228                                                       ! Initial mass required for an iceberg of each class 
    229       rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    230                                                       ! Proportion of calving mass to apportion to each class   
    231       rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    232                                                       ! Ratio between effective and real iceberg mass (non-dim) 
    233                                                       ! i.e. number of icebergs represented at a point          
    234       rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    235                                                       ! thickness of newly calved bergs (m) 
    236       rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    237       rn_rho_bergs             = 850.                 ! Density of icebergs 
    238       rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
    239       ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
    240       rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    241       rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    242       ln_passive_mode          = .false.              ! iceberg - ocean decoupling    
    243       nn_test_icebergs         =   8                  ! Create test icebergs of this class (-1 = no) 
    244                                                       ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    245       !rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    246       rn_test_box              = -180.0,  180.0,  70.0,  90.0     ! 
    247       rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
     206   cn_dir = './'  !  root directory for the location of the runoff files 
     207   !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     208   !           !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     209   !           !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     210   sn_rnf      = 'runoff-icb_DaiTrenberth_Depoorter.nc',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     211   sn_i_rnf    = 'runoff-icb_DaiTrenberth_Depoorter.nc',        -1         , 'Icb_flux',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     212   sn_cnf      = 'runoff-icb_DaiTrenberth_Depoorter.nc',         0         , 'socoeff' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
     213   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     214   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     215   sn_dep_rnf  = 'runoffs_eORCA1.1_depths.nc'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
     216/ 
     217!----------------------------------------------------------------------- 
     218&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     219!----------------------------------------------------------------------- 
     220   ! 
     221   ! ---------------- ice shelf load ------------------------------- 
     222   ! 
     223   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     224      rn_isfload_T = -1.9 
     225      rn_isfload_S = 34.4 
     226   ! 
     227   ! ---------------- ice shelf melt formulation ------------------------------- 
     228   ! 
     229   ln_isf = .true.           ! activate ice shelf module 
     230      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     231      cn_isfdir   = './'         ! directory for all ice shelf input file 
     232      ! 
     233      ! ---------------- cavities opened ------------------------------- 
     234      ! 
     235      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     236         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     237         !                       ! spe = fwfisf is read from a forcing field 
     238         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     239         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     240         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     241         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     242         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     243         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     244         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     245         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     246         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     247         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     248         ! 
     249         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     250         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     251         ! 
     252         !* 'spe' and 'oasis' case 
     253         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     254         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     255         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     256         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     257      ! 
     258      ! ---------------- cavities parametrised ------------------------------- 
     259      ! 
     260      ln_isfpar_mlt = .true.   ! ice shelf melting parametrised 
     261         cn_isfpar_mlt = 'oasis'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     262         !                      ! spe   = fwfisf is read from a forcing field 
     263         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     264         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     265         !* bg03 case 
     266         rn_isfpar_bg03_gt0 = 1.0e-4 ! gamma coeficient used in bg03 paper [m/s] 
     267         ! 
     268         ! 
     269         !* all cases 
     270         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     271         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     272         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     273         sn_isfpar_zmax = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -12      ,'sodepmax_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     274         sn_isfpar_zmin = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -12      ,'sodepmin_isf' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     275         !* 'spe' and 'oasis' case 
     276         sn_isfpar_fwf = 'runoff-icb_DaiTrenberth_Depoorter.nc' ,   -12      ,'sornfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     277         !* 'bg03' case 
     278         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     279      ! 
     280      ! ---------------- ice sheet coupling ------------------------------- 
     281      ! 
     282      ln_isfcpl = .false. 
     283         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     284         ln_isfcpl_cons = .false. 
     285/ 
     286!----------------------------------------------------------------------- 
     287&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
     288!----------------------------------------------------------------------- 
     289/ 
     290!----------------------------------------------------------------------- 
     291&namberg       !   iceberg parameters                                   (default: OFF) 
     292!----------------------------------------------------------------------- 
     293/ 
     294!!====================================================================== 
     295!!               ***  Lateral boundary condition  ***                 !! 
     296!!                                                                    !! 
     297!!   namlbc        lateral momentum boundary condition                  (default: NO selection) 
     298!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif") 
     299!!   nam_tide      Tidal forcing                                        (default: OFF) 
     300!!   nambdy        Unstructured open boundaries                         (default: OFF) 
     301!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy) 
     302!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF) 
     303!!====================================================================== 
     304! 
     305!----------------------------------------------------------------------- 
     306&namlbc        !   lateral momentum boundary condition                  (default: NO selection) 
     307!----------------------------------------------------------------------- 
     308   rn_shlat    =    0.     !  no slip 
     309/ 
     310!----------------------------------------------------------------------- 
     311&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     312!----------------------------------------------------------------------- 
     313/ 
     314!!====================================================================== 
     315!!                ***  Top/Bottom boundary condition  ***             !! 
     316!!                                                                    !! 
     317!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
     318!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
     319!!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     320!!   nambbc        bottom temperature boundary condition                (default: OFF) 
     321!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
     322!!====================================================================== 
     323! 
     324!----------------------------------------------------------------------- 
     325&namdrg        !   top/bottom drag coefficient                          (default: NO selection) 
     326!----------------------------------------------------------------------- 
     327   ln_non_lin  = .true.   !  non-linear  drag: Cd = Cd0 |U| 
     328   ln_loglayer = .false.   !  logarithmic drag: Cd = vkarmn/log(z/z0) |U| 
     329   ! 
     330   ln_drgimp   = .true.    !  implicit top/bottom friction flag 
     331/ 
     332!----------------------------------------------------------------------- 
     333&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     334!----------------------------------------------------------------------- 
     335   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
     336   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
     337   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     338   rn_z0       =  3.e-3    !  roughness [m] (ln_loglayer=T) 
     339   ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant 
     340      rn_boost =  50.         !  local boost factor  [-] 
     341/ 
     342!----------------------------------------------------------------------- 
     343&nambbc        !   bottom temperature boundary condition                (default: OFF) 
     344!----------------------------------------------------------------------- 
     345   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
     346      nn_geoflx   =    2      !  geothermal heat flux: = 2 read variable flux [mW/m2] 
    248347 
    249 !              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    250 !              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    251       sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , '' 
    252     
    253       cn_dir = './'  
    254 / 
    255 !----------------------------------------------------------------------- 
    256 &namlbc        !   lateral momentum boundary condition 
    257 !----------------------------------------------------------------------- 
    258    rn_shlat    =    0.0    !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    259 / 
    260 !----------------------------------------------------------------------- 
    261 &namcla        !   cross land advection 
    262 !----------------------------------------------------------------------- 
    263 / 
    264 !----------------------------------------------------------------------- 
    265 &nambfr        !   bottom friction 
    266 !----------------------------------------------------------------------- 
    267    nn_bfr      =    2      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    268 / 
    269 !----------------------------------------------------------------------- 
    270 &nambbc        !   bottom temperature boundary condition 
    271 !----------------------------------------------------------------------- 
    272    sn_qgh      ='Goutorbe_ghflux.nc',  -12.  , 'gh_flux'    ,   .false.     , .true. , 'yearly'  , 'weights_Goutorbe1_2_eorca1_bilinear.nc'       , ''       , '' 
    273    ! 
    274    cn_dir      = './'      !  root directory for the location of the runoff files 
    275    nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
    276 / 
    277 !----------------------------------------------------------------------- 
    278 &nambbl        !   bottom boundary layer scheme 
    279 !----------------------------------------------------------------------- 
    280 / 
    281 !----------------------------------------------------------------------- 
    282 &nameos        !   ocean physical parameters 
    283 !----------------------------------------------------------------------- 
    284 / 
    285 !----------------------------------------------------------------------- 
    286 &namtra_adv    !   advection scheme for tracer 
    287 !----------------------------------------------------------------------- 
    288    ln_traadv_tvd    =  .true.   !  TVD scheme 
    289    ln_traadv_ubs    =  .false.  !  UBS scheme 
    290 / 
    291 !----------------------------------------------------------------------- 
    292 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    293 !----------------------------------------------------------------------- 
    294 / 
    295 !---------------------------------------------------------------------------------- 
    296 &namtra_ldf    !   lateral diffusion scheme for tracers 
    297 !---------------------------------------------------------------------------------- 
    298    ln_traldf_grif   =  .false.   !  griffies skew flux formulation       (require "key_ldfslp") 
    299    ln_traldf_gdia   =  .false.   !  griffies operator strfn diagnostics  (require "key_ldfslp") 
    300    ln_botmix_grif   =  .false.   !  griffies operator with lateral mixing on bottom (require "key_ldfslp") 
    301    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    302    rn_aeiv_0        =  1000.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
    303 / 
    304 !----------------------------------------------------------------------- 
    305 &namtra_dmp    !   tracer: T & S newtonian damping 
    306 !----------------------------------------------------------------------- 
    307    ln_tradmp   =  .false.  !  add a damping termn (T) or not (F) 
    308 / 
    309 !----------------------------------------------------------------------- 
    310 &namdyn_adv    !   formulation of the momentum advection 
    311 !----------------------------------------------------------------------- 
    312    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    313    ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    314    ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    315    nn_dynkeg     = 1       !  scheme for grad(KE): =0  C2  ; =1  Hollingsworth correction 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 &nam_vvl    !   vertical coordinate options 
    319 !----------------------------------------------------------------------- 
    320 / 
    321 !----------------------------------------------------------------------- 
    322 &namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    323 !----------------------------------------------------------------------- 
    324 / 
    325 !----------------------------------------------------------------------- 
    326 &namdyn_hpg    !   Hydrostatic pressure gradient option 
    327 !----------------------------------------------------------------------- 
    328    ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
    329    ln_hpg_sco  = .true.    !  s-coordinate (standard jacobian formulation) 
    330    !ln_hpg_isf  = .true.    !  s-coordinate (sco ) adapted to isf 
    331    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    332                            !           centered      time scheme  (F) 
    333 / 
    334 !----------------------------------------------------------------------- 
    335 &namdyn_ldf    !   lateral diffusion on momentum 
    336 !----------------------------------------------------------------------- 
    337    rn_ahm_0_lap     = 20000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    338 / 
    339 !----------------------------------------------------------------------- 
    340 &namzdf        !   vertical physics 
    341 !----------------------------------------------------------------------- 
    342 / 
    343 !----------------------------------------------------------------------- 
    344 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    345 !----------------------------------------------------------------------- 
    346   nn_etau     = 0         !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    347                           !    = 0 no penetration 
    348                           !    = 1 add a tke source below the ML 
    349                           !    = 2 add a tke source just at the base of the ML 
    350                           !    = 3 as = 1 applied on HF part of the stress    ("key_oasis3") 
    351   nn_mxl0     = 2         ! type of scaling under sea-ice 
    352                           !    = 0 no scaling under sea-ice 
    353                           !    = 1 scaling with constant sea-ice thickness 
    354                           !    = 2  scaling with mean sea-ice thickness 
    355                           !    = 3  scaling with maximum sea-ice thickness 
    356   rn_hice    = 10.        ! max constant ice thickness value when scaling under sea-ice ( nn_mxl0=1) 
    357   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
    358   rn_lc      =  0.20      !  coef. associated to Langmuir cells 
    359 / 
    360 !----------------------------------------------------------------------- 
    361 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    362 !----------------------------------------------------------------------- 
    363 / 
    364 !----------------------------------------------------------------------- 
    365 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    366 !----------------------------------------------------------------------- 
    367 / 
    368 !----------------------------------------------------------------------- 
    369 &namzdf_tmx_new    !   new tidal mixing parameterization                ("key_zdftmx_new") 
    370 !----------------------------------------------------------------------- 
    371    nn_zpyc     = 2         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
    372    ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
     348   cn_dir = './'  !  root directory for the geothermal data location 
     349   !___________!____________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     350   !           !  file name         ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     351   !           !                    !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     352   sn_qgh      ='Lucazeau_ghflux.nc',        -12.       , 'gh_flux' ,   .false.   , .true. , 'yearly'  , 'weights_Lucazeau1_2_bilinear.nc'       , ''       , '' 
     353/ 
     354!----------------------------------------------------------------------- 
     355&nambbl        !   bottom boundary layer scheme                         (default: OFF) 
     356!----------------------------------------------------------------------- 
     357   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
     358      nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     359      nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     360      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     361      rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     362/ 
     363!!====================================================================== 
     364!!                        Tracer (T-S) namelists                      !! 
     365!!                                                                    !! 
     366!!   nameos        equation of state                                    (default: NO selection) 
     367!!   namtra_adv    advection scheme                                     (default: NO selection) 
     368!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection) 
     369!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF) 
     370!!   namtra_eiv    eddy induced velocity param.                         (default: OFF) 
     371!!   namtra_dmp    T & S newtonian damping                              (default: OFF) 
     372!!====================================================================== 
     373! 
     374!----------------------------------------------------------------------- 
     375&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
     376!----------------------------------------------------------------------- 
     377   ln_teos10   = .true.         !  = Use TEOS-10 
     378/ 
     379!----------------------------------------------------------------------- 
     380&namtra_adv    !   advection scheme for tracer                          (default: NO selection) 
     381!----------------------------------------------------------------------- 
     382   ln_traadv_fct = .true.     !  FCT scheme 
     383      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     384      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     385/ 
     386!----------------------------------------------------------------------- 
     387&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection) 
     388!----------------------------------------------------------------------- 
     389   ln_traldf_lap   = .true.    !    laplacian operator 
     390   ln_traldf_iso   = .true.    !  iso-neutral (Standard operator) 
     391   ln_traldf_msc   = .true.    !  Method of Stabilizing Correction      (both operators) 
     392   !                       !  Coefficients: 
     393   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     394      !                             !   = 20     aht = 1/2  Ud. max(e1,e2) 
     395      rn_Ud        = 0.018          !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     396      rn_Ld        = 100.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     397/ 
     398!----------------------------------------------------------------------- 
     399&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)        (default: OFF) 
     400!----------------------------------------------------------------------- 
     401   ln_mle      = .true.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     402/ 
     403!----------------------------------------------------------------------- 
     404&namtra_eiv    !   eddy induced velocity param.                         (default: OFF) 
     405!----------------------------------------------------------------------- 
     406   ln_ldfeiv   = .true.    ! use eddy induced velocity parameterization 
     407      !                        !  Coefficients: 
     408      nn_aei_ijk_t  = 21          ! space/time variation of the eiv coeficient 
     409      !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     410      !                           !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
     411      rn_Ue        = 0.018             !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     412      rn_Le        = 100.e+3           !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     413      ! 
     414      ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     415/ 
     416!----------------------------------------------------------------------- 
     417&namtra_dmp    !   tracer: T & S newtonian damping                      (default: OFF) 
     418!----------------------------------------------------------------------- 
     419   ln_tradmp   =  .false.   !  add a damping term (using resto.nc coef.) 
     420      nn_zdmp  =    0         !  vertical shape =0    damping throughout the water column 
     421/ 
     422!!====================================================================== 
     423!!                      ***  Dynamics namelists  ***                  !! 
     424!!                                                                    !! 
     425!!   nam_vvl       vertical coordinate options                          (default: z-star) 
     426!!   namdyn_adv    formulation of the momentum advection                (default: NO selection) 
     427!!   namdyn_vor    advection scheme                                     (default: NO selection) 
     428!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection) 
     429!!   namdyn_spg    surface pressure gradient                            (default: NO selection) 
     430!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection) 
     431!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only) 
     432!!====================================================================== 
     433! 
     434!----------------------------------------------------------------------- 
     435&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
     436!----------------------------------------------------------------------- 
     437   ln_dynadv_vec = .true.  !  vector form - 2nd centered scheme 
     438     nn_dynkeg     = 1        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction 
     439/ 
     440!----------------------------------------------------------------------- 
     441&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection) 
     442!----------------------------------------------------------------------- 
     443   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     444      nn_e3f_typ = 0          ! =0   e3f = mean masked e3t divided by 4 
     445/ 
     446!----------------------------------------------------------------------- 
     447&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection) 
     448!----------------------------------------------------------------------- 
     449   ln_hpg_sco  = .true.   !  s-coordinate (standard jacobian formulation) 
     450/ 
     451!----------------------------------------------------------------------- 
     452&namdyn_spg    !   surface pressure gradient                            (default: NO selection) 
     453!----------------------------------------------------------------------- 
     454   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     455/ 
     456!----------------------------------------------------------------------- 
     457&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection) 
     458!----------------------------------------------------------------------- 
     459   ln_dynldf_lap =  .true.     !    laplacian operator 
     460   ln_dynldf_lev =  .true.     !  iso-level 
     461   nn_ahm_ijk_t  = -30           !  space/time variation of eddy coefficient : 
     462      !                             !  =-30  read in eddy_viscosity_3D.nc file 
     463/ 
     464!!====================================================================== 
     465!!                     vertical physics namelists                     !! 
     466!!                                                                    !! 
     467!!    namzdf        vertical physics manager                            (default: NO selection) 
     468!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T) 
     469!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T) 
     470!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T) 
     471!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T) 
     472!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T) 
     473!!====================================================================== 
     474! 
     475!----------------------------------------------------------------------- 
     476&namzdf        !   vertical physics manager                             (default: NO selection) 
     477!----------------------------------------------------------------------- 
     478   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     479   ln_zdfevd   = .true.       !  Enhanced Vertical Diffusion scheme 
     480      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
     481      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
     482   ln_zdfddm   = .true.    ! double diffusive mixing 
     483      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     484      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     485   ln_zdfiwm   = .true.       ! internal wave-induced mixing            (T =>   fill namzdf_iwm) 
     486   !                       !  Coefficients 
     487   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     488   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     489   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     490   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     491/ 
     492!----------------------------------------------------------------------- 
     493&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     494!----------------------------------------------------------------------- 
     495   nn_mxl      =   2       !  mixing length: = 2 first vertical derivative of mixing length bounded by 1 
     496   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     497      nn_mxlice   = 2         ! type of scaling under sea-ice 
     498                              !    = 2  scaling with mean sea-ice thickness 
     499   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs 
     500                               !        = 0 none ; = 1 add a tke source below the ML 
     501   ln_lc      = .true.     !  Langmuir cell parameterisation (Axell 2002) 
     502   rn_lc      =  0.20      !  coef. associated to Langmuir cells 
     503   nn_eice     =   3       !  attenutaion of langmuir & surface wave breaking under ice 
     504   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     505/ 
     506!----------------------------------------------------------------------- 
     507&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T) 
     508!----------------------------------------------------------------------- 
     509   ln_mevar    = .false.    !  variable (T) or constant (F) mixing efficiency 
    373510   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    374 / 
    375 !----------------------------------------------------------------------- 
    376 &namsol        !   elliptic solver / island / free surface 
    377 !----------------------------------------------------------------------- 
    378 / 
    379 !----------------------------------------------------------------------- 
    380 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    381 !----------------------------------------------------------------------- 
    382    ln_nnogather=  .true.   ! 
    383    jpni        =   22      !  jpni   number of processors following i (set automatically if < 1) 
    384    jpnj        =   22      !  jpnj   number of processors following j (set automatically if < 1) 
    385    jpnij       =  360      !  jpnij  number of local domains (set automatically if < 1), 360 for eORCA1/IPSLCM6-LR 
    386 / 
    387 !----------------------------------------------------------------------- 
    388 &namctl        !   Control prints & Benchmark 
    389 !-----------------------------------------------------------------------  
    390 / 
    391 !----------------------------------------------------------------------- 
    392 &namptr       !   Poleward Transport Diagnostic 
    393 !----------------------------------------------------------------------- 
    394    ln_diaptr  = .true.      !  Poleward heat and salt transport (T) or not (F) 
    395    ln_subbas  = .true.      !  Atlantic/Pacific/Indian basins computation (T) or not 
    396                             !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    397 / 
    398 !----------------------------------------------------------------------- 
    399 &namhsb       !  Heat and salt budgets 
    400 !----------------------------------------------------------------------- 
    401    ln_diahsb  = .true.   
    402 / 
    403 !----------------------------------------------------------------------- 
    404 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    405 !----------------------------------------------------------------------- 
    406 / 
    407 !----------------------------------------------------------------------- 
    408 &nam_vvl    !   vertical coordinate options 
    409 !----------------------------------------------------------------------- 
    410 / 
    411 !----------------------------------------------------------------------- 
    412 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    413 !----------------------------------------------------------------------- 
    414 / 
    415 !----------------------------------------------------------------------- 
    416 &namtrd        !   diagnostics on dynamics and/or tracer trends 
    417 !              !       and/or mixed-layer trends and/or barotropic vorticity 
     511 
     512   cn_dir      = './'      !  root directory for the iwm data location 
     513   !___________!____________________!___________________!_____________!_____________!________!___________!__________________!__________!_______________! 
     514   !           !  file name         ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     515   !           !                    !  (if <0  months)  !   name      !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     516   sn_mpb      = 'zdfiwm_forcing.nc',      -12.         , 'power_bot' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     517   sn_mpc      = 'zdfiwm_forcing.nc',      -12.         , 'power_cri' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     518   sn_mpn      = 'zdfiwm_forcing.nc',      -12.         , 'power_nsq' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     519   sn_mps      = 'zdfiwm_forcing.nc',      -12.         , 'power_sho' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     520   sn_dsb      = 'zdfiwm_forcing.nc',      -12.         , 'scale_bot' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     521   sn_dsc      = 'zdfiwm_forcing.nc',      -12.         , 'scale_cri' ,    .false.  , .true. , 'yearly' , '' , ''  , '' 
     522/ 
     523!!====================================================================== 
     524!!                  ***  Diagnostics namelists  ***                   !! 
     525!!                                                                    !! 
     526!!   namtrd       dynamics and/or tracer trends                         (default: OFF) 
     527!!   namptr       Poleward Transport Diagnostics                        (default: OFF) 
     528!!   namhsb       Heat and salt budgets                                 (default: OFF) 
     529!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
     530!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
     531!!   namflo       float parameters                                      (default: OFF) 
     532!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
     533!!   nam_diadct   transports through some sections                      (default: OFF) 
     534!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     535!!   nam_dia25h   25h Mean Output                                       (default: OFF) 
     536!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
     537!!====================================================================== 
     538!----------------------------------------------------------------------- 
     539&namtrd        !   trend diagnostics                                    (default: OFF) 
    418540!----------------------------------------------------------------------- 
    419541   ln_tra_trd  = .true.   ! (T) 3D tracer trend output 
    420542/ 
    421 !----------------------------------------------------------------------- 
    422 &namsto       ! Stochastic parametrization of EOS 
    423 !----------------------------------------------------------------------- 
    424 / 
     543! 
     544!!====================================================================== 
     545!!               ***  Observation & Assimilation  ***                 !! 
     546!!                                                                    !! 
     547!!   namobs       observation and model comparison                      (default: OFF) 
     548!!   nam_asminc   assimilation increments                               ('key_asminc') 
     549!!====================================================================== 
     550! 
     551!!====================================================================== 
     552!!                  ***  Miscellaneous namelists  ***                 !! 
     553!!                                                                    !! 
     554!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi") 
     555!!   namctl            Control prints                                   (default: OFF) 
     556!!   namsto            Stochastic parametrization of EOS                (default: OFF) 
     557!!====================================================================== 
     558! 
     559!----------------------------------------------------------------------- 
     560&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
     561!----------------------------------------------------------------------- 
     562   nn_hls      =   1       !  halo width (applies to both rows and columns) 
     563/ 
     564!----------------------------------------------------------------------- 
     565&namctl        !   Control prints                                       (default: OFF) 
     566!----------------------------------------------------------------------- 
     567   sn_cfctl%l_runstat = .true.    ! switches and which areas produce reports with the proc integer settings. 
     568   sn_cfctl%l_trcstat = .true.   ! The default settings for the proc integers should ensure 
     569/ 
     570!----------------------------------------------------------------------- 
     571&namsto        ! Stochastic parametrization of EOS                      (default: OFF) 
     572!----------------------------------------------------------------------- 
     573/ 
     574!----------------------------------------------------------------------- 
     575&namhsb        !  Heat and salt budgets                                 (default: OFF) 
     576!----------------------------------------------------------------------- 
     577   ln_diahsb   = .true.   !  output the heat and salt budgets (T) or not (F) 
     578/ 

  • CONFIG/UNIFORM/v7/IPSLCM7/GENERAL/PARAM/namelist_top_ORCA1_cfg

    r5479 r6346  
    33!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!----------------------------------------------------------------------- 
    5 &namtrc_run     !   run information 
     5&namtrc_run      !   run information 
    66!----------------------------------------------------------------------- 
    77   ln_top_euler  = .true.    !  use Euler time-stepping for TOP 
    8    ln_rsttr      = _AUTOBLOCKER_   !  AUTO - start from a restart file (T) or not (F) 
    9    nn_rsttr      = _AUTOBLOCKER_   !  AUTO - restart control = 0 initial time step is not compared to the restart file value 
     8   ln_rsttr      = _AUTO_   !  AUTO - start from a restart file (T) or not (F) 
     9   nn_rsttr      = _AUTO_   !  AUTO - restart control = 0 initial time step is not compared to the restart file value 
    1010                             !                  = 1 do not use the value in the restart file 
    1111                             !                  = 2 calendar parameters read in the restart file 
    … …  
    1414/ 
    1515!----------------------------------------------------------------------- 
    16 &namtrc     !   tracers definition 
     16&namtrc          !   tracers definition 
    1717!----------------------------------------------------------------------- 
     18   jp_bgc        =  24 
     19! 
     20   ln_pisces     =  .true. 
     21   ln_my_trc     =  .false. 
     22   ln_age        =  _AUTO_ 
     23   ln_cfc11      =  _AUTO_  
     24   ln_cfc12      =  _AUTO_ 
     25   ln_c14        =  .false. 
     26 
     27! 
    1828   ln_trcdta     =  .true.  !  Initialisation from data input file (T) or not (F) 
    19    ln_trcdmp_clo =  .true.  !  restoring on closed seas (T) or not (F)  
    20  
    21  
    22 !                !    name   !           title of the field              ! initial data ! initial data ! save   ! 
    23 !                !           !                                           !  units       ! from file    ! or not !  
    24 !                !           !                                           !              ! or not       !        ! 
    25    sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true.     ,  .true. 
    26    sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true. 
    27    sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true.     ,  .true. 
    28    sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    29    sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true.     ,  .true. 
    30    sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    31    sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    32    sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false.    ,  .true. 
    33    sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    34    sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .true.     ,  .true. 
    35    sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    36    sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    37    sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    38    sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true.     ,  .true. 
    39    sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    40    sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    41    sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    42    sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false.    ,  .true. 
    43    sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false.    ,  .true. 
    44    sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false.    ,  .true. 
    45    sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    46    sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    47    sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    48    sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false.    ,  .true. 
     29   ln_trcbc      =  .true.  !  Enables Boundary conditions 
     30   ln_trcdmp_clo =  .true.  !  damping term (T) or not (F) on closed seas 
     31   ln_trcais     =  .false.  !  Antarctic Ice Sheet nutrient supply 
     32!                !           !                                          !             ! 
     33!                !    name   !           title of the field              !   units     ! init    ! sbc    ! cbc    !  obc    !  ais 
     34   sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' , .true.  , .false., .true. , .false. , .false. 
     35   sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   , .true.  , .false., .true. , .false. , .false. 
     36   sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' , .true.  , .false., .false., .false. , .false. 
     37   sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     38   sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' , .true.  , .true. , .true. , .false. , .false. 
     39   sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     40   sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' , .true.  , .true. , .true. , .false. , .false. 
     41   sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     42   sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     43   sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' , .true.  , .false., .true. , .false. , .false. 
     44   sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     45   sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     46   sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     47   sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' , .true.  , .true. , .true. , .false. , .false. 
     48   sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     49   sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     50   sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     51   sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     52   sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     53   sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     54   sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     55   sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     56   sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' , .true.  , .true. , .true. , .false. , .false. 
     57   sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' , .false. , .false., .false., .false. , .false. 
     58/ 
     59!----------------------------------------------------------------------- 
     60&namage          !   AGE  
     61!----------------------------------------------------------------------- 
    4962/ 
    5063!----------------------------------------------------------------------- 
    5164&namtrc_dta      !    Initialisation from data input file 
    5265!----------------------------------------------------------------------- 
    53 !                !  file name                          ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    54 !                !                                     !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    55    sn_trcdta(1)  = 'DIC_GLODAP_annual_eORCA_R1.nc'     ,        -12        ,  'DIC'     ,    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
    56    sn_trcdta(2)  = 'Alkalini_GLODAP_annual_eORCA_R1.nc',        -12        ,  'Alkalini',    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
    57    sn_trcdta(3)  = 'O2_WOA2009_monthly_eORCA_R1.nc'    ,        -1         ,  'O2'      ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    58    sn_trcdta(5)  = 'PO4_WOA2009_monthly_eORCA_R1.nc'   ,        -1         ,  'PO4'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    59    sn_trcdta(7)  = 'Si_WOA2009_monthly_eORCA_R1.nc'    ,        -1         ,  'Si'      ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    60    sn_trcdta(10) = 'DOC_PISCES_monthly_eORCA_R1.nc'    ,        -1         ,  'DOC'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    61    sn_trcdta(14) = 'Fer_PISCES_monthly_eORCA_R1.nc'    ,        -1         ,  'Fer'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    62    sn_trcdta(23) = 'NO3_WOA2009_monthly_eORCA_R1.nc'   ,        -1         ,  'NO3'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    63    rn_trfac(1)   =   1.028e-06  !  multiplicative factor 
    64    rn_trfac(2)   =   1.028e-06  !  -      -      -     - 
    65    rn_trfac(3)   =  44.6e-06  !  -      -      -     - 
    66    rn_trfac(5)   = 122.0e-06  !  -      -      -     - 
    67    rn_trfac(7)   =   1.0e-06  !  -      -      -     - 
    68    rn_trfac(10)  =   1.0e-06  !  -      -      -     - 
    69    rn_trfac(14)  =   1.0e-06  !  -      -      -     - 
    70    rn_trfac(23)  =   7.6e-06  !  -      -      -     - 
     66!                !  file name          ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     67!                !                     !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     68   sn_trcdta(1)  = 'data_DIC_nomask.nc',        -12        ,  'PiDIC'   ,    .false.   , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     69   sn_trcdta(2)  = 'data_ALK_nomask.nc',        -12        ,  'TALK'    ,    .false.   , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     70   sn_trcdta(3)  = 'data_OXY_nomask.nc',        -12        ,  'O2'      ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     71   sn_trcdta(5)  = 'data_PO4_nomask.nc',        -12        ,  'PO4'     ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     72   sn_trcdta(7)  = 'data_SIL_nomask.nc',        -12        ,  'Si'      ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     73   sn_trcdta(10) = 'data_DOC_nomask.nc',        -1         ,  'DOC'     ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     74   sn_trcdta(14) = 'data_FER_nomask.nc',        -12        ,  'Fer'     ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     75   sn_trcdta(23) = 'data_NO3_nomask.nc',        -12        ,  'NO3'     ,    .true.    , .true. , 'yearly'  , 'weights_3D_r360x180_bilin.nc'       , ''   , '' 
     76   rn_trfac(1)   =   1.028e-06   !  multiplicative factor 
     77   rn_trfac(2)   =   1.028e-06   !  -      -      -     - 
     78   rn_trfac(3)   =  44.6e-06     !  -      -      -     - 
     79   rn_trfac(5)   = 117.0e-06     !  -      -      -     - 
     80   rn_trfac(7)   =   1.0e-06     !  -      -      -     - 
     81   rn_trfac(10)  =   1.0e-06     !  -      -      -     - 
     82   rn_trfac(14)  =   1.0e-06     !  -      -      -     - 
     83   rn_trfac(23)  =   7.3125e-06  !  -      -      -     - 
    7184/ 
    7285!----------------------------------------------------------------------- 
    73 &namtrc_adv    !   advection scheme for passive tracer  
     86&namtrc_adv      !   advection scheme for passive tracer                (default: NO selection) 
    7487!----------------------------------------------------------------------- 
    75    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    76    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     88   ln_trcadv_mus =  .true.  !  MUSCL scheme 
     89      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    7790/ 
    7891!----------------------------------------------------------------------- 
    79 &namtrc_ldf    !   lateral diffusion scheme for passive tracer  
     92&namtrc_ldf      !   lateral diffusion scheme for passive tracer        (default: NO selection) 
    8093!----------------------------------------------------------------------- 
    81    rn_fact_lap      =     15.    !     enhanced zonal eddy diffusivity 
     94   ln_trcldf_tra   = .true.      !  use active tracer setting 
     95   rn_fact_lap     =     15.    !     enhanced zonal eddy diffusivity 
    8296/ 
    8397!----------------------------------------------------------------------- 
    84 &namtrc_zdf        !   vertical physics 
     98&namtrc_rad      !  treatment of negative concentrations  
    8599!----------------------------------------------------------------------- 
    86100/ 
    87101!----------------------------------------------------------------------- 
    88 &namtrc_rad        !  treatment of negative concentrations  
     102&namtrc_snk      !  sedimentation of particles 
    89103!----------------------------------------------------------------------- 
    90104/ 
    91105!----------------------------------------------------------------------- 
    92 &namtrc_dmp    !   passive tracer newtonian damping 
     106&namtrc_dcy      !  Diurnal cycle 
     107!----------------------------------------------------------------------- 
     108   ln_trcdc2dm   =  .true.   !  Diurnal cycle for TOP 
     109/ 
     110!----------------------------------------------------------------------- 
     111&namtrc_dmp      !   passive tracer newtonian damping    
    93112!----------------------------------------------------------------------- 
    94113/ 
    95114!----------------------------------------------------------------------- 
    96 &namtrc     !   tracers definition 
     115&namtrc_ice      !    Representation of sea ice growth & melt effects 
    97116!----------------------------------------------------------------------- 
    98117/ 
    99118!----------------------------------------------------------------------- 
    100 &namtrc_ice       !    Representation of sea ice growth & melt effects 
     119&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc') 
     120!---------------------------------------------------------------------- 
     121/ 
     122!---------------------------------------------------------------------- 
     123&namtrc_bc       !   data for boundary conditions 
     124!----------------------------------------------------------------------- 
     125!                !  file name  ! frequency (hours) ! variable      ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     126!                !             !  (if <0  months)  !   name        !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     127   sn_trcsbc(5)  = 'dustdep'   ,       -1          , 'dustpo4'     ,  .true.      , .true. , 'yearly'  , 'weights_2D_r360x180_bilin.nc'       , ''    , '' 
     128   sn_trcsbc(7)  = 'dustdep'   ,       -1          , 'dustsi'      ,  .true.      , .true. , 'yearly'  , 'weights_2D_r360x180_bilin.nc'       , ''    , '' 
     129   sn_trcsbc(14) = 'dustdep'   ,       -1          , 'dustfer'     ,  .true.      , .true. , 'yearly'  , 'weights_2D_r360x180_bilin.nc'       , ''    , '' 
     130   sn_trcsbc(23) = 'nitdep'    ,       -12         , 'ndep2'       ,  .true.      , .true. , 'yearly'  , 'weights_2D_r360x180_bilin.nc'       , ''    , '' 
     131   rn_trsfac(5)  = 7.9258065e-02    !  (  0.021 / 31. * 117 ) 
     132   rn_trsfac(7)  = 3.1316726e-01    !  ( 8.8   / 28.1 ) 
     133   rn_trsfac(14) = 6.2667860e-04    !  (  0.035 / 55.85 ) 
     134   rn_trsfac(23) = 5.2232143e-01    !  ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.3125/14 ) 
     135   rn_sbc_time   =  1.          !  Time scaling factor for SBC and CBC data (seconds in a day) 
     136   ! 
     137   sn_trccbc(1)  = 'river.orca' ,    -1            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     138   sn_trccbc(2)  = 'river.orca' ,    -1            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     139   sn_trccbc(5)  = 'river.orca' ,    -1            , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     140   sn_trccbc(7)  = 'river.orca' ,    -1            , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     141   sn_trccbc(10) = 'river.orca' ,    -1            , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     142   sn_trccbc(14) = 'river.orca' ,    -1            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     143   sn_trccbc(23) = 'river.orca' ,    -1            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     144   rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
     145   rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
     146   rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
     147   rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
     148   rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 
     149   rn_trcfac(14) = 4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
     150   rn_trcfac(23) = 5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
     151   rn_cbc_time   = 3.1536e+7      !  Time scaling factor for CBC data (seconds in a year) 
     152/ 
     153!---------------------------------------------------------------------- 
     154&namtrc_bdy      !   Setup of tracer boundary conditions 
    101155!----------------------------------------------------------------------- 
    102156/ 
    103157!----------------------------------------------------------------------- 
    104 &namtrc_trd       !   diagnostics on tracer trends        ('key_trdtrc') 
    105 !                          or mixed-layer trends          ('key_trdmld_trc') 
    106 !---------------------------------------------------------------------- 
     158&namtrc_ais      !  Representation of Antarctic Ice Sheet tracers supply 
     159!----------------------------------------------------------------------- 
     160   rn_trafac(14) =  4.476e-07   !  (  0.5e-3 / 55.85 * 0.05 ) 
     161! 
     162   nn_ais_tr     =  0      !  tracer concentration in iceberg and ice shelf 
     163                           !    = 0 (null concentrations) 
     164                           !    = 1 prescribed concentrations 
     165   rn_icbdep     =  120.   ! Mean underwater depth of iceberg (m) 
    107166/ 
    108 !----------------------------------------------------------------------- 
    109 &namtrc_dia       !   parameters for passive tracer additional diagnostics 
    110 !---------------------------------------------------------------------- 
    111 / 
    112 !---------------------------------------------------------------------- 
    113 &namtrc_bc        !   data for boundary conditions 
    114 !----------------------------------------------------------------------- 
    115 / 
Note: See TracChangeset for help on using the changeset viewer.