Changeset 1353 for palm/trunk/SOURCE/data_output_dvrp.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_dvrp.f90
r1347 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 75 75 slicer_range_limits_dvrp(1,islice_dvrp) ) 76 76 77 scale = MODULO( 180.0 + 180.0* scale, 360.0_wp )77 scale = MODULO( 180.0_wp + 180.0_wp * scale, 360.0_wp ) 78 78 79 79 color = (/ scale, 0.5_wp, 1.0_wp, 0.0_wp /) … … 231 231 !-- Definition of characteristics of particle material 232 232 ! tmp_r = 0.1; tmp_g = 0.7; tmp_b = 0.1; tmp_t = 0.0 233 tmp_r = 0.0 ; tmp_g = 0.0; tmp_b = 0.0; tmp_t = 0.0233 tmp_r = 0.0_wp; tmp_g = 0.0_wp; tmp_b = 0.0_wp; tmp_t = 0.0_wp 234 234 CALL DVRP_MATERIAL_RGB( m-1, 1, tmp_r, tmp_g, tmp_b, tmp_t ) 235 235 … … 284 284 ALLOCATE( psize(dvrp_nop), p_t(dvrp_nop), p_c(dvrp_nop), & 285 285 p_x(dvrp_nop), p_y(dvrp_nop), p_z(dvrp_nop) ) 286 psize = 0.0; p_t = 0; p_c = 0.0; p_x = 0.0; p_y = 0.0 287 p_z = 0.0 286 psize = 0.0_wp; p_t = 0_wp; p_c = 0.0_wp 287 p_x = 0.0_wp; p_y = 0.0_wp 288 p_z = 0.0_wp 288 289 k = 0 289 290 DO n = 1, number_of_particles … … 305 306 p_y(dvrp_not*maximum_number_of_tailpoints), & 306 307 p_z(dvrp_not*maximum_number_of_tailpoints) ) 307 psize = 0.0; p_t = 0; p_c = 0.0; p_x = 0.0; p_y = 0.0 308 p_z = 0.0; 308 psize = 0.0_wp; p_t = 0_wp; p_c = 0.0_wp 309 p_x = 0.0_wp; p_y = 0.0_wp 310 p_z = 0.0_wp 309 311 i = 0 310 312 k = 0 … … 364 366 365 367 366 ELSEIF ( ( mode_dvrp(m)(1:10) == 'isosurface' .OR. &367 mode_dvrp(m)(1:6) == 'slicer' ) &368 ELSEIF ( ( mode_dvrp(m)(1:10) == 'isosurface' .OR. & 369 mode_dvrp(m)(1:6) == 'slicer' ) & 368 370 .AND. output_variable /= ' ' ) THEN 369 371 370 372 ! 371 373 !-- Create an intermediate array, properly dimensioned for plot-output 372 ALLOCATE( local_pf(nxl_dvrp:nxr_dvrp+1,nys_dvrp:nyn_dvrp+1, &374 ALLOCATE( local_pf(nxl_dvrp:nxr_dvrp+1,nys_dvrp:nyn_dvrp+1, & 373 375 nzb:nz_do3d) ) 374 376 … … 390 392 !-- Replace mirrored values at lower surface by real surface 391 393 !-- values 392 IF ( output_variable == 'u_xz' .OR. &394 IF ( output_variable == 'u_xz' .OR. & 393 395 output_variable == 'u_yz' ) THEN 394 IF ( ibc_uv_b == 0 ) local_pf(:,:,nzb) = 0.0 396 IF ( ibc_uv_b == 0 ) local_pf(:,:,nzb) = 0.0_wp 395 397 ENDIF 396 398 … … 407 409 !-- Replace mirrored values at lower surface by real surface 408 410 !-- values 409 IF ( output_variable == 'v_xz' .OR. &411 IF ( output_variable == 'v_xz' .OR. & 410 412 output_variable == 'v_yz' ) THEN 411 IF ( ibc_uv_b == 0 ) local_pf(:,:,nzb) = 0.0 413 IF ( ibc_uv_b == 0 ) local_pf(:,:,nzb) = 0.0_wp 412 414 ENDIF 413 415 … … 453 455 DO j = nys_dvrp, nyn_dvrp+1 454 456 DO k = nzb, nz_do3d 455 local_pf(i,j,k) = pt(k,j,i) + l_d_cp * &457 local_pf(i,j,k) = pt(k,j,i) + l_d_cp * & 456 458 pt_d_t(k) * ql(k,j,i) 457 459 ENDDO … … 470 472 ENDDO 471 473 ELSE 472 message_string = 'if humidity/passive_scalar = ' // &473 'FALSE output of ' // TRIM( output_variable ) // &474 message_string = 'if humidity/passive_scalar = ' // & 475 'FALSE output of ' // TRIM( output_variable ) // & 474 476 'is not provided' 475 477 CALL message( 'data_output_dvrp', 'PA0183',& … … 487 489 ENDDO 488 490 ELSE 489 message_string = 'if cloud_physics = FALSE ' // &490 'output of ' // TRIM( output_variable) // &491 message_string = 'if cloud_physics = FALSE ' // & 492 'output of ' // TRIM( output_variable) // & 491 493 'is not provided' 492 494 CALL message( 'data_output_dvrp', 'PA0184',& … … 529 531 DO j = nys_dvrp, nyn_dvrp+1 530 532 DO k = nzb, nz_do3d 531 local_pf(i,j,k) = 0.0 533 local_pf(i,j,k) = 0.0_wp 532 534 ENDDO 533 535 ENDDO … … 547 549 tmp_g = isosurface_color(2,n_isosurface) 548 550 tmp_b = isosurface_color(3,n_isosurface) 549 tmp_t = 0.0 551 tmp_t = 0.0_wp 550 552 CALL DVRP_MATERIAL_RGB( m-1, 1, tmp_r, tmp_g, tmp_b, tmp_t ) 551 553 552 554 ! 553 555 !-- Compute and plot isosurface in dvr-format 554 CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &556 CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, & 555 557 cyclic_dvrp, cyclic_dvrp, cyclic_dvrp ) 556 558 … … 573 575 tmp_th = threshold(n_isosurface) 574 576 ELSE 575 tmp_th = 1.0 577 tmp_th = 1.0_wp ! nothing is plotted because array values are 0 576 578 ENDIF 577 579 … … 590 592 ! 591 593 !-- Material and color definitions 592 tmp_r = 0.0 ; tmp_g = 0.0; tmp_b = 0.0; tmp_t = 0.0594 tmp_r = 0.0_wp; tmp_g = 0.0_wp; tmp_b = 0.0_wp; tmp_t = 0.0_wp 593 595 CALL DVRP_MATERIAL_RGB( m-1, 1, tmp_r, tmp_g, tmp_b, tmp_t ) 594 596 … … 610 612 CALL user_dvrp_coltab( 'slicer', output_variable ) 611 613 612 CALL DVRP_COLORTABLE_HLS( m-1, 1, interval_values_dvrp, &613 interval_h_dvrp, interval_l_dvrp, &614 CALL DVRP_COLORTABLE_HLS( m-1, 1, interval_values_dvrp, & 615 interval_h_dvrp, interval_l_dvrp, & 614 616 interval_s_dvrp, interval_a_dvrp ) 615 617 616 618 ! 617 619 !-- Compute and plot slicer in dvr-format 618 CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &620 CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, & 619 621 cyclic_dvrp, cyclic_dvrp, cyclic_dvrp ) 620 622 tmp_pos = slicer_position … … 638 640 local_pfi(2,i,j,k) = v(k,j,i) 639 641 local_pfi(3,i,j,k) = w(k,j,i) 640 tmp_norm = SQRT( u(k,j,i) * u(k,j,i) + &641 v(k,j,i) * v(k,j,i) + &642 tmp_norm = SQRT( u(k,j,i) * u(k,j,i) + & 643 v(k,j,i) * v(k,j,i) + & 642 644 w(k,j,i) * w(k,j,i) ) 643 tmp_alpha = ACOS( 0.0 * u(k,j,i) / tmp_norm +&644 0.0 * v(k,j,i) / tmp_norm -&645 1.0 * w(k,j,i) / tmp_norm )646 tmp_alpha_w = tmp_alpha / pi * 180.0 645 tmp_alpha = ACOS( 0.0_wp * u(k,j,i) / tmp_norm + & 646 0.0_wp * v(k,j,i) / tmp_norm - & 647 1.0_wp * w(k,j,i) / tmp_norm ) 648 tmp_alpha_w = tmp_alpha / pi * 180.0_wp 647 649 local_pfi(4,i,j,k) = tmp_alpha_w 648 650 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.