Little_NEMO Peru12 - Etude de sensibilité

dans  par_depth.h90:
- pour 46 niveaux
     &     ppacr =    9.0000000000000_wp ,  &  !: (non dimensional): stretching factor
- pour 75 niveaux
     &     ppacr =    7.0000000000000_wp ,  &  !: (non dimensional): stretching factor 

diff erai4l / erai20c:

-------------------------------------------------------------------------------- diff -b -w erai20c/param_inputs/namelist.base.oce erai4l/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 552,553c552,553 < ln_traldf_lap = .true. ! laplacian operator < ln_traldf_bilap = .false. ! bilaplacian operator --- > ln_traldf_lap = .false. ! laplacian operator > ln_traldf_bilap = .true. ! bilaplacian operator 556,557c556,557 < ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) < ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") --- > ln_traldf_hor = .true. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) > ln_traldf_iso = .false. ! iso-neutral (needs "key_ldfslp") 568c568 < rn_aht_0 = 300. ! horizontal eddy diffusivity for tracers [m2/s] --- > rn_aht_0 = -1.e9 ! horizontal eddy diffusivity for tracers [m2/s] 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -2.5e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai20c/compile/nemo/cpp_PERU12.fcm erai4l/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_ldfslp key_iomput --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

diff erai5f / erai20c:

-------------------------------------------------------------------------------- diff -b -w erai20c/param_inputs/namelist.base.oce erai5f/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s] -------------------------------------------------------------------------------- diff -b -w erai20c/compile/nemo/cpp_PERU12.fcm erai5f/compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_ldfslp key_iomput --- > bld::tool::fppkeys key_peru12 key_depth=46 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput

diff erai6a / erai20c:

-------------------------------------------------------------------------------- diff -b -w erai20c/param_inputs/namelist.base.oce erai6a/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 542c542 < ln_traadv_tvd = .true. ! TVD scheme --- > ln_traadv_tvd = .false. ! TVD scheme 545c545 < ln_traadv_ubs = .false. ! UBS scheme --- > ln_traadv_ubs = .true. ! UBS scheme 552c552 < ln_traldf_lap = .true. ! laplacian operator --- > ln_traldf_lap = .false. ! laplacian operator 601c601 < ln_dynadv_vec = .true. ! vector form (T) or flux form (F) --- > ln_dynadv_vec = .false. ! vector form (T) or flux form (F) 603c603 < ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme --- > ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme 644c644 < rn_ahm_0_blp = -1.e10 ! horizontal bilaplacian eddy viscosity [m4/s] --- > rn_ahm_0_blp = -1.e9 ! horizontal bilaplacian eddy viscosity [m4/s]

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt=169220 max abs(U): 21.95 , i j k: 13 376 1

diff erai6a / erai20c:

diff 19930101_19931231_1_PLANTE/namelist 19930101_19931231_3_PLANTE_sans_bilap_dyn/namelist 636c636 < ln_dynldf_bilap = .true. ! bilaplacian operator --- > ln_dynldf_bilap = .false. ! bilaplacian operator

===>>> : E R R O R (output.abort) =========== stpctl: the zonal velocity is larger than 20 m/s ====== kt=168314 max abs(U): 197.7 , i j k: 9 376 1