Changeset 3780 for palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/src/create_kpp_module.C

Timestamp:

Mar 5, 2019 11:19:45 AM (6 years ago)

Author:

forkel

Message:

removed read from unit 10 in chemistry_model_mod.f90, added get_mechanismname

File:

• palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/src/create_kpp_module.C (modified) (4 diffs)

palm/trunk/UTIL/chemistry/gasphase_preproc/kpp4palm/src/create_kpp_module.C

@@ -17 +17 @@
 //-----------------
 //$Id: create_kpp_module.C 3453 2018-10-30 13:21:51Z forkel $
+// Added create_set_cs and cs_mech and get_mechanismname in module_header (05.03.2019, forkel)
+//
 // exclude kco_compress from handling by global_variables2vector (30.10.2018, forkel)
 //
@@ -191 +193 @@
 
    generate_module_header();
+
+// create_set_cs
+   create_set_cs();
 
 // Create kpp_integrate subroutine (chem_gasphase_integrate) for skalar and vector mode
@@ -458 +463 @@
 }
 
+void create_kpp_module::create_set_cs() {
+   fortran_file          kppi;          
+   vector<Vvar>::iterator               iv;
+   string                               xline;
+     
+   string                          buf;
+   ifstream                        in;
+   program_line                    line;
+
+   kppi.set_name("get_mechanismname");
+   kppi.add_line("SUBROUTINE get_mechanismname                                        ");
+   kppi.add_line("                                                                    ");
+   kppi.add_line("  IMPLICIT NONE                                                     ");
+// Read mechanism from set_cm
+// Tis got an own own subroutine to aviod being called at each timestep
+
+   in.open("set_cm");
+   if( !in ) {
+      cout << "cannot open " << endl; my_abort("set_cm");
+   }
+
+   while ( 1 ) {
+     getline (in, buf);
+     if( in.eof() ) break;
+     if( in.bad() ) my_abort("ERROR_READ_4");
+     line.set_line(buf);
+     kppi.add_line(line);
+   }
+   in.close();
+
+   kppi.add_line("                                                                    ");
+   kppi.add_line("  return                                                            ");
+   kppi.add_line("END SUBROUTINE get_mechanismname                                    ");
+   kppi.add_line("                                                                    ");
+   kpp_subroutines.push_back(kppi);
+
+   return;
+}
+
+
 void create_kpp_module::create_kpp_integrate() {
    fortran_file          kppi;
    vector<Vvar>::iterator               iv;
-   string                               line;
+   string                               xline;
+
 
    kppi.set_name("chem_gasphase_integrate");
@@ -598 +644 @@
    kppi.add_line("                                                                    ");
    kppi.add_line("  return                                                            ");
-   kppi.add_line("END SUBROUTINE chem_gasphase_integrate                                        ");
+   kppi.add_line("END SUBROUTINE chem_gasphase_integrate                              ");
 
 //   e5_subroutines.push_back(kppi);