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

Timestamp:

Oct 17, 2019 11:29:38 AM (5 years ago)

Author:

schwenkel

Message:

Introducing module interface for boundary conditions and move module specific boundary conditions into their modules

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Moving module specific boundary conditions from time_integration to module
+! 
+! 4230 2019-09-11 13:58:14Z suehring
 ! Bugfix, initialize mean profiles also in restart runs. Also initialize 
 ! array used for Runge-Kutta tendecies in restart runs.  
@@ -380 +383 @@
     END INTERFACE chem_boundary_conds
 
+    INTERFACE chem_boundary_conditions
+       MODULE PROCEDURE chem_boundary_conditions
+    END INTERFACE chem_boundary_conditions
+
     INTERFACE chem_check_data_output
        MODULE PROCEDURE chem_check_data_output
@@ -537 +544 @@
 
 
-    PUBLIC chem_3d_data_averaging, chem_boundary_conds,                       &
+    PUBLIC chem_3d_data_averaging, chem_boundary_conds, chem_boundary_conditions, &
             chem_boundary_conds_decycle, chem_check_data_output,              &
          chem_check_data_output_pr, chem_check_parameters,                    &
@@ -799 +806 @@
  END SUBROUTINE chem_boundary_conds
 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Subroutine for boundary conditions
+!------------------------------------------------------------------------------!
+ SUBROUTINE chem_boundary_conditions
+
+    IMPLICIT NONE
+
+    INTEGER(iwp) ::  lsp             !<
+    INTEGER(iwp) ::  lsp_usr         !<
+
+!
+!--       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
+
+
+ END SUBROUTINE chem_boundary_conditions
 
 !------------------------------------------------------------------------------!