Ignore:
Timestamp:
Feb 27, 2009 11:42:39 AM (13 years ago)
Author:
raasch
Message:

further changes for dvrp clipping

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r242 r246  
    7474    USE cloud_parameters
    7575    USE constants
     76    USE control_parameters
    7677    USE cpulog
    7778    USE DVRP
     
    8384    USE particle_attributes
    8485    USE pegrid
    85     USE control_parameters
    8686
    8787    IMPLICIT NONE
     
    165165!--       If clipping is active and if this subdomain is clipped, find out the
    166166!--       number of particles and tails to be plotted; otherwise, all
    167 !--       particles/tails are plotted
     167!--       particles/tails are plotted. dvrp_mask is used to mark the partikles.
    168168          IF ( .NOT. use_particle_tails )  THEN
    169169             ALLOCATE( dvrp_mask(number_of_particles) )
     
    299299!
    300300!--       Create an intermediate array, properly dimensioned for plot-output
    301           ALLOCATE( local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) )
     301          ALLOCATE( local_pf(nxl_dvrp:nxr_dvrp+1,nys_dvrp:nyn_dvrp+1, &
     302                             nzb:nz_do3d) )
    302303
    303304!
    304305!--       Move original array to intermediate array
    305           SELECT CASE ( output_variable )
    306 
    307              CASE ( 'u', 'u_xy', 'u_xz', 'u_yz' )
    308                 DO  i = nxl, nxr+1
    309                    DO  j = nys, nyn+1
    310                       DO  k = nzb, nz_do3d
    311                          local_pf(i,j,k) = u(k,j,i)
    312                       ENDDO
    313                    ENDDO
    314                 ENDDO
    315 !
    316 !--             Replace mirrored values at lower surface by real surface values
    317                 IF ( output_variable == 'u_xz'  .OR. &
    318                      output_variable == 'u_yz' )  THEN
    319                    IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
    320                 ENDIF
    321 
    322 
    323              CASE ( 'v', 'v_xy', 'v_xz', 'v_yz' )
    324                 DO  i = nxl, nxr+1
    325                    DO  j = nys, nyn+1
    326                       DO  k = nzb, nz_do3d
    327                          local_pf(i,j,k) = v(k,j,i)
    328                       ENDDO
    329                    ENDDO
    330                 ENDDO
    331 !
    332 !--             Replace mirrored values at lower surface by real surface values
    333                 IF ( output_variable == 'v_xz'  .OR. &
    334                      output_variable == 'v_yz' )  THEN
    335                    IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
    336                 ENDIF
    337 
    338              CASE ( 'w', 'w_xy', 'w_xz', 'w_yz' )
    339                 DO  i = nxl, nxr+1
    340                    DO  j = nys, nyn+1
    341                       DO  k = nzb, nz_do3d
    342                          local_pf(i,j,k) = w(k,j,i)
    343                       ENDDO
    344                    ENDDO
    345                 ENDDO
     306          IF ( dvrp_overlap )  THEN
     307
     308             SELECT CASE ( output_variable )
     309
     310                CASE ( 'u', 'u_xy', 'u_xz', 'u_yz' )
     311                   DO  i = nxl_dvrp, nxr_dvrp+1
     312                      DO  j = nys_dvrp, nyn_dvrp+1
     313                         DO  k = nzb, nz_do3d
     314                            local_pf(i,j,k) = u(k,j,i)
     315                         ENDDO
     316                      ENDDO
     317                   ENDDO
     318!
     319!--                Replace mirrored values at lower surface by real surface
     320!--                values
     321                   IF ( output_variable == 'u_xz'  .OR. &
     322                        output_variable == 'u_yz' )  THEN
     323                      IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
     324                   ENDIF
     325
     326
     327                CASE ( 'v', 'v_xy', 'v_xz', 'v_yz' )
     328                   DO  i = nxl_dvrp, nxr_dvrp+1
     329                      DO  j = nys_dvrp, nyn_dvrp+1
     330                         DO  k = nzb, nz_do3d
     331                            local_pf(i,j,k) = v(k,j,i)
     332                         ENDDO
     333                      ENDDO
     334                   ENDDO
     335!
     336!--                Replace mirrored values at lower surface by real surface
     337!--                values
     338                   IF ( output_variable == 'v_xz'  .OR. &
     339                        output_variable == 'v_yz' )  THEN
     340                      IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
     341                   ENDIF
     342
     343                CASE ( 'w', 'w_xy', 'w_xz', 'w_yz' )
     344                   DO  i = nxl_dvrp, nxr_dvrp+1
     345                      DO  j = nys_dvrp, nyn_dvrp+1
     346                         DO  k = nzb, nz_do3d
     347                            local_pf(i,j,k) = w(k,j,i)
     348                         ENDDO
     349                      ENDDO
     350                   ENDDO
    346351! Averaging for Langmuir circulation
    347 !                DO  k = nzb, nz_do3d
    348 !                   DO  j = nys+1, nyn
    349 !                      DO  i = nxl, nxr+1
    350 !                         local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
    351 !                                           0.50 * local_pf(i,j,k)   + &
    352 !                                           0.25 * local_pf(i,j+1,k)
     352!                   DO  k = nzb, nz_do3d
     353!                      DO  j = nys_dvrp+1, nyn_dvrp
     354!                         DO  i = nxl_dvrp, nxr_dvrp+1
     355!                            local_pf(i,j,k) = 0.25 * local_pf(i,j-1,k) + &
     356!                                              0.50 * local_pf(i,j,k)   + &
     357!                                              0.25 * local_pf(i,j+1,k)
     358!                         ENDDO
    353359!                      ENDDO
    354360!                   ENDDO
    355 !                ENDDO
    356 
    357              CASE ( 'p', 'p_xy', 'p_xz', 'p_yz' )
    358                 DO  i = nxl, nxr+1
    359                    DO  j = nys, nyn+1
    360                       DO  k = nzb, nz_do3d
    361                          local_pf(i,j,k) = p(k,j,i)
    362                       ENDDO
    363                    ENDDO
    364                 ENDDO
    365 
    366              CASE ( 'pt', 'pt_xy', 'pt_xz', 'pt_yz' )
    367                 IF ( .NOT. cloud_physics ) THEN
    368                    DO  i = nxl, nxr+1
    369                       DO  j = nys, nyn+1
     361
     362                CASE ( 'p', 'p_xy', 'p_xz', 'p_yz' )
     363                   DO  i = nxl_dvrp, nxr_dvrp+1
     364                      DO  j = nys_dvrp, nyn_dvrp+1
    370365                         DO  k = nzb, nz_do3d
    371                             local_pf(i,j,k) = pt(k,j,i)
    372                          ENDDO
    373                       ENDDO
    374                    ENDDO
    375                 ELSE
    376                    DO  i = nxl, nxr+1
    377                       DO  j = nys, nyn+1
    378                          DO  k = nzb, nz_do3d
    379                             local_pf(i,j,k) = pt(k,j,i) + l_d_cp * pt_d_t(k) * &
    380                                                           ql(k,j,i)
    381                          ENDDO
    382                       ENDDO
    383                    ENDDO
    384                 ENDIF
    385 
    386              CASE ( 'q', 'q_xy', 'q_xz', 'q_yz' )
    387                 IF ( humidity  .OR.  passive_scalar )  THEN
    388                    DO  i = nxl, nxr+1
    389                       DO  j = nys, nyn+1
    390                          DO  k = nzb, nz_do3d
    391                             local_pf(i,j,k) = q(k,j,i)
    392                          ENDDO
    393                       ENDDO
    394                    ENDDO           
    395                 ELSE
    396                    IF ( myid == 0 )  THEN
    397                       PRINT*, '+++ data_output_dvrp: if humidity/passive_scalar = ', &
    398                               'FALSE output of ', output_variable,            &
    399                               'is not provided'
     366                            local_pf(i,j,k) = p(k,j,i)
     367                         ENDDO
     368                      ENDDO
     369                   ENDDO
     370
     371                CASE ( 'pt', 'pt_xy', 'pt_xz', 'pt_yz' )
     372                   IF ( .NOT. cloud_physics ) THEN
     373                      DO  i = nxl_dvrp, nxr_dvrp+1
     374                         DO  j = nys_dvrp, nyn_dvrp+1
     375                            DO  k = nzb, nz_do3d
     376                               local_pf(i,j,k) = pt(k,j,i)
     377                            ENDDO
     378                         ENDDO
     379                      ENDDO
     380                   ELSE
     381                      DO  i = nxl_dvrp, nxr_dvrp+1
     382                         DO  j = nys_dvrp, nyn_dvrp+1
     383                            DO  k = nzb, nz_do3d
     384                               local_pf(i,j,k) = pt(k,j,i) + l_d_cp * &
     385                                                 pt_d_t(k) * ql(k,j,i)
     386                            ENDDO
     387                         ENDDO
     388                      ENDDO
    400389                   ENDIF
    401                 ENDIF
     390
     391                CASE ( 'q', 'q_xy', 'q_xz', 'q_yz' )
     392                   IF ( humidity  .OR.  passive_scalar )  THEN
     393                      DO  i = nxl_dvrp, nxr_dvrp+1
     394                         DO  j = nys_dvrp, nyn_dvrp+1
     395                            DO  k = nzb, nz_do3d
     396                               local_pf(i,j,k) = q(k,j,i)
     397                            ENDDO
     398                         ENDDO
     399                      ENDDO           
     400                   ELSE
     401                      IF ( myid == 0 )  THEN
     402                         PRINT*, '+++ data_output_dvrp: if humidity/passive_scalar = ', &
     403                                 'FALSE output of ', output_variable,            &
     404                                 'is not provided'
     405                      ENDIF
     406                   ENDIF
    402407             
    403              CASE ( 'ql', 'ql_xy', 'ql_xz', 'ql_yz' )
    404                 IF ( cloud_physics  .OR.  cloud_droplets )  THEN
    405                    DO  i = nxl, nxr+1
    406                       DO  j = nys, nyn+1
    407                          DO  k = nzb, nz_do3d
    408                             local_pf(i,j,k) = ql(k,j,i)
    409                          ENDDO
    410                       ENDDO
    411                    ENDDO
    412                 ELSE
    413                    IF ( myid == 0 ) THEN
    414                       PRINT*, '+++ data_output_dvrp: if cloud_physics = FALSE ', &
    415                               'output of ', output_variable, 'is not provided'
     408                CASE ( 'ql', 'ql_xy', 'ql_xz', 'ql_yz' )
     409                   IF ( cloud_physics  .OR.  cloud_droplets )  THEN
     410                      DO  i = nxl_dvrp, nxr_dvrp+1
     411                         DO  j = nys_dvrp, nyn_dvrp+1
     412                            DO  k = nzb, nz_do3d
     413                               local_pf(i,j,k) = ql(k,j,i)
     414                            ENDDO
     415                         ENDDO
     416                      ENDDO
     417                   ELSE
     418                      IF ( myid == 0 ) THEN
     419                         PRINT*, '+++ data_output_dvrp: if cloud_physics = FALSE ', &
     420                                 'output of ', output_variable, 'is not provided'
     421                      ENDIF
    416422                   ENDIF
    417                 ENDIF
    418 
    419              CASE ( 'u*_xy' )
    420                 DO  i = nxl, nxr+1
    421                    DO  j = nys, nyn+1
    422                       local_pf(i,j,nzb+1) = us(j,i)
    423                    ENDDO
    424                 ENDDO
    425                 slicer_position = zu(nzb+1)
    426 
    427              CASE ( 't*_xy' )
    428                 DO  i = nxl, nxr+1
    429                    DO  j = nys, nyn+1
    430                       local_pf(i,j,nzb+1) = ts(j,i)
    431                    ENDDO
    432                 ENDDO
    433                 slicer_position = zu(nzb+1)
    434 
    435 
    436              CASE DEFAULT
    437 !
    438 !--             The DEFAULT case is reached either if output_variable contains
    439 !--             unsupported variable or if the user has coded a special case in
    440 !--             the user interface. There, the subroutine user_data_output_dvrp
    441 !--             checks which of these two conditions applies.
    442                 CALL user_data_output_dvrp( output_variable, local_pf )
    443 
    444 
    445           END SELECT
    446 
     423
     424                CASE ( 'u*_xy' )
     425                   DO  i = nxl_dvrp, nxr_dvrp+1
     426                      DO  j = nys_dvrp, nyn_dvrp+1
     427                         local_pf(i,j,nzb+1) = us(j,i)
     428                      ENDDO
     429                   ENDDO
     430                   slicer_position = zu(nzb+1)
     431
     432                CASE ( 't*_xy' )
     433                   DO  i = nxl_dvrp, nxr_dvrp+1
     434                      DO  j = nys_dvrp, nyn_dvrp+1
     435                         local_pf(i,j,nzb+1) = ts(j,i)
     436                      ENDDO
     437                   ENDDO
     438                   slicer_position = zu(nzb+1)
     439
     440
     441                CASE DEFAULT
     442!
     443!--                The DEFAULT case is reached either if output_variable
     444!--                contains unsupported variable or if the user has coded a
     445!--                special case in the user interface. There, the subroutine
     446!--                user_data_output_dvrp checks which of these two conditions
     447!--                applies.
     448                   CALL user_data_output_dvrp( output_variable, local_pf )
     449
     450
     451             END SELECT
     452
     453          ENDIF
    447454
    448455          IF ( mode_dvrp(m)(1:10) == 'isosurface' )  THEN
     
    489496             CALL DVRP_THRESHOLD( -(m-1)-1, tmp_c_alpha )
    490497
    491              tmp_th = threshold(tv)
     498             IF ( dvrp_overlap )  THEN
     499                tmp_th = threshold(tv)
     500             ELSE
     501                tmp_th = 1.0   ! nothing is plotted because array values are 0
     502             ENDIF
     503
    492504             CALL DVRP_THRESHOLD( m-1, tmp_th )
    493505             CALL DVRP_VISUALIZE( m-1, 21, dvrp_filecount )
Note: See TracChangeset for help on using the changeset viewer.