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

Timestamp:

Oct 23, 2019 3:18:57 PM (5 years ago)

Author:

schwenkel

Message:

further modularization of boundary conditions: moving boundary conditions to their modules

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Further modularization of boundary conditions: moved boundary conditions to
+! respective modules
+!
+! 4268 2019-10-17 11:29:38Z schwenkel
 ! Moving module specific boundary conditions from time_integration to module
 ! 
@@ -819 +823 @@
 
 !
-!--       Boundary conditions for prognostic quantitites of other modules:
-!--       Here, only decycling is carried out
-
-          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_p,                          &
-                                             chem_species(lsp)%conc_pr_init )
-                ENDIF
-                lsp_usr = lsp_usr + 1
-             ENDDO
-          ENDDO
+!-- Top/bottom boundary conditions for chemical species
+    CALL chem_boundary_conds( 'set_bc_bottomtop' )
+!
+!-- Lateral boundary conditions for chemical species
+    CALL chem_boundary_conds( 'set_bc_lateral' )
+
+!
+!--  Boundary conditions for prognostic quantitites of other modules:
+!--  Here, only decycling is carried out
+     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_p,                          &
+                                        chem_species(lsp)%conc_pr_init )
+           ENDIF
+           lsp_usr = lsp_usr + 1
+        ENDDO
+     ENDDO