Changeset 430


Ignore:
Timestamp:
Jan 29, 2010 9:14:09 AM (15 years ago)
Author:
raasch
Message:

bugfix concerning gradient_level indices in header

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r421 r430  
    3333Errors:
    3434------
     35Bugfix: index problem concerning gradient_level indices removed (header)
     36
    3537Dimension of array stat in cascade change to prevent type problems with
    3638mpi2 libraries (poisfft_hybrid)
     
    5456Bugfix: calculation of collision_efficiency (collision_efficiency)
    5557
    56 advec_particles, calc_precipitation, collision_efficiency, disturb_field, modules, poisfft_hybrid, sum_up_3d_data, time_integration
     58advec_particles, calc_precipitation, collision_efficiency, disturb_field, header, modules, poisfft_hybrid, sum_up_3d_data, time_integration
    5759
    5860
  • palm/trunk/SOURCE/header.f90

    r411 r430  
    55! -----------------
    66! +large scale vertical motion (subsidence/ascent)
     7! Bugfix: index problem concerning gradient_level indices removed
    78!
    89!
     
    12281229    DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
    12291230     
     1231       print*, 'i=',i,' level_u=',ug_vertical_gradient_level_ind(i)
    12301232       WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
    12311233       ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
     
    12401242       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
    12411243
    1242        i = i + 1
     1244       IF ( i == 10 )  THEN
     1245          EXIT
     1246       ELSE
     1247          i = i + 1
     1248       ENDIF
     1249
    12431250    ENDDO
    12441251
     
    12551262    DO  WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
    12561263
     1264       print*, 'i=',i,' level_v=',vg_vertical_gradient_level_ind(i)
    12571265       WRITE (coor_chr,'(F6.2,1X)')  vg(vg_vertical_gradient_level_ind(i))
    12581266       vgcomponent = TRIM( vgcomponent ) // '  ' // TRIM( coor_chr )
     
    12671275       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
    12681276
    1269        i = i + 1
     1277       IF ( i == 10 )  THEN
     1278          EXIT
     1279       ELSE
     1280          i = i + 1
     1281       ENDIF
     1282 
    12701283    ENDDO
    12711284
     
    12831296    DO  WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
    12841297
     1298       print*, 'i=',i,' level_pt=',pt_vertical_gradient_level_ind(i)
    12851299       WRITE (coor_chr,'(F7.2)')  pt_init(pt_vertical_gradient_level_ind(i))
    12861300       temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
     
    12951309       coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
    12961310
    1297        i = i + 1
     1311       IF ( i == 10 )  THEN
     1312          EXIT
     1313       ELSE
     1314          i = i + 1
     1315       ENDIF
     1316
    12981317    ENDDO
    12991318
     
    13241343          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
    13251344
    1326           i = i + 1
     1345          IF ( i == 10 )  THEN
     1346             EXIT
     1347          ELSE
     1348             i = i + 1
     1349          ENDIF
     1350
    13271351       ENDDO
    13281352
     
    13591383          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
    13601384
    1361           i = i + 1
     1385          IF ( i == 10 )  THEN
     1386             EXIT
     1387          ELSE
     1388             i = i + 1
     1389          ENDIF
     1390
    13621391       ENDDO
    13631392
     
    13901419          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
    13911420
    1392           i = i + 1
     1421          IF ( i == 10 )  THEN
     1422             EXIT
     1423          ELSE
     1424             i = i + 1
     1425          ENDIF
     1426
    13931427       ENDDO
    13941428
Note: See TracChangeset for help on using the changeset viewer.