From 7fdcc5687efd5a9cd30814cd239fc3d457e05f35 Mon Sep 17 00:00:00 2001 From: "L. Diana Amorim" Date: Wed, 18 Dec 2019 16:28:41 -0800 Subject: Added print line to see absent interactions The search for ionization should only happen for particles other than elementary --- Python/pywarpx/picmi.py | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) (limited to 'Python/pywarpx/picmi.py') diff --git a/Python/pywarpx/picmi.py b/Python/pywarpx/picmi.py index 13774bff2..12accdee7 100644 --- a/Python/pywarpx/picmi.py +++ b/Python/pywarpx/picmi.py @@ -55,15 +55,15 @@ class Species(picmistandard.PICMI_Species): self.element = element if self.mass is None: self.mass = element.mass*periodictable.constants.atomic_mass_constant - - if 'ionization' in self.interactions: - if self.interactions[1]=='ADK': - species.do_field_ionization=1 - species.physical_element=self.particle_type - species.ionization_product_species = self.interactions[2] - species.ionization_initial_level = self.charge_state - else: - raise Exception('WarpX only supports the ionization model of ADK') + print(self.name,self.interactions) + if 'ionization' in self.interactions: + if self.interactions[1]=='ADK': + species.do_field_ionization=1 + species.physical_element=self.particle_type + species.ionization_product_species = self.interactions[2] + species.ionization_initial_level = self.charge_state + else: + raise Exception('WarpX only supports the ionization model of ADK') def initialize_inputs(self, layout, initialize_self_fields=False): self.species_number = pywarpx.particles.nspecies -- cgit v1.2.3