diff --git a/options/SPECIES/SPECIES_050 b/options/SPECIES/SPECIES_050
new file mode 100644
index 0000000000000000000000000000000000000000..b327d2fa1bc2053ab7205796ee825e9e4c1de907
--- /dev/null
+++ b/options/SPECIES/SPECIES_050
@@ -0,0 +1,27 @@
+&SPECIES_PARAMS
+ PSPECIES="MP",               ! Tracer name
+ PDECAY=-9.9,                 ! Species half life
+ PWETA_GAS=-0.9E-9,           ! Below-cloud scavenging (gases) - A (weta_gas)
+ PWETB_GAS=-9.9,              ! Below-cloud scavenging (gases) - B (wetb_gas)
+ PCRAIN_AERO=1.0,             ! Below-cloud scavenging (particles) - Crain (crain_aero)
+ PCSNOW_AERO=1.0,             ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
+ PCCN_AERO=0.05,               ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
+ PIN_AERO=0.15,                ! In-cloud scavenging (particles) - INeff (in_aero)
+ PDENSITY=1234.0,             ! Dry deposition (particles) - rho
+ PDQUER=1.0E-06,             ! Dry deposition (particles) - dquer (equivalent diameter for shape)
+ PDSIGMA=3.3,                 ! Dry deposition (particles) - dsig
+ PNDIA=1,                     ! Dry deposition (particles) - ndia
+ PDRYVEL=-9.9,                ! Alternative: dry deposition velocity
+ PRELDIFF=-9.9,               ! Dry deposition (gases) - D
+ PHENRY=-0.9E-9,              ! Dry deposition (gases) - Henrys const.
+ PF0=-9,                      ! Dry deposition (gases) - f0 (reactivity)
+ PWEIGHTMOLAR=-9.9,           ! molweight
+ POHCCONST=-0.9E-9,           ! OH Reaction rate - C [cm^3/molecule/sec]
+ POHDCONST=-9.9,              ! OH Reaction rate - D [K]
+ POHNCONST=2.0,               ! OH Reaction rate - C [cm^3/molecule/sec]
+ PSHAPE=1,                    ! 0 for sphere, 1 any shape (defined by axes PLA,PIA,PSA), 2-cylinder, 3-cube, 4-tetrahedron, 5-octahedron, 6-ellipsoid
+ PASPECTRATIO=50.0            ! Aspect ratio of cylinders: works for PSHAPE=2 only
+ PLA=1.8E-06,                 ! Longest axis in meters (Bagheri & Bonadonna 2016): only for PSHAPE=1
+ PIA=0.6E-06,                 ! Intermediate axis in meters: only for PSHAPE=1
+ PSA=0.6E-06,                 ! Smallest axis in meters: only for PSHAPE=1
+ PORIENT=0,                   !  0 for random orientation of particles, 1 for horizontal, 2 for an average between random and horizontal
\ No newline at end of file