Changeset 430 for palm/trunk

Timestamp:

Jan 29, 2010 9:14:09 AM (14 years ago)

Author:

raasch

Message:

bugfix concerning gradient_level indices in header

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (2 diffs)
• header.f90 (modified) (10 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -33 +33 @@
 Errors:
 ------
+Bugfix: index problem concerning gradient_level indices removed (header)
+
 Dimension of array stat in cascade change to prevent type problems with
 mpi2 libraries (poisfft_hybrid)
@@ -54 +56 @@
 Bugfix: calculation of collision_efficiency (collision_efficiency)
 
-advec_particles, calc_precipitation, collision_efficiency, disturb_field, modules, poisfft_hybrid, sum_up_3d_data, time_integration
+advec_particles, calc_precipitation, collision_efficiency, disturb_field, header, modules, poisfft_hybrid, sum_up_3d_data, time_integration
 
 

palm/trunk/SOURCE/header.f90

@@ -5 +5 @@
 ! -----------------
 ! +large scale vertical motion (subsidence/ascent)
+! Bugfix: index problem concerning gradient_level indices removed
 ! 
 !
@@ -1228 +1229 @@
     DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
       
+       print*, 'i=',i,' level_u=',ug_vertical_gradient_level_ind(i)
        WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
        ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
@@ -1240 +1242 @@
        coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
 
-       i = i + 1
+       IF ( i == 10 )  THEN
+          EXIT
+       ELSE
+          i = i + 1
+       ENDIF
+
     ENDDO
 
@@ -1255 +1262 @@
     DO  WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
 
+       print*, 'i=',i,' level_v=',vg_vertical_gradient_level_ind(i)
        WRITE (coor_chr,'(F6.2,1X)')  vg(vg_vertical_gradient_level_ind(i))
        vgcomponent = TRIM( vgcomponent ) // '  ' // TRIM( coor_chr )
@@ -1267 +1275 @@
        coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
 
-       i = i + 1 
+       IF ( i == 10 )  THEN
+          EXIT
+       ELSE
+          i = i + 1
+       ENDIF
+ 
     ENDDO
 
@@ -1283 +1296 @@
     DO  WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
 
+       print*, 'i=',i,' level_pt=',pt_vertical_gradient_level_ind(i)
        WRITE (coor_chr,'(F7.2)')  pt_init(pt_vertical_gradient_level_ind(i))
        temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
@@ -1295 +1309 @@
        coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
 
-       i = i + 1
+       IF ( i == 10 )  THEN
+          EXIT
+       ELSE
+          i = i + 1
+       ENDIF
+
     ENDDO
 
@@ -1324 +1343 @@
           coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
 
-          i = i + 1
+          IF ( i == 10 )  THEN
+             EXIT
+          ELSE
+             i = i + 1
+          ENDIF
+
        ENDDO
 
@@ -1359 +1383 @@
           coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
 
-          i = i + 1
+          IF ( i == 10 )  THEN
+             EXIT
+          ELSE
+             i = i + 1
+          ENDIF
+
        ENDDO
 
@@ -1390 +1419 @@
           coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
 
-          i = i + 1
+          IF ( i == 10 )  THEN
+             EXIT
+          ELSE
+             i = i + 1
+          ENDIF
+
        ENDDO