!$POLYPHEMUS/bin/download_dataset preprocessing/mozart4
### 1/2 Update the local cache... cd /libre_RAID1/dores/data/cache/preprocessing_mozart4 sending incremental file list sent 207 bytes received 14 bytes 442.00 bytes/sec total size is 601,025,916 speedup is 2,719,574.28 None ### 2/2 Use the local cache to create the requested files... cd /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example sending incremental file list ./ sent 214 bytes received 21 bytes 470.00 bytes/sec total size is 601,025,916 speedup is 2,557,557.09 None
from test_utils.notebook.env import *
Polyphemus is at /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1 Atmopy is at /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/include/atmopy This notebook is at /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example
!cd $POLYPHEMUS/preprocessing/mozart4/example; ls . ./input ./raw_data/mozart4/
.: data input mozart4.cfg raw_data example_mozart4.ipynb levels.dat mozart4-ic-beta.cfg general.cfg mozart4-bc-beta.cfg polyphemus ./input: molecular_weights_mozart4_aer.txt speciation_mozart4_aer.dat molecular_weights_mozart4_du.txt speciation_mozart4_du.dat molecular_weights_mozart4_gas_sorted.txt speciation_mozart4_gas.dat molecular_weights_mozart4_gas.txt speciation_mozart4_ss.dat molecular_weights_mozart4_ss.txt ./raw_data/mozart4/: mozart4geos5-20081231-20090101.nc mozart4-Gravelines-201308.nc
!cd $POLYPHEMUS/preprocessing/mozart4/ && scons -j 8
assert _exit_code == 0
scons: Reading SConscript files ... Checking for C library blas... (cached) yes Checking for C library atlas... (cached) no Checking for C library lapack... (cached) yes Checking for C library g2c... (cached) no Checking for C library gfortran... (cached) yes Checking for C library gslcblas... (cached) yes Checking for C library blitz... (cached) yes Checking for C library newran... (cached) no Checking for C library netcdf_c++... (cached) yes Checking for C library netcdf... (cached) yes Checking for C library netcdf_c++... (cached) yes Checking for C library ifcore... (cached) no scons: done reading SConscript files. scons: Building targets ... scons: `.' is up to date. scons: done building targets.
for chem in "gas", "aer", "DU", "SS":
result = "data/ic-" + chem
!rm -fr $result
!mkdir -p $result
!../ic-mozart4-$chem general.cfg mozart4.cfg 20081231_06:00
assert _exit_code == 0
assert_no_difference(compare_dir(result))
==> GET BOUNDARY CONDITIONS FOR GAS <===
Get MOZART4-GEOS5 data for 2008-12-31_06
Read configuration... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute IC for gas species...
<< From MOZART4 : Nt = 12 Nz = 57 Ny = 21 Nx = 21
>> to POLAIR : Nt = 1 Nz = 9 Ny = 50 Nx = 65
AACD ... compute with :
\=> CH3COOH (1) : read, arrange, convert (mw = 60), interpolate
(min: 0.00011326) (mean: 0.270595) (max: 4.89675)
ALD2 ... compute with :
\=> CH3CHO (1) : read, arrange, convert (mw = 44), interpolate
\=> GLYALD (1) : read, arrange, convert (mw = 60), interpolate
(min: 0.00792833) (mean: 0.722762) (max: 8.91801)
CH4 ... compute with :
\=> CH4 (1) : read, arrange, convert (mw = 16), interpolate
(min: 390.425) (mean: 1036.99) (max: 1373.32)
CO ... compute with :
\=> CO (1) : read, arrange, convert (mw = 28), interpolate
(min: 28.9818) (mean: 160.304) (max: 854.075)
CRES ... compute with :
\=> CRESOL (1) : read, arrange, convert (mw = 108), interpolate
(min: 0.00715653) (mean: 0.398086) (max: 3.7472)
ETH ... compute with :
\=> C2H4 (1) : read, arrange, convert (mw = 28), interpolate
(min: 0.000438194) (mean: 0.388121) (max: 4.691)
ETHA ... compute with :
\=> C2H6 (1) : read, arrange, convert (mw = 30), interpolate
(min: 0.188036) (mean: 1.57391) (max: 9.6277)
ETOH ... compute with :
\=> C2H5OH (1) : read, arrange, convert (mw = 46), interpolate
(min: 0.00146884) (mean: 0.287846) (max: 3.98339)
FORM ... compute with :
\=> CH2O (1) : read, arrange, convert (mw = 30), interpolate
(min: 0.0155643) (mean: 0.572379) (max: 9.98827)
H2O2 ... compute with :
\=> H2O2 (1) : read, arrange, convert (mw = 34), interpolate
(min: 1.33389e-11) (mean: 0.179746) (max: 0.89616)
HNO3 ... compute with :
\=> HNO3 (1) : read, arrange, convert (mw = 63), interpolate
(min: 0) (mean: 1.36509) (max: 18.8538)
HNO4 ... compute with :
\=> HO2NO2 (1) : read, arrange, convert (mw = 79), interpolate
(min: 1.19558e-05) (mean: 0.0442625) (max: 1.27807)
HO2 ... compute with :
\=> HO2 (1) : read, arrange, convert (mw = 33), interpolate
(min: 2.18877e-06) (mean: 0.00031088) (max: 0.0018179)
ISON ... compute with :
\=> ISOPNO3 (1) : read, arrange, convert (mw = 162), interpolate
(min: 2.20611e-13) (mean: 0.000194206) (max: 0.0050937)
ISOP ... compute with :
\=> ISOP (1) : read, arrange, convert (mw = 68), interpolate
(min: 3.8228e-11) (mean: 0.0030632) (max: 0.0576258)
ISPD ... compute with :
\=> ISOPNO3 (1) : read, arrange, convert (mw = 162), interpolate
\=> MACR (1) : read, arrange, convert (mw = 70), interpolate
\=> MPAN (1) : read, arrange, convert (mw = 147), interpolate
\=> MVK (1) : read, arrange, convert (mw = 70), interpolate
(min: 0.00272995) (mean: 0.25758) (max: 2.0011)
MACR ... compute with :
\=> MACR (1) : read, arrange, convert (mw = 70), interpolate
(min: 2.60573e-05) (mean: 0.0690642) (max: 0.799585)
MEO2 ... compute with :
\=> CH3O2 (1) : read, arrange, convert (mw = 47), interpolate
(min: 7.30318e-06) (mean: 0.00159206) (max: 0.00840628)
MEOH ... compute with :
\=> CH3OH (1) : read, arrange, convert (mw = 32), interpolate
(min: 0.105701) (mean: 0.80641) (max: 5.90559)
MEPX ... compute with :
\=> CH3OOH (1) : read, arrange, convert (mw = 48), interpolate
(min: 0.0302843) (mean: 0.250798) (max: 0.852966)
MGLY ... compute with :
\=> CH3COCHO (1) : read, arrange, convert (mw = 72), interpolate
(min: 3.41319e-06) (mean: 0.0617791) (max: 1.1893)
MPAN ... compute with :
\=> MPAN (1) : read, arrange, convert (mw = 147), interpolate
(min: 0.00164876) (mean: 0.111447) (max: 0.470638)
N2O5 ... compute with :
\=> N2O5 (1) : read, arrange, convert (mw = 108), interpolate
(min: 2.25178e-06) (mean: 0.0313484) (max: 0.375367)
NH3 ... compute with :
\=> NH3 (1) : read, arrange, convert (mw = 17), interpolate
(min: 0) (mean: 0.574054) (max: 32.876)
NO ... compute with :
\=> NO (1) : read, arrange, convert (mw = 30), interpolate
(min: 3.95304e-08) (mean: 0.780851) (max: 75.7862)
NO2 ... compute with :
\=> NO2 (1) : read, arrange, convert (mw = 46), interpolate
(min: 0.00706546) (mean: 4.339) (max: 37.4535)
NO3 ... compute with :
\=> NO3 (1) : read, arrange, convert (mw = 62), interpolate
(min: 3.15594e-08) (mean: 0.00828678) (max: 0.070198)
NTR ... compute with :
\=> ONIT (1) : read, arrange, convert (mw = 119), interpolate
\=> ONITR (1) : read, arrange, convert (mw = 147), interpolate
(min: 0.000299448) (mean: 0.0800173) (max: 1.41813)
O3 ... compute with :
\=> O3 (1) : read, arrange, convert (mw = 48), interpolate
(min: 0.00351451) (mean: 56.7972) (max: 90.6539)
OH ... compute with :
\=> OH (1) : read, arrange, convert (mw = 17), interpolate
(min: 6.29982e-09) (mean: 1.77965e-06) (max: 4.17128e-05)
OLE ... compute with :
\=> C3H6 (1) : read, arrange, convert (mw = 42), interpolate
\=> BIGENE (1) : read, arrange, convert (mw = 56), interpolate
(min: 8.01105e-06) (mean: 0.191584) (max: 2.21856)
OPEN ... compute with :
\=> BIGALD (1) : read, arrange, convert (mw = 98), interpolate
(min: 3.05854e-05) (mean: 0.102489) (max: 5.14947)
PACD ... compute with :
\=> CH3COOOH (1) : read, arrange, convert (mw = 76), interpolate
(min: 0.031849) (mean: 0.18121) (max: 0.442409)
PAN ... compute with :
\=> PAN (1) : read, arrange, convert (mw = 121), interpolate
(min: 0.0937012) (mean: 2.16411) (max: 14.3627)
PAR ... compute with :
\=> C3H6 (1) : read, arrange, convert (mw = 42), interpolate
\=> C3H8 (1.5) : read, arrange, convert (mw = 44), interpolate
\=> CH3COCH3 (3) : read, arrange, convert (mw = 58), interpolate
\=> BIGALK (3) : read, arrange, convert (mw = 72), interpolate
\=> BIGENE (1) : read, arrange, convert (mw = 56), interpolate
\=> DMS (2) : read, arrange, convert (mw = 62), interpolate
\=> MEK (4) : read, arrange, convert (mw = 72), interpolate
\=> HYAC (2) : read, arrange, convert (mw = 74), interpolate
(min: 1.63329) (mean: 34.6132) (max: 310.3)
SO2 ... compute with :
\=> SO2 (1) : read, arrange, convert (mw = 64), interpolate
(min: 1.52017e-11) (mean: 5.14065) (max: 114.137)
TERP ... compute with :
\=> C10H16 (1) : read, arrange, convert (mw = 136), interpolate
(min: 5.48538e-21) (mean: 0.0222025) (max: 0.562507)
TOL ... compute with :
\=> TOLUENE (1) : read, arrange, convert (mw = 92), interpolate
(min: 0.00306372) (mean: 1.90691) (max: 24.4759)
==============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/ic-gas
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/ic-gas
==> GET BOUNDARY CONDITIONS FOR AER <===
Get MOZART4-GEOS5 data for 2008-12-31_06
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute IC for aerosols species...
<< From MOZART4 : 8 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
CB1 ... read ... arrange ... convert (mw = 12) ... interpolate
CB2 ... read ... arrange ... convert (mw = 12) ... interpolate
OC1 ... read ... arrange ... convert (mw = 12) ... interpolate
OC2 ... read ... arrange ... convert (mw = 12) ... interpolate
SO4 ... read ... arrange ... convert (mw = 96) ... interpolate
NH4NO3 ... read ... arrange ... convert (mw = 80) ... interpolate
SOA ... read ... arrange ... convert (mw = 144) ... interpolate
NH4 ... read ... arrange ... convert (mw = 18) ... interpolate
>> to POLAIR : Nt = 1 Nz = 9 Ny = 50 Nx = 65
PBC_0 ... compute with : CB1 (1) CB2 (1) ... write
(min: 0.000258914) (mean: 0.0301521) (max: 0.427784)
PBC_1 ... compute with : CB1 (1) CB2 (1) ... write
(min: 0.00184202) (mean: 0.214515) (max: 3.04343)
PBC_2 ... compute with : CB1 (1) CB2 (1) ... write
(min: 0.000697439) (mean: 0.0812217) (max: 1.15233)
PBC_3 ... compute with : CB1 (1) CB2 (1) ... write
(min: 1.09362e-05) (mean: 0.0012736) (max: 0.018069)
PBC_4 ... compute with : CB1 (1) CB2 (1) ... write
(min: 4.56628e-09) (mean: 5.31771e-07) (max: 7.54451e-06)
PBiA2D_0 ... compute with : SOA (1) ... write
(min: 4.70411e-05) (mean: 0.000236434) (max: 0.00118173)
PBiA2D_1 ... compute with : SOA (1) ... write
(min: 0.00154276) (mean: 0.00775412) (max: 0.0387559)
PBiA2D_2 ... compute with : SOA (1) ... write
(min: 0.00397973) (mean: 0.0200026) (max: 0.0999754)
PBiA2D_3 ... compute with : SOA (1) ... write
(min: 0.000926708) (mean: 0.00465775) (max: 0.02328)
PBiA2D_4 ... compute with : SOA (1) ... write
(min: 1.60469e-05) (mean: 8.06536e-05) (max: 0.000403116)
PNH4_0 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
(min: 8.47866e-07) (mean: 0.000196523) (max: 0.00322886)
PNH4_1 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
(min: 0.000443503) (mean: 0.102798) (max: 1.68896)
PNH4_2 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
(min: 0.00809912) (mean: 1.87731) (max: 30.8432)
PNH4_3 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
(min: 0.00795601) (mean: 1.84415) (max: 30.2982)
PNH4_4 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
(min: 0.000418524) (mean: 0.0970085) (max: 1.59383)
PNO3_0 ... compute with : NH4NO3 (0.775) ... write
(min: 0) (mean: 0.000293274) (max: 0.00515794)
PNO3_1 ... compute with : NH4NO3 (0.775) ... write
(min: 0) (mean: 0.153406) (max: 2.69802)
PNO3_2 ... compute with : NH4NO3 (0.775) ... write
(min: 0) (mean: 2.80146) (max: 49.2705)
PNO3_3 ... compute with : NH4NO3 (0.775) ... write
(min: 0) (mean: 2.75194) (max: 48.3999)
PNO3_4 ... compute with : NH4NO3 (0.775) ... write
(min: 0) (mean: 0.144767) (max: 2.54606)
PSO4_0 ... compute with : SO4 (1) ... write
(min: 3.69323e-06) (mean: 7.52836e-05) (max: 0.000623806)
PSO4_1 ... compute with : SO4 (1) ... write
(min: 0.00193186) (mean: 0.0393794) (max: 0.326301)
PSO4_2 ... compute with : SO4 (1) ... write
(min: 0.035279) (mean: 0.719135) (max: 5.95881)
PSO4_3 ... compute with : SO4 (1) ... write
(min: 0.0346556) (mean: 0.70643) (max: 5.85353)
PSO4_4 ... compute with : SO4 (1) ... write
(min: 0.00182305) (mean: 0.0371617) (max: 0.307923)
PSOAhP_0 ... compute with : OC1 (0.79) OC2 (0.79) ... write
(min: 3.87342e-05) (mean: 0.00187562) (max: 0.0288132)
PSOAhP_1 ... compute with : OC1 (0.79) OC2 (0.79) ... write
(min: 0.00127033) (mean: 0.0615131) (max: 0.944957)
PSOAhP_2 ... compute with : OC1 (0.79) OC2 (0.79) ... write
(min: 0.00327695) (mean: 0.15868) (max: 2.43763)
PSOAhP_3 ... compute with : OC1 (0.79) OC2 (0.79) ... write
(min: 0.000763062) (mean: 0.0369502) (max: 0.567619)
PSOAhP_4 ... compute with : OC1 (0.79) OC2 (0.79) ... write
(min: 1.32132e-05) (mean: 0.000639823) (max: 0.00982889)
PSOAlP_0 ... compute with : OC1 (0.46) OC2 (0.46) ... write
(min: 2.25541e-05) (mean: 0.00109214) (max: 0.0167773)
PSOAlP_1 ... compute with : OC1 (0.46) OC2 (0.46) ... write
(min: 0.000739683) (mean: 0.0358177) (max: 0.550228)
PSOAlP_2 ... compute with : OC1 (0.46) OC2 (0.46) ... write
(min: 0.0019081) (mean: 0.0923954) (max: 1.41938)
PSOAlP_3 ... compute with : OC1 (0.46) OC2 (0.46) ... write
(min: 0.000444314) (mean: 0.0215152) (max: 0.330512)
PSOAlP_4 ... compute with : OC1 (0.46) OC2 (0.46) ... write
(min: 7.69375e-06) (mean: 0.000372555) (max: 0.00572315)
PSOAmP_0 ... compute with : OC1 (0.58) OC2 (0.58) ... write
(min: 2.84378e-05) (mean: 0.00137705) (max: 0.021154)
PSOAmP_1 ... compute with : OC1 (0.58) OC2 (0.58) ... write
(min: 0.000932644) (mean: 0.0451615) (max: 0.693766)
PSOAmP_2 ... compute with : OC1 (0.58) OC2 (0.58) ... write
(min: 0.00240586) (mean: 0.116499) (max: 1.78965)
PSOAmP_3 ... compute with : OC1 (0.58) OC2 (0.58) ... write
(min: 0.000560223) (mean: 0.0271276) (max: 0.416733)
PSOAmP_4 ... compute with : OC1 (0.58) OC2 (0.58) ... write
(min: 9.70081e-06) (mean: 0.00046974) (max: 0.00721614)
==============>>> DONE <<<===============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/ic-aer
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/ic-aer
==> GET BOUNDARY CONDITIONS FOR DU <===
Get MOZART4-GEOS5 data for 2008-12-31_06
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute IC for dust species...
<< From MOZART4 : 4 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
DUST1 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST2 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST3 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST4 ... read ... arrange ... convert (mw = 34) ... interpolate
>> to POLAIR : Nt = 1 Nz = 9 Ny = 50 Nx = 65
PMD_0 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
(min: 0) (mean: 0) (max: 0)
PMD_1 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
(min: 0.0284213) (mean: 0.322203) (max: 3.88441)
PMD_2 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
(min: 0.135305) (mean: 1.64307) (max: 20.2395)
PMD_3 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
(min: 0.26607) (mean: 4.25109) (max: 59.5587)
PMD_4 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
(min: 0.00412765) (mean: 0.335116) (max: 10.0888)
=============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/ic-DU
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/ic-DU
==> GET BOUNDARY CONDITIONS FOR SS <===
Get MOZART4-GEOS5 data for 2008-12-31_06
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute IC for sea-salt species...
<< From MOZART4 : 4 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
SA1 ... read ... arrange ... convert (mw = 58) ... interpolate
SA2 ... read ... arrange ... convert (mw = 58) ... interpolate
SA3 ... read ... arrange ... convert (mw = 58) ... interpolate
SA4 ... read ... arrange ... convert (mw = 58) ... interpolate
>> to POLAIR : Nt = 1 Nz = 9 Ny = 50 Nx = 65
PHCL_0 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
(min: 0) (mean: 0) (max: 0)
PHCL_1 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
(min: 0.000136893) (mean: 0.00543397) (max: 0.092896)
PHCL_2 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
(min: 0.00313707) (mean: 0.0598886) (max: 1.42853)
PHCL_3 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
(min: 0.000135435) (mean: 0.28984) (max: 10.1431)
PHCL_4 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
(min: 3.97351e-06) (mean: 0.310872) (max: 14.2021)
PNA_0 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
(min: 0) (mean: 0) (max: 0)
PNA_1 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
(min: 9.12619e-05) (mean: 0.00362265) (max: 0.0619306)
PNA_2 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
(min: 0.00209138) (mean: 0.0399262) (max: 0.952352)
PNA_3 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
(min: 9.02897e-05) (mean: 0.193227) (max: 6.76204)
PNA_4 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
(min: 2.64901e-06) (mean: 0.207251) (max: 9.46809)
=============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/ic-SS
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/ic-SS
for chem in "gas", "aer", "DU", "SS":
result = "data/bc-" + chem
!rm -fr $result
!mkdir -p $result
!../bc-mozart4-$chem general.cfg mozart4.cfg 20081231_06:00 20090103_00:00|tee bc.log
assert _exit_code == 0
assert_no_difference(compare_dir(result))
==> GET BOUNDARY CONDITIONS FOR GAS <===
Get MOZART4-GEOS5 data from 2008-12-31_06 to 2009-01-03_00
Read configuration... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute BC for gas species...
<< From MOZART4 : Nt = 12 Nz = 57 Ny = 21 Nx = 21
>> to POLAIR : Nt = 12 Nz = 9 Ny = 50 Nx = 65
AACD ... compute with :
\=> CH3COOH (1) : read, arrange, convert (mw = 60), interpolate
in x direction (min: 0.000127491) (mean: 0.180436) (max: 2.84054)
in y direction (min: 0.000164181) (mean: 0.131343) (max: 0.910704)
in z direction (min: 0.000402213) (mean: 0.0036205) (max: 0.0234307)
ALD2 ... compute with :
\=> CH3CHO (1) : read, arrange, convert (mw = 44), interpolate
\=> GLYALD (1) : read, arrange, convert (mw = 60), interpolate
in x direction (min: 0.00744348) (mean: 0.533821) (max: 5.171)
in y direction (min: 0.0083073) (mean: 0.333037) (max: 1.46714)
in z direction (min: 0.00975393) (mean: 0.0329469) (max: 0.0701719)
CH4 ... compute with :
\=> CH4 (1) : read, arrange, convert (mw = 16), interpolate
in x direction (min: 420.468) (mean: 1044.1) (max: 1377)
in y direction (min: 427.966) (mean: 1038.4) (max: 1370.14)
in z direction (min: 480.823) (mean: 565.399) (max: 596.948)
CO ... compute with :
\=> CO (1) : read, arrange, convert (mw = 28), interpolate
in x direction (min: 29.3154) (mean: 144.804) (max: 533.665)
in y direction (min: 29.3609) (mean: 126.579) (max: 273.404)
in z direction (min: 41.3955) (mean: 52.6854) (max: 63.0493)
CRES ... compute with :
\=> CRESOL (1) : read, arrange, convert (mw = 108), interpolate
in x direction (min: 0.00578785) (mean: 0.327022) (max: 2.34685)
in y direction (min: 0.00624516) (mean: 0.2179) (max: 0.850514)
in z direction (min: 0.005375) (mean: 0.0412457) (max: 0.0907278)
ETH ... compute with :
\=> C2H4 (1) : read, arrange, convert (mw = 28), interpolate
in x direction (min: 0.000256952) (mean: 0.286541) (max: 3.89275)
in y direction (min: 0.000313174) (mean: 0.198079) (max: 1.57532)
in z direction (min: 0.000149483) (mean: 0.00518636) (max: 0.0172028)
ETHA ... compute with :
\=> C2H6 (1) : read, arrange, convert (mw = 30), interpolate
in x direction (min: 0.185683) (mean: 1.33649) (max: 6.76523)
in y direction (min: 0.189107) (mean: 1.11569) (max: 3.19501)
in z direction (min: 0.228478) (mean: 0.331105) (max: 0.448982)
ETOH ... compute with :
\=> C2H5OH (1) : read, arrange, convert (mw = 46), interpolate
in x direction (min: 0.00095825) (mean: 0.180564) (max: 2.11188)
in y direction (min: 0.00118274) (mean: 0.127573) (max: 0.913957)
in z direction (min: 0.000708625) (mean: 0.0058014) (max: 0.0150276)
FORM ... compute with :
\=> CH2O (1) : read, arrange, convert (mw = 30), interpolate
in x direction (min: 0.0112904) (mean: 0.395759) (max: 4.99235)
in y direction (min: 0.0114844) (mean: 0.305) (max: 1.30055)
in z direction (min: 0.0184515) (mean: 0.0364586) (max: 0.0770876)
H2O2 ... compute with :
\=> H2O2 (1) : read, arrange, convert (mw = 34), interpolate
in x direction (min: 1.13423e-11) (mean: 0.155071) (max: 0.803926)
in y direction (min: 1.39624e-11) (mean: 0.206952) (max: 1.04394)
in z direction (min: 0.0203781) (mean: 0.100002) (max: 0.341412)
HNO3 ... compute with :
\=> HNO3 (1) : read, arrange, convert (mw = 63), interpolate
in x direction (min: 0) (mean: 1.10664) (max: 8.50741)
in y direction (min: 3.76652e-16) (mean: 0.80165) (max: 9.59732)
in z direction (min: 0.00039789) (mean: 0.0310657) (max: 0.0960317)
HNO4 ... compute with :
\=> HO2NO2 (1) : read, arrange, convert (mw = 79), interpolate
in x direction (min: 1.34048e-05) (mean: 0.0371056) (max: 1.15283)
in y direction (min: 1.32752e-05) (mean: 0.0139246) (max: 0.323555)
in z direction (min: 0.000194699) (mean: 0.0221041) (max: 0.0692699)
HO2 ... compute with :
\=> HO2 (1) : read, arrange, convert (mw = 33), interpolate
in x direction (min: 1.63972e-06) (mean: 0.00109551) (max: 0.0141063)
in y direction (min: 1.83617e-06) (mean: 0.00148255) (max: 0.0145221)
in z direction (min: 4.33558e-06) (mean: 0.000887039) (max: 0.00709991)
ISON ... compute with :
\=> ISOPNO3 (1) : read, arrange, convert (mw = 162), interpolate
in x direction (min: 2.98109e-18) (mean: 0.000377292) (max: 0.0247055)
in y direction (min: 3.202e-18) (mean: 0.00061921) (max: 0.0373332)
in z direction (min: 5.54076e-16) (mean: 3.67882e-07) (max: 0.000121314)
ISOP ... compute with :
\=> ISOP (1) : read, arrange, convert (mw = 68), interpolate
in x direction (min: 1.69266e-12) (mean: 0.00694344) (max: 0.343246)
in y direction (min: 1.77613e-12) (mean: 0.00959333) (max: 0.266161)
in z direction (min: 2.05033e-10) (mean: 1.02054e-05) (max: 0.000746733)
ISPD ... compute with :
\=> ISOPNO3 (1) : read, arrange, convert (mw = 162), interpolate
\=> MACR (1) : read, arrange, convert (mw = 70), interpolate
\=> MPAN (1) : read, arrange, convert (mw = 147), interpolate
\=> MVK (1) : read, arrange, convert (mw = 70), interpolate
in x direction (min: 0.00216782) (mean: 0.211557) (max: 1.32316)
in y direction (min: 0.0019952) (mean: 0.169302) (max: 1.10928)
in z direction (min: 0.00516021) (mean: 0.0387227) (max: 0.0936107)
MACR ... compute with :
\=> MACR (1) : read, arrange, convert (mw = 70), interpolate
in x direction (min: 2.2328e-05) (mean: 0.057143) (max: 0.513574)
in y direction (min: 3.05539e-05) (mean: 0.0485049) (max: 0.41352)
in z direction (min: 1.22238e-05) (mean: 0.000633236) (max: 0.00292725)
MEO2 ... compute with :
\=> CH3O2 (1) : read, arrange, convert (mw = 47), interpolate
in x direction (min: 1.30928e-05) (mean: 0.00156976) (max: 0.0142795)
in y direction (min: 1.05811e-05) (mean: 0.00191982) (max: 0.0147608)
in z direction (min: 7.33052e-05) (mean: 0.000635475) (max: 0.00293495)
MEOH ... compute with :
\=> CH3OH (1) : read, arrange, convert (mw = 32), interpolate
in x direction (min: 0.0900427) (mean: 0.580819) (max: 3.89176)
in y direction (min: 0.102747) (mean: 0.464373) (max: 2.22518)
in z direction (min: 0.118977) (mean: 0.177731) (max: 0.301138)
MEPX ... compute with :
\=> CH3OOH (1) : read, arrange, convert (mw = 48), interpolate
in x direction (min: 0.0240251) (mean: 0.24716) (max: 0.878184)
in y direction (min: 0.0320893) (mean: 0.273586) (max: 0.907751)
in z direction (min: 0.0668912) (mean: 0.183146) (max: 0.448084)
MGLY ... compute with :
\=> CH3COCHO (1) : read, arrange, convert (mw = 72), interpolate
in x direction (min: 8.6753e-06) (mean: 0.0398744) (max: 0.653322)
in y direction (min: 1.62737e-05) (mean: 0.020264) (max: 0.126049)
in z direction (min: 4.27464e-06) (mean: 0.000537604) (max: 0.0018139)
MPAN ... compute with :
\=> MPAN (1) : read, arrange, convert (mw = 147), interpolate
in x direction (min: 0.000744089) (mean: 0.0939158) (max: 0.397485)
in y direction (min: 0.00085792) (mean: 0.0742876) (max: 0.317416)
in z direction (min: 0.00429828) (mean: 0.037321) (max: 0.0921996)
N2O5 ... compute with :
\=> N2O5 (1) : read, arrange, convert (mw = 108), interpolate
in x direction (min: 5.36046e-07) (mean: 0.0569153) (max: 0.970275)
in y direction (min: 4.89771e-07) (mean: 0.0404569) (max: 0.503007)
in z direction (min: 0.000102728) (mean: 0.00246071) (max: 0.00830516)
NH3 ... compute with :
\=> NH3 (1) : read, arrange, convert (mw = 17), interpolate
in x direction (min: 0) (mean: 0.226885) (max: 8.56279)
in y direction (min: 0) (mean: 0.0646894) (max: 1.51748)
in z direction (min: 0) (mean: 1.31342e-05) (max: 0.00149458)
NO ... compute with :
\=> NO (1) : read, arrange, convert (mw = 30), interpolate
in x direction (min: 1.76325e-08) (mean: 0.223204) (max: 17.247)
in y direction (min: 1.84778e-08) (mean: 0.0961627) (max: 7.86189)
in z direction (min: 4.71254e-08) (mean: 0.00258312) (max: 0.0297092)
NO2 ... compute with :
\=> NO2 (1) : read, arrange, convert (mw = 46), interpolate
in x direction (min: 0.00258976) (mean: 3.04491) (max: 36.6905)
in y direction (min: 0.00162135) (mean: 1.5734) (max: 26.6557)
in z direction (min: 0.00301395) (mean: 0.0188875) (max: 0.0784962)
NO3 ... compute with :
\=> NO3 (1) : read, arrange, convert (mw = 62), interpolate
in x direction (min: 2.28058e-07) (mean: 0.00680625) (max: 0.138228)
in y direction (min: 2.69778e-07) (mean: 0.00689209) (max: 0.0747368)
in z direction (min: 2.63496e-07) (mean: 0.000119526) (max: 0.000296482)
NTR ... compute with :
\=> ONIT (1) : read, arrange, convert (mw = 119), interpolate
\=> ONITR (1) : read, arrange, convert (mw = 147), interpolate
in x direction (min: 0.000116131) (mean: 0.049911) (max: 0.53825)
in y direction (min: 0.000111241) (mean: 0.0395198) (max: 0.351498)
in z direction (min: 9.81932e-05) (mean: 0.00216205) (max: 0.00846856)
O3 ... compute with :
\=> O3 (1) : read, arrange, convert (mw = 48), interpolate
in x direction (min: 0.0610762) (mean: 61.6215) (max: 103.284)
in y direction (min: 8.08429) (mean: 63.7825) (max: 103.543)
in z direction (min: 36.3972) (mean: 47.6833) (max: 63.5759)
OH ... compute with :
\=> OH (1) : read, arrange, convert (mw = 17), interpolate
in x direction (min: 4.72197e-09) (mean: 5.94987e-06) (max: 0.000125959)
in y direction (min: 3.76897e-09) (mean: 8.33374e-06) (max: 0.000131175)
in z direction (min: 5.13033e-09) (mean: 4.82977e-06) (max: 6.51778e-05)
OLE ... compute with :
\=> C3H6 (1) : read, arrange, convert (mw = 42), interpolate
\=> BIGENE (1) : read, arrange, convert (mw = 56), interpolate
in x direction (min: 2.41043e-06) (mean: 0.131526) (max: 2.29814)
in y direction (min: 4.82226e-06) (mean: 0.136267) (max: 1.42217)
in z direction (min: 2.50549e-06) (mean: 0.000292588) (max: 0.00228821)
OPEN ... compute with :
\=> BIGALD (1) : read, arrange, convert (mw = 98), interpolate
in x direction (min: 1.90243e-05) (mean: 0.0263859) (max: 1.95774)
in y direction (min: 2.59451e-05) (mean: 0.00647664) (max: 0.229751)
in z direction (min: 1.97233e-05) (mean: 0.00013241) (max: 0.000594703)
PACD ... compute with :
\=> CH3COOOH (1) : read, arrange, convert (mw = 76), interpolate
in x direction (min: 0.0213156) (mean: 0.157959) (max: 0.408245)
in y direction (min: 0.0241399) (mean: 0.173751) (max: 0.647759)
in z direction (min: 0.0393123) (mean: 0.0886785) (max: 0.188767)
PAN ... compute with :
\=> PAN (1) : read, arrange, convert (mw = 121), interpolate
in x direction (min: 0.078533) (mean: 1.70621) (max: 14.5611)
in y direction (min: 0.0816936) (mean: 1.07385) (max: 4.27246)
in z direction (min: 0.0963284) (mean: 0.299751) (max: 0.489334)
PAR ... compute with :
\=> C3H6 (1) : read, arrange, convert (mw = 42), interpolate
\=> C3H8 (1.5) : read, arrange, convert (mw = 44), interpolate
\=> CH3COCH3 (3) : read, arrange, convert (mw = 58), interpolate
\=> BIGALK (3) : read, arrange, convert (mw = 72), interpolate
\=> BIGENE (1) : read, arrange, convert (mw = 56), interpolate
\=> DMS (2) : read, arrange, convert (mw = 62), interpolate
\=> MEK (4) : read, arrange, convert (mw = 72), interpolate
\=> HYAC (2) : read, arrange, convert (mw = 74), interpolate
in x direction (min: 1.60657) (mean: 27.123) (max: 185.513)
in y direction (min: 1.6732) (mean: 20.5736) (max: 78.0519)
in z direction (min: 1.97934) (mean: 4.11296) (max: 6.80535)
SO2 ... compute with :
\=> SO2 (1) : read, arrange, convert (mw = 64), interpolate
in x direction (min: 1.565e-11) (mean: 5.90473) (max: 124.115)
in y direction (min: 6.0691e-07) (mean: 1.06741) (max: 12.7289)
in z direction (min: 1.18277e-11) (mean: 0.00215863) (max: 0.0102977)
TERP ... compute with :
\=> C10H16 (1) : read, arrange, convert (mw = 136), interpolate
in x direction (min: 3.41007e-20) (mean: 0.0118665) (max: 0.246605)
in y direction (min: 4.29438e-19) (mean: 0.0146822) (max: 0.512346)
in z direction (min: 3.00884e-14) (mean: 3.1937e-07) (max: 6.09901e-05)
TOL ... compute with :
\=> TOLUENE (1) : read, arrange, convert (mw = 92), interpolate
in x direction (min: 0.00137877) (mean: 1.40692) (max: 15.8782)
in y direction (min: 0.00190924) (mean: 1.00644) (max: 6.86484)
in z direction (min: 0.00109195) (mean: 0.0417381) (max: 0.152854)
==============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/bc-gas
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/bc-gas
==> GET BOUNDARY CONDITIONS FOR AER <===
Get MOZART4-GEOS5 data from 2008-12-31_06 to 2009-01-03_00
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute BC for AER species...
<< From MOZART4 : 8 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
CB1 ... read ... arrange ... convert (mw = 12) ... interpolate
CB2 ... read ... arrange ... convert (mw = 12) ... interpolate
OC1 ... read ... arrange ... convert (mw = 12) ... interpolate
OC2 ... read ... arrange ... convert (mw = 12) ... interpolate
SO4 ... read ... arrange ... convert (mw = 96) ... interpolate
NH4NO3 ... read ... arrange ... convert (mw = 80) ... interpolate
SOA ... read ... arrange ... convert (mw = 144) ... interpolate
NH4 ... read ... arrange ... convert (mw = 18) ... interpolate
>> to POLAIR : Nt = 12 Nz = 9 Ny = 50 Nx = 65
PBC_0 ... compute with : CB1 (1) CB2 (1) ... write
in x direction (min: 0.000162916) (mean: 0.0189037) (max: 0.328643)
in y direction (min: 0.000182913) (mean: 0.00853049) (max: 0.0702749)
in z direction (min: 0.000181435) (mean: 0.000508676) (max: 0.00136556)
PBC_1 ... compute with : CB1 (1) CB2 (1) ... write
in x direction (min: 0.00115905) (mean: 0.134489) (max: 2.3381)
in y direction (min: 0.00130132) (mean: 0.0606896) (max: 0.499964)
in z direction (min: 0.0012908) (mean: 0.0036189) (max: 0.00971517)
PBC_2 ... compute with : CB1 (1) CB2 (1) ... write
in x direction (min: 0.000438849) (mean: 0.050921) (max: 0.885267)
in y direction (min: 0.000492715) (mean: 0.0229788) (max: 0.1893)
in z direction (min: 0.000488732) (mean: 0.00137022) (max: 0.00367843)
PBC_3 ... compute with : CB1 (1) CB2 (1) ... write
in x direction (min: 6.88136e-06) (mean: 0.000798466) (max: 0.0138814)
in y direction (min: 7.726e-06) (mean: 0.000360317) (max: 0.00296831)
in z direction (min: 7.66356e-06) (mean: 2.14855e-05) (max: 5.76795e-05)
PBC_4 ... compute with : CB1 (1) CB2 (1) ... write
in x direction (min: 2.87323e-09) (mean: 3.33389e-07) (max: 5.79603e-06)
in y direction (min: 3.2259e-09) (mean: 1.50446e-07) (max: 1.23939e-06)
in z direction (min: 3.19983e-09) (mean: 8.97115e-09) (max: 2.40834e-08)
PBiA2D_0 ... compute with : SOA (1) ... write
in x direction (min: 2.73835e-05) (mean: 0.000222262) (max: 0.000919446)
in y direction (min: 2.97429e-05) (mean: 0.000209591) (max: 0.000682385)
in z direction (min: 2.84636e-05) (mean: 7.68768e-05) (max: 0.000138905)
PBiA2D_1 ... compute with : SOA (1) ... write
in x direction (min: 0.00089807) (mean: 0.00728928) (max: 0.0301541)
in y direction (min: 0.000975446) (mean: 0.00687374) (max: 0.0223795)
in z direction (min: 0.000933492) (mean: 0.00252126) (max: 0.00455553)
PBiA2D_2 ... compute with : SOA (1) ... write
in x direction (min: 0.00231667) (mean: 0.0188037) (max: 0.0777861)
in y direction (min: 0.00251628) (mean: 0.0177316) (max: 0.0577305)
in z direction (min: 0.00240805) (mean: 0.00650382) (max: 0.0117515)
PBiA2D_3 ... compute with : SOA (1) ... write
in x direction (min: 0.000539454) (mean: 0.00437854) (max: 0.0181131)
in y direction (min: 0.000585933) (mean: 0.00412894) (max: 0.013443)
in z direction (min: 0.000560732) (mean: 0.00151447) (max: 0.00273643)
PBiA2D_4 ... compute with : SOA (1) ... write
in x direction (min: 9.34119e-06) (mean: 7.58189e-05) (max: 0.000313646)
in y direction (min: 1.0146e-05) (mean: 7.14968e-05) (max: 0.000232778)
in z direction (min: 9.70964e-06) (mean: 2.62244e-05) (max: 4.7384e-05)
PNH4_0 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
in x direction (min: 5.41769e-07) (mean: 0.00011741) (max: 0.00184139)
in y direction (min: 5.94088e-07) (mean: 3.04464e-05) (max: 0.000433717)
in z direction (min: 5.63305e-07) (mean: 1.70389e-06) (max: 4.30661e-06)
PNH4_1 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
in x direction (min: 0.000283389) (mean: 0.0614149) (max: 0.9632)
in y direction (min: 0.000310757) (mean: 0.0159259) (max: 0.226869)
in z direction (min: 0.000294655) (mean: 0.000891278) (max: 0.00225271)
PNH4_2 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
in x direction (min: 0.00517517) (mean: 1.12154) (max: 17.5897)
in y direction (min: 0.00567494) (mean: 0.290835) (max: 4.14302)
in z direction (min: 0.00538089) (mean: 0.0162762) (max: 0.0411382)
PNH4_3 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
in x direction (min: 0.00508373) (mean: 1.10172) (max: 17.2789)
in y direction (min: 0.00557467) (mean: 0.285696) (max: 4.06981)
in z direction (min: 0.00528582) (mean: 0.0159887) (max: 0.0404113)
PNH4_4 ... compute with : NH4 (1) NH4NO3 (0.225) ... write
in x direction (min: 0.000267428) (mean: 0.057956) (max: 0.90895)
in y direction (min: 0.000293254) (mean: 0.015029) (max: 0.214091)
in z direction (min: 0.000278059) (mean: 0.000841082) (max: 0.00212583)
PNO3_0 ... compute with : NH4NO3 (0.775) ... write
in x direction (min: 0) (mean: 0.000166776) (max: 0.0028438)
in y direction (min: 0) (mean: 3.36749e-05) (max: 0.000651205)
in z direction (min: 0) (mean: 2.18165e-07) (max: 2.14546e-06)
PNO3_1 ... compute with : NH4NO3 (0.775) ... write
in x direction (min: 0) (mean: 0.0872373) (max: 1.48754)
in y direction (min: 0) (mean: 0.0176147) (max: 0.340633)
in z direction (min: 0) (mean: 0.000114118) (max: 0.00112225)
PNO3_2 ... compute with : NH4NO3 (0.775) ... write
in x direction (min: 0) (mean: 1.59311) (max: 27.165)
in y direction (min: 0) (mean: 0.321674) (max: 6.22054)
in z direction (min: 0) (mean: 0.00208399) (max: 0.0204942)
PNO3_3 ... compute with : NH4NO3 (0.775) ... write
in x direction (min: 0) (mean: 1.56495) (max: 26.685)
in y direction (min: 0) (mean: 0.315991) (max: 6.11063)
in z direction (min: 0) (mean: 0.00204716) (max: 0.0201321)
PNO3_4 ... compute with : NH4NO3 (0.775) ... write
in x direction (min: 0) (mean: 0.0823239) (max: 1.40376)
in y direction (min: 0) (mean: 0.0166226) (max: 0.321448)
in z direction (min: 0) (mean: 0.000107691) (max: 0.00105904)
PSO4_0 ... compute with : SO4 (1) ... write
in x direction (min: 2.02822e-06) (mean: 6.05442e-05) (max: 0.000566023)
in y direction (min: 2.41894e-06) (mean: 3.50357e-05) (max: 0.000173389)
in z direction (min: 2.28883e-06) (mean: 6.06666e-06) (max: 1.30434e-05)
PSO4_1 ... compute with : SO4 (1) ... write
in x direction (min: 0.00106093) (mean: 0.0316696) (max: 0.296076)
in y direction (min: 0.0012653) (mean: 0.0183265) (max: 0.0906966)
in z direction (min: 0.00119725) (mean: 0.00317337) (max: 0.00682274)
PSO4_2 ... compute with : SO4 (1) ... write
in x direction (min: 0.0193743) (mean: 0.578341) (max: 5.40685)
in y direction (min: 0.0231065) (mean: 0.334672) (max: 1.65627)
in z direction (min: 0.0218637) (mean: 0.0579509) (max: 0.124595)
PSO4_3 ... compute with : SO4 (1) ... write
in x direction (min: 0.019032) (mean: 0.568122) (max: 5.31132)
in y direction (min: 0.0226983) (mean: 0.32876) (max: 1.62701)
in z direction (min: 0.0214774) (mean: 0.0569269) (max: 0.122393)
PSO4_4 ... compute with : SO4 (1) ... write
in x direction (min: 0.00100117) (mean: 0.0298858) (max: 0.2794)
in y direction (min: 0.00119404) (mean: 0.0172943) (max: 0.0855884)
in z direction (min: 0.00112981) (mean: 0.00299462) (max: 0.00643846)
PSOAhP_0 ... compute with : OC1 (0.79) OC2 (0.79) ... write
in x direction (min: 3.09713e-05) (mean: 0.00155073) (max: 0.0270875)
in y direction (min: 2.97383e-05) (mean: 0.000524032) (max: 0.00386534)
in z direction (min: 2.25836e-05) (mean: 8.02277e-05) (max: 0.000206219)
PSOAhP_1 ... compute with : OC1 (0.79) OC2 (0.79) ... write
in x direction (min: 0.00101573) (mean: 0.0508575) (max: 0.88836)
in y direction (min: 0.000975295) (mean: 0.0171861) (max: 0.126768)
in z direction (min: 0.000740651) (mean: 0.00263113) (max: 0.00676316)
PSOAhP_2 ... compute with : OC1 (0.79) OC2 (0.79) ... write
in x direction (min: 0.0026202) (mean: 0.131193) (max: 2.29163)
in y direction (min: 0.00251589) (mean: 0.0443336) (max: 0.327012)
in z direction (min: 0.00191059) (mean: 0.00678732) (max: 0.0174464)
PSOAhP_3 ... compute with : OC1 (0.79) OC2 (0.79) ... write
in x direction (min: 0.000610133) (mean: 0.0305491) (max: 0.533622)
in y direction (min: 0.000585842) (mean: 0.0103234) (max: 0.0761471)
in z direction (min: 0.000444896) (mean: 0.00158049) (max: 0.00406251)
PSOAhP_4 ... compute with : OC1 (0.79) OC2 (0.79) ... write
in x direction (min: 1.05651e-05) (mean: 0.00052899) (max: 0.0092402)
in y direction (min: 1.01444e-05) (mean: 0.000178759) (max: 0.00131856)
in z direction (min: 7.70382e-06) (mean: 2.73674e-05) (max: 7.03464e-05)
PSOAlP_0 ... compute with : OC1 (0.46) OC2 (0.46) ... write
in x direction (min: 1.80339e-05) (mean: 0.000902955) (max: 0.0157724)
in y direction (min: 1.73159e-05) (mean: 0.00030513) (max: 0.00225071)
in z direction (min: 1.31499e-05) (mean: 4.67146e-05) (max: 0.000120077)
PSOAlP_1 ... compute with : OC1 (0.46) OC2 (0.46) ... write
in x direction (min: 0.00059144) (mean: 0.0296133) (max: 0.517273)
in y direction (min: 0.000567893) (mean: 0.0100071) (max: 0.0738141)
in z direction (min: 0.000431265) (mean: 0.00153205) (max: 0.00393804)
PSOAlP_2 ... compute with : OC1 (0.46) OC2 (0.46) ... write
in x direction (min: 0.00152569) (mean: 0.0763908) (max: 1.33437)
in y direction (min: 0.00146495) (mean: 0.0258145) (max: 0.190412)
in z direction (min: 0.0011125) (mean: 0.00395213) (max: 0.0101586)
PSOAlP_3 ... compute with : OC1 (0.46) OC2 (0.46) ... write
in x direction (min: 0.000355267) (mean: 0.0177881) (max: 0.310716)
in y direction (min: 0.000341123) (mean: 0.00601108) (max: 0.0443388)
in z direction (min: 0.000259053) (mean: 0.000920278) (max: 0.00236551)
PSOAlP_4 ... compute with : OC1 (0.46) OC2 (0.46) ... write
in x direction (min: 6.15181e-06) (mean: 0.00030802) (max: 0.00538037)
in y direction (min: 5.90689e-06) (mean: 0.000104088) (max: 0.000767771)
in z direction (min: 4.48577e-06) (mean: 1.59355e-05) (max: 4.09612e-05)
PSOAmP_0 ... compute with : OC1 (0.58) OC2 (0.58) ... write
in x direction (min: 2.27384e-05) (mean: 0.00113851) (max: 0.019887)
in y direction (min: 2.18332e-05) (mean: 0.000384732) (max: 0.00283785)
in z direction (min: 1.65804e-05) (mean: 5.89013e-05) (max: 0.000151401)
PSOAmP_1 ... compute with : OC1 (0.58) OC2 (0.58) ... write
in x direction (min: 0.000745729) (mean: 0.0373384) (max: 0.652213)
in y direction (min: 0.00071604) (mean: 0.0126177) (max: 0.09307)
in z direction (min: 0.000543769) (mean: 0.00193173) (max: 0.00496536)
PSOAmP_2 ... compute with : OC1 (0.58) OC2 (0.58) ... write
in x direction (min: 0.00192369) (mean: 0.0963187) (max: 1.68246)
in y direction (min: 0.00184711) (mean: 0.0325487) (max: 0.240085)
in z direction (min: 0.00140272) (mean: 0.00498312) (max: 0.0128087)
PSOAmP_3 ... compute with : OC1 (0.58) OC2 (0.58) ... write
in x direction (min: 0.000447946) (mean: 0.0224286) (max: 0.391773)
in y direction (min: 0.000430112) (mean: 0.00757919) (max: 0.0559055)
in z direction (min: 0.000326632) (mean: 0.00116035) (max: 0.0029826)
PSOAmP_4 ... compute with : OC1 (0.58) OC2 (0.58) ... write
in x direction (min: 7.75663e-06) (mean: 0.000388372) (max: 0.00678394)
in y direction (min: 7.44782e-06) (mean: 0.000131241) (max: 0.000968059)
in z direction (min: 5.65597e-06) (mean: 2.00927e-05) (max: 5.16467e-05)
==============>>> DONE <<<===============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/bc-aer
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/bc-aer
==> GET BOUNDARY CONDITIONS FOR DU <===
Get MOZART4-GEOS5 data from 2008-12-31_06 to 2009-01-03_00
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute BC for dust species...
<< From MOZART4 : 4 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
DUST1 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST2 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST3 ... read ... arrange ... convert (mw = 34) ... interpolate
DUST4 ... read ... arrange ... convert (mw = 34) ... interpolate
>> to POLAIR : Nt = 12 Nz = 9 Ny = 50 Nx = 65
PMD_0 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
in x direction (min: 0) (mean: 0) (max: 0)
in y direction (min: 0) (mean: 0) (max: 0)
in z direction (min: 0) (mean: 0) (max: 0)
PMD_1 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
in x direction (min: 0.0282469) (mean: 0.360582) (max: 3.25539)
in y direction (min: 0.0278979) (mean: 0.732533) (max: 4.22354)
in z direction (min: 0.0301188) (mean: 0.0589696) (max: 0.158664)
PMD_2 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
in x direction (min: 0.134597) (mean: 1.84192) (max: 16.8596)
in y direction (min: 0.132584) (mean: 3.77329) (max: 22.0131)
in z direction (min: 0.147217) (mean: 0.29521) (max: 0.807604)
PMD_3 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
in x direction (min: 0.264265) (mean: 4.78993) (max: 47.6084)
in y direction (min: 0.258852) (mean: 10.3643) (max: 65.0957)
in z direction (min: 0.314208) (mean: 0.701676) (max: 2.1052)
PMD_4 ... compute with : DUST1 (1) DUST2 (1) DUST3 (1) DUST4 (1) ... write
in x direction (min: 0.0040325) (mean: 0.40124) (max: 6.42162)
in y direction (min: 0.0034178) (mean: 1.20437) (max: 11.6496)
in z direction (min: 0.00701471) (mean: 0.0398671) (max: 0.197547)
=============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/bc-DU
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/bc-DU
==> GET BOUNDARY CONDITIONS FOR SS <===
Get MOZART4-GEOS5 data from 2008-12-31_06 to 2009-01-03_00
Read configuration... Compute aerosol size distribution... Get speciation from input to output species... Get molecular weights of input species... Define MOZART4 domain... Define output domain... done
Memory allocation for input fields... for output fields... done
Extract Temperature ... Compute level heights ... Arrange data ... done
Compute BC for sea-salt species...
<< From MOZART4 : 4 specie(s), Nt = 12 Nz = 57 Ny = 21 Nx = 21
SA1 ... read ... arrange ... convert (mw = 58) ... interpolate
SA2 ... read ... arrange ... convert (mw = 58) ... interpolate
SA3 ... read ... arrange ... convert (mw = 58) ... interpolate
SA4 ... read ... arrange ... convert (mw = 58) ... interpolate
>> to POLAIR : Nt = 12 Nz = 9 Ny = 50 Nx = 65
PHCL_0 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
in x direction (min: 0) (mean: 0) (max: 0)
in y direction (min: 0) (mean: 0) (max: 0)
in z direction (min: 0) (mean: 0) (max: 0)
PHCL_1 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
in x direction (min: 0.000211575) (mean: 0.0119106) (max: 0.107979)
in y direction (min: 0.000357228) (mean: 0.00739101) (max: 0.0555785)
in z direction (min: 0.000635796) (mean: 0.00307855) (max: 0.00989984)
PHCL_2 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
in x direction (min: 0.002114) (mean: 0.14854) (max: 1.62535)
in y direction (min: 0.00403497) (mean: 0.0759409) (max: 0.846064)
in z direction (min: 0.0043034) (mean: 0.019742) (max: 0.0773042)
PHCL_3 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
in x direction (min: 0.000117272) (mean: 0.83675) (max: 11.1964)
in y direction (min: 0.000239659) (mean: 0.331549) (max: 5.9283)
in z direction (min: 0.000245623) (mean: 0.0136234) (max: 0.245479)
PHCL_4 ... compute with : SA1 (0.6) SA2 (0.6) SA3 (0.6) SA4 (0.6) ... write
in x direction (min: 3.75748e-06) (mean: 0.937989) (max: 14.8006)
in y direction (min: 1.49721e-05) (mean: 0.335462) (max: 8.27117)
in z direction (min: 6.96381e-06) (mean: 0.00468909) (max: 0.135639)
PNA_0 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
in x direction (min: 0) (mean: 0) (max: 0)
in y direction (min: 0) (mean: 0) (max: 0)
in z direction (min: 0) (mean: 0) (max: 0)
PNA_1 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
in x direction (min: 0.00014105) (mean: 0.00794039) (max: 0.0719862)
in y direction (min: 0.000238152) (mean: 0.00492735) (max: 0.0370523)
in z direction (min: 0.000423864) (mean: 0.00205238) (max: 0.0065999)
PNA_2 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
in x direction (min: 0.00140933) (mean: 0.0990274) (max: 1.08356)
in y direction (min: 0.00268998) (mean: 0.0506272) (max: 0.564043)
in z direction (min: 0.00286893) (mean: 0.0131614) (max: 0.0515361)
PNA_3 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
in x direction (min: 7.81816e-05) (mean: 0.557833) (max: 7.46428)
in y direction (min: 0.000159773) (mean: 0.221033) (max: 3.9522)
in z direction (min: 0.000163748) (mean: 0.00908237) (max: 0.163653)
PNA_4 ... compute with : SA1 (0.4) SA2 (0.4) SA3 (0.4) SA4 (0.4) ... write
in x direction (min: 2.50499e-06) (mean: 0.625324) (max: 9.86707)
in y direction (min: 9.98139e-06) (mean: 0.223641) (max: 5.51411)
in z direction (min: 4.64254e-06) (mean: 0.00312607) (max: 0.0904257)
=============>>> DONE <<<==============
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data/bc-SS
with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data/bc-SS
assert_no_difference(
compare_dir("data", "general.cfg"))
Comparing /libre_RAID1/dores/dev/polyphemus/src/Polyphemus-1.8.1/preprocessing/mozart4/example/data with /cerea_raid/share/results/polyphemus/dev/preprocessing/mozart4/example/data