Changeset 4268 for palm/trunk/SOURCE/time_integration.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/time_integration.f90 (modified) (9 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Removing module specific boundary conditions an put them into their modules
+! 
+! 4227 2019-09-10 18:04:34Z gronemeier
 ! implement new palm_date_time_mod
 ! 
@@ -254 +257 @@
 
     USE module_interface,                                                                          &
-        ONLY:  module_interface_actions, module_interface_swap_timelevel
+        ONLY:  module_interface_actions, module_interface_swap_timelevel,                          &
+               module_interface_boundary_conditions
 
     USE multi_agent_system_mod,                                                                    &
@@ -583 +587 @@
 !--    Execute all other module actions routunes
        CALL module_interface_actions( 'before_timestep' )
-       
+
 !
 !--    Start of intermediate step loop
@@ -750 +754 @@
 !--       velocities at the outflow in case of a non-cyclic lateral wall)
           CALL boundary_conds
-
-!
-!--       Boundary conditions for prognostic quantitites of other modules:
-!--       Here, only decycling is carried out
-          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_p,                          &
-                                                chem_species(lsp)%conc_pr_init )
-                   ENDIF
-                   lsp_usr = lsp_usr + 1
-                ENDDO
-             ENDDO
-
-          ENDIF
-
-          IF ( salsa  .AND.  time_since_reference_point >= skip_time_do_salsa )  THEN
-
-             DO  ib = 1, nbins_aerosol
-                CALL salsa_boundary_conds( aerosol_number(ib)%conc_p, aerosol_number(ib)%init )
-                DO  ic = 1, ncomponents_mass
-                   icc = ( ic - 1 ) * nbins_aerosol + ib
-                   CALL salsa_boundary_conds( aerosol_mass(icc)%conc_p, aerosol_mass(icc)%init )
-                ENDDO
-             ENDDO
-             IF ( .NOT. salsa_gases_from_chem )  THEN
-                DO  ig = 1, ngases_salsa
-                   CALL salsa_boundary_conds( salsa_gas(ig)%conc_p, salsa_gas(ig)%init )
-                ENDDO
-             ENDIF
-
-          ENDIF
-
+!
+!--       Boundary conditions for module-specific variables
+          CALL module_interface_boundary_conditions
 !
 !--       Incrementing timestep counter
@@ -1009 +980 @@
 
              ELSE
-!                
+!
 !--             Mass (volume) flux correction to ensure global mass conservation for child domains.
                 IF ( child_domain )  THEN
@@ -1018 +989 @@
                    ENDIF
                 ENDIF
-                
+
                 CALL pres
 
@@ -1084 +1055 @@
                 CALL lsm_energy_balance( .FALSE., 3 )
                 CALL lsm_soil_model( .FALSE., 3, .TRUE. )
-                
+
 !
 !--             At the end, set boundary conditons for potential temperature
@@ -1091 +1062 @@
                 CALL lsm_boundary_condition
 
-                
+
                 CALL cpu_log( log_point(54), 'land_surface', 'stop' )
              ENDIF
@@ -1099 +1070 @@
              IF (urban_surface) THEN
                 CALL cpu_log( log_point(74), 'urban_surface', 'start' )
-                
+
                 CALL usm_surface_energy_balance( .FALSE. )
                 IF ( usm_material_model )  THEN