In [1]:

    !$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

In [2]:

    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

# Mozart4-geos5 integration test¶

## The ic and bc preprocessing use the following files in config/:¶

- general.cfg: data domain, with the path to the level file.
- levels.dat: vertical level altitudes.
- mozart4.cfg: Mozart4 parameters, input files and output directories.
- input/molecularweights_mozart4{gas,aer,ss,du}.txt: species molecular weights
- input/speciationmozart4{gas,aer,ss,du}.dat: speciation factors

In [3]:

    !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

In [4]:

    !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.

In [5]:

    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

In [6]:

    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

In [7]:

    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