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

Timestamp:

Jan 7, 2019 3:29:59 PM (6 years ago)

Author:

forkel

Message:

Checks added for chemistry mechanism, parameter chem_mechanism added

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Checks added for chemistry mechanism, parameter chem_mechanism added (basit)
+! 
+! 
+! 3646 2018-12-28 17:58:49Z kanani
 ! Bugfix: use time_since_reference_point instead of simulated_time (relevant
 ! when using wall/soil spinup)
@@ -1120 +1124 @@
     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)
 
 !
@@ -1153 +1164 @@
        CALL message( 'chem_check_parameters', 'CM0426', 1, 2, 0, 6, 0 )
     ENDIF
+!-- check for chemical mechanism used 
+    IF (chem_mechanism /= trim(cs_mech) )  THEN
+       message_string = 'Incorrect chemical mechanism selected, please check spelling and/or chem_gasphase_mod'
+       CALL message( 'chem_check_parameters', 'CM0462', 1, 2, 0, 6, 0 )
+    ENDIF
+
 
 !---------------------
 !-- chem_check_parameters is called before the array chem_species is allocated!
 !-- temporary switch of this part of the check
-    RETURN
+!    RETURN                !bK commented
 !---------------------
 
@@ -1173 +1190 @@
        IF ( .not. found ) THEN
           message_string = 'Unused/incorrect input for initial surface value: ' // TRIM(cs_name(lsp_usr))
-          CALL message( 'chem_check_parameters', 'CM0427', 0, 1, 0, 6, 0 )
+          CALL message( 'chem_check_parameters', 'CM0427', 1, 2, 0, 6, 0 )
        ENDIF
        lsp_usr = lsp_usr + 1
@@ -1193 +1210 @@
           message_string = 'Unused/incorrect input of chemical species for surface emission fluxes: '  &
                             // TRIM(surface_csflux_name(lsp_usr))
-          CALL message( 'chem_check_parameters', 'CM0428', 0, 1, 0, 6, 0 )
+          CALL message( 'chem_check_parameters', 'CM0428', 1, 2, 0, 6, 0 )
        ENDIF
        lsp_usr = lsp_usr + 1
@@ -1590 +1607 @@
     CHARACTER (LEN=80)       :: docsflux_chr
     CHARACTER (LEN=80)       :: docsinit_chr
-
-!
+    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)
+
 !-- Write chemistry model  header
     WRITE( io, 1 )
@@ -1658 +1682 @@
 !-- number of variable and fix chemical species and number of reactions
     cs_fixed = nspec - nvar
+
+    WRITE ( io, * ) '   --> Chemical Mechanism        : ', cs_mech 
     WRITE ( io, * ) '   --> Chemical species, variable: ', nvar
     WRITE ( io, * ) '   --> Chemical species, fixed   : ', cs_fixed
@@ -1960 +1986 @@
     REAL(wp), PARAMETER              ::  t_std   = 273.15_wp                 !< standard pressure (Pa)
     REAL(wp), PARAMETER              ::  p_std   = 101325.0_wp               !< standard pressure (Pa)
-    REAL(wp), PARAMETER              ::  vmolcm  = 22.414e3_wp               !< Mole volume (22.414 l) in cm^{-3}
+    REAL(wp), PARAMETER              ::  vmolcm  = 22.414e3_wp               !< Mole volume (22.414 l) in cm^3
     REAL(wp), PARAMETER              ::  xna     = 6.022e23_wp               !< Avogadro number (molecules/mol)
 
@@ -2076 +2102 @@
          chem_debug2,                      &
          chem_gasphase_on,                 &
+         chem_mechanism,                   &          
          cs_heights,                       &
          cs_name,                          &
@@ -2560 +2587 @@
     INTEGER(iwp) ::  m    !<
     INTEGER(iwp) ::  lpr  !< running index chem spcs
+    LOGICAL      ::  ts_calc = .TRUE.
+
     !    REAL(wp),                                                                                      &
     !    DIMENSION(dots_num_palm+1:dots_max) ::                                                         &
@@ -2589 +2618 @@
           ENDDO
        ENDDO
-    ELSEIF ( mode == 'time_series' )  THEN
+    ELSEIF ( mode == 'time_series' .and. ts_calc )  THEN
        CALL location_message( 'Time series not calculated for chemistry', .TRUE. )
+       ts_calc = .FALSE.
     ENDIF