Changeset 4511 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Apr 30, 2020 12:20:40 PM (4 years ago)

Author:

raasch

Message:

chemistry decycling replaced by explicit setting of lateral boundary conditions

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (7 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! chemistry decycling replaced by explicit setting of lateral boundary conditions
+! 
+! 4508 2020-04-24 13:32:20Z raasch
 ! salsa decycling replaced by explicit setting of lateral boundary conditions
 !
@@ -242 +245 @@
 
     USE chem_modules,                                                                              &
-        ONLY:  bc_cs_t_val, chem_species, emissions_anthropogenic, emiss_read_legacy_mode,         &
-               n_matched_vars
-
-#if defined( __parallel )
-    USE chem_modules,                                                                              &
-        ONLY:  cs_name
-#endif
-
-    USE chemistry_model_mod,                                                                       &
-        ONLY:  chem_boundary_conds
+        ONLY:  bc_cs_t_val, chem_species, communicator_chem, emissions_anthropogenic,              &
+               emiss_read_legacy_mode, n_matched_vars
 
     USE control_parameters,                                                                        &
@@ -469 +464 @@
     INTEGER(iwp) ::  icc                 !< additional index for aerosol mass bins
     INTEGER(iwp) ::  ig                  !< index for salsa gases
-    INTEGER(iwp) ::  lsp                 !<
     INTEGER(iwp) ::  mid                 !< masked output running index
-#if defined( __parallel )
-    INTEGER(iwp) ::  lsp_usr             !<
     INTEGER(iwp) ::  n                   !< loop counter for chemistry species
-#endif
 
     REAL(wp) ::  dt_3d_old  !< temporary storage of timestep to be used for
@@ -675 +666 @@
            bc_q_t_val  = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
            IF ( air_chemistry )  THEN
-              DO  lsp = 1, nvar
-                 bc_cs_t_val = (  chem_species(lsp)%conc_pr_init(nzt+1)                            &
-                                - chem_species(lsp)%conc_pr_init(nzt) )                            &
+              DO  n = 1, nvar
+                 bc_cs_t_val = (  chem_species(n)%conc_pr_init(nzt+1)                              &
+                                - chem_species(n)%conc_pr_init(nzt) )                              &
                                / dzu(nzt+1)
               ENDDO
@@ -824 +815 @@
           IF ( passive_scalar )  CALL exchange_horiz( s_p, nbgp )
           IF ( air_chemistry )  THEN
-             DO  lsp = 1, nvar
-                CALL exchange_horiz( chem_species(lsp)%conc_p, nbgp )
+             DO  n = 1, nvar
+                CALL exchange_horiz( chem_species(n)%conc_p, nbgp,                                 &
+                                     alternative_communicator = communicator_chem )
              ENDDO
           ENDIF
@@ -932 +924 @@
                 IF ( air_chemistry )  THEN
                    DO  n = 1, nvar
-                      CALL exchange_horiz( chem_species(n)%conc, nbgp )
+                      CALL exchange_horiz( chem_species(n)%conc, nbgp,                             &
+                                           alternative_communicator = communicator_chem )
                    ENDDO
                 ENDIF
@@ -960 +953 @@
 !--          Set boundary conditions again after interpolation and anterpolation.
              CALL pmci_boundary_conds
-
-!
-!--          Set chemistry boundary conditions (decycling)
-             IF ( air_chemistry )  THEN
-                DO  lsp = 1, nvar
-                   lsp_usr = 1
-                   DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' )
-                      IF ( TRIM( chem_species(lsp)%name ) == TRIM( cs_name(lsp_usr) ) )  THEN
-                         CALL chem_boundary_conds( chem_species(lsp)%conc,                         &
-                                                   chem_species(lsp)%conc_pr_init )
-                      ENDIF
-                      lsp_usr = lsp_usr + 1
-                   ENDDO
-                ENDDO
-             ENDIF
 
              CALL cpu_log( log_point(60), 'nesting', 'stop' )