Table des matières
Little_NEMO Peru12 - Etude de sensibilité - 361 x 505 x 46
Questions en cours... Prochains Runs Inputs
Forçages - blk - ERAI OBC
46 niveaux (run ORCA025-B83) 75 niveaux (run ORCA025-B83) interpolé depuis 46 niveaux avec matlab 75 Niveaux (run ORCA025.L75-MJM95-S) Nemo, XIOS et Scripts de lancement
Scripts: changer de pas de temps en cours d'expérience 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 Détail des runs
run erai64c: 1er run 75 couches (bilap/bilap) basé sur erai41c au lieu de (UBS/UBS) run erai65c: 1er run 75 couches (UBS/UBS) basé sur erai61c (UBS/UBS) run erai66c: run 75 niveaux (UBS/UBS) - remake erai65c sans bug x2 sur UBS traceurs Diagnostiques
Sections
erai2: bilapdyn + laptrac + avec Understress erai3: bilapdyn + laptrac + sans Understress erai4: bilapdyn + bilaptrac erai5 : bilapdyn + UBStrac erai6 : UBSdyn + UBStrac grad SST
erai2: bilapdyn + laptrac + avec Understress erai3: bilapdyn + laptrac + sans Understress erai4: bilapdyn + bilaptrac erai5 : bilapdyn + UBStrac erai6 : UBSdyn + UBStrac ROMS 12°: UBSdyn + UBStrac PDF
erai2: bilapdyn + laptrac + avec Understress erai3: bilapdyn + laptrac + sans Understress erai4: bilapdyn + bilaptrac erai5 : bilapdyn + UBStrac erai6 : UBSdyn + UBStrac ROMS 12°: UBSdyn + UBStrac Reacpitulatif 1990-1994 : erai40c / erai50c / erai61c FAQ - Errors
Pb create new config
Light penetration : verifier rbg versus 2 Bandes
seb solver nn_solv = 1 au lieu de 2=> rn_sor n'intervient pas.
mais nn_solv = 2 plus efficace si on règle rn_sor qui dépend de la config.
Régler par dichotomie (tester 1.8 , 1.9 puis 2)
regarder dans solver.stat et tracer la courbe niteration = f(pdt)
nombre de procs
bottom friction même que trop12?
remettre le restart obc. Ne coûte rien. et démarrage plus efficace (solver moins d'itérations au début)
Pb obc conservation transport si on conserve la bathy mais qu'on interpole les champs... programme idl Seb+Gurvan basé sur les lignes de courant.
fréquence lecture obc + blk
patcher dans little_nemo le bug 2x diffusion sur UBS
mercator extraction gilles et romain (run libre au 1/4 de degré)
On vire le rappel sss
Sn_cnf dans namelist au lieu de sn_cnf: majuscule importe? si oui, quelle conséquence? (à priori fortran => sans importance)
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 figures (demander run à Vincent) avec ou sans restart.obc? dessin matlab des zones (100:300,50:200) ou (50:400,70:end)
commiter sur little_NEMO traadv_ubs.F90
Vincent: pourquoi etat zero levitus au lieu de Drakar?!?
plot temporel du fond pour voir si dérive...
2eme run 75 couches
question du choix de la bathy peru12 par rapport à celle de trop12
Extraction OBC depuis le 75 couches de Molines.
porter le orca75 et les scripts de steph sur curie (pour idl)
Stéphane: quelles modifs pour utiliser tes scripts OBC sur ORCA025.L75-MJM95-S au lieu de ORCA025-B83
3eme run 75 couches
Le même que le 2eme avec:
restart OBC TRUE (discussion Rachid non, Stephane oui, Seb je ne sais pas donc je mets)
-fp model precise
-i4 -r8 -O3 -xAVX -fp-model precise
-fp-model precise (pour reproductibilité avec ifort; a mettre systematiquement)
-xAVX (pour thin nodes)
cf chez Seb: /ccc/work/cont005/ra0542/massons/now/models/little_nemo/ARCH/CNRS
update little_nemo => modif OBC pour remettre div=0
Avec 75 niveaux sur 1992-2000 (autant démarrer en 1990)
Paramétrisation Drakkar 12° - Molines
Kepsilon au lieu de TKE
???
<!--field field_ref="20d" /-->
Disponibles sur reee084@gaya:/ homegaya/rech/eee/reee084/peru12_nemo/inputs/blk/ERAI/MASK (1990-2011)
Reprise des scripts de Stéphane (1990-1995) pour prolonger (1995-2011) : scripts sur reee084@ulam:/workdir/rech/eee/reee084/wd/pre/peru12_nemo/blk/ERAI (explications détaillées dans "README_CHRIS")
extract_ERAI.ksh : extrait les données era-i sur la region Peru.
depuis compte Sophie Cloché : /u/rech/psl/rpsl376/ERAI/NETCDF/GLOBAL_075/4xdaily
IDL: era_inter_6h.pro:
mask
interpolation temporelle (analyses)
derivée des cumuls (forecast)
calcul de wspd10
calcul de q2 = f(d2,t2,msl)
calcul de precip = f(cp,lsp)
nocs_weights_v1.0/data: Calcul du fichier poids pour l'atmosphere (à faire une seule fois par grille)
Disponibles reee084@gaya:/u/rech/eee/reee084/peru12_nemo/inputs/obc/ORCA025-B83
Rreprise des scripts de Stéphane (1990-1995) pour prolonger (1995-2007) :
scripts sur reee084@ulam:/workdir/rech/eee/reee084/wd/pre/peru12_nemo/obc/ORCA025-B83
(explications détaillées dans "p0.sh")
regional_job_peru.ksh qui à partir d'idl et nco va extraire les obc dans le run global ORCA025-B83 (p1.sh) puis les reconstruit dans un deuxième temps en fichiers annuels (p2.sh)
Run ORCA025.L75-MJM95-S
https://forge.ipsl.jussieu.fr/mentor
svn checkout svn+ssh:// cholod @forge.ipsl.jussieu.fr/ipsl/forge/projets/mentor/svn/trunk/little_NEMO/exp/exp02 erai64c
svn checkout svn+ssh://cholod @forge.ipsl.jussieu.fr/ipsl/forge/projets/little_nemo/svn/trunk little_NEMO
svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/trunk XIOS etcdf4_internal 2 >&1 | tee listing_XIOS_prod.txt
Il suffit de mettre date_debut_exp=date_debut_job => nit_0000=1 et tout le reste OK...
On peut laisser time_origine à la même date. (le nombre de pdt n'est plus cohérent avec la date du début d'expérience réel, mais l'indication de la date de début reste valable.
ln_rst=0 (puis remettre normal des que le job est passé pour que le suivant reprenne le processus normal)
LEVITUS=0 poour repartir des restarts même si date_deb_job = date_deb_exp (puis remettre à 1 pour si on recupere les scripts pour une autre expérience, on parte bien de levitus et non du restart indiqué dans getrst_oce.sh)
Expériences suivi:
/ccc/work/cont005/ra0542/hourdinc/wd/peru12_nemo/exp/erai20c, erai4l...
Inputs:
/ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/inputs/Stephane
Execution en cours:
/ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/execute/erai20c, erai4l...
Outputs:
/ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/outputs/erai20c, erai4l...
run
base
Momentum
Tracer
ppacr
pdt
periode
understress
restart
verion
machine
qui?
cpu
advection
diffusion
advection
diffusion
proc
time
ERAI2 Vector form bilaplacian bilaplacian -1e10 TVD scheme laplacian iso-neutral 300 7.000_wp 800 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
800
1990-1994
avec
levitus
little_nemo r48
curie fat nodes
Christophe
64 nemo +1 xios
~6h/an
moyenne < 6 heures / an sur curie Fat nodes 64 nemo.exe + 1 xios_server.exe
19900101_19901231/erai21c_19900101_19901231.jobid_795489.txt:795489 erai21c_1990 96 3 22/10/2012 18:47:41 5:54:10 RUNNING None
19910101_19911231/erai21c_19910101_19911231.jobid_795957.txt:795957 erai21c_1991 96 3 23/10/2012 00:42:37 5:54:26 RUNNING None
19920101_19921231/erai21c_19920101_19921231.jobid_796333.txt:796333 erai21c_1992 96 3 23/10/2012 06:37:24 5:52:27 RUNNING None
19930101_19931231/erai21c_19930101_19931231.jobid_797257.txt:797257 erai21c_1993 96 3 23/10/2012 15:34:22 5:57:30 RUNNING None
19940101_19941231/erai21c_19940101_19941231.jobid_798724.txt:798724 erai21c_1994 96 3 24/10/2012 09:10:38 5:54:00 RUNNING None
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
pdt
periode
understress
restart
verion
machine
qui?
cpu
advection
diffusion
advection
diffusion
proc
time
ERAI3 ERAI2 Vector form bilaplacian bilaplacian -1e10 TVD scheme laplacian iso-neutral 300 7.000_wp 800 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
800
1990-1994
sans
levitus
nemo_indian
vargas
Stéphane
erai31v en cours... ERAI3B
Vector form
bilaplacian
bilaplacian
-1e10
TVD scheme
laplacian
iso-neutral
300
9.000_wp
800
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
800
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
800
1990-1994
sans
levitus
little_nemo r48
curie fat nodes
Christophe
64 nemo +1 xios
~6h/an
moyenne < 6 heures / an sur curie Fat nodes 64 nemo.exe + 1 xios_server.exe
19900101_19901231/erai31c_19900101_19901231.jobid_798167.txt:798167 erai31c_1990 96 3 23/10/2012 17:59:15 5:59:17 RUNNING None
19910101_19911231/erai31c_19910101_19911231.jobid_798887.txt:798887 erai31c_1991 96 3 24/10/2012 09:10:38 5:57:42 RUNNING None
19920101_19921231/erai31c_19920101_19921231.jobid_800364.txt:800364 erai31c_1992 96 3 24/10/2012 15:13:28 5:57:21 RUNNING None
19930101_19931231/erai31c_19930101_19931231.jobid_801282.txt:801282 erai31c_1993 96 3 24/10/2012 21:14:34 5:52:19 RUNNING None
19940101_19941231/erai31c_19940101_19941231.jobid_801761.txt:801761 erai31c_1994 96 3 25/10/2012 03:10:02 5:52:29 RUNNING None
run
base
Momentum
Tracer
ppacr
pdt
periode
understress
restart
verion
machine
qui?
cpu
advection
diffusion
advection
diffusion
proc
time
erai4l
erai21c
Vector form
bilaplacian
horizontal
-2.5e9
TVD scheme
bilaplacian
horizontal
-1e9
9.000_wp
800
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
800
1990-1996 - plante le 27 Fev 1997
sans levitus
little_nemo r48
curie fat nodes
Christophe
erai41c en cours... erai40c
Vector form
bilaplacian
horizontal
-2.5e9
TVD scheme
bilaplacian
horizontal
-1e9
9.000_wp
400
1990-1994
sans levitus
little_nemo r62
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
pdt
periode
understress
restart
verion
machine
qui?
cpu
advection
diffusion
advection
diffusion
proc
time
erai5f
erai21c
Vector form
bilaplacian
horizontal
-1e9 UBS
none
none
none
9.000_wp
800
1989 -1995
avec
levitus
little_nemo r48
curie fat nodes
erai50c erai31c
Vector form
bilaplacian
horizontal
-1e9 UBS
none
none
none
9.000_wp
800
1990-1994
sans levitus
little_nemo r48
curie fat nodes
Christophe
erai51c erai31c
Vector form
bilaplacian
horizontal
-2.5e9 UBS
none
none
none
9.000_wp
800
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
diff erai50c / erai31c:
diff erai51c / erai50c:
run
base
Momentum
Tracer
ppacr
pdt
periode
understress
restart
verion
machine
qui?
cpu
advection
diffusion
advection
diffusion
proc
time
erai6a erai21c
UBS
bilaplacian
horizontal
-1e10
UBS
none
none
none
9.000_wp
800
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
800
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
800
1990-plante à 1990-01-01 + 70 jours
sans levitus
little_nemo r48
curie fat nodes
Christophe
erai61c erai60c
UBS
none
none
none
UBS
none
none
none
9.000_wp
400
1990-1994
sans levitus
little_nemo r48
curie thin nodes
Christophe
64 nemo +1 xios
~6h/an
erai62c erai61c
UBS
none
none
none
UBS
none
none
none
9.000_wp
720
1990-plante à 1990-01-01 + 83 jours
sans levitus
little_nemo r48
curie thin nodes
Christophe
64 nemo +1 xios
erai63c erai62c
UBS
none
none
none
UBS
none
none
none
9.000_wp
600
1990-plante à 1991-01-01 + 112 jours
sans levitus
little_nemo r48
curie thin nodes
Christophe
64 nemo +1 xios
erai67c erai61c
UBS
none
none
none
UBS
none
none
none
9.000_wp
400
1990-plante le 27 Mai 1997 (el niño)
sans levitus
little_nemo r62
curie thin nodes
Christophe
64 nemo +1 xios
moyenne < 6 heures / an sur curie Thin nodes 64 nemo.exe + 1 xios_server.exe + pdt=600s
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
diff erai61c / erai31c:
--------------------------------------------------------------------------------
diff -b -w erai31c/param_inputs/namelist.base.oce erai61c/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/scripts/namelist_exp.sh erai61c/scripts/namelist_exp.sh
--------------------------------------------------------------------------------
23c23
< export TSP_OCE=800
---
> export TSP_OCE=400
--------------------------------------------------------------------------------
diff -b -w erai31c/compile/nemo/cpp_PERU12.fcm erai61c/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
diff erai62c / erai61c:
--------------------------------------------------------------------------------
diff -b -w erai61c/scripts/namelist_exp.sh erai62c/scripts/namelist_exp.sh
--------------------------------------------------------------------------------
23c23
< export TSP_OCE=400
---
> export TSP_OCE=720
===>>> : E R R O R (output.abort)
===========
stpctl: the zonal velocity is larger than 20 m/s
======
kt= 9972 max abs(U): 42.72 , i j k: 13 379 13
diff erai63c / erai62c:
--------------------------------------------------------------------------------
diff -b -w erai62c/scripts/namelist_exp.sh erai63c/scripts/namelist_exp.sh
--------------------------------------------------------------------------------
23c23
< export TSP_OCE=720
---
> export TSP_OCE=600
===>>> : E R R O R (output.abort)
===========
stpctl: the zonal velocity is larger than 20 m/s
======
kt= 68824 max abs(U): 31.52 , i j k: 11 378 1
erai6a tourne avec les mêmes forçages que erai21c. Ils ne sont donc pas en cause dans le plantage.
OBC issues du run global Drakkar 46 niveaux ORCA025-B83; ensuite interpollées sur 75 niveaux via matlab
bathy/blk/release nemo identiques aux runs 46 niveaux précédents
mld (warning) XIOS r416
version little_nemo entre r48 et r62 rien n'impacte peru12. Après r62, impacte sur les OBC (div=0 remis au goût du jour)
Démarrage de levitus avec 75 niveaux. Comme pour trop12, le modèle explose très vite. Il faut démarrer avec un pdt très faible.
Janvier 1990
19900101_19900105: rdt=60 OK
19900106_19900110: rdt=120 OK
19900111_19900115: rdt=200 OK
19900116_19900130: rdt=400... plante a time.step=98
19900116_19900131: rdt=200 OK 2700s /16j 170s/j
Fevrier 1990
19900201_19900228: rdt=400... plante a time.step= 438 ... ~2 jours
19900201_19900228: rdt=300... plante a time.step= 1147... ~4 jours
19900201_19900228: rdt=200... plante a time.step= 3742... ~8 jours 1/2
19900201_19900228: rdt=160... plante a time.step= 6096... ~11 jours
Année 1990: 64 coeurs nemo + 1 coeur xios: ~9 heures / an
19900101_19900131: rdt=600 OK
19900201_19900228: rdt=600 OK
19900301_19900331: rdt=600 OK
19900401_19900430: rdt=600 OK
19900501_19900531: rdt=600... PLANTE
the zonal velocity is larger than 20 m/s : kt= 20670 max abs(U): 23.71 , i j k: 4 384 62
à time.step= 20670=~143,5 jours depuis le 19900101
20670-17281 = 3389=23,5 depuis le 19900501
19900501_19900531: rdt=400 OK
19900601_19900630: rdt=400 OK
19900701_19900731: rdt=400 OK
19900801_19900831: rdt=400 OK
19900901_19900930: rdt=400 OK
19901001_19901031: rdt=400 OK
19901101_19901130: rdt=400 OK
19901201_19901231: rdt=400 OK
Année 1991: 63 coeurs nemo + 1 coeur xios: ~9 heures / an
Recompilation avec
%FCFLAGS -i4 -r8 -O3 -xAVX -fp-model precise
au lieu de
%FCFLAGS -i4 -r8
63 nemo + 1 xios au lieu de 64 nemo + 1 xios (préciser le découpage dans la namelist)
jpni = 7 ! jpni number of processors following i (set automatically if < 1)
jpnj = 9 ! jpnj number of processors following j (set automatically if < 1)
jpnij = 63 ! jpnij number of local domains (set automatically if < 1)
PENSER à concatener les restarts (car 63 au lieu de 64) et à modifier l'appel dans getrst_oce.sh
rebuild_nemo -t 4 restart_oce_erai65c_19901231 64
19910101_19911231: rdt=400 OK ~9 heures / an
Année 1992
PENSER à remettre l'appel dans getrst_oce.sh aux 63 parties des restarts
19920101_19921231: rdt=400 OK ~9 heures / an
Année 1993
19930101_19931231: rdt=400 OK ~9 heures / an
Année 1994: 126 coeurs nemo + 2 coeur xios: ???? heures / an
modif traadv_ubs.F90 (MY_SRC) (=> 2 fois moins diffusif)
pas de raison de changer le pas de temps par rapport à bilap/bilap => démarre à 800 s
Ajout de
19900101_19900131: pdt=800 plante a ...
19900101=> 19930331: pdt=600 OK
19930401_19930430_A_pdt_600_PLANTE/
19930401 => 19940331: pdt=400 OK
19940401_19940430_A_pdt_400_PLANTE/
Mais démarrage depuis lévitus => JETER cette partie et garder jusqu'à PLANTE avec pdt=400 => pour comparer sans bug traceur x 2 et 2 types d'obc différentes.
run erai67c: run 75 niveaux (UBS/UBS) - remake erai65c sans bug x2 sur UBS traceurs + rustine obc west
Rustine François: interpole la bathy sur la frontière west au lieu de la dupliquer façon Stephane + Charles.
19900101_19900131: pdt=800 plante a ...
stpctl: the zonal velocity is larger than 20 m/s
19910401_19910430_A_pdt_600_PLANTE/
après OK => 19951201_19951231/
run erai40c: run 46 niveaux (Bilap/Bilap) - pour voir si passe le el niño de 1997
stpctl: the zonal velocity is larger than 20 m/s
======
kt=485320 max abs(U): 22.14 , i j k: 11 374 1
Plante à kt=485320 - 479521 = 5779 pdt => 5779 * 400 = 2311600 secondes => 2311600 / 86400 = 26,7 jours => le 26 Mai 1997
86400= 2^7 * 3^3 * 5^2
1990
1991
1992
1993
1994
1991-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
1991-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
1991-1994
erai4l (bilapdyn + bilaptrac + avec Understress)
erai40c (bilapdyn + bilaptrac + sans Understress)
1990
1991
1992
1993
1994
1991-1994
erai5f (bilapdyn + UBStrac + avec Understress)
erai50c (bilapdyn + UBStrac + sans Understress)
erai51c (bilapdyn + UBStrac + sans Understress)
1990
1991
1992
1993
1994
1991-1994
erai61c (UBSdyn + UBStrac + sans Understress)
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)
erai51c (bilapdyn + UBStrac + sans Understress)
1990
1991
1992
1993
erai61c (UBSdyn + UBStrac + sans Understress)
1990
1991
1992
1993
1994
roms12° (UBSdyn + UBStrac + sans Understress?)
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)
erai51c (bilapdyn + UBStrac + sans Understress)
1990
1991
1992
1993
1994
erai61c (UBSdyn + UBStrac + sans Understress)
1990
1991
1992
1993
1994
roms12° (UBSdyn + UBStrac + sans Understress?)
ke_time_serie (1990-1994)
section (1991-1994)
grad SST (1994 - day365)
Variance ssh serie (1991-1994)
Gradient ssh serie (1991-1994)
Annual ssh Variance (1991-1994)
rms W (1991-1994)
mean SST (1993-1994)
mean SST-obs (1993-1994)
mean mldr10_1 (1993-1994)
erai40c (bilapdyn + bilaptrac + sans Understress)
erai50c (bilapdyn + UBStrac + sans Understress) rn_ahm_0_blp=-1e9
erai51c (bilapdyn + UBStrac + sans Understress) - rn_ahm_0_blp=-2.5e9
erai61c (UBSdyn + UBStrac + sans Understress)
erai40c erai51c erai61c
Observations
Surface Mixed Layer Depth (MLD)
./makenemo -n PERU12 -m X64_CURIE -j0 2> &1 | tee listing_PERU12.txt
Creating PERU12/WORK = OPA_SRC NST_SRC for PERU12
mkdir: cannot create directory `PERU12/MY_SRC': No such file or directory
mkdir: cannot create directory `PERU12/WORK': No such file or directory
MY_SRC directory is : PERU12/MY_SRC
ln: target `PERU12/WORK' is not a directory
ln: target `PERU12/WORK' is not a directory
ln: target `PERU12/WORK' is not a directory
Solution:
C'est qu'une version de cette config a déjà eisté et qu'il faut détruire la ligne correspondant à cette config dans NEMOGCM/CONFIG/cfg.txt
