particles can only be made from evolving gases

Nicholaswogan · May 31, 2024 · 023ce56 · 023ce56
1 parent fc66611
commit 023ce56
Showing 1 changed file with 9 additions and 3 deletions.
diff --git a/src/input/photochem_input_read.f90 b/src/input/photochem_input_read.f90
@@ -393,14 +393,20 @@ subroutine get_rxmechanism(mapping, infile, dat, var, err)
         if (dat%particle_formation_method(i) == CondensingParticle) then
           ! if a condensing molecule
           ind = findloc(dat%species_names,dat%particle_gas_phase(i))
-          if (ind(1) /= 0) then
-            dat%particle_gas_phase_ind(i) = ind(1)
-          else
+          if (ind(1) == 0) then
             err = "IOError: particle "//trim(dat%particle_names(i))// &
                   " can not be made from "//trim(dat%particle_gas_phase(i))// &
                   " because "//trim(dat%particle_gas_phase(i))//" is not a gas"// &
                   " in the model."
             return
+          elseif (ind(1) > dat%nq .or. ind(1) < dat%ng_1) then
+            err = "Particle "//trim(dat%particle_names(i))// &
+                  " can not be made from "//trim(dat%particle_gas_phase(i))// &
+                  " because "//trim(dat%particle_gas_phase(i))//" is"// &
+                  " is a particle, short-lived species, or background gas."
+            return
+          else
+            dat%particle_gas_phase_ind(i) = ind(1)
           endif
         endif
       enddo