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

Timestamp:

May 13, 2019 11:04:01 AM (6 years ago)

Author:

suehring

Message:

Updates from chemistriy branched merged into trunk: code cleaning and formatting, code structure optimizations

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! - added "emiss_lod" which serves the same function as "mode_emis"
+!   both will be synchronized with emiss_lod having pirority over
+!   mode_emis (see informational messages)
+! - modified existing error message and introduced new informational messages
+!    - CM0436 - now also applies to invalid LOD settings
+!    - CM0463 - emiss_lod take precedence in case of conflict with mod_emis
+!    - CM0464 - emiss_lod valued assigned based on mode_emis if undefined
+! 
+! 3930 2019-04-24 14:57:18Z forkel
 ! Changed chem_non_transport_physics to chem_non_advective_processes
 ! 
@@ -1682 +1691 @@
 !
 !-- Emission mode info
-    IF ( mode_emis == "DEFAULT" )  THEN 
-       WRITE( io, 5 ) 
-    ELSEIF ( mode_emis == "PARAMETERIZED" )  THEN 
-       WRITE( io, 6 )
-    ELSEIF ( mode_emis == "PRE-PROCESSED" )  THEN 
-       WRITE( io, 7 )
-    ENDIF 
+!-- At the moment the evaluation is done with both emiss_lod and mode_emis
+!-- but once salsa has been migrated to emiss_lod the .OR. mode_emis
+!-- conditions can be removed (ecc 20190513)
+    IF     ( (emiss_lod == 1) .OR. (mode_emis == 'DEFAULT') )        THEN
+        WRITE ( io, 5 )
+    ELSEIF ( (emiss_lod == 0) .OR. (mode_emis == 'PARAMETERIZED') )  THEN 
+        WRITE ( io, 6 )
+    ELSEIF ( (emiss_lod == 2) .OR. (mode_emis == 'PRE-PROCESSED') )  THEN 
+        WRITE ( io, 7 )
+    ENDIF
 !
 !-- Photolysis scheme info
@@ -2209 +2221 @@
          decycle_method,                   &
          deposition_dry,                   &
-         emissions_anthropogenic,          &  
+         emissions_anthropogenic,          &
+         emiss_lod,                        &
          emiss_factor_main,                &
          emiss_factor_side,                &                      
@@ -2271 +2284 @@
 !
 !--    check for emission mode for chem species
-    IF ( (mode_emis /= 'PARAMETERIZED')  .AND. ( mode_emis /= 'DEFAULT' ) .AND. ( mode_emis /= 'PRE-PROCESSED'  ) )  THEN
-       message_string = 'Incorrect mode_emiss  option select. Please check spelling'
-       CALL message( 'chem_check_parameters', 'CM0436', 1, 2, 0, 6, 0 )
+
+    IF ( emiss_lod < 0 )  THEN   !- if LOD not defined in namelist
+       IF ( ( mode_emis /= 'PARAMETERIZED'  )    .AND.      &
+            ( mode_emis /= 'DEFAULT'        )    .AND.      &
+            ( mode_emis /= 'PRE-PROCESSED'  ) )  THEN
+          message_string = 'Incorrect mode_emiss  option select. Please check spelling'
+          CALL message( 'chem_check_parameters', 'CM0436', 1, 2, 0, 6, 0 )
+       ENDIF
+    ELSE
+       IF ( ( emiss_lod /= 0 )    .AND.         &
+            ( emiss_lod /= 1 )    .AND.         &
+            ( emiss_lod /= 2 ) )  THEN
+          message_string = 'Invalid value for emiss_lod (0, 1, or 2)'
+          CALL message( 'chem_check_parameters', 'CM0436', 1, 2, 0, 6, 0 )
+       ENDIF
+    ENDIF
+
+! for reference (ecc)
+!    IF ( (mode_emis /= 'PARAMETERIZED')  .AND. ( mode_emis /= 'DEFAULT' ) .AND. ( mode_emis /= 'PRE-PROCESSED'  ) )  THEN
+!       message_string = 'Incorrect mode_emiss  option select. Please check spelling'
+!       CALL message( 'chem_check_parameters', 'CM0436', 1, 2, 0, 6, 0 )
+!    ENDIF
+
+!
+!-- conflict resolution for emiss_lod and mode_emis
+!-- 1) if emiss_lod is defined, have mode_emis assume same setting as emiss_lod
+!-- 2) if emiss_lod it not defined, have emiss_lod assuem same setting as mode_emis
+!-- this check is in place to retain backward compatibility with salsa until the
+!-- code is migrated completed to emiss_lod
+    IF  ( emiss_lod >= 0 ) THEN
+       SELECT CASE  ( emiss_lod )
+          CASE (0)  !- parameterized mode
+             mode_emis = 'PARAMETERIZED'
+          CASE (1)  !- default mode
+             mode_emis = 'DEFAULT'
+          CASE (2)  !- preprocessed mode
+             mode_emis = 'PRE-PROCESSED'
+       END SELECT
+       
+       message_string = 'Synchronizing mode_emis to defined emiss_lod'               //  &
+                        CHAR(10)  //  '          '                                   //  &
+                        'NOTE - mode_emis will be depreciated in future releases'    //  &
+                        CHAR(10)  //  '          '                                   //  &
+                        'please use emiss_lod to define emission mode'
+
+       CALL message ( 'parin_chem', 'CM0463', 0, 0, 0, 6, 0 )
+    ELSE  ! if emiss_lod is not set
+       SELECT CASE ( mode_emis )
+          CASE ('PARAMETERIZED')
+             emiss_lod = 0
+          CASE ('DEFAULT')
+             emiss_lod = 1
+          CASE ('PRE-PROCESSED')
+             emiss_lod = 2
+       END SELECT
+
+       message_string = 'emiss_lod undefined.  Using existing mod_emiss setting'     //  &
+                        CHAR(10)  //  '          '                                   //  &
+                        'NOTE - mode_emis will be depreciated in future releases'    //  &
+                        CHAR(10)  //  '          '                                   //  &
+                        '       please use emiss_lod to define emission mode'
+
+       CALL message ( 'parin_chem', 'CM0464', 0, 0, 0, 6, 0 )
     ENDIF