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 5782 for branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2015-10-08T15:51:36+02:00 (9 years ago)
Author:
cetlod
Message:

Final step of improvements/simplifications of ADV & LDF momentum trends. The branch is now phased with the trunk at revision 5721 and is ready to be used

Location:
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r5778 r5782  
    44!----------------------------------------------------------------------- 
    55&namrun        !   parameters of the run 
     6 nn_it000=1 
    67!----------------------------------------------------------------------- 
    78   cn_exp      = "Agulhas" !  experience name  
    8    nn_itend    =   10950   !  last  time step 
     9nn_itend=150 
    910   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    1011   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    11    ln_clobber  = .true.    !  clobber (overwrite) an existing file 
     12ln_clobber=.true. 
    1213/ 
    1314!----------------------------------------------------------------------- 
    … …  
    152153&namtra_ldf    !   lateral diffusion scheme for tracers 
    153154!----------------------------------------------------------------------- 
    154    !                       !  Operator type: 
    155    !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
    156    ln_traldf_lap   =  .true.   !    laplacian operator 
    157    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    158    !                       !  Direction of action: 
    159    ln_traldf_lev   =  .false.  !  iso-level 
    160    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    161    ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
    162    ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    163    ! 
    164    !                       !  iso-neutral options:         
    165    ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
    166    rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    167    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    168    rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    169    ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    170    ! 
    171    !                       !  Coefficients: 
    172    nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    173    !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    174    !                                !   =  0           constant  
    175    !                                !   = 10 F(k)      =ldf_c1d  
    176    !                                !   = 20 F(i,j)    =ldf_c2d  
    177    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    178    !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    179    !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    180155   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    181    rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    182156/ 
    183157!---------------------------------------------------------------------------------- 
    … …  
    207181&namdyn_ldf    !   lateral diffusion on momentum 
    208182!----------------------------------------------------------------------- 
    209    !                       !  Type of the operator :  
    210    ln_dynldf_lap    =  .false.  !  laplacian operator          
    211    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    212    rn_ahm_0_blp     =  8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     183   !                       !  Type of the operator : 
     184   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F 
     185   ln_dynldf_lap =  .false.    !    laplacian operator 
     186   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     187   !                       !  Direction of action  : 
     188   ln_dynldf_lev =  .true.     !  iso-level 
     189   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     190   ln_dynldf_iso =  .false.    !  iso-neutral 
     191   !                       !  Coefficient 
     192   nn_ahm_ijk_t  =   0         !  space/time variation of eddy coef 
     193   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     194   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     195   !                                !  =  0  constant 
     196   !                                !  = 10  F(k)=c1d 
     197   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     198   !                                !  = 30  F(i,j,k)=c2d*c1d 
     199   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     200   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     201   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     202   rn_bhm_0      = 8.5e+11      !  horizontal bilaplacian eddy viscosity [m4/s] 
    213203/ 
    214204!----------------------------------------------------------------------- 
    … …  
    243233!----------------------------------------------------------------------- 
    244234&namctl        !   Control prints & Benchmark 
    245 !----------------------------------------------------------------------- 
    246235/ 
    247236!----------------------------------------------------------------------- 

  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r5778 r5782  
    66&namrun        !   parameters of the run 
    77!----------------------------------------------------------------------- 
    8    nn_no       =       0   !  job number (no more used...) 
    9    cn_exp      =  "ORCA2"  !  experience name 
    10    nn_it000    =       1   !  first time step 
    11    nn_itend    =    1050   !  last  time step (std 5475) 
    12    nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    13    nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    14    ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    15    nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
    16    nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    17                            !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    18                            !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    19                            !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    20    cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
    21    cn_ocerst_indir = "."       !  directory from which to read input ocean restarts 
    22    cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    23    cn_ocerst_outdir = "."      !  directory in which to write output ocean restarts 
    24    nn_istate   =       0   !  output the initial state (1) or not (0) 
    25    ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    26    nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
    27    nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 
    28    nn_write    =    105   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    29    ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    30    ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    31    ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard 
    32    ln_clobber  = .true.   !  clobber (overwrite) an existing file 
    33    nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    348/ 
    359!----------------------------------------------------------------------- 
    … …  
    137111&namtra_adv    !   advection scheme for tracer 
    138112!----------------------------------------------------------------------- 
    139    ln_traadv_cen =  .false.   !  2nd order centered scheme 
    140       nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
    141       nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    142113   ln_traadv_fct =  .true.    !  FCT scheme 
    143114      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
    144115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    145       nn_fct_zts =  0               !  >=1, 2nd order FCT scheme with vertical sub-timestepping 
     116      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
    146117      !                             !        (number of sub-timestep = nn_fct_zts) 
    147    ln_traadv_mus  =  .false.  !  MUSCL scheme 
    148       ln_mus_ups  =  .false.         !  use upstream scheme near river mouths 
    149    ln_traadv_ubs  =  .false.  !  UBS scheme 
    150       nn_ubs_v    =  2               !  =2  , vertical 2nd order FCT 
    151    ln_traadv_qck  =  .false.  !  QUICKEST scheme 
    152 / 
    153 !----------------------------------------------------------------------- 
    154 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    155 !----------------------------------------------------------------------- 
    156118/ 
    157119!----------------------------------------------------------------------- 
    … …  
    168130   ln_traldf_lev   =  .false.  !  iso-level 
    169131   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    170    ln_traldf_iso   =  .true.   !  iso-neutral 
    171    ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     132   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     133   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    172134   ! 
    173135   !                       !  iso-neutral options:         
    174136   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
    175137   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    176    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     138   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    177139   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    178    ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     140   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    179141   ! 
    180142   !                       !  Coefficients: 
    181    nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
     143   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
    182144   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    183    !                                !   =  0           constant  
    184    !                                !   = 10 F(k)      =ldf_c1d  
    185    !                                !   = 20 F(i,j)    =ldf_c2d  
     145   !                                !   =  0           constant 
     146   !                                !   = 10 F(k)      =ldf_c1d 
     147   !                                !   = 20 F(i,j)    =ldf_c2d 
    186148   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    187149   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    … …  
    198160   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    199161   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    200    !                                !   =  0           constant  
    201    !                                !   = 10 F(k)      =ldf_c1d  
    202    !                                !   = 20 F(i,j)    =ldf_c2d  
     162   !                                !   =  0           constant 
     163   !                                !   = 10 F(k)      =ldf_c1d 
     164   !                                !   = 20 F(i,j)    =ldf_c2d 
    203165   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    204166   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    205 /!----------------------------------------------------------------------- 
     167/ 
     168!----------------------------------------------------------------------- 
    206169&namtra_dmp    !   tracer: T & S newtonian damping 
    207170!----------------------------------------------------------------------- 
    … …  
    218181   ln_dynvor_mix = .false. !  mixed scheme 
    219182   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    220       nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    221    ln_dynvor_msk = .true. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
     183      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     184   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    222185/ 
    223186!----------------------------------------------------------------------- 
    … …  
    289252!----------------------------------------------------------------------- 
    290253/ 
    291 !----------------------------------------------------------------------- 
    292 &namobs       !  observation usage                                      ('key_diaobs') 
    293 !----------------------------------------------------------------------- 
    294 / 
    295 !----------------------------------------------------------------------- 
    296 &nam_asminc   !   assimilation increments                               ('key_asminc') 
    297 !----------------------------------------------------------------------- 
    298 / 

  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r5778 r5782  
    208208!----------------------------------------------------------------------- 
    209209   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    210    ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
     210   ln_hpg_zps  = .false.    !  z-coordinate - partial steps (interpolation) 
    211211   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    212212   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    213    ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
     213   ln_hpg_prj  = .true.   !  s-coordinate (Pressure Jacobian scheme) 
    214214   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    215215                                 !           centered      time scheme  (F) 

  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r5778 r5782  
    106106   ! 
    107107   !                       !  Coefficients: 
    108    nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     108   nn_aht_ijk_t    = 0        !  space/time variation of eddy coef 
    109109   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    110110   !                                !   =  0           constant  
    … …  
    120120&namtra_ldfeiv !   eddy induced velocity param. 
    121121!---------------------------------------------------------------------------------- 
    122    ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
     122   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     123   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     124   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     125   nn_aei_ijk_t  = 0      ! space/time variation of the eiv coeficient 
     126   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     127   !                                !   =  0           constant 
     128   !                                !   = 10 F(k)      =ldf_c1d 
     129   !                                !   = 20 F(i,j)    =ldf_c2d 
     130   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     131   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    123132/ 
    124133!----------------------------------------------------------------------- 
    … …  
    141150   sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    142151   sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    143    sn_ahu  = 'dyna_grid_U' ,    120            , 'vozoahtu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    144    sn_ahv  = 'dyna_grid_V' ,    120            , 'vomeahtv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    145    sn_ahw  = 'dyna_grid_W' ,    120            , 'voveahtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    146    sn_eiu  = 'dyna_grid_U' ,    120            , 'vozoaeiu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    147    sn_eiv  = 'dyna_grid_V' ,    120            , 'vomeaeiv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    148    sn_eiw  = 'dyna_grid_W' ,    120            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    149152! 
    150153   cn_dir      = './'      !  root directory for the location of the dynamical files 
    151    ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad") 
    152154   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file 
    153155   ln_dynbbl   =  .true.   !  bbl coef are in files, so read them - requires ("key_trabbl") 
Note: See TracChangeset for help on using the changeset viewer.