added multiline output keyword support; propogated ions into product mixture through legacy cli (#52)

Author: markleader

source/input.f90

@@ -250,6 +250,11 @@ subroutine read_problem(fin, problem, ierr)
                 line = adjustl(line)
                 if (is_empty(line)) cycle
                 if (is_keyword(line)) then
+                    if (dsname == 'outp' .and. line(1:4) == 'outp') then
+                        ! Legacy behavior allows repeated output datasets; merge them.
+                        buffer = buffer // ' ' // trim(line)
+                        cycle
+                    end if
                     backspace(fin)
                     exit
                 else

source/input_test.pf

@@ -120,6 +120,19 @@ contains
         return
     end subroutine
 
+    @test
+    subroutine test_parse_outp_repeated_dataset_keywords
+        type(OutputDataset) :: outp
+
+        outp = parse_outp('outp short outp massf outp siunits outp trace=1e-8 outp transport')
+
+        @assertEqual(.true., outp%mass_fractions)
+        @assertEqual(.true., outp%transport)
+        @assertTrue(allocated(outp%trace))
+        @assertEqual(1.0d-8, outp%trace)
+        @assertEqual(.true., outp%siunit)
+    end subroutine
+
     @test
     subroutine test_parse_reac_custom_species
         type(ReactantInput), allocatable :: reac(:)

source/main.f90

@@ -266,7 +266,7 @@ subroutine run_thermo_problem(prob, thermo, solver, solutions, partials)
         else
             product_names = reactants%get_products(thermo)
         end if
-        products = Mixture(thermo, product_names)
+        products = Mixture(thermo, product_names, ions=prob%problem%include_ions)
 
         ! Get the loop sizes
         num_state1 = 1
@@ -535,7 +535,7 @@ subroutine run_shock_problem(prob, thermo, solver, solutions)
         else
             product_names = reactants%get_products(thermo)
         end if
-        products = Mixture(thermo, product_names)
+        products = Mixture(thermo, product_names, ions=prob%problem%include_ions)
 
         ! Get the problem flags
         if (prob%problem%shk_incident) then
@@ -767,7 +767,7 @@ subroutine run_detonation_problem(prob, thermo, solver, solutions)
         else
             product_names = reactants%get_products(thermo)
         end if
-        products = Mixture(thermo, product_names)
+        products = Mixture(thermo, product_names, ions=prob%problem%include_ions)
 
         ! Get the problem flags
         if (prob%problem%frozen) then