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 5329 for trunk/NEMOGCM/CONFIG – NEMO

Changeset 5329 for trunk/NEMOGCM/CONFIG


Ignore:
Timestamp:
2015-06-01T15:11:25+02:00 (9 years ago)
Author:
pabouttier
Message:

Add stochastic parametrization of EOS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5321 r5329  
    1010!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    1111!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
    12 !!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb) 
     12!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    1313!!             10 - miscellaneous    (namsol, nammpp, namctl) 
    1414!!             11 - Obs & Assim      (namobs, nam_asminc) 
     
    5050!!                      ***  Domain namelists  *** 
    5151!!====================================================================== 
    52 !!   namcfg       parameters of the configuration       
     52!!   namcfg       parameters of the configuration 
    5353!!   namzgr       vertical coordinate 
    5454!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     
    5858! 
    5959!----------------------------------------------------------------------- 
    60 &namcfg     !   parameters of the configuration       
     60&namcfg     !   parameters of the configuration 
    6161!----------------------------------------------------------------------- 
    6262   cp_cfg      =  "default"            !  name of the configuration 
     
    7272   jperio      =       0               !  lateral cond. type (between 0 and 6) 
    7373                                       !  = 0 closed                 ;   = 1 cyclic East-West 
    74                                        !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot  
     74                                       !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    7575                                       !  = 4 cyclic East-West AND North fold T-point pivot 
    7676                                       !  = 5 North fold F-point pivot 
    7777                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    78    ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present  
     78   ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present 
    7979                                       !  in netcdf input files, as the start j-row for reading 
    8080/ 
     
    101101                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    102102   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    103    rn_bb       =    0.8    !  stretching with SH94 s-sigma    
     103   rn_bb       =    0.8    !  stretching with SH94 s-sigma 
    104104                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
    105105   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     
    110110   rn_zb_b     =   -0.2    !  offset for calculating Zb 
    111111                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
    112    rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
     112   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1) 
    113113/ 
    114114!----------------------------------------------------------------------- 
     
    164164   nn_baro       =    30               !  Number of iterations of barotropic mode 
    165165                                       !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
    166    rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T  
     166   rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T 
    167167   nn_bt_flt     =    1                !  Time filter choice 
    168168                                       !  = 0 None 
    169169                                       !  = 1 Boxcar over   nn_baro barotropic steps 
    170                                        !  = 2 Boxcar over 2*nn_baro     "        "   
     170                                       !  = 2 Boxcar over 2*nn_baro     "        " 
    171171/ 
    172172!----------------------------------------------------------------------- 
     
    245245   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
    246246   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
    247                            !  0 =no isf                  1 = presence of ISF  
    248                            !  2 = bg03 parametrisation   3 = rnf file for isf    
     247                           !  0 =no isf                  1 = presence of ISF 
     248                           !  2 = bg03 parametrisation   3 = rnf file for isf 
    249249                           !  4 = ISF fwf specified 
    250250                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
     
    277277&namsbc_flx    !   surface boundary condition : flux formulation 
    278278!----------------------------------------------------------------------- 
    279 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !  
     279!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    280280!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    281281   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     
    320320   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    321321   rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
    322    rn_zu       = 10.        !  Wind vector reference height (m)                  
     322   rn_zu       = 10.        !  Wind vector reference height (m) 
    323323   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    324324   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    325    rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     325   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity 
    326326                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    327327/ 
     
    373373!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    374374!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    375    sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''       
     375   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    376376   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    377377   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     
    422422/ 
    423423!----------------------------------------------------------------------- 
    424 &namsbc_isf    !  Top boundary layer (ISF)  
     424&namsbc_isf    !  Top boundary layer (ISF) 
    425425!----------------------------------------------------------------------- 
    426426!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     
    499499                                                      ! Initial mass required for an iceberg of each class 
    500500      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    501                                                       ! Proportion of calving mass to apportion to each class   
     501                                                      ! Proportion of calving mass to apportion to each class 
    502502      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    503503                                                      ! Ratio between effective and real iceberg mass (non-dim) 
    504                                                       ! i.e. number of icebergs represented at a point          
     504                                                      ! i.e. number of icebergs represented at a point 
    505505      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    506506                                                      ! thickness of newly calved bergs (m) 
     
    511511      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    512512      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    513       ln_passive_mode          = .false.              ! iceberg - ocean decoupling    
     513      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
    514514      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
    515515                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    516516      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    517       rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
     517      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    518518 
    519519!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    520520!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    521521      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , '' 
    522     
    523       cn_dir = './'  
     522 
     523      cn_dir = './' 
    524524/ 
    525525 
     
    597597                                          !  = 2, use tidal harmonic forcing data from files 
    598598                                          !  = 3, use external data AND tidal harmonic forcing 
    599     cn_dyn3d      =  'none'               !   
     599    cn_dyn3d      =  'none'               ! 
    600600    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
    601601                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    602     cn_tra        =  'none'               !  
     602    cn_tra        =  'none'               ! 
    603603    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    604604                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    605     cn_ice_lim      =  'none'             !   
     605    cn_ice_lim      =  'none'             ! 
    606606    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
    607607                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     
    612612    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    613613    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
    614     rn_time_dmp   =  1.                   ! Damping time scale in days  
     614    rn_time_dmp   =  1.                   ! Damping time scale in days 
    615615    rn_time_dmp_out =  1.                 ! Outflow damping time scale 
    616616    nn_rimwidth   = 10                    !  width of the relaxation zone 
     
    665665   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    666666   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    667    rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
     667   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T 
    668668   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    669669   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     
    711711!----------------------------------------------------------------------- 
    712712   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency 
    713                                  !  =-1, TEOS-10  
    714                                  !  = 0, EOS-80  
     713                                 !  =-1, TEOS-10 
     714                                 !  = 0, EOS-80 
    715715                                 !  = 1, S-EOS   (simplified eos) 
    716716   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     
    811811&nam_vvl    !   vertical coordinate options 
    812812!----------------------------------------------------------------------- 
    813    ln_vvl_zstar  = .true.           !  zstar vertical coordinate                    
     813   ln_vvl_zstar  = .true.           !  zstar vertical coordinate 
    814814   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
    815815   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
     
    996996!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
    997997!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
     998!!   namsto            Stochastic parametrization of EOS 
    998999!!====================================================================== 
    9991000! 
     
    10541055   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    10551056/ 
     1057!----------------------------------------------------------------------- 
     1058&namsto       ! Stochastic parametrization of EOS 
     1059!----------------------------------------------------------------------- 
     1060   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T) 
     1061   ln_rstseed = .true.           ! read seed of RNG from restart file 
     1062   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input) 
     1063   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output) 
     1064 
     1065   ln_sto_eos = .false.          ! stochastic equation of state 
     1066   nn_sto_eos = 1                ! number of independent random walks 
     1067   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points) 
     1068   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points) 
     1069   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps) 
     1070   nn_eos_ord  = 1               ! order of autoregressive processes 
     1071   nn_eos_flt  = 0               ! passes of Laplacian filter 
     1072   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0) 
     1073/ 
    10561074 
    10571075!!====================================================================== 
     
    11651183   ln_sst     = .false.     ! Logical switch for SST observations 
    11661184   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    1167    ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
     1185   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations 
    11681186 
    11691187   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
     
    11921210   sstfbfiles = 'sst_01.nc' 
    11931211                           !     seaicefiles             Sea Ice input observation file names 
    1194    seaicefiles = 'seaice_01.nc'   
     1212   seaicefiles = 'seaice_01.nc' 
    11951213                           !     velavcurfiles           Vel. cur. daily av. input file name 
    11961214                           !     velhvcurfiles           Vel. cur. high freq. input file name 
Note: See TracChangeset for help on using the changeset viewer.