Changeset 333 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Jun 9, 2009 8:22:04 AM (15 years ago)

Author:

heinze

Message:

bugfix in init_3_model concerning ngp_2dh_s_inner

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Current revisions:
 ! -----------------
+! bugfix: avoid that ngp_2dh_s_inner becomes zero
 ! initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill', now
 ! independent of turbulent_inflow
@@ -1473 +1474 @@
 !-- buoyancy, etc. A zero value will occur for cases where all grid points of
 !-- the respective subdomain lie below the surface topography
-    ngp_2dh_outer = MAX( 1, ngp_2dh_outer(:,:) ) 
-    ngp_3d_inner  = MAX( 1, ngp_3d_inner(:)    )
+    ngp_2dh_outer   = MAX( 1, ngp_2dh_outer(:,:)   ) 
+    ngp_3d_inner    = MAX( 1, ngp_3d_inner(:)      )
+    ngp_2dh_s_inner = MAX( 1, ngp_2dh_s_inner(:,:) ) 
 
     DEALLOCATE( ngp_2dh_l, ngp_2dh_outer_l, ngp_3d_inner_l )