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
erai6a tourne avec les mêmes forçages que erai21c. Ils ne sont donc pas en cause dans le plantage.
Normalement UBS est plus robuste que les autres shémas puisque la diffusion est proportionnelle à la dynamique.
erai6a_2 plante plus vite que erai6a, ce qui est cohérent puisqu'on ajoutaient (inutilement) dans erai6a un opérateur de diffusion à UBS.
erai60c: Sans understress => plus de meso => explose plus vite?
stpctl: the zonal velocity is larger than 20 m/s ====== kt= 20670 max abs(U): 23.71 , i j k: 4 384 62 nn_it000 = 17281 ( le 01 05 1990) kt - 17281 = 20670 - 17281 = 3389 3389 * 600s = 2033400 / 86400 = 23 jours
stpctl: the zonal velocity is larger than 20 m/s ====== kt=173537 max abs(U): 28.44 , i j k: 14 377 5 nn_it000 = 170641 ( le 01 04 1993) kt - 170641 = 173537 - 170641 = 2896 2896 * 600s = 1737600 / 86400 = 20 jours
stpctl: the zonal velocity is larger than 20 m/s ====== kt= 81128 max abs(U): 23.04 , i j k: 5 376 1 nn_it000 = 78841 ( le 01 04 1994) kt - 78841 = 81128 - 78841 = 2287 2287 * 400s = 914800 / 86400 = 10 jours
stpctl: the zonal velocity is larger than 20 m/s ====== kt= 69835 max abs(U): 28.39 , i j k: 8 379 1 nn_it000 = 65521 ( le 01 04 1990) kt - 65521 = 69835 - 65521 = 4314 4314 * 600s = 2.588.400 / 86400 = 29,95 jours
stpctl: the zonal velocity is larger than 20 m/s ====== kt=485320 max abs(U): 22.14 , i j k: 11 374 1 nn_it000 = 479521 ( le 01 05 1997) kt - 479521 = 485320 - 479521 = 5799 5799 * 400s = 2.319.600 / 86400 = 26 jours
stpctl: the zonal velocity is larger than 20 m/s ====== kt=700768 max abs(U): 21.85 , i j k: 6 380 1 nn_it000 = 696817 ( le 01 02 2000) kt - 696817 = 700768 - 696817 = 3951 3951 * 400s = 1.580.400 / 86400 = 18,29 jours
-------------------------------------------------------------------------------- diff -b -w erai69c/compile/nemo/arch-X64_CURIE.fcm erai70a//compile/nemo/arch-X64_CURIE.fcm -------------------------------------------------------------------------------- 16,17c16,17 < %NCDF_INC -I/usr/local/netcdf-4.2_hdf5/include < %NCDF_LIB -L/usr/local/netcdf-4.2_hdf5/lib -L/usr/local/hdf5-1.8.8/lib -lnetcdff -lnetcdf -lhdf5_hl -lhdf5 -lhdf5 -lz -lcurl -lstdc++ --- > %NCDF_INC -I/usr/local/netcdf-4.2_hdf5_parallel/include > %NCDF_LIB -L/usr/local/netcdf-4.2_hdf5_parallel/lib -L/usr/local/hdf5-1.8.9_parallel/lib -lnetcdff -lnetcdf -lhdf5_hl -lhdf5 -lhdf5 -lz -lcurl -lstdc++
-------------------------------------------------------------------------------- diff -b -w erai69c/param_inputs/namelist.base.oce erai70a/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 41c41 < nn_write = 108 ! frequency of write in the output file (modulo referenced to nn_it000) --- > nn_write = -1 ! frequency of write in the output file (modulo referenced to nn_it000) 127c127 < nn_fsbc = <fsbc> ! frequency of surface boundary condition computation --- > nn_fsbc = 1 ! frequency of surface boundary condition computation 141c141 < ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) --- > ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 251,252c251,252 < ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration < ln_qsr_2bd = .false. ! 2 bands light penetration --- > ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration > ln_qsr_2bd = .true. ! 2 bands light penetration 272,273c272,273 < ln_rnf_mouth = .true. ! specific treatment at rivers mouths < rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used --- > ln_rnf_mouth = .false. ! specific treatment at rivers mouths > rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used 499c499 < ln_bfrimp = .false. ! implicit bottom friction (requires ln_zdfexp = .false. if true) --- > ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 513c513 < nn_bbl_ldf = 0 ! diffusive bbl (=1) or not (=0) --- > nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) 568c568 < rn_aht_0 = 500. ! horizontal eddy diffusivity for tracers [m2/s] --- > rn_aht_0 = 300. ! horizontal eddy diffusivity for tracers [m2/s] 621c621 < ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) --- > ln_dynhpg_imp = .true. ! time stepping: semi-implicit time scheme (T) 661,662c661,662 < rn_avm0 = 1e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") < rn_avt0 = 1e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") --- > rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") > rn_avt0 = 1.e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 692c692 < rn_ebb = 60. ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) --- > rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 758c758 < ln_tmx_itf = .true. ! ITF specific parameterisation --- > ln_tmx_itf = .false. ! ITF specific parameterisation
------------------------------------------------------------------------ r73 | smasson | 2013-02-13 17:03:28 +0100 (Wed, 13 Feb 2013) | 1 line peru12 nested into trop025 ------------------------------------------------------------------------ r72 | smasson | 2013-02-12 23:04:18 +0100 (Tue, 12 Feb 2013) | 1 line put bugfixes from NEMO v3.4.1 ------------------------------------------------------------------------ r71 | smasson | 2013-02-12 10:05:47 +0100 (Tue, 12 Feb 2013) | 1 line add timing of nemo_init and nemo_finalize ------------------------------------------------------------------------ r70 | smasson | 2013-02-12 10:02:33 +0100 (Tue, 12 Feb 2013) | 1 line update arch-X64_CURIE ------------------------------------------------------------------------ r69 | smasson | 2013-01-21 21:39:42 +0100 (Mon, 21 Jan 2013) | 1 line bugfix for oa3mct + agrif + xios in serveur mode ------------------------------------------------------------------------ r68 | smasson | 2013-01-04 21:34:53 +0100 (Fri, 04 Jan 2013) | 1 line minor change in step and finaize timing ------------------------------------------------------------------------ r67 | smasson | 2013-01-04 21:22:11 +0100 (Fri, 04 Jan 2013) | 1 line minor bugfixes related to previous changeset ------------------------------------------------------------------------ r66 | smasson | 2013-01-03 10:25:56 +0100 (Thu, 03 Jan 2013) | 1 line minor bugfixes related to previous changeset ------------------------------------------------------------------------ r65 | smasson | 2013-01-02 18:50:00 +0100 (Wed, 02 Jan 2013) | 1 line bugfix for agrf, obc and xios ------------------------------------------------------------------------ r64 | smasson | 2012-12-27 21:55:06 +0100 (Thu, 27 Dec 2012) | 1 line add back missing part of NEMO 3.4: div must be zero behind the open boundary ------------------------------------------------------------------------ r63 | smasson | 2012-12-27 21:41:35 +0100 (Thu, 27 Dec 2012) | 1 line correct pb with agrif and obc, following changset 3692 of NEMO
-------------------------------------------------------------------------------- diff -b -w -r erai69c/compile/nemo/MY_SRC/dommsk.F90 erai70a//compile/nemo/MY_SRC/dommsk.F90 -------------------------------------------------------------------------------- 182c182 < ispo = 20 ! nombre de points de grille --- > ispo = 21 ! nombre de points de grille
!----------------------------------------------------------------------- &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") !----------------------------------------------------------------------- rn_emin = 1.e-6 ! minimum value of e [m2/s2] rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) rn_clim_galp = 0.53 ! galperin limit ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux rn_charn = 70000. ! Charnock constant for wb induced roughness length nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
-------------------------------------------------------------------------------- diff -b -w erai70a//compile/nemo/cpp_PERU12.fcm erai71a//compile/nemo/cpp_PERU12.fcm -------------------------------------------------------------------------------- 1c1 < bld::tool::fppkeys key_peru12 key_depth=75 key_dynspg_flt key_zdftke key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput key_diahth --- > bld::tool::fppkeys key_peru12 key_depth=75 key_dynspg_flt key_zdfgls key_dtatem key_dtasal key_traldf_c2d key_dynldf_c2d key_mpp_mpi key_obc key_iomput key_diahth
-------------------------------------------------------------------------------- diff -b -w erai71a/param_inputs/namelist.base.oce erai71b/param_inputs/namelist.base.oce -------------------------------------------------------------------------------- 665c665 < ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) --- > ln_zdfevd = .false. ! enhanced vertical diffusion (evd) (T) or not (F)
!----------------------------------------------------------------------- &namzdf ! vertical physics !----------------------------------------------------------------------- rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") rn_avt0 = 1.e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) nn_havtb = 1 ! horizontal shape for avtb (=1) or not (=0) ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) nn_evdm = 1 ! evd apply on tracer (=0) or on tracer and momentum (=1) rn_avevd = 10. ! evd mixing coefficient [m2/s] ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) nn_npc = 0 ! frequency of application of npc nn_npcp = 365 ! npc control print frequency ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T
!----------------------------------------------------------------------- &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") !----------------------------------------------------------------------- rn_emin = 1.e-6 ! minimum value of e [m2/s2] rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) rn_clim_galp = 0.53 ! galperin limit ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux rn_charn = 70000. ! Charnock constant for wb induced roughness length nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
Par contre je suis pas sur que la namelist standart de GLS soit OK. Voila notre bloc (et pas oublie ln_zdfevd=.false.): !----------------------------------------------------------------------- &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") !----------------------------------------------------------------------- rn_emin = 1.e-6 ! minimum value of e [m2/s2] rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) rn_clim_galp = 0.267 ! galperin limit ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux rn_charn = 70000. ! Charnock constant for wb induced roughness length rn_hsro = 0.02 ! Minimum surface roughness rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met=2) nn_z0_met = 2 ! Method for surface roughness computation (0/1/2) nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum) nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
|
|
diags de comparaison Véra |
On retrouve bien des sections comparables pour erai69c (avant phasage) et erai70b (après phasage + retour à RGB). C'est bien le passage de RGB au 2 bands lors du phasage qui réchauffait la surface à la côte de 1° (19° au lieu de 18°)
 |
 |
Globalement pas beaucoup d'impact. Un peu plus de méso peut-être sur le snapshot de grad sst.
 |
diags de comparaison Véra |
Pas d'impact. Runs strictement identiques.
la condition div=0 intervient sur le calcul de la div et donc dans
le bilap, mais pas de div dans UBS donc sans impact sur UBS! Merci Gurvan.
wps sur curie: /ccc/work/cont005/ra0542/hourdinc/wd/pre/peru12_wrf/make_new_grid/wps
Francois avait généré sa grille a partir d'une latitude de reference differente
sans doute une valeur de la grille trop12 avant que Seb ne la modifie.
< ref_lat = -10.44172,
---
> ref_lat = -10.4359,
(Valeur à modifier dans namelist.wps)
Seb avait du augmenter la precision du dx et dy pour retouver la periodicite parfaite pour le bouclage en longitude de la grille trop12
dx = 9264.7894,
dy = 9264.7894,
sur curie: /ccc/work/cont005/ra0542/hourdinc/wd/pre/peru12_now/make_lake_mask
create a lake file and add lake file to geo_em.d01.nc:
sur curie: /ccc/work/cont005/ra0542/hourdinc/wd/pre/peru12_now/oasis
Les fichiers sont crees au format oasis3 dans dir = '/ccc/scratch/cont005/ra0542/hourdinc/peru12_now/inputs/' Il faut ensuite les convertir au format oasis3_mct 1/ Generation des fichiers poids --------------------------------- idl IDL> retall IDL> .r build_peru12_mozaic_oasis IDL> build_peru12_mozaic_oasis, 'peru12' 2/ Build restart oasis ---------------------- Creation des fichiers rst oasis a partir des restarts nemo (pas besoin de restarts wrf car pour la partie atmosphere, on genere des flux a zero pour le restart oasis) dans build_restart_oasis.pro: changer le chemin du restart nemo en fonction de la date voulue penser avant a le reconstruire le restart nemo cd /ccc/scratch/cont005/ra0542/hourdinc/peru12_nemo/restarts/erai70a/20061231 rebuild_nemo -t4 restart_oce_erai70a_20061231 126 idl IDL> retall IDL> .r build_restart_oasis IDL> build_restart_oasis, 'peru12' 3/ oasis3 to oasis3_mct ---------------------- IDL> retall IDL> .r build_oasis2oasismct IDL> build_oasis2oasismct, 'peru12'
Dans build_peru12_mozaic_oasis on traite les points bads.
;------------------------------------------------- ; mozaic oce -> atm (sst + courants de surface) ;------------------------------------------------- bads = where( amask - omask EQ 1, cntbad ) ; keep only the inner domain ;------------------------------------------------- ; mozaic atm -> oce (flux, stress?... ;------------------------------------------------- bads = where( omask - amask EQ 1, cntbad ) ; keep only the inner domain
Au merge des namelistes nemo et wrf, on est parti pour la compilation par erreur du run 71a (test GLS) au lieu de 70a => GLS qui explose dans le forcé en quelques pas de temps faisait également exploser le couplé
 |
fichier: /ccc/scratch/cont005/ra0542/hourdinc/peru12_now/execute/cpl00/20060101_20060101_F_nemo_output_1pdt/auto_T_2D.nc |
Pas vu tout de suite car différences des clés et non des namelistes
=> Toujours comparer les ocean.output
 |
"taum" dans les sorties nemo de cpl00 fichier: /ccc/scratch/cont005/ra0542/hourdinc/peru12_now/execute/cpl00/20071127_20071231/output |
Pour trouver la raison, on a fait de nombreux diagnostiques avec idl à l'intérieur du programme "build_peru12_mozaic_oasis.pro" de génération des fichiers poids oasis (on a commencé par vérifier que le traitement des points bads étaient corrects). Il suffit de relancer le script IDL /ccc/work/cont005/ra0542/hourdinc/wd/pre/peru12_now/oasis/build_peru12_mozaic_oasis.pro pour voir les diags.
Puis on a plotté ensemble les bads et les points "bizarres" vus dans le couplé. C'est points n'étaient pas du tout corrélés :
 |
Le traitement des points bads n'étant à priori pas la cause, on a fait tourné le même run en forcé (à partir des mêmes scripts exactement). On a retrouvé les mêmes points "bizarres". En revanche dans les runs de test de sensibilité peru12_wrf de François, les points "bizarres" n'apparaissaient pas. C'est comme ça qu'on s'est aperçu que:
Raison: Les fichiers inputs wrf n’étant pas encore prêt pour la nouvelle grille, nous avons tourné avec les fichiers d’input de l’ancienne grille pour tester la mise en place du couplé. Les 2 grilles wrf nouvelle et ancienne étant quasi identiques, il y a juste quelques points à la côte qui changent (terre/mer). Cpl00 a donc tourné sur la nouvelle grille avec des inputs ancienne grille, générant ainsi des problèmes sur ces quelques points où wrf calculait un stress sur un point (terre/mer).
 |
"ust" dans wrf de cpl00 et cpl01 fichier: /ccc/scratch/cont005/ra0542/hourdinc/peru12_now/execute/cpl00/20070101_20070130/output |
de 20000101 a 20090630
Les fichiers quikscat sont déjà traités par Seb avec IDL pour combler les trous des données quikscat
On les recupère sur son Mac "visiteur@arete:/Volumes/data/DATAin/wind/orgsrc"
(Sur ciclad.ipsl.jussieu.fr)
On étend avec IDL les valeurs d'océan sur la terre pour que les fichiers puissent être utilisés par NEMO/TOOLS/WEIGHTS pour générer les fichiers de poids pour le on the fly.
Sur ciclad : loceanfs/pulsation/cholod/wd/data/quikscat/1_fill_land.sh
NEMO/TOOLS/WEIGHTS fait une interpolation bicubique et a besoin de 4 voisins => il lui faut pour les points le long de la côte des valeurs sur la terre. Parfois ces valeurs existent, mais si ce sont par exemple des vents de terre, il faut les écraser par des vents de mer!
Pour utiliser les sorties de peru12_wrf pour forcer peru12_nemo, pas besoin d'étendre les valeurs sur la côte car la grille est la même pour les forçages (peru12_wrf) que pour nemo (peru12_nemo)! Inutile donc de faire du "on the fly"
(Sur curie.ccc.cea.fr)
cf maketools WEIGHTS
les namelistes sur curie: /ccc/work/cont005/ra0542/hourdinc/wd/pre/peru12_nemo/weights
(Sur curie.ccc.cea.fr)
Pour phaser, on modifie la bathy de nemo uniquement. On ne touche pas à wrf!
suite:
 |
| erai70e | erai70c / erai70e |
 |
 |
| JFM | AMJ | JAS | OND | annual |
| erai70c (2000-2007) | ||||
 |
 |
 |
 |
 |
| erai70e (2000-2007) | ||||
 |
 |
 |
 |
 |
| erai70e - erai70c (2000-2007) | ||||
 |
 |
 |
 |
 |
