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)
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
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------