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

Timestamp:

May 24, 2019 11:32:38 AM (5 years ago)

Author:

suehring

Message:

chemistry: perform basic checks only when anthropenic emissions are switched on; virtual flights: allow arbitrary start/end positions also in return mode; bugfix in 2d data output

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! in subroutine chem_parin check emiss_lod / mod_emis only
+! when emissions_anthropogenic is activated in namelist (E.C. Chan)
+! 
+! 3968 2019-05-13 11:04:01Z suehring
 ! - added "emiss_lod" which serves the same function as "mode_emis"
 !   both will be synchronized with emiss_lod having pirority over
@@ -2282 +2286 @@
 
  20 CONTINUE
+
+!
+!-- synchronize emiss_lod and mod_emis only if emissions_anthropogenic
+!-- is activated in the namelist.  Otherwise their values are "don't care"
+    IF ( emissions_anthropogenic )  THEN
 !
 !--    check for emission mode for chem species
 
-    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 )
+       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
-    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
@@ -2313 +2323 @@
 !-- 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
+!-- note that 
+
+       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
+          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
+
+    ENDIF  ! if emissions_anthropengic
 
     t_steps = my_steps