Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (23 months ago)
Author:
motisi
Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:
1 edited

Legend:

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

    r4331 r4346  
    2525! -----------------
    2626! $Id$
     27! Introduction of wall_flags_total_0, which currently sets bits based on static
     28! topography information used in wall_flags_static_0
     29!
     30! 4331 2019-12-10 18:25:02Z suehring
    2731! - Modularize 2-m potential temperature output
    2832! - New output for 10-m wind speed
     
    114118               nzb,                                                            &
    115119               nzt,                                                            &
    116                wall_flags_static_0
     120               wall_flags_total_0
    117121
    118122    USE kinds
     
    704708             DO  k = nzb_do, nzt_do
    705709                local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i),                &
    706                                          REAL( fill_value, KIND = wp ),        &
    707                                          BTEST( wall_flags_static_0(k,j,i), flag_nr ) )
     710                                     REAL( fill_value, KIND = wp ),            &
     711                                     BTEST( wall_flags_total_0(k,j,i), flag_nr ) )
    708712             ENDDO
    709713          ENDDO
     
    807811             DO  k = nzb_do, nzt_do
    808812                local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i),                &
    809                                          REAL( fill_value, KIND = wp ),        &
    810                                          BTEST( wall_flags_static_0(k,j,i), flag_nr ) )
     813                                     REAL( fill_value, KIND = wp ),            &
     814                                     BTEST( wall_flags_total_0(k,j,i), flag_nr ) )
    811815             ENDDO
    812816          ENDDO
     
    924928                im = mask_i(mid,i)
    925929                jm = mask_j(mid,j)
    926                 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), &
     930                ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), &
    927931                              DIM = 1 ) - 1
    928932                DO  k = 1, mask_size_l(mid,3)
     
    930934!
    931935!--                Set value if not in building
    932                    IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     936                   IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    933937                      local_pf(i,j,k) = fill_value
    934938                   ELSE
     
    10711075                          - (   u(k,j+1,i) + u(k,j+1,i+1)                      &
    10721076                              - u(k,j-1,i) - u(k,j-1,i+1) ) * ddy     )**2  )  &
    1073                        * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0) )
     1077                       * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 0) )
    10741078                   ENDDO
    10751079                ENDDO
     
    10821086                   DO  k = nzb+1, nzt
    10831087                      uu(k,j,i) = u(k,j,i) * u(k,j,i)                          &
    1084                        * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 1) )
     1088                       * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 1) )
    10851089                   ENDDO
    10861090                ENDDO
     
    10931097                   DO  k = nzb+1, nzt
    10941098                      vv(k,j,i) = v(k,j,i) * v(k,j,i)                          &
    1095                        * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 2) )
     1099                       * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 2) )
    10961100                   ENDDO
    10971101                ENDDO
     
    11041108                   DO  k = nzb+1, nzt-1
    11051109                      ww(k,j,i) = w(k,j,i) * w(k,j,i)                          &
    1106                        * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 3) )
     1110                       * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), 3) )
    11071111                   ENDDO
    11081112                ENDDO
Note: See TracChangeset for help on using the changeset viewer.