Table des matières
Little_NEMO Peru12 - Etude de sensibilité
Questions en cours... Description des runs
Où
Christophe sur Curie (06 09 2012) - run "erai2..." : Remake du run ERAI2 de Stephane (AVEC understress) : => bug ppacr
Erreur ppacr=7 au lieu de 9 dans la définition des 46 niveaux verticaux (06 09 2012) - run "erai3..." : Remake du run ERAI3 de Stephane (SANS understress) : => run contrôle passage nemo_indian -> little_nemo + tests multi machines (10 09 2012) - run "erai4" = bilplacien au lieu de laplacien sur les traceurs (10 09 2012) - run "erai5" = UBS sur les traceurs (12 09 2012) - run "erai6" = UBS sur traceurs et Momentum
Plantage avec UBS sur le Momentum Diagnostiques
Sn_cnf dans namelist au lieu de sn_cnf: majuscule importe? si oui, quelle conséquence?
ln_tmx_itf = .true. ! ITF specific parameterisation
=> mettre a false (mer indonesiennes) Conséquence? Verifier que bug fixe TKE OK dans little_NEMO!! Si oui, quel impact sur la verticale?
Comparaison des sections aux observations
Utiliser nemo.exe compilé avec la clé key_ldfslp alors qur pas iso-neutral: Source de bug possible ;-( ? Pas de raison,
mais compte tenu du passé... comparaison roms12 => refaire les pdf à échelle fixe! problème des fenêtres java qui disparaissent en arrière plan. Différence de section pour un même run entre saison à partir de 3 x 1 mois et saison à partir de 18 fois 5d avec ou sans restart.obc?
Expériences suivi:
run
base
Momentum
Tracer
ppacr
periode
understress
restart
verion
machine
qui?
advection
diffusion
advection
diffusion
ERAI2 Vector form bilaplacian bilaplacian -1e10 TVD scheme laplacian iso-neutral 300 7.000_wp 1990-1994 avec levitus nemo_indian vargas Stéphane
erai21c ERAI2
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
1990-1994
avec
levitus
little_nemo r48
curie fat nodes
Christophe
Stéphane a laissé la valeur de ppacr = 7.0000000000000_wp des 75 niveaux dans la définition des 46 niveaux des simus ERAI[1,2,3]
Du coup les niveaux verticaux des OBCs et de levitus issus des simus à 46 niveaux (ppacr=9) ne sont pas compatibles avec les 46 niveaux (ppacr=7) de nos runs.
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
run
base
Momentum
Tracer
ppacr
periode
understress
restart
verion
machine
qui?
advection
diffusion
advection
diffusion
ERAI3 ERAI2 Vector form bilaplacian bilaplacian -1e10 TVD scheme laplacian iso-neutral 300 7.000_wp 1994 sans 1994 de ERAI2 nemo_indian vargas Stéphane
ERAI3B (ERAI4)ERAI3
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
1990-1994
sans
levitus
nemo_indian
vargas
Stéphane
erai31v à faire ERAI3B
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
1990-1994
sans
levitus
little_nemo r48
vargas
Christophe
erai31j à faire ERAI3B
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
1990-1994
sans
levitus
little_nemo r48
jade
Christophe
erai31c ERAI3B
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
1990-1994
sans
levitus
little_nemo r48
curie
Christophe
run
base
Momentum
Tracer
ppacr
periode
understress
restart
verion
machine
qui?
advection
diffusion
advection
diffusion
erai4l
erai21c
Vector form
bilaplacian
horizontal
-2.5e9
TVD scheme
bilaplacian
horizontal
-1e9
9.000_wp
1989 -1995
avec
levitus
little_nemo r48
curie fat nodes
Christophe
erai40c erai31c
Vector form
bilaplacian
horizontal
-2.5e9
TVD scheme
bilaplacian
horizontal
-1e9
9.000_wp
1990-1994
sans levitus
little_nemo r48
curie fat nodes
Christophe
diff erai4l / erai21c:
--------------------------------------------------------------------------------
diff -b -w erai21c/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 erai21c/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
run
base
Momentum
Tracer
ppacr
periode
understress
restart
verion
machine
qui?
advection
diffusion
advection
diffusion
erai5f
erai21c
Vector form
bilaplacian
horizontal
-1e10
UBS
none
none
none
9.000_wp
1989 -1995
avec
levitus
little_nemo r48
curie fat nodes
Christophe
erai50c erai31c
Vector form
bilaplacian
horizontal
-1e10
UBS
none
none
none
9.000_wp
1990-1994
sans levitus
little_nemo r48
curie fat nodes
Christophe
diff erai5f / erai21c:
--------------------------------------------------------------------------------
diff -b -w erai21c/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 erai21c/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
run
base
Momentum
Tracer
ppacr
periode
understress
restart
verion
machine
qui?
advection
diffusion
advection
diffusion
erai6a erai21c
UBS
bilaplacian
horizontal
-1e10
UBS
none
none
none
9.000_wp
1989- plante à 1993-01-01 + 106 jours
avec
levitus
little_nemo r48
curie fat nodes
Christophe
erai6a_2 erai6a
UBS
none
none
none
UBS
none
none
none
9.000_wp
1993-plante à 1993-01-01 + 98 jours
avec
1993-01-01
little_nemo r48
curie fat nodes
Christophe
erai60c erai31c
UBS
none
none
none
UBS
none
none
none
9.000_wp
1990-plante à 1990-01-01 + 70 jours
sans levitus
little_nemo r48
curie fat nodes
Christophe
erai6a tourne avec les mêmes forçages que erai21c. Ils ne sont donc pas en cause dans le plantage.
diff erai6a / erai21c:
--------------------------------------------------------------------------------
diff -b -w erai21c/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 / erai21c:
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
diff erai60c / erai31c:
--------------------------------------------------------------------------------
diff -b -w erai31c/param_inputs/namelist.base.oce erai60c/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
636c636
< ln_dynldf_bilap = .true. ! bilaplacian operator
---
> ln_dynldf_bilap = .false. ! bilaplacian 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 erai31c/compile/nemo/cpp_PERU12.fcm erai60c/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_iomput key_ldfslp
---
> 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
===>>> : E R R O R (output.abort)
===========
stpctl: the zonal velocity is larger than 20 m/s
======
kt= 8532 max abs(U): 22.51 , i j k: 18 377 9
1990
1991
1992
1993
1994
ERAI2 (ppacr=7.000) (bilapdyn + laptrac + avec Understress) + nemo_indian + vargas
erai21c (bilapdyn + laptrac + avec Understress) + little_nemo + curie
1990
1991
1992
1993
1994
ERAI3 (ppacr=7.000) (bilapdyn + laptrac + sans Understress) + nemo_indian + vargas
ERAI3B (bilapdyn + laptrac + sans Understress) + nemo_indian + vargas
erai31c (bilapdyn + laptrac + sans Understress) + little_nemo + curie
1990
1991
1992
1993
1994
erai4l (bilapdyn + bilaptrac + avec Understress)
erai40c (bilapdyn + bilaptrac + sans Understress)
1990
1991
1992
1993
1994
erai5f (bilapdyn + UBStrac + avec Understress)
erai50c (bilapdyn + UBStrac + sans Understress)
