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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11844 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs – NEMO

Ignore:
Timestamp:
2019-11-01T21:54:15+01:00 (4 years ago)
Author:
mathiot
Message:

rm useless lbclnk + fix reproducibility WED025 + add isf debug option

Location:
NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED/namelist_ref

    r11835 r11844  
    455455         !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    456456         !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    457          rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    458          rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
     457         rn_gammat0  = 1.4e-2    ! gammat coefficient used in blk formula 
     458         rn_gammas0  = 4.e-4    ! gammas coefficient used in blk formula 
    459459         ! 
    460460         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    … …  
    465465         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    466466         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    467          sn_isfcav_fwf = 'isfmlt_cav',      -12       , 'fwflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     467         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    468468      ! 
    469469      ! ---------------- cavities parametrised ------------------------------- 
    … …  
    482482         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    483483         !* 'spe' and 'oasis' case 
    484          sn_isfpar_fwf = 'isfmlt_par' ,      -12       ,'fwfisf'   ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     484         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    485485         !* 'bg03' case 
    486          sn_isfpar_Leff = 'isfmlt_par',       0        ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     486         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    487487      ! 
    488488      ! ---------------- ice sheet coupling ------------------------------- 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/WED025/EXPREF/file_def_nemo-ice.xml

    r11487 r11844  
    1717       <!-- ice mask --> 
    1818       <field field_ref="icemask"          name="simsk"   /> 
    19        <field field_ref="icemask05"        name="simsk05" /> 
    2019       <field field_ref="icemask15"        name="simsk15" /> 
    2120        
    … …  
    2827       <field field_ref="iceapnd"          name="siapnd" /> 
    2928       <field field_ref="icevpnd"          name="sivpnd" /> 
    30        <field field_ref="sst_m"            name="sst_m"  /> 
    31        <field field_ref="sss_m"            name="sss_m"  /> 
    3229        
    3330       <!-- heat --> 
    … …  
    7976       <field field_ref="vfxsnw"           name="vfxsnw" /> 
    8077 
    81        <!-- categories --> 
    82        <field field_ref="icemask_cat"      name="simskcat"/> 
    83        <field field_ref="snwthic_cat"      name="snthicat"/> 
    84        <field field_ref="iceconc_cat"      name="siconcat"/> 
    85        <field field_ref="icethic_cat"      name="sithicat"/> 
    86        <field field_ref="icesalt_cat"      name="sisalcat"/> 
    87        <field field_ref="icetemp_cat"      name="sitemcat"/> 
    88        <field field_ref="snwtemp_cat"      name="sntemcat"/> 
    89        <field field_ref="icettop_cat"      name="sitopcat"/> 
    90  
    9178     </file> 
    9279      

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/WED025/EXPREF/file_def_nemo-oce.xml

    r11495 r11844  
    99    --> 
    1010    
    11     <file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1mo" min_digits="4"> 
     11    <file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="5d" min_digits="4"> 
    1212     
    1313      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE.">  <!-- 5d files -->    
    … …  
    2020          <field field_ref="sss"          name="sos"   /> 
    2121          <field field_ref="ssh"          name="zos"   /> 
    22           <field field_ref="sst"          name="tosstd"       long_name="sea surface temperature standard deviation"         operation="instant" freq_op="5d" > sqrt( @sst2 - @sst * @sst ) </field> 
    23           <field field_ref="ssh"          name="zosstd"       long_name="sea surface height above geoid standard deviation"  operation="instant" freq_op="5d" > sqrt( @ssh2 - @ssh * @ssh ) </field> 
    24           <field field_ref="sst"          name="sstdcy"       long_name="amplitude of sst diurnal cycle"                     operation="average" freq_op="1d" > @sstmax - @sstmin </field> 
    2522 
    26           <field field_ref="mldkz5"   /> 
    2723          <field field_ref="mldr10_1" /> 
    28           <field field_ref="mldr10_1"     name="mldr10_1dcy"  long_name="amplitude of mldr10_1 diurnal cycle"                operation="average" freq_op="1d" > @mldr10_1max - @mldr10_1min </field> 
    2924 
    3025          <field field_ref="sbt"                          /> 
    … …  
    4439          <field field_ref="isfgammat"    name="sogammat_cav"  /> 
    4540          <field field_ref="isfgammas"    name="sogammas_cav"  /> 
     41          <field field_ref="isfustar"     name="ustar_cav"  /> 
    4642 
    4743          <!-- ice shelf tbl properties  --> 
    … …  
    5349 
    5450          <!-- ice shelf fluxes  --> 
    55           <field field_ref="fwfisf_cav"    name="sowflisf_cav" /> 
    56           <field field_ref="qoceisf_cav"   name="qoceisf_cav"  /> 
    57           <field field_ref="qlatisf_cav"   name="qlatisf_cav"  /> 
    58           <field field_ref="qhcisf_cav"    name="qhcisf_cav"   /> 
     51          <field field_ref="fwfisf_cav"      name="sowflisf_cav" /> 
     52          <field field_ref="qoceisf_cav"     name="qoceisf_cav"  /> 
     53          <field field_ref="qlatisf_cav"     name="qlatisf_cav"  /> 
     54          <field field_ref="qhcisf_cav"      name="qhcisf_cav"   /> 
     55          <field field_ref="qconisf"         name="qhcond_cav"   /> 
     56          <field field_ref="fwfisf3d_cav"    name="sowflisf3d_cav" /> 
    5957           
    60           <field field_ref="fwfisf_par"    name="sowflisf_par" /> 
    61           <field field_ref="qoceisf_par"   name="qoceisf_par"  /> 
    62           <field field_ref="qlatisf_par"   name="qlatisf_par"  /> 
    63           <field field_ref="qhcisf_par"    name="qhcisf_par"   /> 
     58          <field field_ref="fwfisf_par"      name="sowflisf_par" /> 
     59          <field field_ref="qoceisf_par"     name="qoceisf_par"  /> 
     60          <field field_ref="qlatisf_par"     name="qlatisf_par"  /> 
     61          <field field_ref="qhcisf_par"      name="qhcisf_par"   /> 
     62          <field field_ref="fwfisf3d_par"    name="sowflisf3d_par" /> 
    6463 
    6564   </file> 
    … …  
    8685   </file> 
    8786 
    88         <file id="file15" name_suffix="_scalar" description="scalar variables" > 
    89           <!-- global drifts (conservation checks) --> 
    90           <field field_ref="bgtemper"     grid_ref="grid_1point"   name="bgtemper"    /> 
    91           <field field_ref="bgsaline"     grid_ref="grid_1point"   name="bgsaline"    /> 
    92           <field field_ref="bgheatco"     grid_ref="grid_1point"   name="bgheatco"    /> 
    93           <field field_ref="bgheatfx"     grid_ref="grid_1point"   name="bgheatfx"    /> 
    94           <field field_ref="bgsaltco"     grid_ref="grid_1point"   name="bgsaltco"    /> 
    95           <field field_ref="bgvolssh"     grid_ref="grid_1point"   name="bgvolssh"    /> 
    96           <field field_ref="bgvole3t"     grid_ref="grid_1point"   name="bgvole3t"    /> 
    97  
    98           <!-- global surface forcings  --> 
    99           <field field_ref="bgfrcvol"     grid_ref="grid_1point"   name="bgfrcvol"    /> 
    100           <field field_ref="bgfrctem"     grid_ref="grid_1point"   name="bgfrctem"    /> 
    101           <field field_ref="bgfrchfx"     grid_ref="grid_1point"   name="bgfrchfx"    /> 
    102           <field field_ref="bgfrcsal"     grid_ref="grid_1point"   name="bgfrcsal"    /> 
    103         </file> 
    104  
    10587      </file_group>  
    10688 

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/WED025/EXPREF/namelist_cfg

    r11494 r11844  
    211211!----------------------------------------------------------------------- 
    212212   ! 
    213    ! ---------------- general ------------------------------- 
     213   ! ---------------- ice shelf load ------------------------------- 
     214   ! 
     215   cn_isfload = 'isomip'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     216   ! 
     217   ! ---------------- ice shelf melt formulation ------------------------------- 
    214218   ! 
    215219   ln_isf = .true.            ! activate ice shelf module 
    216    cn_isfdir = './'           ! directory for all ice shelf input file 
    217    ! 
    218    ! ---------------- cavities opened ------------------------------- 
    219    ! 
    220    cn_isfload = 'isomip' 
    221    ! 
    222    ln_isfcav_mlt = .true.     ! ice shelf melting into the cavity 
    223       cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
    224       !                       ! spe = fwfisf is read from a forcing field 
    225       !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    226       !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    227       !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
    228       !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
    229       cn_gammablk = 'ad15'    ! scheme to compute gammat/s (spe,ad15,hj99) 
    230       !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    231       !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    232       rn_gammat0  = 1.4e-2    ! gammat coefficient used in blk formula 
    233       rn_gammas0  = 4.e-4    ! gammas coefficient used in blk formula 
     220      cn_isfdir = './'           ! directory for all ice shelf input file 
    234221      ! 
    235       rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    236       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     222      ! ---------------- cavities opened ------------------------------- 
    237223      ! 
    238       !* 'spe' and 'oasis' case 
    239       !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    240       !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    241       !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    242       sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    243    ! 
    244    ! ---------------- cavities parametrised ------------------------------- 
    245    ! 
    246    ln_isfpar_mlt = .true.    ! ice shelf melting parametrised 
    247       cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
    248       !                      ! spe   = fwfisf is read from a forcing field 
    249       !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
    250       !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     224      ln_isfcav_mlt = .true.     ! ice shelf melting into the cavity 
     225         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     226         !                       ! spe = fwfisf is read from a forcing field 
     227         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     228         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     229         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     230         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     231         cn_gammablk = 'ad15'    ! scheme to compute gammat/s (spe,ad15,hj99) 
     232         !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     233         !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     234         rn_gammat0  = 1.4e-2    ! gammat coefficient used in blk formula 
     235         rn_gammas0  = 4.e-4    ! gammas coefficient used in blk formula 
     236         ! 
     237         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     238         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     239         ! 
     240         !* 'spe' and 'oasis' case 
     241         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     242         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     243         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     244         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    251245      ! 
    252       !* all cases 
    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_isfpar_zmax = 'isfmlt_par',      -12.      , 'sozisfmax' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    257       sn_isfpar_zmin = 'isfmlt_par',      -12.      , 'sozisfmin' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    258       !* 'spe' and 'oasis' case 
    259       sn_isfpar_fwf = 'isfmlt_par' ,      -12.      , 'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    260       !* 'bg03' case 
    261       sn_isfpar_Leff = 'isfmlt_par',       0.       , 'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    262    ! 
    263    ! ---------------- ice sheet coupling ------------------------------- 
    264    ! 
    265    ln_isfcpl = .false. 
    266       nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    267       ln_isfcpl_cons = .false. 
     246      ! ---------------- cavities parametrised ------------------------------- 
     247      ! 
     248      ln_isfpar_mlt = .true.    ! ice shelf melting parametrised 
     249         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     250         !                      ! spe   = fwfisf is read from a forcing field 
     251         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     252         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     253         ! 
     254         !* all cases 
     255         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     256         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     257         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     258         sn_isfpar_zmax = 'isfmlt_par',      -12.      , 'sozisfmax' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     259         sn_isfpar_zmin = 'isfmlt_par',      -12.      , 'sozisfmin' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     260         !* 'spe' and 'oasis' case 
     261         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      , 'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     262         !* 'bg03' case 
     263         sn_isfpar_Leff = 'isfmlt_par',       0.       , 'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     264      ! 
     265      ! ---------------- ice sheet coupling ------------------------------- 
     266      ! 
     267      ln_isfcpl = .false. 
     268         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     269         ln_isfcpl_cons = .false. 
    268270/ 
    269271!----------------------------------------------------------------------- 
    … …  
    595597&namhsb        !  Heat and salt budgets                                 (default: OFF) 
    596598!----------------------------------------------------------------------- 
    597    ln_diahsb  = .true.    !  check the heat and salt budgets (T) or not (F) 
    598599/ 
    599600!----------------------------------------------------------------------- 
    … …  
    664665&namctl        !   Control prints                                       (default: OFF) 
    665666!----------------------------------------------------------------------- 
     667   ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
     668     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
     669       sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
     670       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
     671       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
     672       sn_cfctl%l_layout  = .FALSE. ! 
     673       sn_cfctl%l_mppout  = .FALSE. ! 
     674       sn_cfctl%l_mpptop  = .FALSE. ! 
     675       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
     676       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
     677       sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
     678       sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
     679   nn_print    =    0      !  level of print (0 no extra print) 
     680   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
     681   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
     682   nn_jctls    =    0      !  start j indice of control               over a subdomain) 
     683   nn_jctle    =    0      !  end   j indice of control 
     684   nn_isplt    =    1      !  number of processors in i-direction 
     685   nn_jsplt    =    1      !  number of processors in j-direction 
     686   ln_timing   = .true.    !  timing by routine write out in timing.output file 
     687   ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
    666688/ 
    667689!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.