Changeset 3780 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

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

Author:

forkel

Message:

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

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Removed READ from unit 10, added CALL get_mechanismname
+! 
+! 
+! 3767 2019-02-27 08:18:02Z raasch
 ! unused variable for file index removed from rrd-subroutines parameter list
 ! 
@@ -278 +282 @@
 
     USE control_parameters,                                                                        &
-         ONLY:  bc_lr_cyc, bc_ns_cyc, dt_3d, humidity, initializing_actions, message_string,        &
-         omega, tsc, intermediate_timestep_count, intermediate_timestep_count_max,           &
+         ONLY:  bc_lr_cyc, bc_ns_cyc, dt_3d, humidity, initializing_actions, message_string,       &
+         omega, tsc, intermediate_timestep_count, intermediate_timestep_count_max,                 &
          max_pr_user, timestep_scheme, use_prescribed_profile_data, ws_scheme_sca         
 
@@ -286 +290 @@
 
     USE chem_gasphase_mod,                                                                         &
-         ONLY:  nspec, spc_names, nkppctrl, nmaxfixsteps, t_steps, chem_gasphase_integrate,         &
-         vl_dim, nvar, nreact,  atol, rtol, nphot, phot_names
+         ONLY:  atol, chem_gasphase_integrate, get_mechanismname, nkppctrl, nmaxfixsteps,          &
+         nphot, nreact, nspec, nvar, phot_names, rtol, spc_names, t_steps, vl_dim
 
     USE chem_modules
@@ -1149 +1153 @@
     INTEGER (iwp) ::  lsp_usr      !< running index for user defined chem spcs 
     INTEGER (iwp) ::  lsp          !< running index for chem spcs.
-    CHARACTER (LEN=30)       ::  cs_mech,a1,b1,string
-
-
-    OPEN (10,FILE="chem_gasphase_mod.f90")   !get the chem_mechanism name from the file.
-    READ (10, 100) a1,b1,string
-    cs_mech = trim(string(16:))
- 100    FORMAT(a)
-        CLOSE(10)
-
 !
 !-- check for chemical reactions status
@@ -1190 +1185 @@
     ENDIF
 !-- check for chemical mechanism used 
+    CALL get_mechanismname
+    WRITE(06,*) 'cs_mech ',cs_mech
+    WRITE(06,*) 'chem_mechanism ', chem_mechanism
     IF (chem_mechanism /= trim(cs_mech) )  THEN
-       message_string = 'Incorrect chemical mechanism selected, please check spelling and/or chem_gasphase_mod'
+       message_string = 'Incorrect chemistry mechanism selected, check spelling in namelist and/or chem_gasphase_mod'
        CALL message( 'chem_check_parameters', 'CM0462', 1, 2, 0, 6, 0 )
     ENDIF