Changeset 631

Timestamp:

Dec 13, 2010 3:06:48 PM (13 years ago)

Author:

maronga

Message:

bugfix: type conversion (init_3d_model, error on ibm machines)

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: type conversion of '1' to 64bit for the MAX function (ngp_3d_inner)
 !
 ! Former revisions:
@@ -1594 +1594 @@
 !-- 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_s_inner = MAX( 1, ngp_2dh_s_inner(:,:) ) 
+    ngp_2dh_outer   = MAX( 1, ngp_2dh_outer(:,:)   )
+    ngp_3d_inner    = MAX( INT(1, KIND = SELECTED_INT_KIND( 18 )),            &
+                           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, ngp_3d_inner_tmp )