Changeset 8fcfd08 in flexpart.git
- Timestamp:
- Apr 3, 2018, 2:48:23 PM (6 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- a9cf4b1
- Parents:
- b1e0742 (diff), 5b7ec80 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent. - Files:
-
- 2 added
- 33 deleted
- 33 edited
-
options/COMMAND (modified) (3 diffs)
-
options/RELEASES (modified) (1 diff)
-
options/SPECIES/SPECIES.README (modified) (1 diff)
-
options/SPECIES/SPECIES_001 (deleted)
-
options/SPECIES/SPECIES_001.oldformat (deleted)
-
options/SPECIES/SPECIES_002 (modified) (1 diff)
-
options/SPECIES/SPECIES_002.oldformat (deleted)
-
options/SPECIES/SPECIES_003 (modified) (1 diff)
-
options/SPECIES/SPECIES_003.oldformat (deleted)
-
options/SPECIES/SPECIES_004 (modified) (1 diff)
-
options/SPECIES/SPECIES_004.oldformat (deleted)
-
options/SPECIES/SPECIES_005 (modified) (1 diff)
-
options/SPECIES/SPECIES_005.oldformat (deleted)
-
options/SPECIES/SPECIES_006 (modified) (1 diff)
-
options/SPECIES/SPECIES_006.oldformat (deleted)
-
options/SPECIES/SPECIES_007 (modified) (1 diff)
-
options/SPECIES/SPECIES_007.oldformat (deleted)
-
options/SPECIES/SPECIES_008 (deleted)
-
options/SPECIES/SPECIES_008.oldformat (deleted)
-
options/SPECIES/SPECIES_009 (modified) (1 diff)
-
options/SPECIES/SPECIES_009.oldformat (deleted)
-
options/SPECIES/SPECIES_010 (modified) (1 diff)
-
options/SPECIES/SPECIES_010.oldformat (deleted)
-
options/SPECIES/SPECIES_011 (modified) (1 diff)
-
options/SPECIES/SPECIES_011.oldformat (deleted)
-
options/SPECIES/SPECIES_012 (modified) (1 diff)
-
options/SPECIES/SPECIES_012.oldformat (deleted)
-
options/SPECIES/SPECIES_013 (modified) (1 diff)
-
options/SPECIES/SPECIES_013.oldformat (deleted)
-
options/SPECIES/SPECIES_014 (modified) (1 diff)
-
options/SPECIES/SPECIES_014.oldformat (deleted)
-
options/SPECIES/SPECIES_015 (modified) (1 diff)
-
options/SPECIES/SPECIES_015.oldformat (deleted)
-
options/SPECIES/SPECIES_016 (modified) (1 diff)
-
options/SPECIES/SPECIES_016.oldformat (deleted)
-
options/SPECIES/SPECIES_017 (modified) (1 diff)
-
options/SPECIES/SPECIES_017.oldformat (deleted)
-
options/SPECIES/SPECIES_018 (modified) (1 diff)
-
options/SPECIES/SPECIES_018.oldformat (deleted)
-
options/SPECIES/SPECIES_019 (modified) (1 diff)
-
options/SPECIES/SPECIES_019.oldformat (deleted)
-
options/SPECIES/SPECIES_020 (modified) (1 diff)
-
options/SPECIES/SPECIES_020.oldformat (deleted)
-
options/SPECIES/SPECIES_021 (modified) (1 diff)
-
options/SPECIES/SPECIES_021.oldformat (deleted)
-
options/SPECIES/SPECIES_022 (modified) (1 diff)
-
options/SPECIES/SPECIES_022.oldformat (deleted)
-
options/SPECIES/SPECIES_023 (modified) (1 diff)
-
options/SPECIES/SPECIES_023.oldformat (deleted)
-
options/SPECIES/SPECIES_024 (modified) (1 diff)
-
options/SPECIES/SPECIES_024.oldformat (deleted)
-
options/SPECIES/SPECIES_025 (modified) (1 diff)
-
options/SPECIES/SPECIES_025.oldformat (deleted)
-
options/SPECIES/SPECIES_026 (modified) (1 diff)
-
options/SPECIES/SPECIES_026.oldformat (deleted)
-
options/SPECIES/SPECIES_027 (modified) (1 diff)
-
options/SPECIES/SPECIES_027.oldformat (deleted)
-
options/SPECIES/SPECIES_028.oldformat (deleted)
-
options/SPECIES/SPECIES_031 (modified) (1 diff)
-
options/SPECIES/SPECIES_031.oldformat (deleted)
-
options/SPECIES/SPECIES_034 (modified) (1 diff)
-
options/SPECIES/SPECIES_034.oldformat (deleted)
-
options/SPECIES/SPECIES_040 (added)
-
options/SPECIES/spec_overview (deleted)
-
options/SPECIES/specoverview.f90 (added)
-
src/concoutput_mpi.f90 (modified) (1 diff)
-
src/readOHfield.f90 (modified) (3 diffs)
-
src/timemanager_mpi.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
options/COMMAND
rf28aa0a r5b7ec80 1 1 &COMMAND 2 2 LDIRECT= 1, 3 IBDATE= 201 40907,3 IBDATE= 20170102, 4 4 IBTIME= 060000, 5 IEDATE= 201 40907,5 IEDATE= 20170102, 6 6 IETIME= 120000, 7 7 LOUTSTEP= 3600, … … 12 12 CTL= -5.0000000 , 13 13 IFINE= 4, 14 IOUT= 9,14 IOUT= 2, 15 15 IPOUT= 0, 16 LSUBGRID= 1,17 LCONVECTION= 1,16 LSUBGRID= 0, 17 LCONVECTION= 0, 18 18 LAGESPECTRA= 0, 19 19 IPIN= 0, … … 24 24 IND_RECEPTOR= 1, 25 25 MQUASILAG= 0, 26 NESTED_OUTPUT= 1,26 NESTED_OUTPUT= 0, 27 27 LINIT_COND= 0, 28 28 SURF_ONLY= 0, -
options/RELEASES
r79abee9 r5b7ec80 1 1 &RELEASES_CTRL 2 2 NSPEC= 1, 3 SPECNUM_REL= 1,3 SPECNUM_REL= 24, 4 4 / 5 5 &RELEASE 6 IDATE1= 201 40907,7 ITIME1= 0 70000,8 IDATE2= 201 40907,9 ITIME2= 0 70000,10 LON1= -75.255,11 LON2= -75.255,12 LAT1= 38.005,13 LAT2= 38.005,14 Z1= 2500.000000,15 Z2= 2500.000000,6 IDATE1= 20170102, 7 ITIME1= 090000, 8 IDATE2= 20170102, 9 ITIME2= 090000, 10 LON1= 0.000 , 11 LON2= 0.000 , 12 LAT1= 0.000 , 13 LAT2= 0.000 , 14 Z1= 50.000 , 15 Z2= 50.000 , 16 16 ZKIND= 1, 17 17 MASS= 1.0000E8 , 18 PARTS= 10000 018 PARTS= 10000 19 19 COMMENT="TEST1", 20 20 / -
options/SPECIES/SPECIES.README
rb836152 r8d70e43 1 1 Specifies all physico-chemical properties for the species. 2 --------------------------------------------------------------------------------------------------------------------------------------------------------| 3 Specie | Radioact. | WetDep(below-Cl)| WetDep(in-Cl)| DryDepo(gas) | DryDepo(particles) Alt. | | OH Reaction | 4 Name | Halflife[s]| A B | Ai Bi | D H f0 | rho dquer dsig vd | molweight| C** D[T] N*** | 5 --------------------------------------------------------------------------------------------------------------------------------------------------------| 6 1 TRACER | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 350.0 | -9.9E-09 -9.9 -9.9 | 7 2 O3 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | 1.5 1.0e-02 1.0 | -9.9E09 -9.99 | 48.00 | -9.9E-09 -9.9 -9.9 | 8 3 NO | -999.9 | 8.0E-06 0.62 | 1.0 1.0 | 1.2 3.0e-03 0.0 | -9.9E09 -9.99 | 30.00 | -9.9E-09 -9.9 -9.9 | 9 4 NO2 | -999.9 | 1.0E-05 0.62 | 1.0 1.0 | 1.6 1.0e-02 0.1 | -9.9E09 -9.99 | 46.00 | -9.9E-09 -9.9 -9.9 | 10 5 HNO3 | -999.9 | 5.0E-05 0.62 | 1.0 1.0 | 1.9 1.0e+14 0.0 | -9.9E09 -9.99 | 63.00 | -9.9E-09 -9.9 -9.9 | 11 6 HNO2 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | 1.6 1.0e+05 0.1 | -9.9E09 -9.99 | 47.00 | -9.9E-09 -9.9 -9.9 | 12 7 H2O2 | -999.9 | 1.0E-04 0.62 | 1.0 1.0 | 1.4 1.0e+05 1.0 | -9.9E09 -9.99 | 34.00 | -9.9E-09 -9.9 -9.9 | 13 8 SO2var | -999.9 | -9.9 -9.9 | -9.9 -9.9 | 2.0 1.0e+05 0.0 | -9.9E09 -9.99 | 64.00 | -9.9E-09 -9.9 -9.9 | 14 9 HCHO | -999.9 | -9.9 -9.9 | -9.9 -9.9 | 1.3 6.0e+03 0.0 | -9.9E09 -9.99 | 30.00 | -9.9E-09 -9.9 -9.9 | 15 10 PAN | -999.9 | -9.9 -9.9 | -9.9 -9.9 | 2.6 3.6e+00 0.1 | -9.9E09 -9.99 | 121.00 | -9.9E-09 -9.9 -9.9 | 16 11 NH3 | -999.9 | 9.9E-05 0.62 | 1.0 1.0 | 1.1 2.0e+14 0.0 | -9.9E09 -9.99 | 17.00 | -9.9E-09 -9.9 -9.9 | 17 12 SO4-aero | -999.9 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.0E03 4.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 18 13 NO3-aero | -999.9 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.0E03 4.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 19 14 I2-131 | 691200.0 | 8.0E-05 0.62 | 1.0 1.0 | 2.7 1.0e+05 0.1 | -9.9E09 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 20 15 I-131 | 691200.0 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.5E03 6.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 21 16 Cs-137 | -999.9 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.5E03 6.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 22 17 Y-91 | 5037120.0 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.5E03 6.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 23 18 Ru-106 |31536000.0 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.5E03 6.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 24 19 Kr-85 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 25 20 Sr-90 | -999.9 | 1.0 1.0 | 0.9 0.1 | -9.9 | 2.5E03 6.0E-7 3.0E-1 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 26 21 Xe-133 | 453168.0 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | -9.99 | -9.9E-09 -9.9 -9.9 | 27 22 CO | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 28.00 | ? ? ? | 28 23 SO2dryOH | | | | | | | ? ? ? | 29 24 AIRTRACER | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 29.00 | -9.9E-09 -9.9 -9.9 | 30 25 AEROTRACER| -999.9 | 1.0 1.0 | 0.9 0.1 | -9.9 | 1.4E+03 2.5E-07 1.25 -9.99 | 29.00 | -9.9E-09 -9.9 -9.9 | 31 26 CH4 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 16.00 | 9.65E-20 1082.0 2.58 | 32 27 C2H6 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 30.00 | 15.2E-18 498.0 2.0 | 33 28 C3H8 | -999.9 | -9.9 -9.9 | -9.9 -9.9 | -9.9 | -9.9E09 -9.99 | 44.00 | 15.5E-18 61.0 2.0 | 34 x PCB28 | -999.9 | 5E-06 0.62 | -9.9 -9.9 | 3.8 3.1e-00 0.0 | -9.9E09 -9.99 | 257.54 | ? ? 2.0 | 35 x PCB101 | -999.9 | 5E-06 0.62 | -9.9 -9.9 | 4.3 3.2e-00 0.0 | -9.9E09 -9.99 | 326.43 | ? ? 2.0 | 36 x PCB180 | -999.9 | 5E-06 0.62 | -9.9 -9.9 | 4.7 1.7e+01 0.0 | -9.9E09 -9.99 | 395.32 | ? ? 2.0 | 37 --------------------------------------------------------------------------------------------------------------------------------------------------------| 38 ** unit [cm^3/molec/s] (in FLEXPART version 9.2 and below this had unit [cm3/s], note the unit is now changed!) 39 *** no unit 2 3 To generate an overview table compile the specoverview.f90 program 4 5 gfortran specoverview.f90 -o specoverview 6 ./specoverview 40 7 41 8 UPDATE FOR FLEXPART VERSION 10.0 DEV: -
options/SPECIES/SPECIES_002
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="O3 ", 3 PRELDIFF= 1.50000000 , 4 PHENRY= 1.0E-02, 5 PF0= 1.00000000 , 6 PWEIGHTMOLAR= 48.0000000 , 2 PSPECIES="Ozone", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=-9.9, ! Dry deposition (particles) - dquer 12 PDSIGMA=-9.9, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.9, ! Alternative: dry deposition velocity 14 PRELDIFF=1.6, ! Dry deposition (gases) - D 15 PHENRY=1.0E-2, ! Dry deposition (gases) - Henrys const. 16 PF0=1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=48.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 7 21 / -
options/SPECIES/SPECIES_003
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="NO ", 3 PWETA_GAS= 8.0E-06, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.2 , 6 PHENRY= 3.0E-03, 7 PWEIGHTMOLAR= 30.0000000 , 2 PSPECIES="NO", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=8.0E-6, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.2, ! Dry deposition (gases) - D 15 PHENRY=3.0E-03, ! Dry deposition (gases) - Henrys const. 16 PF0=0.0, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=30.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9 ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 8 21 / -
options/SPECIES/SPECIES_004
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="NO2 ", 3 PWETA_GAS= 1.0E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.6 , 6 PHENRY= 1.0E-02, 7 PF0= 0.1 , 8 PWEIGHTMOLAR= 46.0000000 , 2 PSPECIES="NO2", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=1.0E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.6, ! Dry deposition (gases) - D 15 PHENRY=1.0E-02, ! Dry deposition (gases) - Henrys const. 16 PF0=0.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=46.0, ! molweight 18 POHCCONST=-0.9E9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 9 21 / -
options/SPECIES/SPECIES_005
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="HNO3 ", 3 PWETA_GAS= 5.0E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.9 , 6 PHENRY= 1.00000000E+14, 7 PWEIGHTMOLAR= 63.0000000 , 2 PSPECIES="HNO3", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=5.0E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.9, ! Dry deposition (gases) - D 15 PHENRY=1E14, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=63.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 8 21 / -
options/SPECIES/SPECIES_006
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="HNO2 ", 3 PRELDIFF= 1.6 , 4 PHENRY= 100000.000 , 5 PF0= 0.1 , 6 PWEIGHTMOLAR= 47.0000000 , 2 PSPECIES="HNO2", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.6, ! Dry deposition (gases) - D 15 PHENRY=100000.0, ! Dry deposition (gases) - Henrys const. 16 PF0=.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=47.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 7 21 / -
options/SPECIES/SPECIES_007
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="H2O2 ", 3 PWETA_GAS= 1.0E-04, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.4 , 6 PHENRY= 100000.000 , 7 PF0= 1.00000000 , 8 PWEIGHTMOLAR= 34.0000000 , 2 PSPECIES="H2O2", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=1.0E-04, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.4, ! Dry deposition (gases) - D 15 PHENRY=100000.0, ! Dry deposition (gases) - Henrys const. 16 PF0=1.0, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=34.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 9 21 / -
options/SPECIES/SPECIES_009
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="HCHO ", 3 PRELDIFF= 1.3 , 4 PHENRY= 6000.00000 , 5 PWEIGHTMOLAR= 30.0000000 , 2 PSPECIES="HCHO", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1,3, ! Dry deposition (gases) - D 15 PHENRY=6000.0, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=30.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 6 21 / -
options/SPECIES/SPECIES_010
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="PAN ", 3 PRELDIFF= 2.6 , 4 PHENRY= 3.6 , 5 PF0= 0.1 , 6 PWEIGHTMOLAR= 121.000000 , 2 PSPECIES="PAN", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=2.6, ! Dry deposition (gases) - D 15 PHENRY=3.6, ! Dry deposition (gases) - Henrys const. 16 PF0=0.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=121.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 7 21 / -
options/SPECIES/SPECIES_011
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="NH3 ", 3 PWETA_GAS= 9.9E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.1 , 6 PHENRY= 2.0E+14, 7 PWEIGHTMOLAR= 17.0000000 , 2 PSPECIES="NH3", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=9.9E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=1.1, ! Dry deposition (gases) - D 15 PHENRY=2.0E14, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=17.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 8 21 / -
options/SPECIES/SPECIES_012
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="SO4-aero ", 3 PCRAIN_AERO=1.0 , 4 PCSNOW_AERO=1.0 , 5 PCCN_AERO=0.9 , 6 PIN_AERO=0.1 , 7 PDENSITY=2000.00000 , 8 PDQUER=4.0E-07 , 9 PDSIGMA=0.3 , 2 PSPECIES="SO4-aero", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=2000.0, ! Dry deposition (particles) - rho 11 PDQUER=4.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E+9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=0.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 10 21 / -
options/SPECIES/SPECIES_013
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="NO3-aero ", 3 PCRAIN_AERO= 1.0 , 4 PCSNOW_AERO= 1.0 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.1 , 7 PDENSITY= 2000.00000 , 8 PDQUER= 4.0E-07, 9 PDSIGMA= 0.3 , 2 PSPECIES="NO3-aero", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=4.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E+9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=2000.0, ! molweight 18 POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 10 21 / -
options/SPECIES/SPECIES_014
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="I2-131 ", 3 PDECAY= 691200.000 , 4 PWETA_GAS= 8.0E-05, 5 PWETB_GAS= 0.62 , 6 PRELDIFF= 2.70 , 7 PHENRY= 100000.000 , 8 PF0= 0.1 , 2 PSPECIES="I2-131", ! Tracer name 3 PDECAY=691200.0, ! Species half life 4 PWETA_GAS=8.0E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=2.7, ! Dry deposition (gases) - D 15 PHENRY=100000.0, ! Dry deposition (gases) - Henrys const. 16 PF0=0.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=-9.9, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 9 21 / -
options/SPECIES/SPECIES_015
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="I-131 ", 3 PDECAY= 691200.000 , 4 PCRAIN_AERO= 1.0 , 5 PCSNOW_AERO= 1.0 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PDENSITY= 2500.00000 , 9 PDQUER= 6.0E-07, 10 PDSIGMA= 0.3 , 2 PSPECIES=""I-131", ! Tracer name 3 PDECAY=691200.0, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=2500.0, ! Dry deposition (particles) - rho 11 PDQUER=6.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-9.9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=350.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 11 21 / -
options/SPECIES/SPECIES_016
re8929a8 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Cs-137 ", 3 PCRAIN_AERO= 1.0 , 4 PCSNOW_AERO= 1.0 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.9 , 7 PDENSITY= 1900.0 , 8 PDQUER= 0.4E-06, 9 PDSIGMA= 0.3 , 2 PSPECIES="Cs-137", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=1900.0, ! Dry deposition (particles) - rho 11 PDQUER=0.4E-06, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=-9.9, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 10 21 / -
options/SPECIES/SPECIES_017
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Y-91 ", 3 PDECAY= 5037120.00 , 4 PCRAIN_AERO= 1.0 , 5 PCSNOW_AERO= 1.0 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PDENSITY= 2500.00000 , 9 PDQUER= 6.0E-07, 10 PDSIGMA= 0.3 , 2 PSPECIES="Y-91", ! Tracer name 3 PDECAY=5037120.0, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=2500.0, ! Dry deposition (particles) - rho 11 PDQUER=6.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=350.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 11 21 / -
options/SPECIES/SPECIES_018
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Ru-106 ", 3 PDECAY= 31536000.0 , 4 PCRAIN_AERO= 1.0 , 5 PCSNOW_AERO= 1.0 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PDENSITY= 2500.00000 , 9 PDQUER= 6.0E-07, 10 PDSIGMA= 0.3 , 2 PSPECIES=""Ru-106", ! Tracer name 3 PDECAY=31536000.0, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=2500.0, ! Dry deposition (particles) - rho 11 PDQUER=6.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=350.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 11 21 / -
options/SPECIES/SPECIES_019
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Kr-85 ", 2 PSPECIES="KR-85", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=-9.9, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 3 21 / -
options/SPECIES/SPECIES_020
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Sr-90 ", 3 PCRAIN_AERO= 1.0 , 4 PCSNOW_AERO= 1.0 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.1 , 7 PDENSITY= 2500.00000 , 8 PDQUER= 6.0E-07, 9 PDSIGMA= 0.3 , 2 PSPECIES="Sr-90", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=2500.0, ! Dry deposition (particles) - rho 11 PDQUER=6.0E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.3, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=-9.9, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 10 21 / -
options/SPECIES/SPECIES_021
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="Xe-133 ", 3 PDECAY= 453168.000 , 2 PSPECIES="Xe-133", ! Tracer name 3 PDECAY=453168.0, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=-9.9, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 4 21 / -
options/SPECIES/SPECIES_022
r341f4b7 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="CO ", 3 PWEIGHTMOLAR= 28.0000000 , 2 PSPECIES="CO", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=28.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 4 21 / -
options/SPECIES/SPECIES_023
r341f4b7 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="SO2dryOH ", 3 PRELDIFF= 2.00000000 , 4 PHENRY= 100000.000 , 5 PWEIGHTMOLAR= 64.0000000 , 2 PSPECIES="SO2", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E-9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=2.0, ! Dry deposition (gases) - D 15 PHENRY=100000.0, ! Dry deposition (gases) - Henrys const. 16 PF0=0, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=64.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 6 21 / -
options/SPECIES/SPECIES_024
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="AIRTRACER ", 3 PWEIGHTMOLAR= 29.0000000 , 2 PSPECIES="AIRTRACER", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=29.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 4 21 / -
options/SPECIES/SPECIES_025
r341f4b7 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="AERO-TRACE ", 3 PCRAIN_AERO= 1.0 , 4 PCSNOW_AERO= 1.0 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.1 , 7 PDENSITY= 1400.00000 , 8 PDQUER= 2.5E-07, 9 PDSIGMA= 1.25000000 , 10 PWEIGHTMOLAR= 29.0000000 , 2 PSPECIES="AERO-TRACE", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=1400.0, ! Dry deposition (particles) - rho 11 PDQUER=2.5E-07, ! Dry deposition (particles) - dquer 12 PDSIGMA=1.25, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=29.0, ! molweight 18 POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=-9.9, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 11 21 / -
options/SPECIES/SPECIES_026
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="CH4 ", 3 PWEIGHTMOLAR= 16.0000000 , 4 POHCCONST= 9.65E-20, 5 POHDCONST= 1082.00000 , 6 POHNCONST= 2.58 , 2 PSPECIES="CH4", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=16.0, ! molweight 18 POHCCONST=9.65E-20, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=1082.0, ! OH Reaction rate - D [K] 20 POHNCONST=2.58, ! OH Reaction rate - N (no unit) 7 21 / -
options/SPECIES/SPECIES_027
re4f26d1 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="C2H6 ", 3 PWEIGHTMOLAR= 30.0000000 , 4 POHDCONST= 498.000000 , 5 POHNCONST= 2.00000000 , 2 PSPECIES="C2H6", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=-9.9, ! Dry deposition (gases) - D 15 PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const. 16 PF0=-9, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=30.0, ! molweight 18 POHCCONST=15.2E-18, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=498.0, ! OH Reaction rate - D [K] 20 POHNCONST=2.0, ! OH Reaction rate - N (no unit) 6 21 / -
options/SPECIES/SPECIES_031
r35fa90d r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="PCB28 ", 3 PWETA_GAS= 3.5E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 4.3 , 6 PHENRY= 3.1 , 7 PF0= 0.1 , 8 PWEIGHTMOLAR= 257.54 , 9 POHCCONST= 1.07E-11, 10 POHDCONST= 1203.00000 , 11 POHNCONST= 0.00000000 , 2 PSPECIES="PCB28", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=3.5E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=4.3, ! Dry deposition (gases) - D 15 PHENRY=3.1, ! Dry deposition (gases) - Henrys const. 16 PF0=0.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=257.54, ! molweight 18 POHCCONST=1.07E-11, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=1203.0, ! OH Reaction rate - D [K] 20 POHNCONST=0.0, ! OH Reaction rate - N (no unit) 12 21 / -
options/SPECIES/SPECIES_034
r341f4b7 r8d70e43 1 1 &SPECIES_PARAMS 2 PSPECIES="G-HCH ", 3 PWETA_GAS= 3.5E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 4.3 , 6 PHENRY= 422.0 , 7 PF0= 0.10 , 8 PWEIGHTMOLAR= 290.83 , 9 POHCCONST= 6.21E-11, 10 POHDCONST= 1203.00000 , 11 POHNCONST= 2.00000000 , 2 PSPECIES="G-HCH", ! Tracer name 3 PDECAY=-9.9, ! Species half life 4 PWETA_GAS=3.5E-05, ! Below-cloud scavenging (gases) - A (weta_gas) 5 PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) 6 PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero) 7 PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) 8 PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) 9 PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero) 10 PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho 11 PDQUER=0.0, ! Dry deposition (particles) - dquer 12 PDSIGMA=0.0, ! Dry deposition (particles) - dsig 13 PDRYVEL=-9.99, ! Alternative: dry deposition velocity 14 PRELDIFF=4.3, ! Dry deposition (gases) - D 15 PHENRY=422.0, ! Dry deposition (gases) - Henrys const. 16 PF0=0.1, ! Dry deposition (gases) - f0 (reactivity) 17 PWEIGHTMOLAR=290.83, ! molweight 18 POHCCONST=6.21E-11, ! OH Reaction rate - C [cm^3/molecule/sec] 19 POHDCONST=1203.0, ! OH Reaction rate - D [K] 20 POHNCONST=2.00, ! OH Reaction rate - N (no unit) 12 21 / -
src/concoutput_mpi.f90
r6a678e3 rb1e0742 113 113 integer :: ierr 114 114 character(LEN=100) :: dates_char 115 ! character :: dates_char116 115 117 116 ! Measure execution time -
src/readOHfield.f90
rd7aab4b rb1e0742 22 22 subroutine readOHfield 23 23 24 !***************************************************************************** 25 ! * 26 ! Reads the OH field into memory * 27 ! * 28 ! AUTHOR: R.L. Thompson, Nov 2014 * 29 ! * 30 !***************************************************************************** 31 ! * 32 ! Variables: * 33 ! * 34 ! path(numpath) contains the path names * 35 ! lonOH(nxOH) longitude of OH fields * 36 ! latOH(nyOH) latitude of OH fields * 37 ! altOH(nzOH) altitude of OH fields * 38 ! etaOH(nzOH) eta-levels of OH fields * 39 ! OH_field(nxOH,nyOH,nzOH,m) OH concentration (molecules/cm3) * 40 ! * 41 ! * 42 !***************************************************************************** 24 !***************************************************************************** 25 ! * 26 ! Reads the OH field into memory * 27 ! * 28 ! AUTHOR: R.L. Thompson, Nov 2014 * 29 ! * 30 ! UPDATES: * 31 ! 03/2018 SEC: Converted original netCDF files to binary format * 32 !***************************************************************************** 33 ! * 34 ! Variables: * 35 ! * 36 ! path(numpath) contains the path names * 37 ! lonOH(nxOH) longitude of OH fields * 38 ! latOH(nyOH) latitude of OH fields * 39 ! altOH(nzOH) altitude of OH fields * 40 ! etaOH(nzOH) eta-levels of OH fields * 41 ! OH_field(nxOH,nyOH,nzOH,m) OH concentration (molecules/cm3) * 42 ! * 43 ! * 44 !***************************************************************************** 43 45 44 46 use oh_mod … … 48 50 implicit none 49 51 50 include 'netcdf.inc'51 52 52 character(len=150) :: thefile 53 53 character(len=2) :: mm 54 54 integer :: nid,ierr,xid,yid,zid,vid,m 55 integer :: i,j,k,l 55 56 real, dimension(:), allocatable :: etaOH 56 57 … … 61 62 real, parameter :: scalehgt=7000. ! scale height in metres 62 63 63 ! Read OH fields and level heights64 !********************************65 64 66 do m=1,12 67 68 ! open netcdf file 69 write(mm,fmt='(i2.2)') m 70 ! thefile=trim(path(1))//'OH_FIELDS/'//'geos-chem.OH.2005'//mm//'01.nc' 71 thefile=trim(ohfields_path)//'OH_FIELDS/'//'geos-chem.OH.2005'//mm//'01.nc' 72 ierr=nf_open(trim(thefile),NF_NOWRITE,nid) 73 if(ierr.ne.0) then 74 write (*,*) 'The OH field at: '//thefile// ' could not be opened' 75 write (*,*) 'please copy the OH fields there, or change the path in the' 76 write (*,*) 'COMMAND namelist!' 77 write(*,*) nf_strerror(ierr) 78 stop 79 endif 65 open(unitOH,file=trim(ohfields_path) & 66 //'OH_FIELDS/OH_variables.bin',status='old', & 67 form='UNFORMATTED', iostat=ierr, convert='little_endian') 80 68 81 ! inquire about variables82 ierr=nf_inq_dimid(nid,'Lon-000',xid)83 if(ierr.ne.0) then84 write(*,*) nf_strerror(ierr)85 stop86 endif87 ierr=nf_inq_dimid(nid,'Lat-000',yid)88 if(ierr.ne.0) then89 write(*,*) nf_strerror(ierr)90 stop91 endif92 ierr=nf_inq_dimid(nid,'Alt-000',zid)93 if(ierr.ne.0) then94 write(*,*) nf_strerror(ierr)95 stop96 endif97 98 if(m.eq.1) then99 100 ! read dimension sizes101 ierr=nf_inq_dimlen(nid,xid,nxOH)102 if(ierr.ne.0) then103 write(*,*) nf_strerror(ierr)104 stop105 endif106 ierr=nf_inq_dimlen(nid,yid,nyOH)107 if(ierr.ne.0) then108 write(*,*) nf_strerror(ierr)109 stop110 endif111 ierr=nf_inq_dimlen(nid,zid,nzOH)112 if(ierr.ne.0) then113 write(*,*) nf_strerror(ierr)114 stop115 endif116 117 ! allocate variables118 allocate(lonOH(nxOH))119 allocate(latOH(nyOH))120 allocate(etaOH(nzOH))121 allocate(altOH(nzOH))122 allocate(OH_field(nxOH,nyOH,nzOH,12))123 allocate(OH_hourly(nxOH,nyOH,nzOH,2))124 125 ! read dimension variables126 ierr=nf_inq_varid(nid,'LON',xid)127 ierr=nf_get_var_real(nid,xid,lonOH)128 if(ierr.ne.0) then129 write(*,*) nf_strerror(ierr)130 stop131 endif132 ierr=nf_inq_varid(nid,'LAT',yid)133 ierr=nf_get_var_real(nid,yid,latOH)134 if(ierr.ne.0) then135 write(*,*) nf_strerror(ierr)136 stop137 endif138 ierr=nf_inq_varid(nid,'ETAC',zid)139 ierr=nf_get_var_real(nid,zid,etaOH)140 if(ierr.ne.0) then141 write(*,*) nf_strerror(ierr)142 stop143 endif144 145 ! convert eta-level to altitude (assume surface pressure of 1010 hPa)146 altOH=log(1010./(etaOH*1010.))*scalehgt147 148 endif ! m.eq.1149 150 ! read OH_field151 ierr=nf_inq_varid(nid,'CHEM-L_S__OH',vid)152 ierr=nf_get_var_real(nid,vid,OH_field(:,:,:,m))153 if(ierr.ne.0) then154 write(*,*) nf_strerror(ierr)155 stop156 endif157 158 ierr=nf_close(nid)159 160 end do161 162 deallocate(etaOH)163 164 ! Read J(O1D) photolysis rates165 !********************************166 167 ! open netcdf file168 ! thefile=trim(path(1))//'OH_FIELDS/jrate_average.nc'169 thefile=trim(ohfields_path)//'OH_FIELDS/jrate_average.nc'170 ierr=nf_open(trim(thefile),NF_NOWRITE,nid)171 69 if(ierr.ne.0) then 172 write(*,*) nf_strerror(ierr)70 write(*,*) 'Cannot read binary OH fields in ',trim(ohfields_path)//'OH_FIELDS/OH_variables.bin' 173 71 stop 174 72 endif 175 73 176 ! read dimension variables 177 ierr=nf_inq_varid(nid,'longitude',xid) 178 ierr=nf_get_var_real(nid,xid,lonjr) 179 if(ierr.ne.0) then 180 write(*,*) nf_strerror(ierr) 181 stop 182 endif 183 ierr=nf_inq_varid(nid,'latitude',yid) 184 ierr=nf_get_var_real(nid,yid,latjr) 185 if(ierr.ne.0) then 186 write(*,*) nf_strerror(ierr) 187 stop 188 endif 74 read(unitOH) nxOH 75 read(unitOH) nyOH 76 read(unitOH) nzOH 77 write(*,*) nxOH,nyOH,nzOH 189 78 190 ! read jrate_average 191 ierr=nf_inq_varid(nid,'jrate',vid)192 ierr=nf_get_var_real(nid,vid,jrate_average)193 if(ierr.ne.0) then194 write(*,*) nf_strerror(ierr)195 stop196 endif79 ! allocate variables 80 allocate(lonOH(nxOH)) 81 allocate(latOH(nyOH)) 82 allocate(etaOH(nzOH)) 83 allocate(altOH(nzOH)) 84 allocate(OH_field(nxOH,nyOH,nzOH,12)) 85 allocate(OH_hourly(nxOH,nyOH,nzOH,2)) 197 86 198 ierr=nf_close(nid) 87 read(unitOH) (lonjr(i),i=1,360) 88 read(unitOH) (latjr(i),i=1,180) 89 read(unitOH) (((jrate_average(i,j,k),i=1,360),j=1,180),k=1,12) 90 read(unitOH) (lonOH(i),i=1,nxOH) 91 read(unitOH) (latOH(i),i=1,nyOH) 92 read(unitOH) (lonOH(i),i=1,nxOH) 199 93 200 return 94 read(unitOH) (altOH(i),i=1,nzOH) 95 read(unitOH) ((((OH_field(i,j,k,l),i=1,nxOH),j=1,nyOH),k=1,nzOH),l=1,12) 96 read(unitOH) ((((OH_hourly(i,j,k,l),i=1,nxOH),j=1,nyOH),k=1,nzOH),l=1,2) 97 write(*,*) 'nzOH: ',nzOH,(altOH(i),i=1,nzOH) 201 98 202 99 end subroutine readOHfield -
src/timemanager_mpi.f90
rd8eed02 rb1e0742 808 808 call initial_cond_calc(itime+lsynctime,j) 809 809 itra1(j)=-999999999 810 !print*, 'terminated particle ',j,'for age' 810 if (verbosity.gt.0) then 811 print*, 'terminated particle ',j,'for age' 812 endif 811 813 endif 812 814 endif … … 819 821 820 822 821 ! Added by mc: counter of "unstable" particle velocity during a time scale822 ! 823 823 ! Counter of "unstable" particle velocity during a time scale 824 ! of maximumtl=20 minutes (defined in com_mod) 825 !************************************************************ 824 826 total_nan_intl=0 825 827 i_nan=i_nan+1 ! added by mc to count nan during a time of maxtl (i.e. maximum tl fixed here to 20 minutes, see com_mod)
Note: See TracChangeset
for help on using the changeset viewer.
