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 10911 for NEMO/trunk/cfgs – NEMO

Changeset 10911 for NEMO/trunk/cfgs


Ignore:
Timestamp:
2019-04-29T18:11:12+02:00 (5 years ago)
Author:
clem
Message:

Major change: the advection scheme UMx has been revisited to clean all the unphysical values which occured. Minor change: the ref config SPITZ12 has been slightly modified to test more parameterizations that are available in the code (melt ponds and landfast)

Location:
NEMO/trunk/cfgs
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/cfgs/ORCA2_ICE_PISCES/EXPREF/file_def_nemo-ice.xml

    r9943 r10911  
    7979       <field field_ref="vfxsnw"           name="vfxsnw" /> 
    8080        
    81        <!-- diag error for negative ice volume after advection --> 
    82        <field field_ref="iceneg_pres"      name="sineg_pres" /> 
    83        <field field_ref="iceneg_volu"      name="sineg_volu" /> 
    84        <field field_ref="iceneg_hfx"       name="sineg_hfx"  /> 
    85  
    8681       <!-- categories --> 
    8782       <field field_ref="icemask_cat"      name="simskcat"/> 
  • NEMO/trunk/cfgs/SHARED/field_def_nemo-ice.xml

    r10578 r10911  
    160160          <field id="hfxdhc"       long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" /> 
    161161 
    162      <!-- diagnostics of the negative values resulting from the advection scheme --> 
    163      <field id="iceneg_pres"  long_name="Fraction of time steps with negative sea ice volume"                    unit=""     /> 
    164      <field id="iceneg_volu"  long_name="Negative sea ice volume per area arising from advection"                unit="m"    /> 
    165           <field id="iceneg_hfx"   long_name="Negative sea ice heat content (eq. heat flux) arising from advection"   unit="W/m2" /> 
    166  
    167162     <!-- sbcssm variables --> 
    168163          <field id="sst_m"    unit="degC" /> 
  • NEMO/trunk/cfgs/SHARED/namelist_ice_ref

    r10611 r10911  
    154154&namthd_do      !   Ice growth in open water 
    155155!------------------------------------------------------------------------------ 
    156    rn_hinew         =   0.1           !  thickness for new ice formation in open water (m), must be larger than rn_hnewice 
     156   rn_hinew         =   0.1           !  thickness for new ice formation in open water (m), must be larger than rn_himin 
    157157   ln_frazil        = .false.         !  Frazil ice parameterization (ice collection as a function of wind) 
    158158      rn_maxfraz    =   1.0           !     maximum fraction of frazil ice collecting at the ice base 
  • NEMO/trunk/cfgs/SPITZ12/EXPREF/file_def_nemo-ice.xml

    r9943 r10911  
    7979       <field field_ref="vfxsnw"           name="vfxsnw" /> 
    8080 
    81        <!-- diag error for negative ice volume after advection --> 
    82        <field field_ref="iceneg_pres"      name="sineg_pres" /> 
    83        <field field_ref="iceneg_volu"      name="sineg_volu" /> 
    84        <field field_ref="iceneg_hfx"       name="sineg_hfx"  /> 
    85         
    8681       <!-- categories --> 
    8782       <field field_ref="icemask_cat"      name="simskcat"/> 
  • NEMO/trunk/cfgs/SPITZ12/EXPREF/namelist_cfg

    r10075 r10911  
    8080                     ! Sea-ice : 
    8181   nn_ice      = 2         !  SI3 
     82   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
     83      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    8284                     ! Misc. options of sbc : 
    8385   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr ) 
     
    9092      ! 
    9193      ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    92       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
     94      ln_Cd_L15   = .true.    !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    9395      ! 
    9496   cn_dir = './'  !  root directory for the bulk data location 
     
    9698   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !          weights filename            ! rotation ! land/sea mask ! 
    9799   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    98 !! ERAI 
    99 !   sn_wndi     = 'u10_era_spitz'               ,         6         , 'u10'     ,   .true.    , .false. , 'yearly'  , 'weights_bicub', 'Uwnd' , '' 
    100 !   sn_wndj     = 'v10_era_spitz'               ,         6         , 'v10'     ,   .true.    , .false. , 'yearly'  , 'weights_bicub', 'Vwnd' , '' 
    101 !   sn_qsr      = 'ssrd_era_spitz'              ,         6         , 'ssrd'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    102 !   sn_qlw      = 'strd_era_spitz'              ,         6         , 'strd'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    103 !   sn_tair     = 't2_era_spitz'                ,         6         , 't2'      ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    104 !   sn_humi     = 'humi_era_spitz'              ,         6         , 'humi'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    105 !   sn_prec     = 'precip_era_spitz'            ,         6         , 'precip'  ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    106 !   sn_snow     = 'snow_era_spitz'              ,         6         , 'snow'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin', '' , '' 
    107 !   sn_tdif     = 'taudif'                      ,         6         , 'taudif'  ,   .true.    , .false. , 'yearly'  , ''             , '' , '' 
    108 !      rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    109 !! MAR 
    110100   sn_wndi     = 'MARv3.6-9km-Svalbard-2hourly_spitz' ,   2 ,  'u10'     ,   .true.    , .false. , 'yearly'  , 'weights_bicub', 'Uwnd' , '' 
    111101   sn_wndj     = 'MARv3.6-9km-Svalbard-2hourly_spitz' ,   2 ,  'v10'     ,   .true.    , .false. , 'yearly'  , 'weights_bicub', 'Vwnd' , '' 
     
    129119   !                       ! type of penetration                        (default: NO selection) 
    130120   ln_qsr_rgb  = .true.       !  RGB light penetration (Red-Green-Blue) 
    131 / 
    132 !----------------------------------------------------------------------- 
    133 &namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
    134 !----------------------------------------------------------------------- 
    135 / 
    136 !----------------------------------------------------------------------- 
    137 &namsbc_rnf    !   runoffs                                              (ln_rnf =T) 
    138 !----------------------------------------------------------------------- 
    139121/ 
    140122!!====================================================================== 
     
    236218!----------------------------------------------------------------------- 
    237219   ln_loglayer = .true.   !  logarithmic drag: Cd = vkarmn/log(z/z0) |U| 
     220   ln_drgimp   = .true.   !  implicit top/bottom friction flag 
    238221/ 
    239222!----------------------------------------------------------------------- 
    240223&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
    241224!----------------------------------------------------------------------- 
    242    rn_Cd0      =  2.5e-3   !!1.e-3    !  drag coefficient [-] 
    243    rn_Cdmax    =  0.01     !!0.1      !  drag value maximum [-] (logarithmic drag) 
    244    rn_ke0      =  0.       !!2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     225   rn_Cd0      =  2.5e-3   !  drag coefficient [-] 
     226   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
     227   rn_ke0      =  0.       !  background kinetic energy  [m2/s2] (non-linear cases) 
    245228   rn_z0       =  3.e-3    !  roughness [m] (ln_loglayer=T) 
    246    ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant 
    247       rn_boost =  50.         !  local boost factor  [-] 
    248 / 
    249 !----------------------------------------------------------------------- 
    250 &nambbc        !   bottom temperature boundary condition                (default: OFF) 
    251 !----------------------------------------------------------------------- 
    252    ln_trabbc   = .false.    !  Apply a geothermal heating at the ocean bottom 
    253229/ 
    254230!----------------------------------------------------------------------- 
     
    258234      nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    259235      nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
    260       rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    261       rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
    262236/ 
    263237!!====================================================================== 
     
    293267   !                       !  Coefficients: 
    294268   nn_aht_ijk_t    = 31        !  space/time variation of eddy coefficient: 
    295       !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     269   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    296270/ 
    297271!!====================================================================== 
     
    320294&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection) 
    321295!----------------------------------------------------------------------- 
    322    ln_dynvor_enT = .true. !  energy conserving scheme (T-point) 
     296   ln_dynvor_eeT = .true.  !  energy conserving scheme (een using e3t) 
    323297/ 
    324298!----------------------------------------------------------------------- 
     
    352326&namzdf        !   vertical physics manager                             (default: NO selection) 
    353327!----------------------------------------------------------------------- 
     328   !                       ! adaptive-implicit vertical advection 
     329   ln_zad_Aimp = .true.      !  Courant number dependent scheme (Shchepetkin 2015) 
    354330   !                       ! type of vertical closure (required) 
    355331   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
    356332   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     333   !                       ! convection 
     334   ln_zdfevd   = .true.       !  enhanced vertical diffusion 
     335   ! 
    357336   ln_zdfddm   = .true.    ! double diffusive mixing 
     337   ! 
    358338   !                       !  Coefficients 
    359339   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     
    384364/ 
    385365!----------------------------------------------------------------------- 
    386 &nam_diaharm   !   Harmonic analysis of tidal constituents              ('key_diaharm') 
    387 !----------------------------------------------------------------------- 
    388 / 
    389 !----------------------------------------------------------------------- 
    390366&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
    391367!----------------------------------------------------------------------- 
  • NEMO/trunk/cfgs/SPITZ12/EXPREF/namelist_ice_cfg

    r10535 r10911  
    3535&namdyn         !   Ice dynamics 
    3636!------------------------------------------------------------------------------ 
     37   ln_landfast_L16  = .true.          !  landfast: parameterization from Lemieux 2016 
    3738/ 
    3839!------------------------------------------------------------------------------ 
     
    4344&namdyn_rhg     !   Ice rheology 
    4445!------------------------------------------------------------------------------ 
     46   ln_rhg_EVP       = .true.          !  EVP rheology 
     47      ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    4548/ 
    4649!------------------------------------------------------------------------------ 
     
    6770&namthd_do      !   Ice growth in open water 
    6871!------------------------------------------------------------------------------ 
    69    rn_hinew         =   0.02          !  thickness for new ice formation in open water (m), must be larger than rn_hnewice 
     72   rn_hinew         =   0.02          !  thickness for new ice formation in open water (m), must be larger than rn_himin 
    7073   ln_frazil        = .true.          !  Frazil ice parameterization (ice collection as a function of wind) 
    7174/ 
     
    7780&namthd_pnd     !   Melt ponds 
    7881!------------------------------------------------------------------------------ 
    79    ln_pnd_H12       = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
    80    ln_pnd_CST       = .false.         !  activate constant melt ponds 
    81    ln_pnd_alb       = .false.         !  melt ponds affect albedo or not 
     82   ln_pnd_H12       = .true.          !  activate evolutive melt ponds (from Holland et al 2012) 
     83   ln_pnd_alb       = .true.          !  melt ponds affect albedo or not 
    8284/ 
    8385 
Note: See TracChangeset for help on using the changeset viewer.