Changeset 2773 for palm/trunk/SOURCE/parin.f90

Timestamp:

Jan 30, 2018 2:12:54 PM (6 years ago)

Author:

suehring

Message:

Nesting for chemical species implemented; Bugfix passive scalar boundary condition after anterpolation; Timeseries output of surface temperature; Enable initialization of 3D topography (was commented out so far)

File:

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

palm/trunk/SOURCE/parin.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Nesting for chemical species implemented
+! 
+! 2766 2018-01-22 17:17:47Z kanani
 ! Removed preprocessor directive __chem
 ! 
@@ -617 +620 @@
              ENDIF
              IF ( nest_domain )  THEN
-                bc_uv_t = 'nested'
-                bc_pt_t = 'nested'
-                bc_q_t  = 'nested'
-                bc_s_t  = 'nested'
-                bc_p_t  = 'neumann'
-             ENDIF
-          ELSE
-          
+                bc_uv_t  = 'nested'
+                bc_pt_t  = 'nested'
+                bc_q_t   = 'nested'
+                bc_s_t   = 'nested'
+                bc_cs_t  = 'nested'
+                bc_p_t   = 'neumann'  
+             ENDIF        
 !
 !--       For other nesting modes only set boundary conditions for 
 !--       nested domains.
              IF ( nest_domain )  THEN
-                bc_lr   = 'nested'
-                bc_ns   = 'nested'
-                bc_uv_t = 'nested'
-                bc_pt_t = 'nested'
-                bc_q_t  = 'nested'
-                bc_s_t  = 'nested'
-                bc_p_t  = 'neumann'
+                bc_lr    = 'nested'
+                bc_ns    = 'nested'
+                bc_uv_t  = 'nested'
+                bc_pt_t  = 'nested'
+                bc_q_t   = 'nested'
+                bc_s_t   = 'nested'
+                bc_cs_t = 'nested'
+                bc_p_t   = 'neumann'
              ENDIF
           ENDIF
 
           IF ( forcing )  THEN
-             bc_lr   = 'forcing'
-             bc_ns   = 'forcing'
-             bc_uv_t = 'forcing'
-             bc_pt_t = 'forcing'
-             bc_q_t  = 'forcing'
-             bc_s_t  = 'forcing'  ! scalar boundary condition is not clear
-             bc_p_t  = 'neumann'
-
+             bc_lr    = 'forcing'
+             bc_ns    = 'forcing'
+             bc_uv_t  = 'forcing'
+             bc_pt_t  = 'forcing'
+             bc_q_t   = 'forcing'
+             bc_s_t   = 'forcing'  ! scalar boundary condition is not clear
+             bc_cs_t  = 'forcing'  ! same for chemical species
+             bc_p_t   = 'neumann'
           ENDIF