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

Timestamp:

Apr 3, 2020 9:38:20 AM (4 years ago)

Author:

raasch

Message:

bugfix for subroutine calls that contain the decycle_chem and decycle_salsa switches as arguments

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! bugfix for subroutine calls that contain the decycle_chem switches as arguments
+! 
+! 4481 2020-03-31 18:55:54Z maronga
 ! use statement for exchange horiz added,
 ! bugfix for call of exchange horiz 2d
@@ -2037 +2040 @@
 !--    appear under shear-free stable conditions. 
        CALL ws_init_flags_scalar(                                              &
-                  bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  decycle_chem_lr, &
-                  bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  decycle_chem_ns, &
-                  bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  decycle_chem_lr, &
-                  bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  decycle_chem_ns, &
+                  bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  ( decycle_chem_lr .AND. nxl == 0  ), &
+                  bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  ( decycle_chem_ns .AND. nyn == ny ), &
+                  bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  ( decycle_chem_lr .AND. nxr == nx ), &
+                  bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  ( decycle_chem_ns .AND. nys == 0  ), &
                   cs_advc_flags_s, .TRUE. )
     ENDIF
@@ -2829 +2832 @@
           IF ( ws_scheme_sca )  THEN
              CALL advec_s_ws( cs_advc_flags_s, chem_species(ilsp)%conc, 'kc',                      &
-                              bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  decycle_chem_lr,         &
-                              bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  decycle_chem_ns,         &
-                              bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  decycle_chem_lr,         &
-                              bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  decycle_chem_ns )
+                  bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  ( decycle_chem_lr .AND. nxl == 0  ), &
+                  bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  ( decycle_chem_ns .AND. nyn == ny ), &
+                  bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  ( decycle_chem_lr .AND. nxr == nx ), &
+                  bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  ( decycle_chem_ns .AND. nys == 0  ) )
           ELSE
              CALL advec_s_pw( chem_species(ilsp)%conc )
@@ -2939 +2942 @@
        IF ( timestep_scheme(1:5) == 'runge' )  THEN
           IF ( ws_scheme_sca )  THEN
-             CALL advec_s_ws( cs_advc_flags_s,                                                     &
-                              i,                                                                   &
-                              j,                                                                   &
-                              chem_species(ilsp)%conc,                                             &
-                              'kc',                                                                &
+             CALL advec_s_ws( cs_advc_flags_s, i, j, chem_species(ilsp)%conc, 'kc',                &
                               chem_species(ilsp)%flux_s_cs,                                        &
                               chem_species(ilsp)%diss_s_cs,                                        &
@@ -2950 +2949 @@
                               i_omp_start,                                                         &
                               tn,                                                                  &
-                              bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  decycle_chem_lr,         &
-                              bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  decycle_chem_ns,         &
-                              bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  decycle_chem_lr,         &
-                              bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  decycle_chem_ns,         &
-                              monotonic_limiter_z )
+                  bc_dirichlet_l  .OR.  bc_radiation_l  .OR.  ( decycle_chem_lr .AND. nxl == 0  ), &
+                  bc_dirichlet_n  .OR.  bc_radiation_n  .OR.  ( decycle_chem_ns .AND. nyn == ny ), &
+                  bc_dirichlet_r  .OR.  bc_radiation_r  .OR.  ( decycle_chem_lr .AND. nxr == nx ), &
+                  bc_dirichlet_s  .OR.  bc_radiation_s  .OR.  ( decycle_chem_ns .AND. nys == 0  ), &
+                  monotonic_limiter_z )
           ELSE
              CALL advec_s_pw( i, j, chem_species(ilsp)%conc )