Changeset 4273 for palm/trunk/SOURCE/salsa_mod.f90

Timestamp:

Oct 24, 2019 1:40:54 PM (4 years ago)

Author:

monakurppa

Message:

Add logical switched nesting_chem and nesting_offline_chem

File:

• palm/trunk/SOURCE/salsa_mod.f90 (modified) (15 diffs)

palm/trunk/SOURCE/salsa_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! - Rename nest_salsa to nesting_salsa
+! - Correct some errors in boundary condition flags
+! - Add a check for not trying to output gas concentrations in salsa if the
+!   chemistry module is applied
+! - Set the default value of nesting_salsa and nesting_offline_salsa to .TRUE.
+! 
+! 4272 2019-10-23 15:18:57Z schwenkel
 ! Further modularization of boundary conditions: moved boundary conditions to
 ! respective modules
@@ -454 +461 @@
     LOGICAL ::  include_emission        = .FALSE.  !< Include or not emissions
     LOGICAL ::  feedback_to_palm        = .FALSE.  !< Allow feedback due to condensation of H2O
-    LOGICAL ::  nest_salsa              = .FALSE.  !< Apply nesting for salsa
-    LOGICAL ::  nesting_offline_salsa   = .FALSE.  !< Apply offline nesting for salsa
+    LOGICAL ::  nesting_salsa           = .TRUE.   !< Apply nesting for salsa
+    LOGICAL ::  nesting_offline_salsa   = .TRUE.   !< Apply offline nesting for salsa
     LOGICAL ::  no_insoluble            = .FALSE.  !< Exclude insoluble chemical components
     LOGICAL ::  read_restart_data_salsa = .FALSE.  !< Read restart data for salsa
@@ -652 +659 @@
        CHARACTER(LEN=15) ::  char_t = 'ls_forcing_top_'   !< leading substring at top boundary
 
-       CHARACTER(LEN=5), DIMENSION(1:ngases_salsa) ::  gas_name = (/'h2so4','hno3 ','nh3  ','ocnv ','ocsv '/)
+       CHARACTER(LEN=5), DIMENSION(1:ngases_salsa) ::  gas_name = (/'H2SO4','HNO3 ','NH3  ','OCNV ','OCSV '/)
 
        CHARACTER(LEN=25),  DIMENSION(:), ALLOCATABLE ::  cc_name    !< chemical component name
@@ -973 +980 @@
            init_aerosol_type,     &
            init_gases_type,       &
-           nest_salsa,            &
+           nesting_salsa,         &
            nesting_offline_salsa, &
            salsa_gases_from_chem, &
@@ -1037 +1044 @@
                                      n_lognorm,                                &
                                      nbin,                                     &
-                                     nest_salsa,                               &
+                                     nesting_salsa,                            &
                                      nesting_offline_salsa,                    &
                                      nf2a,                                     &
@@ -1102 +1109 @@
 
     USE control_parameters,                                                                        &
-        ONLY:  humidity, initializing_actions
+        ONLY:  child_domain, humidity, initializing_actions, nesting_offline
 
     IMPLICIT NONE
@@ -1111 +1118 @@
        WRITE( message_string, * ) 'salsa = ', salsa, ' is not allowed with humidity = ', humidity
        CALL message( 'salsa_check_parameters', 'PA0594', 1, 2, 0, 6, 0 )
+    ENDIF
+!
+!-- For nested runs, explicitly set nesting boundary conditions.
+    IF ( nesting_salsa  .AND. child_domain )  bc_salsa_t = 'nested'
+!
+!-- Set boundary conditions also in case the model is offline-nested in larger-scale models.
+    IF ( nesting_offline )  THEN
+       IF ( nesting_offline_salsa )  THEN
+          bc_salsa_t = 'nesting_offline'
+       ELSE
+          bc_salsa_t = 'neumann'
+       ENDIF
     ENDIF
 !
@@ -1128 +1147 @@
     ELSEIF ( bc_salsa_t == 'neumann' )  THEN
        ibc_salsa_t = 1
-    ELSEIF ( bc_salsa_t == 'nested' )  THEN
+    ELSEIF ( bc_salsa_t == 'initial_gradient' )  THEN
        ibc_salsa_t = 2
+    ELSEIF ( bc_salsa_t == 'nested'  .OR.  bc_salsa_t == 'nesting_offline' )  THEN
+       ibc_salsa_t = 3
     ELSE
        message_string = 'unknown boundary condition: bc_salsa_t = "' // TRIM( bc_salsa_t ) // '"'
        CALL message( 'salsa_check_parameters', 'PA0596', 1, 2, 0, 6, 0 )
-    ENDIF
-!
-!-- If nest_salsa = .F., set top boundary to dirichlet
-    IF ( .NOT. nest_salsa  .AND.  ibc_salsa_t == 2  )  THEN
-       ibc_salsa_t = 0
-       bc_salsa_t = 'dirichlet'
     ENDIF
 !
@@ -1246 +1261 @@
        WRITE( io, 19 )
     ENDIF
-    IF ( nest_salsa )  WRITE( io, 20 )  nest_salsa
-    WRITE( io, 21 ) salsa_emission_mode
+    IF ( nesting_salsa )  WRITE( io, 20 )  nesting_salsa
+    IF ( nesting_offline_salsa )  WRITE( io, 21 )  nesting_offline_salsa
+    WRITE( io, 22 ) salsa_emission_mode
     IF ( salsa_emission_mode == 'uniform' )  THEN
-       WRITE( io, 22 ) surface_aerosol_flux, aerosol_flux_dpg, aerosol_flux_sigmag,                &
+       WRITE( io, 23 ) surface_aerosol_flux, aerosol_flux_dpg, aerosol_flux_sigmag,                &
                        aerosol_flux_mass_fracs_a
     ENDIF
     IF ( SUM( aerosol_flux_mass_fracs_b ) > 0.0_wp  .OR. salsa_emission_mode == 'read_from_file' ) &
     THEN
-       WRITE( io, 23 )
+       WRITE( io, 24 )
     ENDIF
 
@@ -1297 +1313 @@
 19   FORMAT (/'      Size distribution read from a file.')
 20   FORMAT (/'   Nesting for salsa variables: ', L1 )
-21   FORMAT (/'   Emissions: salsa_emission_mode = ', A )
-22   FORMAT (/'      surface_aerosol_flux = ', ES12.4E3, ' #/m**2/s', /                            &
+21   FORMAT (/'   Offline nesting for salsa variables: ', L1 )
+22   FORMAT (/'   Emissions: salsa_emission_mode = ', A )
+23   FORMAT (/'      surface_aerosol_flux = ', ES12.4E3, ' #/m**2/s', /                            &
               '      aerosol_flux_dpg     =  ', 7(F7.3), ' (m)', /                                 &
               '      aerosol_flux_sigmag  =  ', 7(F7.2), /                                         &
               '      aerosol_mass_fracs_a =  ', 7(ES12.4E3) )
-23   FORMAT (/'      (currently all emissions are soluble!)')
+24   FORMAT (/'      (currently all emissions are soluble!)')
 
  END SUBROUTINE salsa_header
@@ -2209 +2226 @@
 !
 !--       Read vertical profiles of gases:
-          CALL get_variable( id_dyn, 'init_atmosphere_h2so4', salsa_gas(1)%init(nzb+1:nzt) )
-          CALL get_variable( id_dyn, 'init_atmosphere_hno3',  salsa_gas(2)%init(nzb+1:nzt) )
-          CALL get_variable( id_dyn, 'init_atmosphere_nh3',   salsa_gas(3)%init(nzb+1:nzt) )
-          CALL get_variable( id_dyn, 'init_atmosphere_ocnv',  salsa_gas(4)%init(nzb+1:nzt) )
-          CALL get_variable( id_dyn, 'init_atmosphere_ocsv',  salsa_gas(5)%init(nzb+1:nzt) )
+          CALL get_variable( id_dyn, 'init_atmosphere_H2SO4', salsa_gas(1)%init(nzb+1:nzt) )
+          CALL get_variable( id_dyn, 'init_atmosphere_HNO3',  salsa_gas(2)%init(nzb+1:nzt) )
+          CALL get_variable( id_dyn, 'init_atmosphere_NH3',   salsa_gas(3)%init(nzb+1:nzt) )
+          CALL get_variable( id_dyn, 'init_atmosphere_OCNV',  salsa_gas(4)%init(nzb+1:nzt) )
+          CALL get_variable( id_dyn, 'init_atmosphere_OCSV',  salsa_gas(5)%init(nzb+1:nzt) )
 !
 !--       Set Neumann top and surface boundary condition for initial + initialise concentrations
@@ -8146 +8163 @@
 
     IF ( time_since_reference_point >= skip_time_do_salsa )  THEN
-
 !
 !--    Surface conditions:
@@ -8247 +8263 @@
           ENDIF
 
-       ELSEIF ( ibc_salsa_t == 2 )  THEN   ! nested
+       ELSEIF ( ibc_salsa_t == 2 )  THEN   ! Initial gradient
 
           DO  ib = 1, nbins_aerosol
@@ -8358 +8374 @@
     flag = 0.0_wp
 !
-!-- Skip input if forcing from larger-scale models is applied.
+!-- Skip input if forcing from a larger-scale models is applied.
     IF ( nesting_offline  .AND.  nesting_offline_salsa )  RETURN
 !
@@ -9778 +9794 @@
 
           CASE ( 'g_H2SO4', 'g_HNO3', 'g_NH3', 'g_OCNV',  'g_OCSV' )
-             IF (  salsa_gases_from_chem )  THEN
+             IF (  air_chemistry )  THEN
                 message_string = 'gases are imported from the chemistry module and thus output '// &
                                  'of "' // TRIM( var ) // '" is not allowed'