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 13540 for NEMO/branches/2020/r12377_ticket2386/cfgs – NEMO

Ignore:
Timestamp:
2020-09-29T12:41:06+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2386: update to latest trunk

Location:
NEMO/branches/2020/r12377_ticket2386
Files:
43 edited
1 copied

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/r12377_ticket2386

    • Property svn:externals

      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13507        sette 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/1_context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg

    r12511 r13540  
    8181   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    8282                     ! Sea-ice : 
    83    nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    84                            !          except in AGRIF zoom where it has to be specified 
     83   nn_ice      = 2         !  =0 no ice boundary condition 
     84      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     85      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8586                     ! Misc. options of sbc :  
    8687   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    … …  
    9495   !                    !  bulk algorithm : 
    9596   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    96    ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    97    ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    98    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    99       ! 
    100       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    101       rn_zu       = 10.       !  Wind vector reference height (m) 
    102       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    103       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    104       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    105       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    106       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    107       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    108       ln_skin_cs = .false.  !  use the cool-skin parameterization (only available in ECMWF and COARE algorithms) !LB 
    109       ln_skin_wl = .false.  !  use the warm-layer        "               "                    " 
    110       ! 
    111       ln_humi_sph = .true.     !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    112       ln_humi_dpt = .false.    !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    113       ln_humi_rlh = .false.    !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    11497   ! 
    11598   cn_dir = './'  !  root directory for the bulk data location 
    … …  
    192175!!                                                                    !! 
    193176!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    194 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    195 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     177!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     178!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    196179!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    197180!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    370353&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    371354!----------------------------------------------------------------------- 
    372       rn_eice     =   0       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
    373355/ 
    374356!!====================================================================== 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/2_context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg

    r12511 r13540  
    7878   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    7979                     ! Sea-ice : 
    80    nn_ice      = 2         !  =0   Use SI3 model 
     80   nn_ice      = 2         !  =0 no ice boundary condition 
     81      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     82      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8183                     ! Misc. options of sbc :  
    8284   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    … …  
    9092   !                    !  bulk algorithm : 
    9193   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    92    ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    93    ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    94    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    95       ! 
    96       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    97       rn_zu       = 10.       !  Wind vector reference height (m) 
    98       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    99       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    100       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    101       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    102       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    103       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    104       ln_skin_cs = .false.  !  use the cool-skin parameterization (only available in ECMWF and COARE algorithms) !LB 
    105       ln_skin_wl = .false.  !  use the warm-layer        "               "                    " 
    106       ! 
    107       ln_humi_sph = .true.     !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    108       ln_humi_dpt = .false.    !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    109       ln_humi_rlh = .false.    !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    11094   ! 
    11195   cn_dir      = './'      !  root directory for the bulk data location 
    … …  
    174158!----------------------------------------------------------------------- 
    175159   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    176    rn_sponge_tra = 1440.   !  coefficient for tracer   sponge layer [m2/s] 
    177    rn_sponge_dyn = 1440.   !  coefficient for dynamics sponge layer [m2/s] 
    178160   ln_chk_bathy  = .true.  !  =T  check the parent bathymetry 
    179161/ 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/3_context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg

    r12511 r13540  
    7878   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    7979                     ! Sea-ice : 
    80    nn_ice      = 2         !  =0   Use SI3 model 
     80   nn_ice      = 2         !  =0 no ice boundary condition 
     81      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     82      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8183                     ! Misc. options of sbc :  
    8284   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    … …  
    156158!----------------------------------------------------------------------- 
    157159   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    158    rn_sponge_tra =  480.   !  coefficient for tracer   sponge layer [m2/s] 
    159    rn_sponge_dyn =  480.   !  coefficient for dynamics sponge layer [m2/s] 
    160160   ln_chk_bathy  = .true.  !  =T  check the parent bathymetry 
    161161/ 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/AGRIF_FixedGrids.in

    r9770 r13540  
    112 
    2 42 82 49 91 1 1 1 
    3 122 153 110 143 4 4 4 
     241 81 49 91 1 1 1 
     3121 152 110 143 4 4 4 
    440 
    551 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg

    r12511 r13540  
    3333/ 
    3434!----------------------------------------------------------------------- 
    35 &namcfg        !   parameters of the configuration                      (default: user defined GYRE) 
     35&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
    3636!----------------------------------------------------------------------- 
    3737   ln_read_cfg = .true.    !  (=T) read the domain configuration file 
    … …  
    4242/ 
    4343!----------------------------------------------------------------------- 
    44 &namtsd        !    Temperature & Salinity Data                         (default: OFF) 
     44&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF) 
    4545!----------------------------------------------------------------------- 
    4646   !                       ! =T  read T-S fields for: 
    … …  
    6363!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
    6464!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only) 
     65!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   ) 
    6566!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
     67!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    6668!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
     69!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    6770!!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    6871!!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    69 !!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    70 !!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    7172!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    7273!!   namberg         iceberg floats                                     (ln_icebergs=T) 
    … …  
    7475! 
    7576!----------------------------------------------------------------------- 
    76 &namsbc        !   Surface Boundary Condition (surface module) 
     77&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
    7778!----------------------------------------------------------------------- 
    7879   nn_fsbc     = 1         !  frequency of SBC module call 
    … …  
    8182   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    8283                     ! Sea-ice : 
    83    nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    84                            !          except in AGRIF zoom where it has to be specified 
     84   nn_ice      = 2         !  =0 no ice boundary condition 
     85      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     86      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8587                     ! Misc. options of sbc :  
    8688   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
     89   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    8790   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
    88    ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    8991   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
    9092/ 
    9193!----------------------------------------------------------------------- 
    92 &namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk =T) 
     94&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
    9395!----------------------------------------------------------------------- 
    9496   !                    !  bulk algorithm : 
    95    ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    96    ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    97    ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    98    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    99       ! 
    100       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    101       rn_zu       = 10.       !  Wind vector reference height (m) 
    102       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    103       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    104       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    105       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    106       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    107       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    108       ln_skin_cs = .false.  !  use the cool-skin parameterization (only available in ECMWF and COARE algorithms) !LB 
    109       ln_skin_wl = .false.  !  use the warm-layer        "               "                    " 
    110       ! 
    111       ln_humi_sph = .true.     !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    112       ln_humi_dpt = .false.    !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    113       ln_humi_rlh = .false.    !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    114    ! 
     97   ln_NCAR    = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
    11598   cn_dir = './'  !  root directory for the bulk data location 
    11699   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    117100   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !          weights filename            ! rotation ! land/sea mask ! 
    118101   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    119    sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
    120    sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
    121    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    122    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    123    sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    124    sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    125    sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    126    sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    127    sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     102   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Uwnd'   , '' 
     103   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Vwnd'   , '' 
     104   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     105   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     106   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     107   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     108   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     109   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     110   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
    128111/ 
    129112!----------------------------------------------------------------------- 
    … …  
    142125/ 
    143126!----------------------------------------------------------------------- 
     127&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
     128!----------------------------------------------------------------------- 
     129   nn_sssr     =     2     !  add a damping term to the surface freshwater flux 
     130      rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
     131      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
     132      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
     133/ 
     134!----------------------------------------------------------------------- 
    144135&namsbc_rnf    !   runoffs                                              (ln_rnf =T) 
    145136!----------------------------------------------------------------------- 
    … …  
    147138      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
    148139      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
     140   rn_rfact    =   1.e0    !  multiplicative factor for runoff 
    149141 
    150142   cn_dir = './'  !  root directory for the location of the runoff files 
    … …  
    159151/ 
    160152!----------------------------------------------------------------------- 
    161 &namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
    162 !----------------------------------------------------------------------- 
    163    nn_sssr     =     2     !  add a damping term to the surface freshwater flux 
    164       rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    165       ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    166       rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    167 / 
    168 !----------------------------------------------------------------------- 
    169 &namberg       !   iceberg parameters                                   (default: No iceberg) 
     153&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
     154!----------------------------------------------------------------------- 
     155/ 
     156!----------------------------------------------------------------------- 
     157&namberg       !   iceberg parameters                                   (default: OFF) 
    170158!----------------------------------------------------------------------- 
    171159   !  iceberg floats are not currently available with AGRIF 
    … …  
    176164!!                                                                    !! 
    177165!!   namlbc        lateral momentum boundary condition                  (default: NO selection) 
    178 !!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     166!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif") 
    179167!!   nam_tide      Tidal forcing                                        (default: OFF) 
    180168!!   nambdy        Unstructured open boundaries                         (default: OFF) 
    … …  
    192180!!                                                                    !! 
    193181!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    194 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    195 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     182!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     183!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    196184!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    197185!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    229217/ 
    230218!!====================================================================== 
    231 !!                        Tracer (T & S) namelists                    !! 
     219!!                        Tracer (T-S) namelists                      !! 
    232220!!                                                                    !! 
    233221!!   nameos        equation of state                                    (default: NO selection) 
    … …  
    250238      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
    251239      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    252 / 
    253 !----------------------------------------------------------------------- 
    254 &namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)        (default: OFF) 
    255 !----------------------------------------------------------------------- 
    256    ln_mle      = .true.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    257240/ 
    258241!----------------------------------------------------------------------- 
    … …  
    272255/ 
    273256!----------------------------------------------------------------------- 
    274 &namtra_eiv !   eddy induced velocity param.                            (default: OFF) 
    275 !----------------------------------------------------------------------- 
    276    ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     257&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)        (default: OFF) 
     258!----------------------------------------------------------------------- 
     259   ln_mle      = .true.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     260/ 
     261!----------------------------------------------------------------------- 
     262&namtra_eiv    !   eddy induced velocity param.                         (default: OFF) 
     263!----------------------------------------------------------------------- 
     264   ln_ldfeiv   = .true.    ! use eddy induced velocity parameterization 
    277265      !                        !  Coefficients: 
    278266      nn_aei_ijk_t  = 21          ! space/time variation of the eiv coeficient 
    … …  
    303291! 
    304292!----------------------------------------------------------------------- 
    305 &nam_vvl       !   vertical coordinate options                          (default: z-star) 
    306 !----------------------------------------------------------------------- 
    307 / 
    308 !----------------------------------------------------------------------- 
    309293&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
    310294!----------------------------------------------------------------------- 
    … …  
    351335! 
    352336!----------------------------------------------------------------------- 
    353 &namzdf        !   vertical physics                                     (default: NO selection) 
     337&namzdf        !   vertical physics manager                             (default: NO selection) 
    354338!----------------------------------------------------------------------- 
    355339   !                       ! type of vertical closure 
    … …  
    370354&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    371355!----------------------------------------------------------------------- 
    372       rn_eice     =   0       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
    373356/ 
    374357!!====================================================================== 
    … …  
    401384!!====================================================================== 
    402385! 
    403 ! 
    404386!----------------------------------------------------------------------- 
    405387&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AMM12/EXPREF/context_nemo.xml

    r12377 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/AMM12/EXPREF/namelist_cfg

    r12511 r13540  
    212212!!                                                                    !! 
    213213!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    214 !!   namdrg_top    top    friction                                      (ln_OFF =F & ln_isfcav=T) 
    215 !!   namdrg_bot    bottom friction                                      (ln_OFF =F) 
     214!!   namdrg_top    top    friction                                      (ln_drg_OFF =F & ln_isfcav=T) 
     215!!   namdrg_bot    bottom friction                                      (ln_drg_OFF =F) 
    216216!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    217217!!   nambbl        bottom boundary layer scheme                         (default: OFF) 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/C1D_PAPA/EXPREF/file_def_nemo-oce.xml

    r9799 r13540  
    5353        <file id="file4" name_suffix="_grid_W" description="ocean W grid variables" > 
    5454          <field field_ref="e3w" /> 
    55           <field field_ref="woce"         name="wo"       /> 
    5655          <field field_ref="avt"          name="difvho"   /> 
    5756        </file> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/C1D_PAPA/EXPREF/namelist_cfg

    r12511 r13540  
    4949&namdom        !   time and space domain 
    5050!----------------------------------------------------------------------- 
     51   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
     52   ! 
    5153   rn_Dt      =  360.     !  time step for the dynamics and tracer 
    5254/ 
    … …  
    256258!!                                                                    !! 
    257259!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    258 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    259 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     260!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     261!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    260262!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    261263!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    268270/ 
    269271!----------------------------------------------------------------------- 
    270 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
    271 !----------------------------------------------------------------------- 
    272 / 
    273 !----------------------------------------------------------------------- 
    274 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     272&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
     273!----------------------------------------------------------------------- 
     274/ 
     275!----------------------------------------------------------------------- 
     276&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    275277!----------------------------------------------------------------------- 
    276278/ 
    … …  
    358360&namdyn_spg    !   surface pressure gradient                            (default: NO selection) 
    359361!----------------------------------------------------------------------- 
    360    ln_dynspg_ts   = .true.   ! split-explicit free surface 
    361       ln_bt_fw      = .false.     ! Forward integration of barotropic Eqs. 
    362       ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    363362/ 
    364363!----------------------------------------------------------------------- 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90

    r12377 r13540  
    3030   PUBLIC   usr_def_zgr        ! called by domzgr.F90 
    3131 
     32   !! * Substitutions 
     33#  include "do_loop_substitute.h90" 
    3234   !!---------------------------------------------------------------------- 
    3335   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    … …  
    157159         pe3vw(:,:,jk) = pe3w_1d (jk) 
    158160      END DO 
    159       DO jj = 1, jpj                      ! bottom scale factors and depth at T- and W-points 
    160          DO ji = 1, jpi 
    161             ik = k_bot(ji,jj) 
    162             pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 
    163             pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) 
    164             pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  )  
    165             ! 
    166             pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp 
    167             pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp 
    168             pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) 
    169          END DO 
    170       END DO          
     161      ! bottom scale factors and depth at T- and W-points 
     162      DO_2D( 1, 1, 1, 1 ) 
     163         ik = k_bot(ji,jj) 
     164         pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 
     165         pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) 
     166         pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  )  
     167         ! 
     168         pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp 
     169         pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp 
     170         pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) 
     171      END_2D         
    171172      !                                   ! bottom scale factors and depth at  U-, V-, UW and VW-points 
    172173      !                                   ! usually Computed as the minimum of neighbooring scale factors 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/GYRE_BFM/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/GYRE_BFM/EXPREF/namelist_cfg

    r12511 r13540  
    101101!!                                                                    !! 
    102102!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    103 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    104 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     103!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     104!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    105105!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    106106!!   nambbl        bottom boundary layer scheme                         (default: OFF) 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/GYRE_PISCES/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/GYRE_PISCES/EXPREF/namelist_cfg

    r12511 r13540  
    9999!!                                                                    !! 
    100100!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    101 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    102 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     101!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     102!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    103103!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    104104!!   nambbl        bottom boundary layer scheme                         (default: OFF) 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_ABL/EXPREF/file_def_nemo-oce.xml

    r12063 r13540  
    5656     <field field_ref="t_abl" /> 
    5757     <field field_ref="q_abl" /> 
     58     <field field_ref="uvz1_abl" /> 
     59     <field field_ref="tz1_abl" /> 
     60     <field field_ref="qz1_abl" /> 
     61     <field field_ref="uvz1_dta" /> 
     62     <field field_ref="tz1_dta" /> 
     63     <field field_ref="qz1_dta" /> 
    5864     <field field_ref="pblh" /> 
    5965   </file> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg

    r12511 r13540  
    8484   ln_abl      = .true.    !  ABL  formulation                          (T => fill namsbc_abl ) 
    8585                     ! Sea-ice : 
    86    nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    87                            !          except in AGRIF zoom where it has to be specified 
     86   nn_ice      = 2         !  =0 no ice boundary condition 
     87      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     88      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8889                     ! Misc. options of sbc :  
    8990   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    … …  
    109110   !                    !  bulk algorithm : 
    110111   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    111    ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    112    ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    113    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    114       rn_zqt      = 10.     !  Air temperature & humidity reference height (m) 
    115       rn_zu       = 10.     !  Wind vector reference height (m) 
    116       ! 
    117       ! Skin is ONLY available in ECMWF and COARE algorithms: 
    118       ln_skin_cs = .false.  !  use the cool-skin parameterization => set nn_fsbc=1 and ln_dm2dc=.true.! 
    119       ln_skin_wl = .false.  !  use the warm-layer        "        => set nn_fsbc=1 and ln_dm2dc=.true.! 
    120       ! 
    121       ln_humi_sph = .true.  !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    122       ln_humi_dpt = .false. !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    123       ln_humi_rlh = .false. !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    124112   ! 
    125113   cn_dir = './'  !  root directory for the bulk data location 
    … …  
    131119   sn_tair     = 'tair_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'tair'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bilinear' , ''    , '' 
    132120   sn_humi     = 'humi_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'humi'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bilinear' , ''    , '' 
    133    sn_hpgi     = 'uhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'uhpg'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic'  , 'UG'  , '' 
    134    sn_hpgj     = 'vhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'vhpg'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic'  , 'VG'  , '' 
    135  
    136121   sn_qsr      = 'ncar_rad.15JUNE2009_fill'                    , 24., 'SWDN_MOD', .false. , .true.  ,  'yearly' , 'weights_core_orca2_bilinear_noc.nc'      , ''    , '' 
    137122   sn_qlw      = 'ncar_rad.15JUNE2009_fill'                    , 24., 'LWDN_MOD', .false. , .true.  ,  'yearly' , 'weights_core_orca2_bilinear_noc.nc'      , ''    , '' 
    … …  
    139124   sn_snow     = 'ncar_precip.15JUNE2009_fill'                 , -1., 'SNOW'    , .false. , .true.  ,  'yearly' , 'weights_core_orca2_bilinear_noc.nc'      , ''    , '' 
    140125   sn_slp      = 'slp.15JUNE2009_fill'                         ,  6., 'SLP'     , .false. , .true.  ,  'yearly' , 'weights_core_orca2_bilinear_noc.nc'      , ''    , '' 
     126   sn_hpgi     = 'uhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'uhpg'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic'  , 'UG'  , '' 
     127   sn_hpgj     = 'vhpg_drwnlnd_ERAI_L25Z10_GLOBAL_F128R_ana1d',  24., 'vhpg'    , .false. , .false. , 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic'  , 'VG'  , '' 
    141128/ 
    142129 
    … …  
    230217!!                                                                    !! 
    231218!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    232 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    233 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     219!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     220!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    234221!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    235222!!   nambbl        bottom boundary layer scheme                         (default: OFF) 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_PISCES/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg

    r12511 r13540  
    2828&namdom        !   time and space domain 
    2929!----------------------------------------------------------------------- 
    30    rn_Dt      = 5400.     !  time step for the dynamics and tracer 
     30   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    3131/ 
    3232!----------------------------------------------------------------------- 
    … …  
    8080   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    8181                     ! Sea-ice : 
    82    nn_ice      = 2         !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    83                            !          except in AGRIF zoom where it has to be specified 
     82   nn_ice      = 2         !  =0 no ice boundary condition 
     83      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     84      !                    !  =2 or 3 for SI3 and CICE, respectively 
    8485                     ! Misc. options of sbc :  
    8586   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    … …  
    109110   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !          weights filename            ! rotation ! land/sea mask ! 
    110111   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    111    sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
    112    sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
    113    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    114    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    115    sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    116    sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    117    sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    118    sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    119    sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     112   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Uwnd'   , '' 
     113   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Vwnd'   , '' 
     114   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     115   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     116   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     117   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     118   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     119   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     120   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     121/ 
     122!----------------------------------------------------------------------- 
     123&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T) 
     124!----------------------------------------------------------------------- 
    120125/ 
    121126!----------------------------------------------------------------------- 
    … …  
    202207!!                                                                    !! 
    203208!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    204 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    205 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     209!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     210!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    206211!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    207212!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    373378                               !        = 2 add a tke source just at the base of the ML 
    374379                               !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
    375       rn_eice     =   0       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
    376380/ 
    377381!----------------------------------------------------------------------- 
    … …  
    381385   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    382386   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
     387 
     388   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     389   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     390   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     391   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     392   sn_mpb      = 'int_wave_mix'      , -12.       , 'mixing_power_bot'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     393   sn_mpp      = 'int_wave_mix'      , -12.       , 'mixing_power_pyc'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     394   sn_mpc      = 'int_wave_mix'      , -12.       , 'mixing_power_cri'   , .false.  , .true. , 'yearly' , '' , ''  , '' 
     395   sn_dsb      = 'int_wave_mix'      , -12.       , 'decay_scale_bot'    , .false.  , .true. , 'yearly' , '' , ''  , '' 
     396   sn_dsc      = 'int_wave_mix'      , -12.       , 'decay_scale_cri'    , .false.  , .true. , 'yearly' , '' , ''  , '' 
    383397/ 
    384398!!====================================================================== 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_top_cfg

    r12377 r13540  
    2020! 
    2121   ln_trcdta     =  .true.  !  Initialisation from data input file (T) or not (F) 
    22    ln_trcbc      =  .true.  !  Enables Boundary conditions 
     22   ln_trcbc      =  .false. !  Enables Boundary conditions 
    2323!                !           !                                           !             !         ! 
    2424!                !    name   !           title of the field              !   units     ! init    ! sbc    ! cbc    !  obc  !  

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_OFF_PISCES/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_cfg

    r12511 r13540  
    190190!!                                                                    !! 
    191191!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    192 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    193 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     192!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     193!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    194194!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    195195!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    201201/ 
    202202!----------------------------------------------------------------------- 
    203 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
    204 !----------------------------------------------------------------------- 
    205 / 
    206 !----------------------------------------------------------------------- 
    207 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     203&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
     204!----------------------------------------------------------------------- 
     205/ 
     206!----------------------------------------------------------------------- 
     207&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    208208!----------------------------------------------------------------------- 
    209209/ 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg

    r12377 r13540  
    2020! 
    2121   ln_trcdta     =  .true.   !  Initialisation from data input file (T) or not (F) 
    22    ln_trcbc      =  .true.   !  Enables Boundary conditions 
     22   ln_trcbc      =  .false.  !  Enables Boundary conditions 
    2323!                !           !                                           !             !         ! 
    2424!                !    name   !           title of the field              !   units     ! init    ! sbc    ! cbc    !  obc  !  

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_OFF_TRC/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg

    r12511 r13540  
    188188!!                                                                    !! 
    189189!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    190 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    191 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     190!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     191!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    192192!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    193193!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    199199/ 
    200200!----------------------------------------------------------------------- 
    201 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
    202 !----------------------------------------------------------------------- 
    203 / 
    204 !----------------------------------------------------------------------- 
    205 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     201&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
     202!----------------------------------------------------------------------- 
     203/ 
     204!----------------------------------------------------------------------- 
     205&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    206206!----------------------------------------------------------------------- 
    207207/ 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_SAS_ICE/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg

    r12377 r13540  
    5959   nn_fsbc     = 1         !  frequency of SBC module call 
    6060   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    61    nn_ice      = 2         !  =2  sea-ice model                         ("key_SI3" or "key_cice") 
     61   nn_ice      = 2         !  =0 no ice boundary condition 
     62      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     63      !                    !  =2 or 3 for SI3 and CICE, respectively 
    6264/ 
    6365!----------------------------------------------------------------------- 
    … …  
    6668   !                    !  bulk algorithm : 
    6769   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    68    ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    69    ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
    70    ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    71       ! 
    72       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    73       rn_zu       = 10.       !  Wind vector reference height (m) 
    74       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    75       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    76       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    77       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    78       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    79       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    80       ln_skin_cs = .false.  !  use the cool-skin parameterization (only available in ECMWF and COARE algorithms) !LB 
    81       ln_skin_wl = .false.  !  use the warm-layer        "               "                    " 
    82       ! 
    83       ln_humi_sph = .true.     !  humidity specified below in "sn_humi" is specific humidity     [kg/kg] if .true. 
    84       ln_humi_dpt = .false.    !  humidity specified below in "sn_humi" is dew-point temperature   [K]   if .true. 
    85       ln_humi_rlh = .false.    !  humidity specified below in "sn_humi" is relative humidity       [%]   if .true. 
    8670   ! 
    8771   cn_dir      = './'      !  root directory for the bulk data location 
    … …  
    8973   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
    9074   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    91    sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
    92    sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
    93    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    94    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    95    sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    96    sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    97    sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    98    sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    99    sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    100 / 
    101 !----------------------------------------------------------------------- 
    102 &namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    103 !----------------------------------------------------------------------- 
     75   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Uwnd'   , '' 
     76   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bicub'  , 'Vwnd'   , '' 
     77   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     78   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     79   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     80   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     81   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     82   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
     83   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core2_orca2_bilin' , ''       , '' 
    10484/ 
    10585!----------------------------------------------------------------------- 
    … …  
    142122!!                                                                    !! 
    143123!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    144 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    145 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     124!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     125!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    146126!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    147127!!   nambbl        bottom boundary layer scheme                         (default: OFF) 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SHARED/field_def_nemo-ice.xml

    r12377 r13540  
    4949          <field id="icehpnd"      long_name="melt pond depth"                                         standard_name="sea_ice_meltpond_depth"                    unit="m" />  
    5050          <field id="icevpnd"      long_name="melt pond volume"                                        standard_name="sea_ice_meltpond_volume"                   unit="m" />  
     51          <field id="icehlid"      long_name="melt pond lid depth"                                     standard_name="sea_ice_meltpondlid_depth"                 unit="m" />  
     52          <field id="icevlid"      long_name="melt pond lid volume"                                    standard_name="sea_ice_meltpondlid_volume"                unit="m" />  
    5153      
    5254     <!-- heat --> 
    … …  
    8183          <field id="icediv"       long_name="Divergence of the sea-ice velocity field"                standard_name="divergence_of_sea_ice_velocity"            unit="s-1"  /> 
    8284          <field id="iceshe"       long_name="Maximum shear of sea-ice velocity field"                 standard_name="maximum_shear_of_sea_ice_velocity"         unit="s-1"  /> 
    83       
     85          <field id="beta_evp"     long_name="Relaxation parameter of ice rheology (beta)"             standard_name="relaxation_parameter_of_ice_rheology"      unit=""  />    
     86  
    8487     <!-- surface heat fluxes --> 
    8588          <field id="qt_ice"       long_name="total heat flux at ice surface"                          standard_name="surface_downward_heat_flux_in_air"         unit="W/m2" /> 
    … …  
    173176          <field id="frq_m"    unit="-"    /> 
    174177 
     178          <!-- rheology convergence tests --> 
     179          <field id="uice_cvg"   long_name="sea ice velocity convergence"      standard_name="sea_ice_velocity_convergence"      unit="m/s" /> 
     180 
    175181     <!-- ================= --> 
    176182          <!-- Add-ons for SIMIP --> 
    … …  
    211217          <field id="dmisum"       long_name="sea-ice mass change through surface melting"             standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"                       unit="kg/m2/s" /> 
    212218          <field id="dmibom"       long_name="sea-ice mass change through bottom melting"              standard_name="tendency_of_sea_ice_amount_due_to_basal_melting"                         unit="kg/m2/s" /> 
     219          <field id="dmilam"       long_name="sea-ice mass change through lateral melting"             standard_name="tendency_of_sea_ice_amount_due_to_lateral_melting"                       unit="kg/m2/s" /> 
    213220          <field id="dmsspr"       long_name="snow mass change through snow fall"                      standard_name="snowfall_flux"                                                           unit="kg/m2/s" /> 
    214221          <field id="dmsmel"       long_name="snow mass change through melt"                           standard_name="surface_snow_melt_flux"                                                  unit="kg/m2/s" /> 
    … …  
    289296          <field id="iceapnd_cat"  long_name="Ice melt pond concentration per category"          unit=""        />  
    290297          <field id="icehpnd_cat"  long_name="Ice melt pond thickness per category"              unit="m"       detect_missing_value="true" />  
     298          <field id="icehlid_cat"  long_name="Ice melt pond lid thickness per category"          unit="m"       detect_missing_value="true" />  
    291299          <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category"               unit=""        />  
     300          <field id="iceaepnd_cat" long_name="Ice melt pond effective fraction per category"     unit=""        />  
    292301          <field id="icemask_cat"  long_name="Fraction of time step with sea ice (per category)" unit=""        /> 
    293302          <field id="iceage_cat"   long_name="Ice age per category"                              unit="days"    detect_missing_value="true" /> 
    … …  
    300309          <field id="snwthic_cat_cmip"     long_name="Snow thickness in thickness categories"          standard_name="snow_thickness_over_categories"        detect_missing_value="true" unit="m"  > snwthic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
    301310          <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%"  > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field> 
     311 
     312          <!-- heat diffusion convergence tests --> 
     313          <field id="tice_cvgerr" long_name="sea ice temperature convergence error"      standard_name="sea_ice_temperature_convergence_err" unit="K" /> 
     314          <field id="tice_cvgstp" long_name="sea ice temperature convergence iterations" standard_name="sea_ice_temperature_convergence_stp" unit=""  /> 
    302315 
    303316   </field_group> <!-- SBC_3D --> 
    … …  
    560573          <field field_ref="dmisum"           name="sidmassmelttop"   /> 
    561574          <field field_ref="dmibom"           name="sidmassmeltbot"   /> 
     575          <field field_ref="dmilam"           name="sidmassmeltlat"   /> 
    562576          <field field_ref="dmsspr"           name="sndmasssnf"       /> 
    563577          <field field_ref="dmsmel"           name="sndmassmelt"      /> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SHARED/field_def_nemo-oce.xml

    r12377 r13540  
    1 <?xml version="1.0"?>  
     1<?xml version="1.0"?> 
    22    <!-- $id$ --> 
    33 
    … …  
    1616                         Configuration of multiple-linear-regression analysis (diamlr) 
    1717       ===================================================================================================== 
    18         
     18 
    1919       This field group configures diamlr for tidal harmonic analysis of field 
    2020       ssh: in addition to a regressor for fitting the mean value (diamlr_r101), 
    … …  
    7373 
    7474    </field_group> 
    75      
    76     <!--  
     75 
     76    <!-- 
    7777============================================================================================================ 
    7878=                                  definition of all existing variables                                    = 
    … …  
    101101    </field_group> 
    102102 
    103     <!--  
     103    <!-- 
    104104============================================================================================================ 
    105105                                  Physical ocean model variables 
    … …  
    108108 
    109109      <!-- T grid --> 
    110        
     110 
    111111      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    112112        <field id="e3t"          long_name="T-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_3D" /> 
    113113        <field id="e3ts"         long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_SFC"/> 
    114114        <field id="e3t_0"        long_name="Initial T-cell thickness"            standard_name="ref_cell_thickness"    unit="m"   grid_ref="grid_T_3D" /> 
    115         <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/>  
     115        <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/> 
    116116        <field id="e3t_300"      field_ref="e3t"                grid_ref="grid_T_zoom_300"       detect_missing_value="true" /> 
    117117        <field id="e3t_vsum300"  field_ref="e3t_300"            grid_ref="grid_T_vsum"   detect_missing_value="true" /> 
    118118   <field id="masscello"    long_name="Sea Water Mass per unit area"   standard_name="sea_water_mass_per_unit_area"   unit="kg/m2"   grid_ref="grid_T_3D"/> 
    119         <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D"/>  
     119        <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D"/> 
    120120        <field id="toce"         long_name="temperature"                         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    121121        <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
    … …  
    129129        <!-- AGRIF sponge --> 
    130130        <field id="agrif_spt"         long_name=" AGRIF t-sponge coefficient"   unit=" " /> 
     131    
     132   <!-- additions to diawri.F90 --> 
     133        <field id="socegrad"    long_name="module of salinity gradient"              unit="psu/m"   grid_ref="grid_T_3D"/> 
     134        <field id="socegrad2"   long_name="square of module of salinity gradient"    unit="psu2/m2" grid_ref="grid_T_3D"/> 
     135        <field id="ke"          long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" /> 
     136        <field id="ke_int"      long_name="vertical integration of kinetic energy"   unit="m3/s2"   /> 
     137        <field id="relvor"      long_name="relative vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
     138        <field id="absvor"      long_name="absolute vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
     139        <field id="potvor"      long_name="potential vorticity"                      unit="s-1"    grid_ref="grid_T_3D"/> 
     140        <field id="salt2c"      long_name="Salt content vertically integrated"       unit="1e-3*kg/m2" /> 
    131141 
    132142        <!-- t-eddy viscosity coefficients (ldfdyn) --> 
    … …  
    146156        <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                                 unit="degC"     /> 
    147157   <field id="temp_3m"      long_name="temperature at 3m"                                                                                      unit="degC"     /> 
    148          
     158 
    149159        <field id="sss"          long_name="sea surface salinity"                               standard_name="sea_surface_salinity"                unit="1e-3"     /> 
    150160        <field id="sss2"         long_name="square of sea surface salinity"                                                                         unit="1e-6"      > sss * sss </field > 
    … …  
    152162        <field id="sssmin"       long_name="min of sea surface salinity"      field_ref="sss"   operation="minimum"                                                 /> 
    153163        <field id="sbs"          long_name="sea bottom salinity"                                                                                    unit="0.001"    /> 
    154         <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" />  
    155  
    156         <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     />  
     164        <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" /> 
     165 
     166        <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     /> 
    157167 
    158168         <!-- Case EOS = TEOS-10 : output potential temperature --> 
    … …  
    177187        <field id="alpha"        long_name="thermal expansion"                                                         unit="degC-1" grid_ref="grid_T_3D" /> 
    178188        <field id="beta"         long_name="haline contraction"                                                        unit="1e3"    grid_ref="grid_T_3D" /> 
    179         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                                           unit="s-1"    grid_ref="grid_T_3D" /> 
    180189        <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" /> 
    181190 
    182191        <!-- Energy - horizontal divergence --> 
    183         <field id="eken"         long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" /> 
    184192        <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" /> 
    185193 
    … …  
    295303        <field id="us_y"        long_name="j component of Stokes drift"                      unit="m/s"     /> 
    296304      </field_group> 
    297        
     305 
    298306      <!-- SBC --> 
    299307      <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 
    … …  
    311319          <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   /> 
    312320          <field id="wclosea"      long_name="closed sea empmr correction"          standard_name="closea_empmr"                                                         unit="kg/m2/s"   /> 
    313       
     321 
    314322          <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           /> 
    315323          <field id="qns"          long_name="non solar Downward Heat Flux"                                                                                               unit="W/m2"                           /> 
    … …  
    321329          <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           /> 
    322330          <field id="wspd"         long_name="wind speed module"                     standard_name="wind_speed"                                                           unit="m/s"                            /> 
    323            
     331 
    324332          <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    325333          <field id="ssh_ib"       long_name="Inverse barometer sea surface height"  standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency"   unit="m"        /> 
    … …  
    369377          <field id="taum_oce"     long_name="wind stress module over open ocean"           standard_name="magnitude_of_surface_downward_stress"               unit="N/m2"  /> 
    370378 
     379          <!-- variables computed by the bulk parameterization algorithms (ln_blk) --> 
     380          <field id="Cd_oce"      long_name="Drag coefficient over open ocean"              standard_name="drag_coefficient_water"                unit=""  /> 
     381          <field id="Ce_oce"      long_name="Evaporaion coefficient over open ocean"        standard_name="evap_coefficient_water"                unit=""  /> 
     382          <field id="Ch_oce"      long_name="Sensible heat coefficient over open ocean"     standard_name="sensible_heat_coefficient_water"       unit=""  /> 
     383          <field id="theta_zt"    long_name="Potential air temperature at z=zt"             standard_name="potential_air_temperature_at_zt"       unit="degC" /> 
     384          <field id="q_zt"        long_name="Specific air humidity at z=zt"                 standard_name="specific_air_humidity_at_zt"           unit="kg/kg" /> 
     385          <field id="theta_zu"    long_name="Potential air temperature at z=zu"             standard_name="potential_air_temperature_at_zu"       unit="degC" /> 
     386          <field id="q_zu"        long_name="Specific air humidity at z=zu"                 standard_name="specific_air_humidity_at_zu"           unit="kg/kg" /> 
     387          <field id="ssq"         long_name="Saturation specific humidity of air at z=0"    standard_name="surface_air_saturation_spec_humidity"  unit="kg/kg" /> 
     388          <field id="wspd_blk"    long_name="Bulk wind speed at z=zu"                       standard_name="bulk_wind_speed_at_zu"                 unit="m/s"   /> 
     389          <!-- ln_blk + key_si3 --> 
     390          <field id="Cd_ice"      long_name="Drag coefficient over ice"                     standard_name="drag_coefficient_ice"                 unit=""  /> 
     391          <field id="Ce_ice"      long_name="Evaporaion coefficient over ice"               standard_name="evap_coefficient_ice"                 unit=""  /> 
     392          <field id="Ch_ice"      long_name="Sensible heat coefficient over ice"            standard_name="sensible_heat_coefficient_ice"        unit=""  /> 
     393 
    371394          <!-- available key_oasis3 --> 
    372395          <field id="snow_ao_cea"  long_name="Snow over ice-free ocean (cell average)"   standard_name="snowfall_flux"                             unit="kg/m2/s"  /> 
    … …  
    405428          <!-- ice field (nn_ice=1)  --> 
    406429          <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            /> 
    407            
     430 
    408431          <!-- dilution --> 
    409432          <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                                                                                              unit="kg*degC/m2/s" /> 
    410           <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" />         
     433          <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" /> 
    411434          <field id="rnf_x_sst"    long_name="Runoff term on SST"                                                                                                              unit="kg*degC/m2/s" /> 
    412435          <field id="rnf_x_sss"    long_name="Runoff term on SSS"                                                                                                              unit="kg*1e-3/m2/s" /> 
    413       
     436 
    414437     <!-- sbcssm variables --> 
    415438          <field id="sst_m"    unit="degC" /> 
    … …  
    422445 
    423446   </field_group> 
    424     
     447 
    425448 
    426449      </field_group> <!-- SBC --> 
    427        
     450 
    428451      <!-- ABL --> 
    429452      <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 
    … …  
    440463          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
    441464          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
    442           <field id="coeft"      long_name="ABL nudging coefficient"       standard_name="coeft"          unit=""         /> 
     465          <field id="u_geo"      long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     466          <field id="v_geo"      long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
    443467          <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    /> 
    444468          <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     /> 
    445469          <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     /> 
    446           <field id="mxl_abl"    long_name="ABL mixing length"             standard_name="abl_mxl"        unit="m"        /> 
     470          <field id="mxlm_abl"   long_name="ABL master mixing length"      standard_name="abl_mxlm"       unit="m"        /> 
     471          <field id="mxld_abl"   long_name="ABL dissipative mixing length" standard_name="abl_mxld"       unit="m"        /> 
    447472   </field_group> 
    448473 
    … …  
    456481          <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
    457482          <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
    458           <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field>  
     483          <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 
    459484          <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
    460485          <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
    … …  
    462487          <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
    463488          <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
    464           <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field>  
     489          <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 
    465490   </field_group> 
    466491 
    467492      </field_group> <!-- ABL --> 
    468493 
    469        
     494 
    470495      <!-- U grid --> 
    471        
     496 
    472497      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    473498        <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               /> 
    … …  
    478503      <field id="uoce_e3u"      long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
    479504      <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical"  field_ref="uoce_e3u"    unit="m3/s"       grid_ref="grid_U_vsum"/> 
    480         <field id="uocetr_vsum"   long_name="ocean transport along i-axis  summed on the vertical"         field_ref="e2u"       unit="m3/s"> this * uoce_e3u_vsum  </field>  
     505        <field id="uocetr_vsum"   long_name="ocean transport along i-axis  summed on the vertical"         field_ref="e2u"       unit="m3/s"> this * uoce_e3u_vsum  </field> 
    481506 
    482507        <field id="uocetr_vsum_op"    long_name="ocean current along i-axis * e3u * e2u summed on the vertical"  read_access="true"  freq_op="1mo"    field_ref="e2u"       unit="m3/s"> @uocetr_vsum </field> 
    483         <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo"  unit="m3/s" />  
    484         <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 
     508        <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo"  unit="m3/s" /> 
     509        <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rho0 </field> 
    485510 
    486511 
    … …  
    534559        <field id="udiff_salttr"  long_name="ocean diffusion salt transport along i-axis"                     standard_name="ocean_salt_x_transport_due_to_diffusion"        unit="1e-3*kg/s"                 /> 
    535560      </field_group> 
    536        
     561 
    537562      <!-- V grid --> 
    538        
     563 
    539564      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    540565        <field id="e1v"          long_name="V-cell width in longitudinal direction"                 standard_name="cell_width"                  unit="m"                              /> 
    … …  
    593618        <field id="vdiff_salttr"  long_name="ocean diffusion salt transport along j-axis"   standard_name="ocean_salt_y_transport_due_to_diffusion"         unit="1e-3*kg/s"                 /> 
    594619      </field_group> 
    595        
     620 
    596621      <!-- W grid --> 
    597        
     622 
    598623      <field_group id="grid_W" grid_ref="grid_W_3D"> 
    599624        <field id="e3w"          long_name="W-cell thickness"                              standard_name="cell_thickness"                         unit="m"    /> 
    600625   <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
    601    <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field>   
     626   <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field> 
    602627        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                                                   unit="m3/s" /> 
    603628 
    … …  
    609634 
    610635   <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    611         <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field>      
     636        <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field> 
    612637        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"         standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    613638        <field id="avm"          long_name="vertical eddy viscosity"                        standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
    614         <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field>  
     639        <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field> 
    615640 
    616641        <!-- avs: /= avt with ln_zdfddm=T --> 
    617642        <field id="avs"          long_name="salt vertical eddy diffusivity"                 standard_name="ocean_vertical_salt_diffusivity"       unit="m2/s" /> 
    618         <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field>  
     643        <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field> 
    619644   <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    620645 
    621646        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    622647        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    623         <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field>  
     648        <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field> 
    624649   <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    625650 
    … …  
    634659        <field id="wstokes"      long_name="Stokes Drift vertical velocity"                 standard_name="upward_StokesDrift_velocity"   unit="m/s" /> 
    635660 
    636         <!-- variables available with diaar5 -->    
     661        <!-- variables available with diaar5 --> 
    637662        <field id="w_masstr"     long_name="vertical mass transport"                        standard_name="upward_ocean_mass_transport"             unit="kg/s"   /> 
    638663        <field id="w_masstr2"    long_name="square of vertical mass transport"              standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    639664 
     665        <!-- EOS --> 
     666        <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" /> 
     667 
    640668      </field_group> 
    641        
     669 
    642670      <!-- F grid --> 
    643671      <!-- AGRIF sponge --> 
    … …  
    683711         <field id="uocetr_vsum_section"  long_name="Total 2D transport in i-direction"               field_ref="uoce_e3u_ave_vsum"    grid_ref="grid_U_scalar"  detect_missing_value="true"> this * e2u </field> 
    684712         <field id="uocetr_strait"        long_name="Total transport across lines in i-direction"     field_ref="uocetr_vsum_section"  grid_ref="grid_U_4strait" /> 
    685          <field id="u_masstr_strait"      long_name="Sea water transport across line in i-direction"  field_ref="uocetr_strait"        grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rau0 </field> 
     713         <field id="u_masstr_strait"      long_name="Sea water transport across line in i-direction"  field_ref="uocetr_strait"        grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rho0 </field> 
    686714 
    687715         <field id="voce_e3v_ave"         long_name="Monthly average of v*e3v"                        field_ref="voce_e3v"                    freq_op="1mo"   freq_offset="_reset_" > @voce_e3v </field> 
    … …  
    689717         <field id="vocetr_vsum_section"  long_name="Total 2D transport of in j-direction"            field_ref="voce_e3v_ave_vsum"    grid_ref="grid_V_scalar"  detect_missing_value="true"> this * e1v </field> 
    690718         <field id="vocetr_strait"        long_name="Total transport across lines in j-direction"     field_ref="vocetr_vsum_section"  grid_ref="grid_V_4strait"  /> 
    691          <field id="v_masstr_strait"      long_name="Sea water transport across line in j-direction"  field_ref="vocetr_strait"        grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rau0 </field> 
     719         <field id="v_masstr_strait"      long_name="Sea water transport across line in j-direction"  field_ref="vocetr_strait"        grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rho0 </field> 
    692720 
    693721         <field id="masstr_strait"        long_name="Sea water transport across line"                                                  grid_ref="grid_4strait"  > u_masstr_strait + v_masstr_strait </field> 
    694722      </field_group> 
    695723 
    696        
    697724      <!-- variables available with ln_floats --> 
    698725 
    … …  
    709736      <!-- variables available with iceberg trajectories --> 
    710737 
    711       <field_group id="icbvar" domain_ref="grid_T"  >  
     738      <field_group id="icbvar" domain_ref="grid_T"  > 
    712739        <field id="berg_melt"          long_name="icb melt rate of icebergs"                       unit="kg/m2/s"                    /> 
    713740        <field id="berg_melt_hcflx"    long_name="icb heat flux to ocean due to melting heat content"   unit="J/m2/s"                /> 
    … …  
    727754      </field_group> 
    728755 
    729       <!-- Poleward transport : ptr -->      
    730       <field_group id="diaptr" >   
     756      <!-- Poleward transport : ptr --> 
     757      <field_group id="diaptr" > 
    731758        <field id="zomsf"         long_name="Overturning Stream-Function : All basins"                     unit="Sv"         grid_ref="grid_znl_W_3D" /> 
    732759        <field id="zotem"         long_name="Zonal Mean Temperature : All basins"                          unit="degree_C"   grid_ref="grid_znl_T_3D" /> 
    … …  
    736763        <field id="sopstove"      long_name="Overturning Salt Transport: All basins"                       unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    737764        <field id="sophtbtr"      long_name="Barotropic Heat Transport: All basins"                        unit="PW"         grid_ref="grid_znl_T_2D" /> 
    738         <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" />  
     765        <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    739766        <field id="sophtadv"      long_name="Advective Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
    740767        <field id="sopstadv"      long_name="Advective Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    … …  
    742769        <field id="sopstldf"      long_name="Diffusive Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    743770        <field id="sophtvtr"      long_name="Heat Transport : All basins"                                  unit="PW"         grid_ref="grid_znl_T_2D" /> 
    744         <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" />   
     771        <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    745772        <field id="sophteiv"      long_name="Heat Transport from mesoscale eddy advection: All basins"     unit="PW"         grid_ref="grid_znl_T_2D" /> 
    746773        <field id="sopsteiv"      long_name="Salt Transport from mesoscale eddy advection : All basins"    unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    … …  
    758785 
    759786 
    760     <!--  
     787    <!-- 
    761788============================================================================================================ 
    762789                  Physical ocean model trend diagnostics : temperature, KE, PE, momentum 
    … …  
    899926     <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                    unit="W/s^3"                        /> 
    900927     <field id="ketrd_tau"     long_name="ke-trend: wind stress "                           unit="W/s^3"   grid_ref="grid_T_2D" /> 
    901      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        />    
    902      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        />    
    903      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        />   
     928     <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        /> 
     929     <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        /> 
     930     <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        /> 
    904931     <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        /> 
    905      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />    
     932     <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        /> 
    906933 
    907934     <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 
    908      <field id="dispkexyfo"    long_name="KE-trend: lateral  mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction"                   unit="W/m^2" grid_ref="grid_T_2D" />    
    909      <field id="dispkevfo"     long_name="KE-trend: vertical mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction"             unit="W/m^2" grid_ref="grid_T_2D" />    
     935     <field id="dispkexyfo"    long_name="KE-trend: lateral  mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction"                   unit="W/m^2" grid_ref="grid_T_2D" /> 
     936     <field id="dispkevfo"     long_name="KE-trend: vertical mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction"             unit="W/m^2" grid_ref="grid_T_2D" /> 
    910937     <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 
    911      <field id="eketrd_eiv"    long_name="EKE-trend due to parameterized eddy advection"   standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection"   unit="W/m^2" grid_ref="grid_T_2D" />    
     938     <field id="eketrd_eiv"    long_name="EKE-trend due to parameterized eddy advection"   standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection"   unit="W/m^2" grid_ref="grid_T_2D" /> 
    912939 
    913940     <!-- variables available with ln_PE_trd --> 
    … …  
    926953     <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"         unit="W/m^3"                        /> 
    927954     <field id="petrd_atf"     long_name="pe-trend: asselin time filter"        unit="W/m^3"                        /> 
    928      <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            />    
    929      <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       />    
    930      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          />    
     955     <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            /> 
     956     <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       /> 
     957     <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          /> 
    931958   </field_group> 
    932959 
    … …  
    945972     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    946973     <field id="utrd_tau"       long_name="i-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_U_2D" /> 
    947      <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    948      <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    949      <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    950      <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        />    
     974     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     975     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     976     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     977     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    951978   </field_group> 
    952979 
    … …  
    965992     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    966993     <field id="vtrd_tau"       long_name="j-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_V_2D" /> 
    967      <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    968      <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    969      <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    970      <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        />    
     994     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     995     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     996     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     997     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    971998   </field_group> 
    972999 
    9731000 
    974     <!--  
     1001    <!-- 
    9751002============================================================================================================ 
    9761003                                        Definitions for iodef_demo.xml 
    … …  
    9901017      <field field_ref="strd_zdfp_li"    name="osaltdiff" /> 
    9911018    </field_group> 
    992      
     1019 
    9931020    <field_group id="mooring" > 
    9941021      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"      /> 
    … …  
    9991026      <field field_ref="avt"          name="difvho"   long_name="ocean_vertical_heat_diffusivity"      /> 
    10001027      <field field_ref="avm"          name="difvmo"   long_name="ocean_vertical_momentum_diffusivity"  /> 
    1001        
     1028 
    10021029      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
    10031030      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
    … …  
    10471074      <field field_ref="BLT"          name="blt"      long_name="Barrier Layer Thickness"                       /> 
    10481075    </field_group> 
    1049      
     1076 
    10501077    <field_group id="groupU" > 
    10511078      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    10521079      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    10531080    </field_group> 
    1054      
     1081 
    10551082    <field_group id="groupV" > 
    10561083      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    10571084      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    10581085    </field_group> 
    1059      
     1086 
    10601087    <field_group id="groupW" > 
    10611088      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  /> 
    … …  
    11001127    </field_group> 
    11011128 
    1102     <!--  
     1129    <!-- 
    11031130============================================================================================================ 
    11041131    --> 
    1105     <!-- output variables for my configuration (example) -->  
    1106      
     1132    <!-- output variables for my configuration (example) --> 
     1133 
    11071134    <field_group id="myvarOCE" > 
    1108       <!-- grid T -->  
     1135      <!-- grid T --> 
    11091136      <field field_ref="e3t"          name="e3t"      long_name="vertical scale factor"           /> 
    11101137      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"         /> 
    11111138      <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"            /> 
    11121139      <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"  /> 
    1113        
    1114       <!-- grid U -->  
     1140 
     1141      <!-- grid U --> 
    11151142      <field field_ref="e3u"          name="e3u"     long_name="vertical scale factor"            /> 
    11161143      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"           /> 
    1117        
    1118       <!-- grid V -->  
     1144 
     1145      <!-- grid V --> 
    11191146      <field field_ref="e3v"          name="e3v"     long_name="vertical scale factor"            /> 
    1120       <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           />      
    1121     </field_group>     
     1147      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           /> 
     1148    </field_group> 
    11221149 
    11231150   </field_definition> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SHARED/namelist_ice_ref

    r12377 r13540  
    4343   ln_cat_usr       = .false.         !  ice categories are defined by rn_catbnd below (m) 
    4444      rn_catbnd     =   0.,0.45,1.1,2.1,3.7,6.0   
    45    rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     45   rn_himin         =   0.1           !  minimum ice thickness (m) allowed 
     46   rn_himax         =  99.0           !  maximum ice thickness (m) allowed 
    4647/ 
    4748!------------------------------------------------------------------------------ 
    … …  
    5657   rn_ishlat        =   2.            !  lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
    5758   ln_landfast_L16  = .false.         !  landfast: parameterization from Lemieux 2016 
    58       rn_depfra     =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast 
     59      rn_lf_depfra  =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast 
    5960                                      !          recommended range: [0.1 ; 0.25] 
    60       rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
    61       rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
    62       rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??] 
     61      rn_lf_bfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
     62      rn_lf_relax   =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
     63      rn_lf_tensile =   0.05          !        isotropic tensile strength [0-0.5??] 
    6364/ 
    6465!------------------------------------------------------------------------------ 
    … …  
    9192!------------------------------------------------------------------------------ 
    9293   ln_rhg_EVP       = .true.          !  EVP rheology 
    93       ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
     94      ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    9495      rn_creepl     =   2.0e-9        !     creep limit [1/s] 
    9596      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
    96       nn_nevp       = 120             !     number of EVP subcycles                              
     97      nn_nevp       = 100             !     number of EVP subcycles                              
    9798      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    98                                       !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     99                                      !        advised value: 1/3 (nn_nevp=100) or 1/9 (nn_nevp=300) 
     100   nn_rhg_chkcvg    =   0             !  check convergence of rheology 
     101                                      !     = 0  no check 
     102                                      !     = 1  check at the main time step (output xml: uice_cvg) 
     103                                      !     = 2  check at both main and rheology time steps (additional output: ice_cvg.nc) 
     104                                      !          this option 2 asks a lot of communications between cpu 
    99105/ 
    100106!------------------------------------------------------------------------------ 
    101107&namdyn_adv     !   Ice advection 
    102108!------------------------------------------------------------------------------ 
    103    ln_adv_Pra       = .true.         !  Advection scheme (Prather) 
    104    ln_adv_UMx       = .false.          !  Advection scheme (Ultimate-Macho) 
     109   ln_adv_Pra       = .true.          !  Advection scheme (Prather) 
     110   ln_adv_UMx       = .false.         !  Advection scheme (Ultimate-Macho) 
    105111      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    106112/ 
    … …  
    109115!------------------------------------------------------------------------------ 
    110116   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-) 
    111    rn_blow_s        =   0.66          !  mesure of snow blowing into the leads 
     117   nn_snwfra        =   2             !  calculate the fraction of ice covered by snow (for zdf and albedo) 
     118                                      !     = 0  fraction = 1 (if snow) or 0 (if no snow) 
     119                                      !     = 1  fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 
     120                                      !     = 2  fraction = hsnw / (hsnw+0.02)    [CICE formulation] 
     121   rn_snwblow       =   0.66          !  mesure of snow blowing into the leads 
    112122                                      !     = 1 => no snow blowing, < 1 => some snow blowing 
    113123   nn_flxdist       =  -1             !  Redistribute heat flux over ice categories 
    … …  
    118128   ln_cndflx        = .false.         !  Use conduction flux as surface boundary conditions (i.e. for Jules coupling) 
    119129      ln_cndemulate = .false.         !     emulate conduction flux (if not provided in the inputs) 
     130   nn_qtrice        =   1             !  Solar flux transmitted thru the surface scattering layer: 
     131                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow)  
     132                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 
    120133/ 
    121134!------------------------------------------------------------------------------ 
    … …  
    126139   ln_icedO         = .true.          !  activate ice growth in open-water (T) or not (F) 
    127140   ln_icedS         = .true.          !  activate brine drainage (T) or not (F) 
     141   ! 
     142   ln_leadhfx       = .true.          !  heat in the leads is used to melt sea-ice before warming the ocean 
    128143/ 
    129144!------------------------------------------------------------------------------ 
    … …  
    135150   rn_cnd_s         =   0.31          !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 
    136151                                      !     Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 
    137    rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice [1/m] 
     152   rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice                     [1/m] 
     153   rn_kappa_s       =  10.0           !  nn_qtrice = 0: radiation attenuation coefficient in snow         [1/m] 
     154   rn_kappa_smlt    =   7.0           !  nn_qtrice = 1: radiation attenuation coefficient in melting snow [1/m] 
     155   rn_kappa_sdry    =  10.0           !                 radiation attenuation coefficient in dry snow     [1/m] 
     156   ln_zdf_chkcvg    = .false.         !  check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp) 
    138157/ 
    139158!------------------------------------------------------------------------------ 
    … …  
    175194&namthd_pnd     !   Melt ponds 
    176195!------------------------------------------------------------------------------ 
    177    ln_pnd           = .false.         !  activate melt ponds or not 
    178      ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
    179      ln_pnd_CST     = .false.         !  activate constant  melt ponds 
    180        rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
    181        rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
    182      ln_pnd_alb     = .false.         !  melt ponds affect albedo or not 
     196   ln_pnd            = .true.         !  activate melt ponds or not 
     197      ln_pnd_LEV     = .true.         !  level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012) 
     198         rn_apnd_min =   0.15         !     minimum ice fraction that contributes to melt pond. range: 0.0 -- 0.15 ?? 
     199         rn_apnd_max =   0.85         !     maximum ice fraction that contributes to melt pond. range: 0.7 -- 0.85 ?? 
     200      ln_pnd_CST     = .false.        !  constant  melt ponds 
     201         rn_apnd     =   0.2          !     prescribed pond fraction, at Tsu=0 degC 
     202         rn_hpnd     =   0.05         !     prescribed pond depth,    at Tsu=0 degC 
     203      ln_pnd_lids    = .true.         !  frozen lids on top of the ponds (only for ln_pnd_LEV) 
     204      ln_pnd_alb     = .true.         !  effect of melt ponds on ice albedo 
    183205/ 
    184206!------------------------------------------------------------------------------ 
    … …  
    186208!------------------------------------------------------------------------------ 
    187209   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
    188    ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F) 
     210   nn_iceini_file   =   0             !     0 = Initialise sea ice based on SSTs 
     211                                      !     1 = Initialise sea ice from single category netcdf file 
     212                                      !     2 = Initialise sea ice from multi category restart file 
    189213   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze) 
    190214   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North 
    … …  
    206230   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North 
    207231   rn_hpd_ini_s     =   0.05          !        "            "             South 
    208    ! -- for ln_iceini_file = T 
    209    sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    210    sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    211    sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    212    sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    213    sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    214    sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    215    sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     232   rn_hld_ini_n     =   0.0           !  initial pond lid depth      (m), North 
     233   rn_hld_ini_s     =   0.0           !        "            "             South 
     234   ! -- for nn_iceini_file = 1 
     235   sn_hti = 'Ice_initialization'    , -12. ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     236   sn_hts = 'Ice_initialization'    , -12. ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     237   sn_ati = 'Ice_initialization'    , -12. ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     238   sn_smi = 'Ice_initialization'    , -12. ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     239   sn_tmi = 'Ice_initialization'    , -12. ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     240   sn_tsu = 'Ice_initialization'    , -12. ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     241   sn_tms = 'NOT USED'              , -12. ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    216242   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 
    217    sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    218    sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     243   sn_apd = 'NOT USED'              , -12. ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     244   sn_hpd = 'NOT USED'              , -12. ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     245   sn_hld = 'NOT USED'              , -12. ,'hld'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    219246   cn_dir='./' 
    220247/ 
    … …  
    238265   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F) 
    239266   ln_icectl        = .false.         !  ice points output for debug (T or F) 
    240    iiceprt          =  10             !  i-index for debug 
    241    jiceprt          =  10             !  j-index for debug 
    242 / 
     267      iiceprt       =  10             !     i-index for debug 
     268      jiceprt       =  10             !     j-index for debug 
     269/ 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SHARED/namelist_pisces_ref

    r12377 r13540  
    352352! 
    353353   cn_dir      = './'      !  root directory for the location of the dynamical files 
    354    ln_ironsed  =  .true.   ! boolean for Fe input from sediments 
    355    ln_ironice  =  .true.   ! boolean for Fe input from sea ice 
    356    ln_hydrofe  =  .true.   ! boolean for from hydrothermal vents 
     354   ln_ironsed  =  .false.   ! boolean for Fe input from sediments 
     355   ln_ironice  =  .false.   ! boolean for Fe input from sea ice 
     356   ln_hydrofe  =  .false.   ! boolean for from hydrothermal vents 
    357357   sedfeinput  =  2.e-9    ! Coastal release of Iron 
    358358   distcoast   =  5.e3     ! Distance off the coast for Iron from sediments 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SHARED/namelist_ref

    r12511 r13540  
    7272   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    7373   ! 
    74    rn_Dt      = 5400.     !  time step for the dynamics and tracer 
     74   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    7575   rn_atfp     =    0.1    !  asselin time filter parameter 
    7676   ! 
    … …  
    217217   nn_ice      = 0         !  =0 no ice boundary condition 
    218218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    219       !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    220       !                    !          except in AGRIF zoom where it has to be specified 
     219      !                    !  =2 or 3 for SI3 and CICE, respectively 
    221220   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    222221      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    … …  
    269268      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
    270269      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
    271       !                     !  - module of the mean stress" data 
     270      ln_crt_fbk = .false.     !  Add surface current feedback to the wind stress (Renault et al. 2020, doi: 10.1029/2019MS001715) 
     271         rn_stau_a = -2.9e-3   !     Alpha from eq. 10: Stau = Alpha * Wnd + Beta 
     272         rn_stau_b =  8.0e-3   !     Beta  
    272273      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
    273274      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
    274       rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity  
    275       !                     !  used to calculate the wind stress 
    276       !                     ! (0. => absolute or 1. => relative winds) 
    277275      ln_skin_cs = .false.  !  use the cool-skin parameterization 
    278276      ln_skin_wl = .false.  !  use the warm-layer parameterization 
    … …  
    281279      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
    282280      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     281      ln_tpot     = .true.  !!GS: compute potential temperature or not 
    283282   ! 
    284283   cn_dir      = './'      !  root directory for the bulk data location 
    … …  
    292291   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    293292   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    294    sn_hpgi     = 'NONE'                       ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
    295    sn_hpgj     = 'NONE'                       ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    296293   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    297294   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    298295   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     296   sn_uoatm    = 'NOT USED'                   ,    6.        , 'UOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Uoceatm', '' 
     297   sn_voatm    = 'NOT USED'                   ,    6.        , 'VOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Voceatm', '' 
     298   sn_cc       = 'NOT USED'                   ,   24.        , 'CC'      ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     299   sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     300   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    299301/ 
    300302!----------------------------------------------------------------------- 
    … …  
    309311   cn_ablrst_outdir = "."             !  directory to write output abl restarts 
    310312 
     313   ln_rstart_abl  = .false. 
    311314   ln_hpgls_frc   = .false. 
    312315   ln_geos_winds  = .false. 
    313    nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring 
     316   ln_smth_pblh   = .false. 
     317   nn_dyn_restore = 0         ! restoring option for dynamical ABL variables: = 0 no restoring 
    314318                              !                                               = 1 equatorial restoring 
    315319                              !                                               = 2 global restoring 
    316    rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
    317    rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
    318    rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
    319    rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     320   rn_ldyn_min   =  4.5       ! dynamics nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     321   rn_ldyn_max   =  1.5       ! dynamics nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
     322   rn_ltra_min   =  4.5       ! tracers  nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     323   rn_ltra_max   =  1.5       ! tracers  nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
    320324   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
    321325                              !                = 1 length-scale based on the distance to the PBL height 
    322326                              !                = 2 Bougeault & Lacarrere 89 length-scale 
    323    rn_Cm         = 0.0667     ! 0.126 in MesoNH 
    324    rn_Ct         = 0.1667     ! 0.143 in MesoNH 
    325    rn_Ce         = 0.4        ! 0.4   in MesoNH 
    326    rn_Ceps       = 0.7        ! 0.85  in MesoNH 
    327    rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented) 
    328    rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities) 
     327                              ! CBR00  ! CCH02  ! MesoNH ! 
     328   rn_Cm          = 0.0667    ! 0.0667 ! 0.1260 ! 0.1260 ! 
     329   rn_Ct          = 0.1667    ! 0.1667 ! 0.1430 ! 0.1430 ! 
     330   rn_Ce          = 0.40      ! 0.40   ! 0.34   ! 0.40   ! 
     331   rn_Ceps        = 0.700     ! 0.700  ! 0.845  ! 0.850  ! 
     332   rn_Ric         = 0.139     ! 0.139  ! 0.143  !   ?    ! Critical Richardson number (to compute PBL height and diffusivities) 
     333   rn_Rod         = 0.15      ! c0 in RMCA17 mixing length formulation (not yet implemented) 
    329334/ 
    330335!----------------------------------------------------------------------- 
    331336&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    332337!----------------------------------------------------------------------- 
    333    nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
    334    ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
    335    !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    336    nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
     338   nn_cplmodel       =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
     339   ln_usecplmask     = .false. !  use a coupling mask file to merge data received from several models 
     340   !                           !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
     341   ln_scale_ice_flux = .false. !  use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 
     342   nn_cats_cpl       =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
    337343   !_____________!__________________________!____________!_____________!______________________!________! 
    338344   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
    … …  
    540546         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    541547         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    542          sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    543          sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     548         sn_isfpar_zmax = 'isfmlt_par',       0.       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     549         sn_isfpar_zmin = 'isfmlt_par',       0.       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    544550         !* 'spe' and 'oasis' case 
    545          sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     551         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    546552         !* 'bg03' case 
    547          sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     553         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    548554      ! 
    549555      ! ---------------- ice sheet coupling ------------------------------- 
    … …  
    639645&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    640646!----------------------------------------------------------------------- 
    641    ln_agrif_2way = .true.  !  activate two way nesting 
    642    ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    643    rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    644    rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
    645    rn_trelax_tra = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
    646    rn_trelax_dyn = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
    647    ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
     647   ln_agrif_2way   = .true.  !  activate two way nesting 
     648   ln_init_chfrpar = .false. !  initialize child grids from parent 
     649   ln_spc_dyn      = .true.  !  use 0 as special value for dynamics 
     650   rn_sponge_tra   = 0.002   !  coefficient for tracer   sponge layer [] 
     651   rn_sponge_dyn   = 0.002   !  coefficient for dynamics sponge layer [] 
     652   rn_trelax_tra   = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
     653   rn_trelax_dyn   = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
     654   ln_chk_bathy    = .false. !  =T  check the parent bathymetry 
    648655/ 
    649656!----------------------------------------------------------------------- 
    … …  
    727734   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    728735   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     736   bn_hil      = 'NOT USED'              ,         24.       , 'sihlid'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    729737   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 
    730738   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice 
    … …  
    733741   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            -- 
    734742   rn_ice_hpnd = 0.05         !       --   pond depth                          -- 
     743   rn_ice_hlid = 0.0          !       --   pond lid depth                      -- 
    735744/ 
    736745!----------------------------------------------------------------------- 
    … …  
    745754!!                                                                    !! 
    746755!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    747 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    748 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     756!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     757!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    749758!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    750759!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    754763&namdrg        !   top/bottom drag coefficient                          (default: NO selection) 
    755764!----------------------------------------------------------------------- 
    756    ln_OFF      = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
     765   ln_drg_OFF  = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
    757766   ln_lin      = .false.   !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top) 
    758767   ln_non_lin  = .false.   !  non-linear  drag: Cd = Cd0 |U| 
    … …  
    760769   ! 
    761770   ln_drgimp   = .true.    !  implicit top/bottom friction flag 
    762 / 
    763 !----------------------------------------------------------------------- 
    764 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
     771      ln_drgice_imp = .true. ! implicit ice-ocean drag 
     772/ 
     773!----------------------------------------------------------------------- 
     774&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
    765775!----------------------------------------------------------------------- 
    766776   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    … …  
    773783/ 
    774784!----------------------------------------------------------------------- 
    775 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     785&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    776786!----------------------------------------------------------------------- 
    777787   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    … …  
    826836                                 ! 
    827837   !                     ! S-EOS coefficients (ln_seos=T): 
    828    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     838   !                             !  rd(T,S,Z)*rho0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    829839   rn_a0       =  1.6550e-1      !  thermal expension coefficient 
    830840   rn_b0       =  7.6554e-1      !  saline  expension coefficient 
    … …  
    9991009      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    10001010         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    1001          nn_e      = 30          ! =F : the number of sub-step in rn_Dt seconds 
     1011         nn_e         = 30          ! =F : the number of sub-step in rn_Dt seconds 
    10021012      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F) 
    10031013/ 
    … …  
    11301140   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it) 
    11311141   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    1132    nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
     1142   nn_mxl      =   3       !  mixing length: = 0 bounded by the distance to surface and bottom 
    11331143   !                       !                 = 1 bounded by the local vertical scale factor 
    11341144   !                       !                 = 2 first vertical derivative of mixing length bounded by 1 
    11351145   !                       !                 = 3 as =2 with distinct dissipative an mixing length scale 
    11361146   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     1147      nn_mxlice    = 2        ! type of scaling under sea-ice 
     1148                              !    = 0 no scaling under sea-ice 
     1149                              !    = 1 scaling with constant sea-ice thickness 
     1150                              !    = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
     1151                              !    = 3 scaling with maximum sea-ice thickness 
     1152      rn_mxlice   = 10.       ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 
    11371153   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
    1138    ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE 
    11391154   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    11401155      rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    … …  
    11471162                              !        = 0  constant 10 m length scale 
    11481163                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1149       rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 
     1164   nn_eice     =   1       !  attenutaion of langmuir & surface wave breaking under ice 
     1165   !                       !           = 0 no impact of ice cover on langmuir & surface wave breaking 
     1166   !                       !           = 1 weigthed by 1-TANH(10*fr_i) 
     1167   !                       !           = 2 weighted by 1-fr_i 
     1168   !                       !           = 3 weighted by 1-MIN(1,4*fr_i)    
    11501169/ 
    11511170!----------------------------------------------------------------------- 
    … …  
    11601179   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    11611180   rn_hsro       =  0.02   !  Minimum surface roughness 
     1181   rn_hsri       =  0.03   !  Ice-ocean roughness 
    11621182   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1) 
    11631183   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3) 
    1164    !                             ! =3 requires ln_wave=T 
     1184   !                       !     = 3 requires ln_wave=T 
     1185   nn_z0_ice     =   1     !  attenutaion of surface wave breaking under ice 
     1186   !                       !           = 0 no impact of ice cover 
     1187   !                       !           = 1 roughness uses rn_hsri and is weigthed by 1-TANH(10*fr_i) 
     1188   !                       !           = 2 roughness uses rn_hsri and is weighted by 1-fr_i 
     1189   !                       !           = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i) 
    11651190   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
    11661191   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
    … …  
    11931218   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    11941219   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1195 / 
    1196  
     1220 
     1221   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     1222   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     1223   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     1224   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     1225   sn_mpb      = 'NOT USED'              , -12.              , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1226   sn_mpp      = 'NOT USED'              , -12.              , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1227   sn_mpc      = 'NOT USED'              , -12.              , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1228   sn_dsb      = 'NOT USED'              , -12.              , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1229   sn_dsc      = 'NOT USED'              , -12.              , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1230/ 
    11971231!!====================================================================== 
    11981232!!                  ***  Diagnostics namelists  ***                   !! 
    … …  
    13821416   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T 
    13831417   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T 
     1418   nn_hls      =   1       !  halo width (applies to both rows and columns) 
    13841419/ 
    13851420!----------------------------------------------------------------------- 
    13861421&namctl        !   Control prints                                       (default: OFF) 
    13871422!----------------------------------------------------------------------- 
    1388    sn_cfctl%l_glochk = .FALSE.    ! Range sanity checks are local (F) or global (T). Set T for debugging only 
    1389    sn_cfctl%l_allon  = .FALSE.    ! IF T activate all options. If F deactivate all unless l_config is T 
    1390      sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    1391        sn_cfctl%l_runstat = .TRUE.  ! switches and which areas produce reports with the proc integer settings. 
    1392        sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    1393        sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    1394        sn_cfctl%l_layout  = .FALSE. ! 
    1395        sn_cfctl%l_prtctl  = .FALSE. ! 
    1396        sn_cfctl%l_prttrc  = .FALSE. ! 
    1397        sn_cfctl%l_oasout  = .FALSE. ! 
    1398        sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    1399        sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
    1400        sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
    1401        sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
    1402    nn_print    =    0      !  level of print (0 no extra print) 
    1403    nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
    1404    nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
    1405    nn_jctls    =    0      !  start j indice of control               over a subdomain) 
    1406    nn_jctle    =    0      !  end   j indice of control 
    1407    nn_isplt    =    1      !  number of processors in i-direction 
    1408    nn_jsplt    =    1      !  number of processors in j-direction 
    1409    ln_timing   = .false.   !  timing by routine write out in timing.output file 
    1410    ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
     1423   sn_cfctl%l_runstat = .TRUE.    ! switches and which areas produce reports with the proc integer settings. 
     1424   sn_cfctl%l_trcstat = .FALSE.   ! The default settings for the proc integers should ensure 
     1425   sn_cfctl%l_oceout  = .FALSE.   ! that  all areas report. 
     1426   sn_cfctl%l_layout  = .FALSE.   ! 
     1427   sn_cfctl%l_prtctl  = .FALSE.   ! 
     1428   sn_cfctl%l_prttrc  = .FALSE.   ! 
     1429   sn_cfctl%l_oasout  = .FALSE.   ! 
     1430   sn_cfctl%procmin   = 0         ! Minimum area number for reporting [default:0] 
     1431   sn_cfctl%procmax   = 1000000   ! Maximum area number for reporting [default:1000000] 
     1432   sn_cfctl%procincr  = 1         ! Increment for optional subsetting of areas [default:1] 
     1433   sn_cfctl%ptimincr  = 1         ! Timestep increment for writing time step progress info 
     1434   nn_ictls    =    0             !  start i indice of control sum (use to compare mono versus 
     1435   nn_ictle    =    0             !  end   i indice of control sum        multi processor runs 
     1436   nn_jctls    =    0             !  start j indice of control               over a subdomain) 
     1437   nn_jctle    =    0             !  end   j indice of control 
     1438   nn_isplt    =    1             !  number of processors in i-direction 
     1439   nn_jsplt    =    1             !  number of processors in j-direction 
     1440   ln_timing   = .false.          !  timing by routine write out in timing.output file 
     1441   ln_diacfl   = .false.          !  CFL diagnostics write out in cfl_diagnostics.ascii 
    14111442/ 
    14121443!----------------------------------------------------------------------- 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SPITZ12/EXPREF/context_nemo.xml

    r12276 r13540  
    1111       <variable id="ref_month" type="int"> 01 </variable> 
    1212       <variable id="ref_day"   type="int"> 01 </variable> 
    13        <variable id="rau0"      type="float" > 1026.0 </variable> 
     13       <variable id="rho0"      type="float" > 1026.0 </variable> 
    1414       <variable id="cpocean"   type="float" > 3991.86795711963 </variable> 
    1515       <variable id="convSpsu"  type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SPITZ12/EXPREF/namelist_cfg

    r12511 r13540  
    205205!!                                                                    !! 
    206206!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    207 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    208 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     207!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     208!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    209209!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    210210!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    216216   ln_loglayer = .true.   !  logarithmic drag: Cd = vkarmn/log(z/z0) |U| 
    217217   ln_drgimp   = .true.   !  implicit top/bottom friction flag 
    218 / 
    219 !----------------------------------------------------------------------- 
    220 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     218      ln_drgice_imp = .true. ! implicit ice-ocean drag 
     219/ 
     220!----------------------------------------------------------------------- 
     221&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    221222!----------------------------------------------------------------------- 
    222223   rn_Cd0      =  2.5e-3   !  drag coefficient [-] 
    … …  
    339340   nn_havtb    =    1         !  horizontal shape for avtb (=1) or not (=0) 
    340341/ 
     342!----------------------------------------------------------------------- 
     343&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     344!----------------------------------------------------------------------- 
     345   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     346      nn_mxlice    = 0        ! type of scaling under sea-ice 
     347                              !    = 0 no scaling under sea-ice 
     348                              !    = 1 scaling with constant sea-ice thickness 
     349                              !    = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
     350                              !    = 3 scaling with maximum sea-ice thickness 
     351   nn_eice     =   0       !  attenutaion of langmuir & surface wave breaking under ice 
     352   !                       !           = 0 no impact of ice cover on langmuir & surface wave breaking 
     353   !                       !           = 1 weigthed by 1-TANH(10*fr_i) 
     354   !                       !           = 2 weighted by 1-fr_i 
     355   !                       !           = 3 weighted by 1-MIN(1,4*fr_i) 
     356/ 
    341357!!====================================================================== 
    342358!!                  ***  Diagnostics namelists  ***                   !! 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/SPITZ12/EXPREF/namelist_ice_cfg

    r11731 r13540  
    5555&namsbc         !   Ice surface boundary conditions 
    5656!------------------------------------------------------------------------------ 
     57   nn_snwfra        =   0             !  calculate the fraction of ice covered by snow (for zdf and albedo) 
     58                                      !     = 0  fraction = 1 (if snow) or 0 (if no snow) 
     59                                      !     = 1  fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 
     60                                      !     = 2  fraction = hsnw / (hsnw+0.02)    [CICE formulation] 
     61   nn_qtrice        =   0             !  Solar flux transmitted thru the surface scattering layer: 
     62                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow) 
     63                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 
    5764/ 
    5865!------------------------------------------------------------------------------ 
    … …  
    8188&namthd_pnd     !   Melt ponds 
    8289!------------------------------------------------------------------------------ 
    83    ln_pnd           = .true.          !  activate melt ponds or not 
    84      ln_pnd_H12     = .true.          !  activate evolutive melt ponds (from Holland et al 2012) 
    85      ln_pnd_alb     = .true.          !  melt ponds affect albedo or not 
     90   ln_pnd           = .false.          !  activate melt ponds or not 
     91     ln_pnd_LEV     = .false.          !  activate level ice melt ponds 
    8692/ 
    8793 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/WED025/EXPREF/context_nemo.xml

    r11487 r13540  
    99    <!-- Year of time origin for NetCDF files; defaults to 1800 --> 
    1010       <variable id="ref_year" type="int"   > 1800 </variable> 
    11        <variable id="rau0"     type="float" > 1026.0 </variable> 
     11       <variable id="rho0"     type="float" > 1026.0 </variable> 
    1212       <variable id="cpocean"  type="float" > 3991.86795711963 </variable> 
    1313       <variable id="convSpsu" type="float" > 0.99530670233846  </variable> 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/WED025/EXPREF/file_def_nemo-ice.xml

    r11844 r13540  
    7878     </file> 
    7979      
    80      <file id="file22" name_suffix="_SBC_scalar" description="scalar variables" enabled=".true." > 
    81        <!-- global contents --> 
    82        <field field_ref="ibgvol_tot"     grid_ref="grid_1point"   name="ibgvol_tot"   /> 
    83        <field field_ref="sbgvol_tot"     grid_ref="grid_1point"   name="sbgvol_tot"   /> 
    84        <field field_ref="ibgarea_tot"    grid_ref="grid_1point"   name="ibgarea_tot"  /> 
    85        <field field_ref="ibgsalt_tot"    grid_ref="grid_1point"   name="ibgsalt_tot"  /> 
    86        <field field_ref="ibgheat_tot"    grid_ref="grid_1point"   name="ibgheat_tot"  /> 
    87        <field field_ref="sbgheat_tot"    grid_ref="grid_1point"   name="sbgheat_tot"  /> 
    88         
    89        <!-- global drifts (conservation checks) --> 
    90        <field field_ref="ibgvolume"      grid_ref="grid_1point"   name="ibgvolume"    /> 
    91        <field field_ref="ibgsaltco"      grid_ref="grid_1point"   name="ibgsaltco"    /> 
    92        <field field_ref="ibgheatco"      grid_ref="grid_1point"   name="ibgheatco"    /> 
    93        <field field_ref="ibgheatfx"      grid_ref="grid_1point"   name="ibgheatfx"    /> 
    94         
    95        <!-- global forcings  --> 
    96        <field field_ref="ibgfrcvoltop"   grid_ref="grid_1point"   name="ibgfrcvoltop" /> 
    97        <field field_ref="ibgfrcvolbot"   grid_ref="grid_1point"   name="ibgfrcvolbot" /> 
    98        <field field_ref="ibgfrctemtop"   grid_ref="grid_1point"   name="ibgfrctemtop" /> 
    99        <field field_ref="ibgfrctembot"   grid_ref="grid_1point"   name="ibgfrctembot" /> 
    100        <field field_ref="ibgfrcsal"      grid_ref="grid_1point"   name="ibgfrcsal"    /> 
    101        <field field_ref="ibgfrchfxtop"   grid_ref="grid_1point"   name="ibgfrchfxtop" /> 
    102        <field field_ref="ibgfrchfxbot"   grid_ref="grid_1point"   name="ibgfrchfxbot" /> 
    103      </file> 
    104       
    10580   </file_group> 
    10681    

  • NEMO/branches/2020/r12377_ticket2386/cfgs/WED025/EXPREF/namelist_cfg

    r12511 r13540  
    55!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 
    66!!                                    namsbc_sas, namtra_qsr, namsbc_rnf, 
    7 !!                                    namsbc_isf, namsbc_iscpl, namsbc_apr,  
     7!!                                    namisf, namsbc_apr,  
    88!!                                    namsbc_ssr, namsbc_wave, namberg) 
    99!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    … …  
    3838   nn_it000    =   1       !  first time step 
    3939   nn_itend    =  26280    !  last  time step (std 5475) 
    40    nn_date0    = 19760301  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
     40   nn_date0    = 20000101  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    4141   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    4242      nn_rstctl   =    2      !  restart control ==> activated only if ln_rstart=T 
    … …  
    6161   ln_tsd_init = .true.          !  ocean initialisation 
    6262   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp) 
    63     
     63 
    6464   cn_dir      = './'      !  root directory for the T-S data location 
    65    !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    66    !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    67    !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    68    sn_tem = 'dta_temp_WED025'            ,         -12       , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    69    sn_sal = 'dta_sal_WED025'             ,         -12       , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     65   !___________!_____________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     66   !           !  file name          ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     67   !           !                     !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     68   sn_tem = 'WED025_init_JRA_200001.nc',       -12       , 'votemper',   .false.   , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     69   sn_sal = 'WED025_init_JRA_200001.nc',       -12       , 'vosaline',   .false.   , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    7070/ 
    7171!----------------------------------------------------------------------- 
    … …  
    116116   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    117117                     ! Sea-ice : 
    118    nn_ice      = 2         !  =0 no ice boundary condition     
     118   nn_ice      = 2         !  =0 no ice boundary condition 
    119119      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    120       !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    121       !                    !          except in AGRIF zoom where it has to be specified 
     120      !                    !  =2 or 3 for SI3 and CICE, respectively 
    122121   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    123122      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    124123                     ! Misc. options of sbc :  
    125124   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    126    ln_dm2dc    = .true.    !  daily mean to diurnal cycle on short wave 
     125   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    127126   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    128127   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
    … …  
    139138!----------------------------------------------------------------------- 
    140139   !                    !  bulk algorithm : 
    141    ln_NCAR     = .true.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     140   ln_NCAR     = .true.     ! "NCAR"      algorithm   (Large and Yeager 2008) 
    142141   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    143    ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    144    ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    145  
     142   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
     143   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
     144   ! 
    146145   cn_dir      = './'      !  root directory for the bulk data location 
    147146   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    148147   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
    149148   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    150    sn_wndi     = 'u10_core'              ,         6         , 'U_10_MOD',   .true.    , .false. , 'yearly'  , 'weights_bicubic_core.nc'           , 'Uwnd'   , '' 
    151    sn_wndj     = 'v10_core'              ,         6         , 'V_10_MOD',   .true.    , .false. , 'yearly'  , 'weights_bicubic_core.nc'           , 'Vwnd'   , '' 
    152    sn_qsr      = 'qsw_core'              ,        24         , 'SWDN_MOD',   .false.   , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    153    sn_qlw      = 'qlw_core'              ,        24         , 'LWDN_MOD',   .false.   , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    154    sn_tair     = 't10_core'              ,         6         , 'T_10_MOD',   .true.    , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    155    sn_humi     = 'q10_core'              ,         6         , 'Q_10_MOD',   .true.    , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    156    sn_prec     = 'precip_core'           ,        -1         , 'TPRECIP',    .true.    , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    157    sn_snow     = 'snow_core'             ,        -1         , 'SNOW'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
    158    sn_slp      = 'slp_core'              ,         6         , 'SLP'     ,   .true.    , .false. , 'yearly'  , 'weights_bilin_core.nc'             , ''       , '' 
     149   sn_wndi     = 'u10_JRA'              ,         3         , 'uas_10m' ,   .true.    , .false. , 'yearly'  , 'weights_bicubic_JRA.nc'           , 'Uwnd'   , '' 
     150   sn_wndj     = 'v10_JRA'              ,         3         , 'vas_10m' ,   .true.    , .false. , 'yearly'  , 'weights_bicubic_JRA.nc'           , 'Vwnd'   , '' 
     151   sn_qsr      = 'rsds_JRA'             ,         3         , 'rsds'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     152   sn_qlw      = 'rlds_JRA'             ,         3         , 'rlds'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     153   sn_tair     = 't10_JRA'              ,         3         , 'tas_10m' ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     154   sn_humi     = 'q10_JRA'              ,         3         , 'huss_10m',   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     155   sn_prec     = 'precip_JRA'           ,         3         , 'prto'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     156   sn_snow     = 'snow_JRA'             ,         3         , 'prsn'    ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
     157   sn_slp      = 'slp_JRA'              ,         3         , 'psl'     ,   .true.    , .false. , 'yearly'  , 'weights_bilin_JRA.nc'             , ''       , '' 
    159158/ 
    160159!----------------------------------------------------------------------- 
    … …  
    201200   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    202201   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    203    sn_rnf      = 'runoff_WED025'         ,  -1               , 'runoff'  ,   .true.    , .false., 'yearly'  , ''               , ''       , '' 
     202   sn_rnf      = 'WED025_icb'         ,  -1               , 'runoff'  ,   .true.    , .false., 'yearly'  , ''               , ''       , '' 
    204203/ 
    205204!----------------------------------------------------------------------- 
    … …  
    221220         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
    222221         !                       ! spe = fwfisf is read from a forcing field 
    223          !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    224          !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     222         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     223         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
    225224         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
    226225         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
    227226         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
    228          !                       ! ad15 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    229          !                       ! hj99 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    230          rn_gammat0  = 1.4e-2    ! gammat coefficient used in blk formula 
    231          rn_gammas0  = 4.e-4    ! gammas coefficient used in blk formula 
     227         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     228         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     229         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     230         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     231         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
    232232         ! 
    233233         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    … …  
    255255         sn_isfpar_zmin = 'isfmlt_par',      -12.      , 'sozisfmin' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    256256         !* 'spe' and 'oasis' case 
    257          sn_isfpar_fwf = 'isfmlt_par' ,      -12.      , 'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     257         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    258258         !* 'bg03' case 
    259          sn_isfpar_Leff = 'isfmlt_par',       0.       , 'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     259         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
    260260      ! 
    261261      ! ---------------- ice sheet coupling ------------------------------- 
    … …  
    300300   ln_tide     = .true.       ! Activate tides 
    301301      ln_tide_pot   = .false.               !  use tidal potential forcing 
    302       clname(1) = 'M2'  !  name of constituent - all tidal components must be set in namelist_cfg 
    303       clname(2) = 'S2' 
    304       clname(3) = 'K1' 
    305       clname(4) = 'O1' 
     302      sn_tide_cnames(1) = 'M2'  !  name of constituent - all tidal components must be set in namelist_cfg 
     303      sn_tide_cnames(2) = 'S2' 
     304      sn_tide_cnames(3) = 'K1' 
     305      sn_tide_cnames(4) = 'O1' 
    306306/ 
    307307!----------------------------------------------------------------------- 
    … …  
    340340   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    341341   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    342    bn_ssh      =    'bdyT_ssh_WED025'    ,         -1        , 'sossheig' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    343    bn_u2d      =    'bdyU_u2d_WED025'    ,         -1        , 'vobtcrtx' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    344    bn_v2d      =    'bdyV_u2d_WED025'    ,         -1        , 'vobtcrty' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    345    bn_u3d      =    'bdyU_u3d_WED025'    ,         -1        , 'vozocrtx' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    346    bn_v3d      =    'bdyV_u3d_WED025'    ,         -1        , 'vomecrty' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    347    bn_tem      =    'bdyT_tra_WED025'    ,         -1        , 'votemper' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    348    bn_sal      =    'bdyT_tra_WED025'    ,         -1        , 'vosaline' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     342   bn_ssh      =    'WED025_bdyT_ssh'    ,         -1        , 'sossheig' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     343   bn_u2d      =    'WED025_bdyU_u2d'    ,         -1        , 'vobtcrtx' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     344   bn_v2d      =    'WED025_bdyV_u2d'    ,         -1        , 'vobtcrty' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     345   bn_u3d      =    'WED025_bdyU_u3d'    ,         -1        , 'vozocrtx' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     346   bn_v3d      =    'WED025_bdyV_u3d'    ,         -1        , 'vomecrty' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     347   bn_tem      =    'WED025_bdyT_tra'    ,         -1        , 'votemper' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     348   bn_sal      =    'WED025_bdyT_tra'    ,         -1        , 'vosaline' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    349349!* for si3 
    350    bn_a_i      =    'bdyT_ice_WED025'    ,         -1        , 'ileadfra' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    351    bn_h_i      =    'bdyT_ice_WED025'    ,         -1        , 'iicethic' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    352    bn_h_s      =    'bdyT_ice_WED025'    ,         -1        , 'isnowthi' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     350   bn_a_i      =    'WED025_bdyT_ice'    ,         -1        , 'ileadfra' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     351   bn_h_i      =    'WED025_bdyT_ice'    ,         -1        , 'iicethic' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
     352   bn_h_s      =    'WED025_bdyT_ice'    ,         -1        , 'isnowthi' ,     .true. , .false., 'yearly'  ,    ''            ,   ''     ,     '' 
    353353/ 
    354354!----------------------------------------------------------------------- 
    355355&nambdy_tide   !  tidal forcing at open boundaries                      (default: OFF) 
    356356!----------------------------------------------------------------------- 
    357    filtide          = 'bdytide_WED025_'         !  file name root of tidal forcing files 
     357   filtide          = 'WED025_bdytide_'         !  file name root of tidal forcing files 
    358358/ 
    359359 
    … …  
    362362!!                                                                    !! 
    363363!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    364 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    365 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     364!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     365!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    366366!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    367367!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    374374/ 
    375375!----------------------------------------------------------------------- 
    376 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
     376&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
    377377!----------------------------------------------------------------------- 
    378378   rn_Cd0      =  2.5e-3    !  drag coefficient [-] 
    379379/ 
    380380!----------------------------------------------------------------------- 
    381 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     381&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    382382!----------------------------------------------------------------------- 
    383383   rn_Cd0      =  2.5e-3    !  drag coefficient [-] 
    … …  
    658658&namctl        !   Control prints                                       (default: OFF) 
    659659!----------------------------------------------------------------------- 
    660    ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
    661      sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    662        sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
    663        sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    664        sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    665        sn_cfctl%l_layout  = .FALSE. ! 
    666        sn_cfctl%l_mppout  = .FALSE. ! 
    667        sn_cfctl%l_mpptop  = .FALSE. ! 
    668        sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    669        sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
    670        sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
    671        sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
    672    nn_print    =    0      !  level of print (0 no extra print) 
    673    nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
    674    nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
    675    nn_jctls    =    0      !  start j indice of control               over a subdomain) 
    676    nn_jctle    =    0      !  end   j indice of control 
    677    nn_isplt    =    1      !  number of processors in i-direction 
    678    nn_jsplt    =    1      !  number of processors in j-direction 
    679    ln_timing   = .true.    !  timing by routine write out in timing.output file 
    680    ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
     660   sn_cfctl%l_runstat = .true.    ! switches and which areas produce reports with the proc integer settings. 
     661   ln_timing   = .true.           !  timing by routine write out in timing.output file 
    681662/ 
    682663!----------------------------------------------------------------------- 

  • NEMO/branches/2020/r12377_ticket2386/cfgs/WED025/EXPREF/namelist_ice_cfg

    r11487 r13540  
    2626&namitd         !   Ice discretization 
    2727!------------------------------------------------------------------------------ 
     28   ln_cat_hfn       = .true.          !  ice categories are defined by a function following rn_himean**(-0.05) 
     29      rn_himean     =   2.0           !  expected domain-average ice thickness (m) 
     30   rn_himin         =   0.01          !  minimum ice thickness (m) used in remapping 
    2831/ 
    2932!------------------------------------------------------------------------------ 
    3033&namdyn         !   Ice dynamics 
    3134!------------------------------------------------------------------------------ 
     35   ln_landfast_L16  = .true.          !  landfast: parameterization from Lemieux 2016 
    3236/ 
    3337!------------------------------------------------------------------------------ 
    … …  
    3842&namdyn_rhg     !   Ice rheology 
    3943!------------------------------------------------------------------------------ 
     44   ln_rhg_EVP       = .true.          !  EVP rheology 
     45      ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    4046/ 
    4147!------------------------------------------------------------------------------ 
    4248&namdyn_adv     !   Ice advection 
    4349!------------------------------------------------------------------------------ 
     50   ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
     51   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
     52      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    4453/ 
    4554!------------------------------------------------------------------------------ 
    4655&namsbc         !   Ice surface boundary conditions 
    4756!------------------------------------------------------------------------------ 
     57   nn_snwfra        =   0             !  calculate the fraction of ice covered by snow (for zdf and albedo) 
     58                                      !     = 0  fraction = 1 (if snow) or 0 (if no snow) 
     59                                      !     = 1  fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 
     60                                      !     = 2  fraction = hsnw / (hsnw+0.02)    [CICE formulation] 
     61   nn_qtrice        =   0             !  Solar flux transmitted thru the surface scattering layer: 
     62                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow) 
     63                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 
    4864/ 
    4965!------------------------------------------------------------------------------ 
    … …  
    6278&namthd_do      !   Ice growth in open water 
    6379!------------------------------------------------------------------------------ 
     80   rn_hinew         =   0.02          !  thickness for new ice formation in open water (m), must be larger than rn_himin 
     81   ln_frazil        = .true.          !  Frazil ice parameterization (ice collection as a function of wind) 
    6482/ 
    6583!------------------------------------------------------------------------------ 
    … …  
    7088&namthd_pnd     !   Melt ponds 
    7189!------------------------------------------------------------------------------ 
     90   ln_pnd            = .false.         !  activate melt ponds or not 
     91      ln_pnd_LEV     = .false.         !  level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012) 
    7292/ 
    7393!------------------------------------------------------------------------------ 
    7494&namini         !   Ice initialization 
    7595!------------------------------------------------------------------------------ 
     96   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
     97   nn_iceini_file   =   1             !     0 = Initialise sea ice based on SSTs 
     98                                      !     1 = Initialise sea ice from single category netcdf file 
     99                                      !     2 = Initialise sea ice from multi category restart file 
     100   ! -- for ln_iceini_file = T 
     101   sn_hti = 'WED025_init_JRA_200001.nc', -12. ,'icethic_cea',  .false.  , .true., 'yearly'  , '' , '', '' 
     102   sn_hts = 'WED025_init_JRA_200001.nc', -12. ,'icesnow_cea',  .false.  , .true., 'yearly'  , '' , '', '' 
     103   sn_ati = 'WED025_init_JRA_200001.nc', -12. ,'ice_cover'  ,  .false.  , .true., 'yearly'  , '' , '', '' 
     104   sn_smi = 'NOT USED'              , -12. ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     105   sn_tmi = 'NOT USED'              , -12. ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     106   sn_tsu = 'NOT USED'              , -12. ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     107   sn_tms = 'NOT USED'              , -12. ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     108   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 
     109   sn_apd = 'NOT USED'              , -12. ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     110   sn_hpd = 'NOT USED'              , -12. ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     111   sn_hld = 'NOT USED'              , -12. ,'hld'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     112   cn_dir='./' 
    76113/ 
    77114!------------------------------------------------------------------------------ 
Note: See TracChangeset for help on using the changeset viewer.