Changeset 3298 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Oct 2, 2018 12:21:11 PM (6 years ago)

Author:

kanani

Message:

Merge chemistry branch at r3297 to trunk

Location:

palm/trunk/SOURCE

Files:

• . (modified) (1 prop)
• prognostic_equations.f90 (modified) (6 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Code added for decycling chemistry (basit)
+!
+! 3294 2018-10-01 02:37:10Z raasch
 ! changes concerning modularization of ocean option
 ! 
@@ -335 +338 @@
 
     USE chemistry_model_mod,                                                   &
-        ONLY:  chem_integrate, chem_prognostic_equations,                      &
-               chem_species, nspec, nvar, spc_names
+        ONLY:  chem_integrate, chem_species,  chem_prognostic_equations,       &
+               nspec, nvar, spc_names, chem_boundary_conds
            
     USE chem_modules,                                                          &
@@ -343 +346 @@
     USE chem_photolysis_mod,                                                   &
         ONLY:  photolysis_control
+
+    USE chem_modules,                                                          &
+        ONLY:  call_chem_at_all_substeps, chem_gasphase_on, cs_name
 
     USE control_parameters,                                                    &
@@ -456 +462 @@
 
     LOGICAL      ::  loop_start          !<
-    INTEGER      ::  n, lsp              !< lsp running index for chem spcs
+    INTEGER(iwp) ::  n
+    INTEGER(iwp) :: lsp
+    INTEGER(iwp) :: lsp_usr              !< lsp running index for chem spcs
 
 
@@ -468 +476 @@
 !-- concentrations of chemical species                                    
     IF ( air_chemistry )  THEN
+       lsp_usr = 1
 !
 !--    If required, calculate photolysis frequencies - 
@@ -492 +501 @@
 !--    Loop over chemical species       
        CALL cpu_log( log_point_s(84), 'chemistry exch-horiz ', 'start' )
-       DO  n = 1, nspec
-          CALL exchange_horiz( chem_species(n)%conc, nbgp )     
+       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), 'chemistry exch-horiz ', 'stop' )