Changeset 13 for palm/trunk
- Timestamp:
- Feb 14, 2007 12:15:07 PM (18 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
advec_particles.f90 (modified) (1 diff)
-
data_output_dvrp.f90 (modified) (15 diffs)
-
init_dvrp.f90 (modified) (20 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_particles.f90
r4 r13 208 208 ! WRITE (9,*) '--- advec_particles: #1' 209 209 ! WRITE (9,*) ' #of p=',number_of_particles,' #of t=',number_of_tails 210 CALL FLUSH_( 9 )210 ! CALL FLUSH_( 9 ) 211 211 ! ENDIF 212 212 ! -
palm/trunk/SOURCE/data_output_dvrp.f90
r4 r13 79 79 80 80 WRITE ( 9, * ) '*** myid=', myid, ' Anfang data_output_dvrp' 81 CALL FLUSH_( 9 ) 81 #if defined( __ibm ) 82 CALL FLUSH_( 9 ) 83 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 84 CALL FLUSH( 9 ) 85 #endif 82 86 CALL cpu_log( log_point(27), 'data_output_dvrp', 'start' ) 83 87 … … 95 99 lock_steering_update = .TRUE. 96 100 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: vor steering_update' 97 CALL FLUSH_( 9 ) 101 #if defined( __ibm ) 102 CALL FLUSH_( 9 ) 103 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 104 CALL FLUSH( 9 ) 105 #endif 98 106 CALL DVRP_STEERING_UPDATE( m-1, data_output_dvrp ) 99 107 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: nach steering_update' 100 CALL FLUSH_( 9 ) 108 #if defined( __ibm ) 109 CALL FLUSH_( 9 ) 110 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 111 CALL FLUSH( 9 ) 112 #endif 101 113 lock_steering_update = .FALSE. 102 114 ENDIF … … 134 146 #if defined( __particles ) 135 147 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: anfang particles' 136 CALL FLUSH_( 9 ) 148 #if defined( __ibm ) 149 CALL FLUSH_( 9 ) 150 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 151 CALL FLUSH( 9 ) 152 #endif 137 153 ! 138 154 !-- DVRP-Calls for plotting particles: … … 163 179 WRITE (9,*) '--- before ALLOCATE simtime=',simulated_time,' #of_tails=', number_of_tails, & 164 180 ' max#of_tp=', maximum_number_of_tailpoints 165 CALL FLUSH_( 9 ) 181 #if defined( __ibm ) 182 CALL FLUSH_( 9 ) 183 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 184 CALL FLUSH( 9 ) 185 #endif 166 186 ALLOCATE( psize(number_of_tails), p_t(number_of_tails), & 167 187 p_c(number_of_tails*maximum_number_of_tailpoints), & … … 170 190 p_z(number_of_tails*maximum_number_of_tailpoints) ) 171 191 WRITE (9,*) '--- after ALLOCATE' 172 CALL FLUSH_( 9 ) 192 #if defined( __ibm ) 193 CALL FLUSH_( 9 ) 194 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 195 CALL FLUSH( 9 ) 196 #endif 173 197 psize = 0.0; p_t = 0; p_c = 0.0; p_x = 0.0; p_y = 0.0 174 198 p_z = 0.0; … … 181 205 IF ( simulated_time > 1338.0 ) THEN 182 206 WRITE (9,*) '--- particle ',n,' tail_id=',nn,' #of_tp=',particles(n)%tailpoints 183 CALL FLUSH_( 9 ) 207 #if defined( __ibm ) 208 CALL FLUSH_( 9 ) 209 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 210 CALL FLUSH( 9 ) 211 #endif 184 212 ENDIF 185 213 DO j = 1, particles(n)%tailpoints … … 195 223 WRITE (9,*) '--- tp= ',i,' x=',p_x(i),' y=',p_y(i), & 196 224 ' z=',p_z(i),' c=',p_c(i) 197 CALL FLUSH_( 9 ) 225 #if defined( __ibm ) 226 CALL FLUSH_( 9 ) 227 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 228 CALL FLUSH( 9 ) 229 #endif 198 230 ENDIF 199 231 ENDDO … … 202 234 IF ( simulated_time > 1338.0 ) THEN 203 235 WRITE (9,*) '--- t= ',k,' psize=',psize(k),' p_t=',p_t(k) 204 CALL FLUSH_( 9 ) 236 #if defined( __ibm ) 237 CALL FLUSH_( 9 ) 238 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 239 CALL FLUSH( 9 ) 240 #endif 205 241 ENDIF 206 242 ENDIF 207 243 ENDDO 208 244 WRITE (9,*) '--- after locally storing the particle attributes' 209 CALL FLUSH_( 9 ) 245 #if defined( __ibm ) 246 CALL FLUSH_( 9 ) 247 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 248 CALL FLUSH( 9 ) 249 #endif 210 250 ENDIF 211 251 … … 231 271 ELSE 232 272 WRITE (9,*) '--- before DVRP_PARTICLES' 233 CALL FLUSH_( 9 ) 273 #if defined( __ibm ) 274 CALL FLUSH_( 9 ) 275 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 276 CALL FLUSH( 9 ) 277 #endif 234 278 CALL DVRP_PARTICLES( m-1, number_of_tails, p_x, p_y, p_z, 15, & 235 279 psize, p_c, p_t ) 236 280 WRITE (9,*) '--- after DVRP_PARTICLES' 237 CALL FLUSH_( 9 ) 281 #if defined( __ibm ) 282 CALL FLUSH_( 9 ) 283 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 284 CALL FLUSH( 9 ) 285 #endif 238 286 ENDIF 239 287 ENDIF … … 241 289 CALL DVRP_VISUALIZE( m-1, 3, dvrp_filecount ) 242 290 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: ende particles' 243 CALL FLUSH_( 9 ) 291 #if defined( __ibm ) 292 CALL FLUSH_( 9 ) 293 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 294 CALL FLUSH( 9 ) 295 #endif 244 296 245 297 DEALLOCATE( psize, p_c, p_t, p_x, p_y, p_z ) … … 392 444 393 445 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: anfang isosurface' 394 CALL FLUSH_( 9 ) 446 #if defined( __ibm ) 447 CALL FLUSH_( 9 ) 448 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 449 CALL FLUSH( 9 ) 450 #endif 395 451 ! 396 452 !-- DVRP-Calls for plotting isosurfaces: … … 413 469 CALL DVRP_VISUALIZE( m-1, 1, dvrp_filecount ) 414 470 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: ende isosurface' 415 CALL FLUSH_( 9 ) 471 #if defined( __ibm ) 472 CALL FLUSH_( 9 ) 473 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 474 CALL FLUSH( 9 ) 475 #endif 416 476 417 477 CALL cpu_log( log_point_s(26), 'dvrp_isosurface', 'stop' ) … … 420 480 421 481 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: anfang slicer' 422 CALL FLUSH_( 9 ) 482 #if defined( __ibm ) 483 CALL FLUSH_( 9 ) 484 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 485 CALL FLUSH( 9 ) 486 #endif 423 487 ! 424 488 !-- DVRP-Calls for plotting slicers: … … 450 514 451 515 WRITE ( 9, * ) '*** myid=', myid, ' data_output_dvrp: ende slicer' 452 CALL FLUSH_( 9 ) 516 #if defined( __ibm ) 517 CALL FLUSH_( 9 ) 518 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 519 CALL FLUSH( 9 ) 520 #endif 453 521 ENDIF 454 522 … … 465 533 CALL cpu_log( log_point(27), 'data_output_dvrp', 'stop' ) 466 534 WRITE ( 9, * ) '*** myid=', myid, ' Ende data_output_dvrp' 467 CALL FLUSH_( 9 ) 535 #if defined( __ibm ) 536 CALL FLUSH_( 9 ) 537 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 538 CALL FLUSH( 9 ) 539 #endif 468 540 469 541 #endif -
palm/trunk/SOURCE/init_dvrp.f90
r4 r13 80 80 ENDIF 81 81 82 IF ( dvrp_file /= 'default' .AND. dvrp_file /= '/dev/null' ) THEN 83 IF ( myid == 0 ) THEN 84 PRINT*, '+++ init_dvrp: dvrp_file="', dvrp_file, '" not allowed' 85 CALL local_stop 82 IF ( dvrp_output /= 'local' ) THEN 83 IF ( dvrp_file /= 'default' .AND. dvrp_file /= '/dev/null' ) THEN 84 IF ( myid == 0 ) THEN 85 PRINT*, '+++ init_dvrp: dvrp_file="', dvrp_file, '" not allowed' 86 CALL local_stop 87 ENDIF 86 88 ENDIF 87 89 ENDIF … … 138 140 139 141 WRITE ( 9, * ) '*** vor dvrp_output_rtsp' 140 CALL FLUSH_( 9 ) 142 #if defined( __ibm ) 143 CALL FLUSH_( 9 ) 144 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 145 CALL FLUSH( 9 ) 146 #endif 147 141 148 dvrp_file = TRIM( mode_dvrp(m) ) // '/camera.dvr' 142 149 dvrp_file_c = dvrp_file … … 145 152 dvrp_file_c ) 146 153 WRITE ( 9, * ) '*** nach dvrp_output_rtsp' 147 CALL FLUSH_( 9 ) 154 #if defined( __ibm ) 155 CALL FLUSH_( 9 ) 156 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 157 CALL FLUSH( 9 ) 158 #endif 148 159 149 160 ELSEIF ( dvrp_output == 'ftp' ) THEN … … 163 174 dvrp_file_local_c = dvrp_file_c 164 175 ENDIF 165 !CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )176 CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c ) 166 177 167 178 ENDIF … … 169 180 CALL DVRP_CAMERA( m-1, center, distance ) 170 181 WRITE ( 9, * ) '*** #1' 171 CALL FLUSH_( 9 ) 182 #if defined( __ibm ) 183 CALL FLUSH_( 9 ) 184 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 185 CALL FLUSH( 9 ) 186 #endif 172 187 173 188 ! … … 182 197 CALL DVRP_EXIT( m-1 ) 183 198 WRITE ( 9, * ) '*** #2' 184 CALL FLUSH_( 9 ) 199 #if defined( __ibm ) 200 CALL FLUSH_( 9 ) 201 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 202 CALL FLUSH( 9 ) 203 #endif 185 204 186 205 … … 199 218 dvrp_file_c ) 200 219 WRITE ( 9, * ) '*** #3' 201 CALL FLUSH_( 9 ) 220 #if defined( __ibm ) 221 CALL FLUSH_( 9 ) 222 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 223 CALL FLUSH( 9 ) 224 #endif 202 225 203 226 ENDIF … … 228 251 CALL DVRP_MATERIAL_RGB( m-1, 1, 0.8, 0.7, 0.6, 0.0 ) 229 252 WRITE ( 9, * ) '*** #4' 230 CALL FLUSH_( 9 ) 253 #if defined( __ibm ) 254 CALL FLUSH_( 9 ) 255 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 256 CALL FLUSH( 9 ) 257 #endif 231 258 232 259 ! … … 242 269 ENDDO 243 270 WRITE ( 9, * ) '*** #4.1' 244 CALL FLUSH_( 9 ) 271 #if defined( __ibm ) 272 CALL FLUSH_( 9 ) 273 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 274 CALL FLUSH( 9 ) 275 #endif 245 276 CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, & 246 277 cyclic_dvrp, cyclic_dvrp, cyclic_dvrp ) 247 278 WRITE ( 9, * ) '*** #4.2' 248 CALL FLUSH_( 9 ) 279 #if defined( __ibm ) 280 CALL FLUSH_( 9 ) 281 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 282 CALL FLUSH( 9 ) 283 #endif 249 284 CALL DVRP_THRESHOLD( m-1, 1.0 ) 250 285 WRITE ( 9, * ) '*** #4.3' 251 CALL FLUSH_( 9 ) 286 #if defined( __ibm ) 287 CALL FLUSH_( 9 ) 288 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 289 CALL FLUSH( 9 ) 290 #endif 252 291 CALL DVRP_VISUALIZE( m-1, 1, 0 ) 253 292 WRITE ( 9, * ) '*** #4.4' 254 CALL FLUSH_( 9 ) 293 #if defined( __ibm ) 294 CALL FLUSH_( 9 ) 295 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 296 CALL FLUSH( 9 ) 297 #endif 255 298 256 299 DEALLOCATE( local_pf ) … … 258 301 CALL DVRP_EXIT( m-1 ) 259 302 WRITE ( 9, * ) '*** #5' 260 CALL FLUSH_( 9 ) 303 #if defined( __ibm ) 304 CALL FLUSH_( 9 ) 305 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 306 CALL FLUSH( 9 ) 307 #endif 261 308 262 309 ! … … 274 321 dvrp_file_c ) 275 322 WRITE ( 9, * ) '*** #6' 276 CALL FLUSH_( 9 ) 323 #if defined( __ibm ) 324 CALL FLUSH_( 9 ) 325 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 326 CALL FLUSH( 9 ) 327 #endif 277 328 278 329 ENDIF … … 303 354 CALL DVRP_MATERIAL_RGB( m-1, 1, 0.0, 0.6, 0.0, 0.0 ) 304 355 WRITE ( 9, * ) '*** #7' 305 CALL FLUSH_( 9 ) 356 #if defined( __ibm ) 357 CALL FLUSH_( 9 ) 358 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 359 CALL FLUSH( 9 ) 360 #endif 306 361 307 362 ! … … 320 375 CALL DVRP_EXIT( m-1 ) 321 376 WRITE ( 9, * ) '*** #8' 322 CALL FLUSH_( 9 ) 377 #if defined( __ibm ) 378 CALL FLUSH_( 9 ) 379 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 380 CALL FLUSH( 9 ) 381 #endif 323 382 324 383 … … 361 420 name_c = TRIM( do3d(0,vn) ) 362 421 WRITE ( 9, * ) '*** #9' 363 CALL FLUSH_( 9 ) 422 #if defined( __ibm ) 423 CALL FLUSH_( 9 ) 424 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 425 CALL FLUSH( 9 ) 426 #endif 364 427 CALL DVRP_STEERING_INIT( m-1, name_c, steering_dvrp(pn)%min, & 365 428 steering_dvrp(pn)%max, threshold(tv) ) 366 429 WRITE ( 9, * ) '*** #10' 367 CALL FLUSH_( 9 ) 430 #if defined( __ibm ) 431 CALL FLUSH_( 9 ) 432 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 433 CALL FLUSH( 9 ) 434 #endif 368 435 369 436 ELSEIF ( mode_dvrp(pn)(1:6) == 'slicer' ) THEN … … 409 476 410 477 WRITE ( 9, * ) '*** #11' 411 CALL FLUSH_( 9 ) 478 #if defined( __ibm ) 479 CALL FLUSH_( 9 ) 480 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 481 CALL FLUSH( 9 ) 482 #endif 412 483 413 484 dvrp_file = TRIM( mode_dvrp(m) ) // '/*****.dvr' … … 417 488 dvrp_file_c ) 418 489 WRITE ( 9, * ) '*** #12' 419 CALL FLUSH_( 9 ) 490 #if defined( __ibm ) 491 CALL FLUSH_( 9 ) 492 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 493 CALL FLUSH( 9 ) 494 #endif 420 495 421 496 ELSEIF ( dvrp_output == 'ftp' ) THEN … … 429 504 430 505 IF ( dvrp_file(1:9) /= '/dev/null' ) THEN 431 dvrp_file_local = TRIM( mode_dvrp(m) ) // ' .camera.dvr'506 dvrp_file_local = TRIM( mode_dvrp(m) ) // '_%05d.dvr' 432 507 dvrp_file_local_c = dvrp_file_local 433 508 ELSE 434 509 dvrp_file_local_c = dvrp_file_c 435 510 ENDIF 436 !CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )511 CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c ) 437 512 438 513 ENDIF … … 469 544 !-- Define the grid used by dvrp 470 545 WRITE ( 9, * ) '*** #13' 471 CALL FLUSH_( 9 ) 546 #if defined( __ibm ) 547 CALL FLUSH_( 9 ) 548 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 549 CALL FLUSH( 9 ) 550 #endif 472 551 473 552 CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, ycoor_dvrp, & 474 553 zcoor_dvrp ) 475 554 WRITE ( 9, * ) '*** #14' 476 CALL FLUSH_( 9 ) 555 #if defined( __ibm ) 556 CALL FLUSH_( 9 ) 557 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 558 CALL FLUSH( 9 ) 559 #endif 477 560 478 561 … … 507 590 !-- Initialize logging of calls by DVRP graphic software 508 591 WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_INIT' 509 CALL FLUSH_( 9 ) 592 #if defined( __ibm ) 593 CALL FLUSH_( 9 ) 594 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 595 CALL FLUSH( 9 ) 596 #endif 510 597 CALL DVRP_LOG_INIT( 'DVRP_LOG' // CHAR( 0 ), 0 ) 511 598 WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_INIT' 512 CALL FLUSH_( 9 ) 599 #if defined( __ibm ) 600 CALL FLUSH_( 9 ) 601 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 602 CALL FLUSH( 9 ) 603 #endif 513 604 514 605 ! 515 606 !-- User-defined logging events: #1 (total time needed by PALM) 516 607 WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_SYMBOL' 517 CALL FLUSH_( 9 ) 608 #if defined( __ibm ) 609 CALL FLUSH_( 9 ) 610 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 611 CALL FLUSH( 9 ) 612 #endif 518 613 CALL DVRP_LOG_SYMBOL( 1, 'PALM_total' // CHAR( 0 ) ) 519 614 WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_SYMBOL' 520 CALL FLUSH_( 9 ) 615 #if defined( __ibm ) 616 CALL FLUSH_( 9 ) 617 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 618 CALL FLUSH( 9 ) 619 #endif 521 620 CALL DVRP_LOG_SYMBOL( 2, 'PALM_timestep' // CHAR( 0 ) ) 522 621 WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_EVENT' 523 CALL FLUSH_( 9 ) 622 #if defined( __ibm ) 623 CALL FLUSH_( 9 ) 624 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 625 CALL FLUSH( 9 ) 626 #endif 524 627 CALL DVRP_LOG_EVENT( 1, 1 ) 525 628 WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_EVENT' 526 CALL FLUSH_( 9 ) 629 #if defined( __ibm ) 630 CALL FLUSH_( 9 ) 631 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 632 CALL FLUSH( 9 ) 633 #endif 527 634 528 635 #if defined( __parallel ) … … 533 640 use_seperate_pe_for_dvrp_output = .TRUE. 534 641 WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_SPLIT' 535 CALL FLUSH_( 9 ) 642 #if defined( __ibm ) 643 CALL FLUSH_( 9 ) 644 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 645 CALL FLUSH( 9 ) 646 #endif 536 647 CALL DVRP_SPLIT( MPI_COMM_WORLD, comm_palm ) 537 648 WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_SPLIT' 538 CALL FLUSH_( 9 ) 649 #if defined( __ibm ) 650 CALL FLUSH_( 9 ) 651 #elif defined( __lcmuk ) || defined( __lctit ) || defined( __nec ) 652 CALL FLUSH( 9 ) 653 #endif 539 654 CALL MPI_COMM_SIZE( comm_palm, numprocs, ierr ) 540 655 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
