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

Timestamp:

Apr 12, 2019 8:47:41 AM (5 years ago)

Author:

schwenkel

Message:

bugfix for chemistry_model_mod via introducing module_interface_exchange_horiz

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -522 +522 @@
        MODULE PROCEDURE chem_non_transport_physics_ij
     END INTERFACE chem_non_transport_physics
+    
+    INTERFACE chem_exchange_horiz
+       MODULE PROCEDURE chem_exchange_horiz
+    END INTERFACE chem_exchange_horiz    
 
     INTERFACE chem_prognostic_equations
@@ -622 +626 @@
          chem_actions, chem_prognostic_equations, chem_rrd_local,             &
          chem_statistics, chem_swap_timelevel, chem_wrd_local, chem_depo,     &
-         chem_non_transport_physics
+         chem_non_transport_physics, chem_exchange_horiz
 
  CONTAINS
@@ -915 +919 @@
     DO  lsp = 1, nspec
 
-       CALL exchange_horiz( chem_species(lsp)%conc_p, nbgp )
+       CALL exchange_horiz( chem_species(lsp)%conc, nbgp )
        lsp_usr = 1
        DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' )
@@ -2460 +2464 @@
 
  END SUBROUTINE chem_non_transport_physics_ij
+ 
+ SUBROUTINE chem_exchange_horiz
+
+
+    INTEGER(iwp) ::  i  !< grid index in x-direction
+    INTEGER(iwp) ::  j  !< grid index in y-direction
+
+    i = i + j
+! !
+! !-- Loop over chemical species       
+!     CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'start' )
+!     DO  lsp = 1, nspec
+!        CALL exchange_horiz( chem_species(lsp)%conc, nbgp )    
+!        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
+!     CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'stop' )
+
+
+ END SUBROUTINE chem_exchange_horiz
+