Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found
Select Git revision

Target

Select target project
  • mesonh/mesonh-code
  • quentin.rodier/mesonh-code-fork
  • benoit.vie/mesonh-code
  • joris.pianezze/mesonh-code
  • 8qq4g5s7/mesonh-code
  • jean.baptiste.filippi/meso-nh-fire-code
  • fdl68d9p/mesonh-code-sophia
  • jean.wurtz/mesonh-code
8 results
Select Git revision
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 1636 additions and 0 deletions
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/CAM1.nam 0 → 100644
View file @ f77924f4
CAM2MESONH
transfer CAMS/ReLACS variables (default values) included aerosols
37 # NUMBER OF OPTIONAL GRIB VARIABLES
(A6,1X,I1,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2)
O3 1 1.00 go3 48.00
CO 1 1.00 co 28.00
H2O2 1 1.00 h2o2 34.00
NO 1 1.00 no 30.00
NO2 1 1.00 no2 46.00
HNO3 1 1.00 hno3 63.00
HNO4 1 1.00 ho2no2 79.00
N2O5 1 1.00 n2o5 108.00
NO3 1 1.00 no3 62.00
HO2 1 1.00 ho2 33.00
SO2 1 1.00 so2 64.10
DMS 1.1.00 dms 62.10
NH3 1 1.00 nh3 17.00
MO2 1 1.00 ch3o2 47.00
ETH 1 1.00 c2h6 30.02
ALKA 4 0.39 c3h8 44.03 0.33 ch3oh 31.01 0.96 c2h5oh 46.02 1.03 par 12.00
ALKE 3 0.88 c2h4 28.00 1.19 c3h6 42.03 1.20 ole 24.00
BIO 2 1.00 c5h8 68.10 0.99 c10h16 136.00
KET 1 0.28 ch3coch3 58.00
CARBO 2 0.96 ch3cocho 72.10 1.07 ispd 70.00
PAN 1 1.00 pan 121.00
ONIT 1 1.00 onit 77.00
HCHO 1 1.00 hcho 30.00
ALD 1 1.00 ald2 24.00
ORA1 1 1.00 hcooh 46.01
OP1 1 1.00 ch3ooh 48.00
OP2 1 1.00 rooh 47.00
SO4I 1 0.20 aermr11 1.00
SO4J 1 0.80 aermr11 1.00
NH3I 1 0.20 nh4 1.00
NH3J 1 0.80 nh4 1.00
NO3I 1 0.20 no3_a 1.00
NO3J 1 0.80 no3_a 1.00
BCI 2 0.20 aermr10 1.00 0.20 aermr09 1.00
OCI 2 0.20 aermr08 1.00 0.20 aermr07 1.00
BCJ 2 0.80 aermr10 1.00 0.80 aermr09 1.00
OCJ 2 0.80 aermr08 1.00 0.80 aermr07 1.00
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/CAM1_relacs.nam 0 → 100644
View file @ f77924f4
CAM2MESONH
transfer CAMS/ReLACS variables (default values) included aerosols
37 # NUMBER OF OPTIONAL GRIB VARIABLES
(A6,1X,I1,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2,1X,F4.2,1X,A8,1X,F6.2)
O3 1 1.00 go3 48.00
CO 1 1.00 co 28.00
H2O2 1 1.00 h2o2 34.00
NO 1 1.00 no 30.00
NO2 1 1.00 no2 46.00
HNO3 1 1.00 hno3 63.00
HNO4 1 1.00 ho2no2 79.00
N2O5 1 1.00 n2o5 108.00
NO3 1 1.00 no3 62.00
HO2 1 1.00 ho2 33.00
SO2 1 1.00 so2 64.10
DMS 1.1.00 dms 62.10
NH3 1 1.00 nh3 17.00
MO2 1 1.00 ch3o2 47.00
ETH 1 1.00 c2h6 30.02
ALKA 4 0.39 c3h8 44.03 0.33 ch3oh 31.01 0.96 c2h5oh 46.02 1.03 par 12.00
ALKE 3 0.88 c2h4 28.00 1.19 c3h6 42.03 1.20 ole 24.00
BIO 2 1.00 c5h8 68.10 0.99 c10h16 136.00
KET 1 0.28 ch3coch3 58.00
CARBO 2 0.96 ch3cocho 72.10 1.07 ispd 70.00
PAN 1 1.00 pan 121.00
ONIT 1 1.00 onit 77.00
HCHO 1 1.00 hcho 30.00
ALD 1 1.00 ald2 24.00
ORA1 1 1.00 hcooh 46.01
OP1 1 1.00 ch3ooh 48.00
OP2 1 1.00 rooh 47.00
SO4I 1 0.20 aermr11 1.00
SO4J 1 0.80 aermr11 1.00
NH3I 1 0.20 nh4 1.00
NH3J 1 0.80 nh4 1.00
NO3I 1 0.20 no3_a 1.00
NO3J 1 0.80 no3_a 1.00
BCI 2 0.20 aermr10 1.00 0.20 aermr09 1.00
OCI 2 0.20 aermr08 1.00 0.20 aermr07 1.00
BCJ 2 0.80 aermr10 1.00 0.80 aermr09 1.00
OCJ 2 0.80 aermr08 1.00 0.80 aermr07 1.00
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/PRE_REAL1.nam 0 → 100644
View file @ f77924f4
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20190328.06.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_06.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_VER_GRID NKMAX=45,
LTHINSHELL=.TRUE.,
YZGRID_TYPE='FUNCTN',
ZDZGRD=60.,
ZDZTOP=1000.,
ZZMAX_STRGRD=500.,
ZSTRGRD=7.,
ZSTRTOP=10. /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/PRE_REAL1.nam_06 0 → 100644
View file @ f77924f4
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20190328.06.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_06.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_VER_GRID NKMAX=45,
LTHINSHELL=.TRUE.,
YZGRID_TYPE='FUNCTN',
ZDZGRD=60.,
ZDZTOP=1000.,
ZZMAX_STRGRD=500.,
ZSTRGRD=7.,
ZSTRTOP=10. /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/PRE_REAL1.nam_12 0 → 100644
View file @ f77924f4
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ecmwf.OD.20190328.12.00',
HATMFILETYPE='GRIBEX',
HPGDFILE='BIOMAIDO_CH_8km.nestne',
HCHEMFILE = 'CAMS_20190328_12.nc',
HCHEMFILETYPE='CAMSEU',
CINIFILE='ECMWF_D1_20190328_12'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_VER_GRID NKMAX=45,
LTHINSHELL=.TRUE.,
YZGRID_TYPE='FUNCTN',
ZDZGRD=60.,
ZDZTOP=1000.,
ZZMAX_STRGRD=500.,
ZSTRGRD=7.,
ZSTRTOP=10. /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/004_prep_real_case_d1/run_prep_real_case.sh 0 → 100755
View file @ f77924f4
#!/bin/sh
#SBATCH -J prep_megan
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 1 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 01:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
ln -sf ../003_prep_nest/BIOM* .
ln -sf ~/SAVE/GRIB_KTEST/CAMS_20190328_* .
ln -sf ~/SAVE/GRIB_KTEST/ecmwf.OD.20190328.* .
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 1"
set -x
set -e
cp PRE_REAL1.nam_06 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cp PRE_REAL1.nam_12 PRE_REAL1.nam
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cd ../005_prep_real_case_d2
sbatch run_prep_real_case.sh
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/005_prep_real_case_d2/PRE_REAL1.nam 0 → 100644
View file @ f77924f4
&NAM_CONFIO LCDF4=.TRUE., LLFIOUT=.FALSE., LLFIREAD=.FALSE.,
NIO_ABORT_LEVEL=0, NGEN_ABORT_LEVEL=0 /
&NAM_REAL_CONF NVERB=10, CEQNSYS='DUR', CPRESOPT='CRESI', LRES=.FALSE., XRES=1.E-07 /
&NAM_FILE_NAMES HATMFILE='ECMWF_D1_20190328_06.spa',
HATMFILETYPE='MESONH',
HPGDFILE='BIOMAIDO_CH_2km.nestne',
CINIFILE='ECMWF_D2_20190328_06'/
&NAM_PREP_TEB LTEB_CANOPY=T /
&NAM_PREP_ISBA LISBA_CANOPY=T /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
&NAM_CONF NHALO=0 /
&NAM_CH_CONF LUSECHEM=T /
&NAM_AERO_CONF LORILAM = .TRUE.,
LSALT = .TRUE.,
LDUST = .TRUE.,
LINITPM = .FALSE.,
CCLOUD = "LIMA",
NMOD_CCN = 3,
NMOD_IFN = 2,
LSLTCAMS = T,
LDSTCAMS = T,
XINIRADIUSI = 0.015,
XINIRADIUSJ = 0.055,
CRGUNIT = "NUMB",
XINISIGI = 1.4,
XINISIGJ = 1.6,
XN0IMIN = 1E2,
XN0JMIN = 1E-2 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/005_prep_real_case_d2/SPAWN1.nam 0 → 100644
View file @ f77924f4
&NAM_LUNIT2_SPA CINIFILE = "ECMWF_D1_20190328_06",
CINIFILEPGD = "BIOMAIDO_CH_8km.nestne",
YDOMAIN = "BIOMAIDO_CH_2km.nestne",
YSPANBR = '' /
&NAM_CONFIO LCDF4 = T, NIO_ABORT_LEVEL=1, NGEN_ABORT_LEVEL=1 /
&NAM_PREP_SURF_ATM NHALO_PREP=0 /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/005_prep_real_case_d2/run_prep_real_case.sh 0 → 100755
View file @ f77924f4
#!/bin/sh
#SBATCH -J prep_d2_megan
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 1 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 00:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
ln -sf ../003_prep_nest/BIOM* .
ln -sf ../004_prep_real_case_d1/ECMWF* .
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
export MPIRUN="Mpirun -np 1"
set -x
set -e
time ${MPIRUN} SPAWNING${XYZ}
time ${MPIRUN} PREP_REAL_CASE${XYZ}
cd ../006_mesonh_283TO3003
sbatch run_mesonh.sh
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/DIAG1.nam 0 → 100644
View file @ f77924f4
&NAM_CONFIO LCDF4 = .TRUE.,
LLFIREAD = .FALSE.,
LLFIOUT = .FALSE. /
&NAM_DIAG_FILE YINIFILE(1) = 'A2830.1.MEGA2.001',
YINIFILEPGD(1) = 'BIOMAIDO_CH_8km.nestne',
YSUFFIX = 'diag' /
&NAM_DIAG LTPZH = .TRUE,
LMSLP = .TRUE.,
LTHW = .TRUE.,
CISO='PREVTK',
LCHEMDIAG = .TRUE. /
&NAM_DIAG_SURFn N2M=2,
LSURF_BUDGET=.TRUE.,
LCOEF=.TRUE. /
&NAM_DIAG_SURF_ATMn LFRAC=.TRUE.,
LDIAG_GRID=.TRUE. /
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/EXSEG1.nam 0 → 100644
View file @ f77924f4
&NAM_CONF CCONF = "START", NVERB=1, NMODEL = 2,
CEXP = "A2830", CSEG = "MEGAD" , NHALO=3,
CSPLIT="BSPLITTING" /
&NAM_CONFIO LCDF4 = .TRUE., LLFIOUT = .FALSE., NIO_ABORT_LEVEL=1, NGEN_ABORT_LEVEL=1 /
&NAM_DYN XSEGLEN = 21600.,
LCORIO = T,
LNUMDIFU = F,
XALKTOP = 0.005,
XALZBOT = 15000. /
&NAM_SERIESn NFREQSERIES = 30 /
&NAM_LUNITn CINIFILE = "ECMWF_D1_20190328_06",
CINIFILEPGD="BIOMAIDO_CH_8km.nestne",
CCPLFILE(1) = "ECMWF_D1_20190328_12" /
&NAM_NESTING NDAD(2) = 1, NDTRATIO(2) = 5, XWAY(2) = 2. /
&NAM_DYNn XTSTEP = 30.,
CPRESOPT = "ZRESI",
NITR = 12,
LHORELAX_UVWTH = F, LHORELAX_RV = F, LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.01 /
&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",NWENO_ORDER=5,CTEMP_SCHEME="RK53",
CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
&NAM_PARAMn CCLOUD = "LIMA", CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "NONE", CSCONV="EDKF" /
&NAM_PARAM_RADn XDTRAD = 600., XDTRAD_CLONLY = 600., NRAD_COLNBR = 400 /
&NAM_PARAM_KAFRn XDTCONV = 360., NICE = 1, LCHTRANS = F,
LREFRESH_ALL = T, LDOWN = T /
&NAM_PARAM_MFSHALLn CMF_UPDRAFT = "EDKF", CMF_CLOUD = "DIRE" /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LSIG_CONV = F, LTURB_FLX=F /
&NAM_CH_MNHCn LUSECHEM = .TRUE.,
LCH_CONV_LINOX = .FALSE.,
LCH_INIT_FIELD = .FALSE.,
LCH_CONV_SCAV = .TRUE.,
CCHEM_INPUT_FILE = "ReLACS_MACCity_aero.nam",
NCH_SUBSTEPS = 2,
LCH_TUV_ONLINE = .FALSE.,
CCH_TUV_LOOKUP = "PHOTO.TUV5",
CCH_TUV_CLOUDS = "CHAN",
XCH_TUV_TUPDATE = 300.,
XCH_TUV_DOBNEW = 250.,
CCH_VEC_METHOD = 'MAX'
NCH_VEC_LENGTH = 50 /
&NAM_CH_SOLVERn CSOLVER = "EXQSSA",
XSLOW = 90.,
XFAST = 0.1,
XDTMIN = 0.1,
XDTMAX = 45.,
XDTFIRST = 10. /
&NAM_BACKUP XBAK_TIME(1,1)=120. ,
XBAK_TIME(1,2)=7200. ,
XBAK_TIME(1,3)=10800. ,
XBAK_TIME(1,4)=14400. ,
XBAK_TIME(1,5)=18000. ,
XBAK_TIME(1,6)=21600. ,
XBAK_TIME(1,7)=25200. ,
XBAK_TIME(1,8)=28800. ,
XBAK_TIME(1,9)=32400. ,
XBAK_TIME(1,10)=36000. ,
XBAK_TIME(1,11)=39600. ,
XBAK_TIME(1,12)=43200. /
&NAM_PARAM_LIMA LACTI = T,
HINI_CCN='AER',
HTYPE_CCN(1) = 'C',
HTYPE_CCN(2) = 'C',
HTYPE_CCN(3) = 'C',
NMOD_CCN = 3, ! Nombre de modes de CCN
LACTIT=T,
LSEDC=T,
LSCAV=T, ! Scavenging
LAERO_MASS=F,
LCCN_HOM=F, ! Concentration en CCN homogène selon z
CCCN_MODES='CAMS_JPP', ! Choix des modes de CCN
XCCN_CONC(1) = 300., ! Concentration du mode 1 de CCN (/cm3)
XCCN_CONC(2) = 140., ! Concentration du mode 2 de CCN (/cm3)
XCCN_CONC(3) = 50., ! Concentration du mode 3 de CCN (/cm3)
XCCN_CONC(4) = 0.01, ! Concentration du mode 4 de CCN (/cm3)
LNUCL = T,
LSEDI = T,
LHHONI = F,
NMOD_IFN = 2, ! Nombre de modes d'IFN
CPRISTINE_ICE_LIMA = 'PLAT',
CHEVRIMED_ICE_LIMA = 'GRAU',
NIND_SPECIE = 1,
LMEYERS = F,
NMOD_IMM = 1, ! Nombre de modes de CCN agissant par IMM
LIFN_HOM = T, ! Concentration en IFN homogène selon z
CIFN_SPECIES = '', !
CINT_MIXING = 'MACC_JPP', ! DM1, DM2, BC, O, MACC, default (Phillips 2008)
XIFN_CONC(1) = 10., ! Concentration du mode 1 d'IFN (/L)
XIFN_CONC(2) = 10. / ! Concentration du mode 2 d'IFN (/L)
&NAM_CH_ORILAM LORILAM = .TRUE.,
LVARSIGI = .FALSE.,
LVARSIGJ = .FALSE.,
LRGFIX = .FALSE.,
LSEDIMAERO = .TRUE.,
LDEPOS_AER(1) = .TRUE.,
LDEPOS_AER(2) = .TRUE.,
XN0IMIN = 5.,
XN0JMIN = 0.001,
XINIRADIUSI= 0.015,
XINIRADIUSJ= 0.055,
XINISIGI = 1.400,
XINISIGJ = 1.600,
CNUCLEATION = "MAATTANEN_NEUTRAL",
LCONDENSATION = .TRUE.,
LCOAGULATION = .TRUE.,
LMODE_MERGING = .TRUE.,
XRADIUS_NUCL = 15E-3,
XSIGMA_NUCL = 1.2 /
&NAM_DIAG_SURFn N2M = 0, LSURF_BUDGET = F /
&NAM_ISBAn CSCOND="NP89", CC1DRY="DEF",
CSOILFRZ="DEF", CDIFSFCOND="DEF",
CSNOWRES="DEF" /
&NAM_SSOn CROUGH="Z01D" /
&NAM_DIAG_ISBAn /
&NAM_DUST LDUST = .FALSE.,
LRGFIX_DST = .FALSE.,
LVARSIG = .FALSE.,
LSEDIMDUST = .TRUE.,
LDEPOS_DST(1) = .TRUE.,
LDEPOS_DST(2) = .TRUE.,
LDEPOS_DST(3) = .TRUE.,
NMODE_DST = 3 /
&NAM_SALT LSALT = .TRUE.,
LRGFIX_SLT = .FALSE.,
LVARSIG_SLT = .FALSE.,
LSEDIMSALT = .TRUE.,
LDEPOS_SLT(1) = .TRUE.,
LDEPOS_SLT(2) = .TRUE.,
CRGUNITS = 'NUMB',
NMODE_SLT = 8 /
&NAM_SURF_SLT CEMISPARAM_SLT='OvB21b' /
&NAM_SURF_DST CEMISPARAM_DST='AMMA', CVERMOD='CMDVER' /
&NAM_SEAFLUXn CSEA_ALB="TA96", LPRECIP = T /
&NAM_CH_SURFn LCH_SURF_EMIS=T, LCH_DMSEMIS=T /
&NAM_CH_CONTROLn CCHEM_SURF_FILE= "ReLACS_MACCity_aero.nam" /
&NAM_CH_ISBAn LCH_BIO_FLUX=T, CCH_DRY_DEP="WES89", CPARAMBVOC="MEGAN",
XDAILYPAR=300., XDAILYTEMP=290. /
&NAM_CH_SEAFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_WATFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_TEBn CCH_DRY_DEP="WES89" /
&NAM_CHS_ORILAM LCH_AERO_FLUX = .TRUE.,
LCO2PM=.FALSE.,
XEMISRADIUSI = 0.015,
XEMISRADIUSJ = 0.055,
XEMISSIGI = 1.40,
XEMISSIGJ = 1.60,
CRGUNIT = "NUMB" /
&NAM_BLANK /
&NAM_NEBn LSIGMAS=T, LSUBG_COND=F /
&NAM_PARAM_ICEn CSUBG_AUCV_RC='NONE'/
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/EXSEG2.nam 0 → 100644
View file @ f77924f4
&NAM_LUNITn CINIFILE = "ECMWF_D2_20190328_06",
CINIFILEPGD="BIOMAIDO_CH_2km.nestne" /
&NAM_DYNn CPRESOPT = "ZRESI",
NITR = 12,
LVE_RELAX = T,
NRIMX = 5, NRIMY = 5, XRIMKMAX = 0.0083,LHORELAX_RH=F /
&NAM_ADVn CUVW_ADV_SCHEME="WENO_K",NWENO_ORDER=5,CTEMP_SCHEME="RK53",
CMET_ADV_SCHEME = "PPM_01", CSV_ADV_SCHEME = "PPM_01" /
&NAM_PARAMn CCLOUD = "LIMA", CTURB = "TKEL", CRAD = "ECMW",
CDCONV = "NONE", CSCONV= "EDKF" /
&NAM_PARAM_MFSHALLn CMF_UPDRAFT='EDKF',CMF_CLOUD='DIRE', LMIXUV=T /
&NAM_PARAM_RADn XDTRAD = 600., XDTRAD_CLONLY = 600., NRAD_COLNBR = 400 /
&NAM_LBCn CLBCX = 2*"OPEN", CLBCY = 2*"OPEN" /
&NAM_TURBn CTURBLEN = "BL89", CTURBDIM = "1DIM", LTURB_FLX=F /
&NAM_CH_MNHCn LUSECHEM = .TRUE.,
LCH_CONV_LINOX = .FALSE.,
LCH_INIT_FIELD = .FALSE.,
LCH_CONV_SCAV = .TRUE.,
CCHEM_INPUT_FILE = "ReLACS_MACCity_aero.nam",
NCH_SUBSTEPS = 2,
LCH_TUV_ONLINE = .FALSE.,
CCH_TUV_LOOKUP = "PHOTO.TUV5",
CCH_TUV_CLOUDS = "CHAN",
XCH_TUV_TUPDATE = 300.,
XCH_TUV_DOBNEW = 250.,
CCH_VEC_METHOD = 'MAX'
NCH_VEC_LENGTH = 50 /
&NAM_CH_SOLVERn CSOLVER = "EXQSSA",
XSLOW = 90.,
XFAST = 0.1,
XDTMIN = 0.1,
XDTMAX = 45.,
XDTFIRST = 10. /
&NAM_DIAG_SURFn N2M = 0, LSURF_BUDGET = F /
&NAM_ISBAn CSCOND="NP89",
CC1DRY="DEF", CSOILFRZ="DEF", CDIFSFCOND="DEF", CSNOWRES="DEF" /
&NAM_SSOn CROUGH="Z01D" /
&NAM_DIAG_ISBAn /
&NAM_SEAFLUXn CSEA_ALB="TA96", LPRECIP = T /
&NAM_CH_SURFn LCH_SURF_EMIS=T /
&NAM_CH_CONTROLn CCHEM_SURF_FILE= "ReLACS_MACCity_aero.nam" /
&NAM_CH_ISBAn LCH_BIO_FLUX=T, CCH_DRY_DEP="WES89", CPARAMBVOC="MEGAN",
XDAILYPAR=150., XDAILYTEMP=290. /
&NAM_CH_SEAFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_WATFLUXn CCH_DRY_DEP="WES89" /
&NAM_CH_TEBn CCH_DRY_DEP="WES89" /
&NAM_NEBn LSIGMAS=T , LSUBG_COND = F /
&NAM_PARAM_ICEn CSUBG_AUCV_RC='NONE'/
\ No newline at end of file
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/ReLACS_MACCity_aero.nam 0 → 100644
View file @ f77924f4
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
FORMPROF
PROFILES LEVELS
12 9
(F7.1,12F7.4)
0.0 0.000 0.660 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
500.0 0.000 0.660 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
1000.0 0.000 0.770 0.400 0.660 0.160 1.000 1.000 1.000 0.500 0.200 1.000 1.000
1500.0 0.000 0.880 0.020 0.330 0.016 1.000 1.000 1.000 0.200 0.200 1.000 1.000
3000.0 0.000 0.660 0.010 0.130 0.003 0.500 1.000 1.000 0.100 0.002 1.000 1.000
4500.0 0.000 0.770 0.003 0.130 0.003 0.500 1.000 1.000 0.050 0.000 1.000 0.010
6000.0 0.000 1.000 0.003 0.130 0.003 1.000 1.000 0.500 0.050 0.000 1.000 0.010
10000.0 0.000 1.000 0.003 0.000 0.003 1.000 1.000 1.000 0.050 0.000 1.000 0.010
30000.0 1.000 1.000 0.003 0.000 0.003 1.000 1.000 1.000 0.050 0.000 1.000 0.010
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
NORMINIT
!!PUT YOUR OWN VALUES HERE!! initial values (units are par per part = MIX)
MIX
3
(1X,a12,1X,E25.8)
'SO2 ' 1.00000000E-12
'NO ' 1.00000000E-12
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
PROFASSO
norm-profiles to be associated
2
(1X,A12,1X,I3)
'SO2 ' 2
'NO ' 7
=====================================================================
*** the following section will be read by ch_init_dep_isban.F90 ***
=====================================================================
SURF_RES
surface resistances (s/m), refer to Seinfeld and Pandis, 1998, p. 975, Tab.19.2
7
(1X,A12,1X,E24.8)
'RCSNOWO3 ' 0.99990000E+04
'RCSNOWSO2 ' 0.99990000E+04
'RCCLAYO3 ' 0.99990000E+04
'RCCLAYSO2 ' 0.99990000E+04
'RCSANDO3 ' 0.99990000E+04
'RCSANDSO2 ' 0.99990000E+04
'LANDREXT ' 0.99990000E+04
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
MASS_MOL
molecular mass (in g/mol) for molecular diffusion
41
(A32,2E15.5)
O3 0.48000E+02
H2O2 0.34000E+02
NO 0.30000E+02
NO2 0.46000E+02
NO3 0.62000E+02
N2O5 0.10800E+03
HONO 0.47000E+02
HNO3 0.63000E+02
HNO4 0.79000E+02
NH3 0.17000E+02
SO2 0.64000E+02
SULF 0.98000E+02
CO 0.28000E+02
OH 0.17000E+02
HO2 0.33000E+02
CH4 0.16000E+02
ETH 0.30000E+02
ALKA 0.50861E+02
ALKE 0.29160E+02
BIO 0.10200E+03
ARO 0.10100E+03
HCHO 0.30000E+02
ALD 0.44000E+02
KET 0.72000E+02
CARBO 0.66900E+02
ONIT 0.11900E+03
PAN 0.12100E+03
OP1 0.48000E+02
OP2 0.63100E+02
ORA1 0.46000E+02
ORA2 0.60000E+02
MO2 0.47000E+02
ALKAP 0.81380E+02
ALKEP 0.83611E+02
BIOP 0.11700E+03
PHO 0.10700E+03
ADD 0.11680E+03
AROP 0.14867E+03
CARBOP 0.85434E+02
OLN 0.13600E+03
XO2 0.10000E+03
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
REA_FACT
reactivity factor with biology, Seinfeld and Pandis, 1998, p. 975, Tab. 19.3
41
(A32,2E15.5)
O3 0.10000E+01
H2O2 0.10000E+01
NO 0.00000E+00
NO2 0.10000E+00
NO3 0.10000E+00
N2O5 0.10000E+00
HONO 0.10000E+00
HNO3 0.00000E+00
HNO4 0.00000E+00
NH3 0.00000E+00
SO2 0.00000E+00
SULF 0.00000E+00
CO 0.00000E+00
OH 0.00000E+00
HO2 0.00000E+00
CH4 0.00000E+00
ETH 0.00000E+00
ALKA 0.00000E+00
ALKE 0.00000E+00
BIO 0.00000E+00
ARO 0.00000E+00
HCHO 0.00000E+00
ALD 0.00000E+00
KET 0.00000E+00
CARBO 0.00000E+00
ONIT 0.00000E+00
PAN 0.10000E+00
OP1 0.30000E+00
OP2 0.10000E+00
ORA1 0.00000E+00
ORA2 0.00000E+00
MO2 0.00000E+00
ALKAP 0.00000E+00
ALKEP 0.00000E+00
BIOP 0.00000E+00
PHO 0.00000E+00
ADD 0.00000E+00
AROP 0.00000E+00
CARBOP 0.00000E+00
OLN 0.00000E+00
XO2 0.00000E+00
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
HENRY_SP
Effective Henrys law factor / exponent, See Leriche et al.2013
35
(A32,2E15.5)
O3 1.03000E-02 -0.28300E+04
H2O2 8.44000E+04 -0.76000E+04
NO 1.92000E-03 -0.17900E+04
NO2 1.20000E-02 -0.25160E+04
NO3 3.80000E-02 -0.87070E+04
N2O5 2.10000E+00 -0.34000E+04
HONO 8.38000E+04 -0.31200E+04
HNO3 1.46000E+13 -0.10500E+05
HNO4 4.78000E+04 -0.69000E+04
NH3 3.24000E+03 0.19000E+03
SO2 5.59000E+04 -0.48950E+04
SULF 6.64000E+14 -0.87000E+04
CO 9.81000E-04 -0.17200E+04
OH 3.90000E+01 -0.00000E+00
HO2 3.49000E+04 -0.00000E+00
CH4 1.41000E-03 -0.20400E+04
ETH 1.88000E-03 -0.28750E+05
ALKA 0.15000E-02 -0.32750E+04
ALKE 0.59600E-02 -0.21700E+04
BIO 0.38500E-01 -0.00000E+00
ARO 0.18000E+00 -0.41000E+04
HCHO 3.23000E+03 -0.71960E+04
ALD 0.12900E+02 -0.58900E+04
KET 0.27800E+02 -0.55300E+04
CARBO 0.36000E+06 -0.75450E+04
ONIT 0.10000E+02 -0.59100E+04
PAN 0.28000E+01 -0.57300E+04
OP1 0.30000E+01 -0.52800E+04
OP2 0.33600E+01 -0.59950E+04
ORA1 5.07000E+06 -0.59500E+04
ORA2 2.66000E+05 -0.62000E+04
MO2 0.24500E+01 -0.23200E+04
ALKAP 0.26600E+03 -0.60000E+04
ALKEP 0.83000E+03 -0.60000E+04
BIOP 0.23110E+04 -0.60000E+04
================================================================
*** the following section will be read by build_emisstabn.F90 ***
*** <called by ch_init_emissionn.F90> ***
*** CON: flux is given in molecules/cm2/s ***
================================================================
EMISUNIT
Emission
CON
================================================================
*** the following section will be read by build_pronoslistn.F90 ***
*** <called by ch_init_emissionn.F90> ***
* Format: Nom espece pronostic Meso-NH - Facteur Aggregation - NOM ESPECE EMISE
* expl: NO (espece Meso-NH) 0.70 (facteur aggregation) NOXAN (espece emise PGD)
================================================================
AGREGATION
EMISSIONS MACC-CITY - schema ReLACS
NO 0.70 NOXAN
NO2 0.30 NOXAN
SO2 1.00 SO2AN
NH3 1.00 NH3AN
CO 1.00 COAN
ETH 1.00 C2H6AN
ALKA 0.39 C3H8AN 1.25 BALKAN 0.33 MTOLAN 1.29 BOLAN
ALKE 0.90 C2H4AN 1.18 C3H6AN 1.20 BENEAN
ARO 1.00 AROAN
HCHO 1.00 CH2OAN
ALD 1.00 BALDAN
KET 0.36 ACETAN 2.35 BONEAN
BCJ 1.0 BCAN 1.0 BCANSEA
OCJ 1.0 OCAN 1.0 OCANSEA
END_AGREGATION
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/ReLACS_MACCity_gas.nam 0 → 100644
View file @ f77924f4
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
FORMPROF
PROFILES LEVELS
12 9
(F7.1,12F7.4)
0.0 0.000 0.660 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
500.0 0.000 0.660 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
1000.0 0.000 0.770 0.400 0.660 0.160 1.000 1.000 1.000 0.500 0.200 1.000 1.000
1500.0 0.000 0.880 0.020 0.330 0.016 1.000 1.000 1.000 0.200 0.200 1.000 1.000
3000.0 0.000 0.660 0.010 0.130 0.003 0.500 1.000 1.000 0.100 0.002 1.000 1.000
4500.0 0.000 0.770 0.003 0.130 0.003 0.500 1.000 1.000 0.050 0.000 1.000 0.010
6000.0 0.000 1.000 0.003 0.130 0.003 1.000 1.000 0.500 0.050 0.000 1.000 0.010
10000.0 0.000 1.000 0.003 0.000 0.003 1.000 1.000 1.000 0.050 0.000 1.000 0.010
30000.0 1.000 1.000 0.003 0.000 0.003 1.000 1.000 1.000 0.050 0.000 1.000 0.010
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
NORMINIT
!!PUT YOUR OWN VALUES HERE!! initial values (units are par per part = MIX)
MIX
3
(1X,a12,1X,E25.8)
'SO2 ' 1.00000000E-12
'NO ' 1.00000000E-12
===================================================================
*** the following section will be read by ch_field_valuen.f90 ***
===================================================================
PROFASSO
norm-profiles to be associated
2
(1X,A12,1X,I3)
'SO2 ' 2
'NO ' 7
=====================================================================
*** the following section will be read by ch_init_dep_isban.F90 ***
=====================================================================
SURF_RES
surface resistances (s/m), refer to Seinfeld and Pandis, 1998, p. 975, Tab.19.2
7
(1X,A12,1X,E24.8)
'RCSNOWO3 ' 0.99990000E+04
'RCSNOWSO2 ' 0.99990000E+04
'RCCLAYO3 ' 0.99990000E+04
'RCCLAYSO2 ' 0.99990000E+04
'RCSANDO3 ' 0.99990000E+04
'RCSANDSO2 ' 0.99990000E+04
'LANDREXT ' 0.99990000E+04
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
MASS_MOL
molecular mass (in g/mol) for molecular diffusion
41
(A32,2E15.5)
O3 0.48000E+02
H2O2 0.34000E+02
NO 0.30000E+02
NO2 0.46000E+02
NO3 0.62000E+02
N2O5 0.10800E+03
HONO 0.47000E+02
HNO3 0.63000E+02
HNO4 0.79000E+02
NH3 0.17000E+02
SO2 0.64000E+02
SULF 0.98000E+02
CO 0.28000E+02
OH 0.17000E+02
HO2 0.33000E+02
CH4 0.16000E+02
ETH 0.30000E+02
ALKA 0.50861E+02
ALKE 0.29160E+02
BIO 0.10200E+03
ARO 0.10100E+03
HCHO 0.30000E+02
ALD 0.44000E+02
KET 0.72000E+02
CARBO 0.66900E+02
ONIT 0.11900E+03
PAN 0.12100E+03
OP1 0.48000E+02
OP2 0.63100E+02
ORA1 0.46000E+02
ORA2 0.60000E+02
MO2 0.47000E+02
ALKAP 0.81380E+02
ALKEP 0.83611E+02
BIOP 0.11700E+03
PHO 0.10700E+03
ADD 0.11680E+03
AROP 0.14867E+03
CARBOP 0.85434E+02
OLN 0.13600E+03
XO2 0.10000E+03
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
REA_FACT
reactivity factor with biology, Seinfeld and Pandis, 1998, p. 975, Tab. 19.3
41
(A32,2E15.5)
O3 0.10000E+01
H2O2 0.10000E+01
NO 0.00000E+00
NO2 0.10000E+00
NO3 0.10000E+00
N2O5 0.10000E+00
HONO 0.10000E+00
HNO3 0.00000E+00
HNO4 0.00000E+00
NH3 0.00000E+00
SO2 0.00000E+00
SULF 0.00000E+00
CO 0.00000E+00
OH 0.00000E+00
HO2 0.00000E+00
CH4 0.00000E+00
ETH 0.00000E+00
ALKA 0.00000E+00
ALKE 0.00000E+00
BIO 0.00000E+00
ARO 0.00000E+00
HCHO 0.00000E+00
ALD 0.00000E+00
KET 0.00000E+00
CARBO 0.00000E+00
ONIT 0.00000E+00
PAN 0.10000E+00
OP1 0.30000E+00
OP2 0.10000E+00
ORA1 0.00000E+00
ORA2 0.00000E+00
MO2 0.00000E+00
ALKAP 0.00000E+00
ALKEP 0.00000E+00
BIOP 0.00000E+00
PHO 0.00000E+00
ADD 0.00000E+00
AROP 0.00000E+00
CARBOP 0.00000E+00
OLN 0.00000E+00
XO2 0.00000E+00
=====================================================================
*** the following section will be read by ch_init_depconst.F90 ***
=====================================================================
HENRY_SP
Effective Henrys law factor / exponent, See Leriche et al.2013
35
(A32,2E15.5)
O3 1.03000E-02 -0.28300E+04
H2O2 8.44000E+04 -0.76000E+04
NO 1.92000E-03 -0.17900E+04
NO2 1.20000E-02 -0.25160E+04
NO3 3.80000E-02 -0.87070E+04
N2O5 2.10000E+00 -0.34000E+04
HONO 8.38000E+04 -0.31200E+04
HNO3 1.46000E+13 -0.10500E+05
HNO4 4.78000E+04 -0.69000E+04
NH3 3.24000E+03 0.19000E+03
SO2 5.59000E+04 -0.48950E+04
SULF 6.64000E+14 -0.87000E+04
CO 9.81000E-04 -0.17200E+04
OH 3.90000E+01 -0.00000E+00
HO2 3.49000E+04 -0.00000E+00
CH4 1.41000E-03 -0.20400E+04
ETH 1.88000E-03 -0.28750E+05
ALKA 0.15000E-02 -0.32750E+04
ALKE 0.59600E-02 -0.21700E+04
BIO 0.38500E-01 -0.00000E+00
ARO 0.18000E+00 -0.41000E+04
HCHO 3.23000E+03 -0.71960E+04
ALD 0.12900E+02 -0.58900E+04
KET 0.27800E+02 -0.55300E+04
CARBO 0.36000E+06 -0.75450E+04
ONIT 0.10000E+02 -0.59100E+04
PAN 0.28000E+01 -0.57300E+04
OP1 0.30000E+01 -0.52800E+04
OP2 0.33600E+01 -0.59950E+04
ORA1 5.07000E+06 -0.59500E+04
ORA2 2.66000E+05 -0.62000E+04
MO2 0.24500E+01 -0.23200E+04
ALKAP 0.26600E+03 -0.60000E+04
ALKEP 0.83000E+03 -0.60000E+04
BIOP 0.23110E+04 -0.60000E+04
================================================================
*** the following section will be read by build_emisstabn.F90 ***
*** <called by ch_init_emissionn.F90> ***
*** CON: flux is given in molecules/cm2/s ***
================================================================
EMISUNIT
Emission
CON
================================================================
*** the following section will be read by build_pronoslistn.F90 ***
*** <called by ch_init_emissionn.F90> ***
* Format: Nom espece pronostic Meso-NH - Facteur Aggregation - NOM ESPECE EMISE
* expl: NO (espece Meso-NH) 0.70 (facteur aggregation) NOXAN (espece emise PGD)
================================================================
AGREGATION
EMISSIONS MACC-CITY - schema ReLACS
NO 0.70 NOXAN
NO2 0.30 NOXAN
SO2 1.00 SO2AN
NH3 1.00 NH3AN
CO 1.00 COAN
ETH 1.00 C2H6AN
ALKA 0.39 C3H8AN 1.25 BALKAN 0.33 MTOLAN 1.29 BOLAN
ALKE 0.90 C2H4AN 1.18 C3H6AN 1.20 BENEAN
ARO 1.00 AROAN
HCHO 1.00 CH2OAN
ALD 1.00 BALDAN
KET 0.36 ACETAN 2.35 BONEAN
END_AGREGATION
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/006_mesonh_283TO3003/run_mesonh.sh 0 → 100755
View file @ f77924f4
#!/bin/sh
#SBATCH -J run_biomaido
#SBATCH -N 1 # nodes number (=NBP)
#SBATCH -n 64 # CPUs number (on all nodes) (=NBP*TPN)
#SBATCH -o run_megan.eo%j #
#SBATCH -e run_megan.eo%j #
#SBATCH -t 04:30:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
. ~/DEV_57/MNH-PHYEX070-b95d84d7/conf/profile_mesonh-LXifort-R8I4-MNH-V5-6-2-ECRAD140-MPIAUTO-O2
ln -sf ../003_prep_nest/BIOM* .
ln -sf ../004_prep_real_case_d1/ECMWF* .
ln -sf ../005_prep_real_case_d2/ECMWF_D2_20190328_06* .
export MPIRUN="Mpirun -np 64"
set -x
set -e
ln -sf ~rodierq/SAVE/GRIB_KTEST/CHARMEX/CHIMIE_FILES/tuv531/DATA* .
time ${MPIRUN} MESONH${XYZ}
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/PYTHON/plot_BIOMAIDO.py 0 → 100644
View file @ f77924f4
#!/usr/bin/env python3
"""
@author: Quentin Rodier
Creation : 07/01/2021
Last modifications
"""
import matplotlib as mpl
mpl.use('Agg')
import cartopy.crs as ccrs
try:
import MNHPy
from MNHPy.read_MNHfile import read_netcdf
from MNHPy.Panel_Plot import PanelPlot
except:
from read_MNHfile import read_netcdf
from Panel_Plot import PanelPlot
import os
Lpath=["../006_mesonh_283TO3003/"]
LnameFiles=['A2830.1.MEGAD.006.nc']
variables=['CO','O3','NO2','NO','BIOP','HCHO','BIO','latitude','longitude',
'ETH', 'CARBO','PHO','SO2','SULF','DMS','CH4','XO2','PAN',
'DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R',
'SLTM01','SLTM02','SLTM03','SLTM04','SLTM05','SLTM06','SLTM07','SLTM08']
Dvar_input = {'f1': variables}
# Read the variables in the files
Dvar = {}
Dvar = read_netcdf(LnameFiles, Dvar_input, path=Lpath, removeHALO=True, get_data_only=True)
os.system('rm -f tempgraph*')
f=Dvar['f1']
nb_subplot=4
lon = [Dvar['f1']['longitude']]*nb_subplot
lat = [Dvar['f1']['latitude']]*nb_subplot
Lxlab = ['longitude']*nb_subplot
Lylab = ['latitude']*nb_subplot
Lprojection = [ccrs.PlateCarree()]*nb_subplot
Llvl = [0]*nb_subplot
################################################################
######### PANEL 1
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['BIOP'],f['BIO'],f['CO'],f['HCHO']]
Ltitle = ['BIOP','BIO','CO','HCHO']
Lminval=[0, 0, 30, 0.001 ]
Lmaxval=[15.0, 1000, 110, 1 ]
Lstep=[0.25, 20, 2, 0.04 ]
Lstepticks=[1, 100, 10, 0.1 ]
Lfacconv=[1E12,1E12,1E9,1E9 ]
Lcbarlabel = ['ppt','ppt','ppb','ppb']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(1,fig)
################################################################
######### PANEL 2
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['O3'],f['NO2'],f['NO'],f['ETH']]
Ltitle = ['O3','NO2','NO','ETH']
Lminval=[15, 0, 30, 125 ]
Lmaxval=[25.0, 1000, 1000, 175 ]
Lstep=[0.1, 20, 20, 0.25 ]
Lstepticks=[1, 100, 100, 5 ]
Lfacconv=[1E9,1E12,1E12,1E12 ]
Lcbarlabel = ['ppb','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(2,fig)
################################################################
######### PANEL 3
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['CARBO'],f['PHO'],f['SO2'],f['SULF']]
Ltitle = ['CARBO','PHO','SO2','SULF']
Lminval=[0, 0, 50, 0 ]
Lmaxval=[400, 0.7, 1000, 1 ]
Lstep=[5, 0.01, 10, 0.01 ]
Lstepticks=[100, 0.1, 100, 0.1 ]
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(3,fig)
################################################################
######### PANEL 4
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['DMS'],f['CH4'],f['XO2'],f['PAN']]
Ltitle = ['DMS','CH4','XO2','PAN']
Lminval=[0, 0, 0.1, 0 ]
Lmaxval=[300, 0.1, 0.6, 500 ]
Lstep=[2, 0.001, 0.01, 5 ]
Lstepticks=[50, 0.01, 0.1, 100 ]
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(4,fig)
################################################################
######### PANEL 5
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['DEAERM31C'],f['DEAERM32C'],f['DEAERM31R'],f['DEAERM32R']]
Ltitle = ['DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R']
Lminval=[0, 0, 0, 0 ]
Lmaxval=[4, 2, 1, 1 ]
Lstep=[0.005, ]*4
Lstepticks=[0.5 ]*4
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(5,fig)
################################################################
######### PANEL 6
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['SLTM01'],f['SLTM02'],f['SLTM03'],f['SLTM04']]
Ltitle = ['SLTM01','SLTM02','SLTM03','SLTM04']
Lminval=[0]*4
Lmaxval=[10]*4
Lstep=[0.01 ]*4
Lstepticks=[1 ]*4
Lfacconv=[1E18,1E18,1E18,1E18 ]
Lcbarlabel = ['1E-6 ppt','1E-6 ppt','1E-6 ppt','1E-6 ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(6,fig)
################################################################
######### PANEL 7
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['SLTM05'],f['SLTM06'],f['SLTM07'],f['SLTM08']]
Ltitle = ['SLTM05','SLTM06','SLTM07','SLTM08']
Lminval=[0 ]*4
Lmaxval=[1 ]*4
Lstep=[0.01 ]*4
Lstepticks=[0.1 ]*4
Lfacconv=[1E18,1E21,1E21,1E21 ]
Lcbarlabel = ['1E-6 ppt','1E-9 ppt','1E-9 ppt','1E-9 ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(7,fig)
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/PYTHON/plot_BIOMAIDO_DMS.py 0 → 100644
View file @ f77924f4
#!/usr/bin/env python3
"""
@author: Quentin Rodier
Creation : 07/01/2021
Last modifications
"""
import matplotlib as mpl
mpl.use('Agg')
import cartopy.crs as ccrs
try:
import MNHPy
from MNHPy.read_MNHfile import read_netcdf
from MNHPy.Panel_Plot import PanelPlot
except:
from read_MNHfile import read_netcdf
from Panel_Plot import PanelPlot
import os
Lpath=["../006_mesonh_283TO3003/"]
LnameFiles=['A2830.1.MEGAD.006.nc']
variables=['CO','O3','NO2','NO','BIOP','HCHO','BIO','latitude','longitude',
'ETH', 'CARBO','PHO','SO2','SULF','DMS','CH4','XO2','PAN',
'DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R',
'SLTM01','SLTM02','SLTM03','SLTM04','SLTM05','SLTM06','SLTM07','SLTM08']
Dvar_input = {'f1': variables}
# Read the variables in the files
Dvar = {}
Dvar = read_netcdf(LnameFiles, Dvar_input, path=Lpath, removeHALO=True, get_data_only=True)
os.system('rm -f tempgraph*')
f=Dvar['f1']
nb_subplot=4
lon = [Dvar['f1']['longitude']]*nb_subplot
lat = [Dvar['f1']['latitude']]*nb_subplot
Lxlab = ['longitude']*nb_subplot
Lylab = ['latitude']*nb_subplot
Lprojection = [ccrs.PlateCarree()]*nb_subplot
Llvl = [0]*nb_subplot
################################################################
######### PANEL 1
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['BIOP'],f['BIO'],f['CO'],f['HCHO']]
Ltitle = ['BIOP','BIO','CO','HCHO']
Lminval=[0, 0, 30, 0.001 ]
Lmaxval=[15.0, 1000, 110, 1 ]
Lstep=[0.25, 20, 2, 0.04 ]
Lstepticks=[1, 100, 10, 0.1 ]
Lfacconv=[1E12,1E12,1E9,1E9 ]
Lcbarlabel = ['ppt','ppt','ppb','ppb']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(1,fig)
################################################################
######### PANEL 2
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['O3'],f['NO2'],f['NO'],f['ETH']]
Ltitle = ['O3','NO2','NO','ETH']
Lminval=[15, 0, 30, 125 ]
Lmaxval=[25.0, 1000, 1000, 175 ]
Lstep=[0.1, 20, 20, 0.25 ]
Lstepticks=[1, 100, 100, 5 ]
Lfacconv=[1E9,1E12,1E12,1E12 ]
Lcbarlabel = ['ppb','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(2,fig)
################################################################
######### PANEL 3
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['CARBO'],f['PHO'],f['SO2'],f['SULF']]
Ltitle = ['CARBO','PHO','SO2','SULF']
Lminval=[0, 0, 50, 0 ]
Lmaxval=[400, 0.7, 1000, 1 ]
Lstep=[5, 0.01, 10, 0.01 ]
Lstepticks=[100, 0.1, 100, 0.1 ]
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(3,fig)
################################################################
######### PANEL 4
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['DMS'],f['CH4'],f['XO2'],f['PAN']]
Ltitle = ['DMS','CH4','XO2','PAN']
Lminval=[0, 0, 0.1, 0 ]
Lmaxval=[300, 0.1, 0.6, 500 ]
Lstep=[2, 0.001, 0.01, 5 ]
Lstepticks=[50, 0.01, 0.1, 100 ]
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(4,fig)
################################################################
######### PANEL 5
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['DEAERM31C'],f['DEAERM32C'],f['DEAERM31R'],f['DEAERM32R']]
Ltitle = ['DEAERM31C','DEAERM32C','DEAERM31R','DEAERM32R']
Lminval=[0, 0, 0, 0 ]
Lmaxval=[10, 40, 1, 1 ]
Lstep=[0.005, ]*4
Lstepticks=[0.5 ]*4
Lfacconv=[1E12,1E12,1E12,1E12 ]
Lcbarlabel = ['ppt','ppt','ppt','ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(5,fig)
################################################################
######### PANEL 6
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['SLTM01'],f['SLTM02'],f['SLTM03'],f['SLTM04']]
Ltitle = ['SLTM01','SLTM02','SLTM03','SLTM04']
Lminval=[0]*4
Lmaxval=[3]*4
Lstep=[0.01 ]*4
Lstepticks=[1 ]*4
Lfacconv=[1E18,1E18,1E18,1E18 ]
Lcbarlabel = ['1E-6 ppt','1E-6 ppt','1E-6 ppt','1E-6 ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(6,fig)
################################################################
######### PANEL 7
###############################################################
Panel = PanelPlot(2,2, [20,20],'KTEST MEGAN Reunion K=0', minmaxpad=1.05)
Lplot=[f['SLTM05'],f['SLTM06'],f['SLTM07'],f['SLTM08']]
Ltitle = ['SLTM05','SLTM06','SLTM07','SLTM08']
Lminval=[0 ]*4
Lmaxval=[1 ]*4
Lstep=[0.01 ]*4
Lstepticks=[0.1 ]*4
Lfacconv=[1E18,1E21,1E21,1E21 ]
Lcbarlabel = ['1E-6 ppt','1E-9 ppt','1E-9 ppt','1E-9 ppt']*nb_subplot
LaddWhite_cm=[True]*len(Lplot)
Lcolormap = ['gist_rainbow_r']*len(Lplot)
fig = Panel.psectionH(lon=lon, lat=lat, Lvar=Lplot, Llevel=Llvl, Lxlab=Lxlab, Lylab=Lylab, Ltitle=Ltitle, Lminval=Lminval, Lmaxval=Lmaxval, Lcarte=[54.5,56,-22,-20],
Lstep=Lstep, Lstepticks=Lstepticks, Lcolormap=Lcolormap, Lcbarlabel=Lcbarlabel, Lproj=Lprojection, Lfacconv=Lfacconv, LaddWhite_cm=LaddWhite_cm)
Panel.save_graph(7,fig)
MY_RUN/INTEGRATION_CASES/HPC/BIOMAIDO_DMS/PYTHON/run_python 0 → 100755
View file @ f77924f4
#!/bin/sh
#SBATCH -J tools
#SBATCH -N 1 # nodes number
#SBATCH -n 1 # CPUs number (on all nodes)
#SBATCH -o tools.eo%j #
#SBATCH -e tools.eo%j #
#SBATCH -t 01:00:00 # time limit
# Echo des commandes
ulimit -c 0
ulimit -s unlimited
# Arrete du job des la premiere erreur
set -e
. ~rodierq/DEV_57/MNH-git_open_source-lfs-0a57991/conf/profile_mesonh-LXifort-R8I4-MNH-V5-7-0-ECRAD140-MPIAUTO-O2
ln -sf ${SRC_MESONH}/src/LIB/Python/* .
module purge
module load python/3.7.6
python3 plot_BIOMAIDO_DMS.py
convert *.png BIOMAIDO_DMS.pdf
MY_RUN/INTEGRATION_CASES/HPC/BLAZE_OROG/01_prep_ideal_case/ConstantSlope.py 0 → 100644
View file @ f77924f4
# -*- coding: utf-8 -*-
__author__ = "Aurélien Costes"
__copyright__ = "copyright 2019, Météo-France Cerfacs"
__credits__ = ["Aurélien Costes"]
__license__ = "CeCILL-C"
__version__ = "0.0.1"
__maintainer__ = "Aurélien Costes"
__email__ = "aurelien.costes@meteo.fr"
__status__ = "Development"
import sys, os
import argparse
import numpy as np
import matplotlib.pyplot as plt
from math import radians, sqrt, pow, ceil, exp, tan
from math import atan, degrees
import csv
# parameters
showfig = False
Savecsv = True
Nx = 200
Ny = 30
dx = 25
dy = 25
minheight = 0
maxheight = 450
# Break slopePoints
# upslope start, downslope start
Offset = [500, 2500]
minzeroslopeend = 500
# up slope angle [deg]
Upslopeangle = 30
# down slope angle [deg]
Downslopeangle = 30
x = np.linspace(dx / 2, Nx * dx - dx / 2, Nx)
y = np.linspace(dy / 2, Ny * dy - dy / 2, Ny)
X, Y = np.meshgrid(x, y)
hmap = np.zeros((Ny, Nx))
upsloperad = radians(Upslopeangle)
downsloperad = radians(Downslopeangle)
z = minheight
for i in range(Nx):
if x[i] <= Offset[0]:
hmap[:, i] = z
continue
if x[i] >= np.max(x) - minzeroslopeend:
hmap[:, i] = z
continue
# downslope
if x[i] >= min(Offset[1], np.max(x) - minzeroslopeend - z / tan(downsloperad)):
z = max(minheight, z - tan(downsloperad) * dx)
hmap[:, i] = z
continue
# upslope
z = min(maxheight, z + tan(upsloperad) * dx)
hmap[:, i] = z
if showfig:
fig = plt.figure()
plt.plot(x, hmap[0, :], ".-", markersize=3)
plt.show()
if Savecsv:
with open("CstSlope.csv", "w") as f:
writer = csv.writer(f)
for j in range(Ny):
row = hmap[j, :]
writer.writerow(row)
MY_RUN/INTEGRATION_CASES/HPC/BLAZE_OROG/01_prep_ideal_case/CstSlope.csv 0 → 100644
View file @ f77924f4
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,14.433756729740644,28.867513459481287,43.30127018922193,57.735026918962575,72.16878364870323,86.60254037844388,101.03629710818453,115.47005383792518,129.90381056766583,144.33756729740648,158.77132402714713,173.20508075688778,187.63883748662843,202.07259421636908,216.50635094610973,230.94010767585038,245.37386440559104,259.80762113533166,274.2413778650723,288.67513459481296,303.1088913245536,317.54264805429426,331.9764047840349,346.41016151377556,360.8439182435162,375.27767497325686,389.7114317029975,404.14518843273817,418.5789451624788,433.01270189221947,447.4464586219601,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,450.0,435.56624327025935,421.1324865405187,406.69872981077805,392.2649730810374,377.83121635129675,363.3974596215561,348.96370289181544,334.5299461620748,320.09618943233414,305.6624327025935,291.22867597285284,276.7949192431122,262.36116251337154,247.9274057836309,233.49364905389024,219.0598923241496,204.62613559440894,190.19237886466829,175.75862213492763,161.32486540518698,146.89110867544633,132.45735194570568,118.02359521596503,103.58983848622438,89.15608175648373,74.72232502674308,60.288568297002435,45.85481156726179,31.421054837521147,16.987298107780504,2.55354137803986,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0

Documentation | Mentions légales | CGU | Accessibilité : non conforme