Ignore:
Timestamp:
May 30, 2017 5:47:52 PM (4 years ago)
Author:
suehring
Message:

Adjustments according new topography and surface-modelling concept implemented

File:
1 edited

Legend:

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

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Adjustments to new topography concept
    2323!
    2424! Former revisions:
     
    105105
    106106       USE indices,                                                            &
    107            ONLY:  nxl, nxr, nyn, nys, nzb, nzb_s_inner, nzt
     107           ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
    108108
    109109       USE kinds
     
    153153             blackbody_emission(nzb) = sigma * temperature**4
    154154
    155              DO  k = nzb_s_inner(j,i)+1, nzt
     155             DO  k = nzb+1, nzt
    156156
    157157                k_help = ( nzt+nzb+1 ) - k
     
    163163
    164164                temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
    165                 blackbody_emission(k) = sigma * temperature**4
     165                blackbody_emission(k) = sigma * temperature**4                 &
     166                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     167                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    166168
    167169             ENDDO
     
    179181             impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
    180182
    181              DO  k = nzb_s_inner(j,i)+1, nzt
     183             DO  k = nzb+1, nzt
    182184!
    183185!--             Save some computational time, but this may cause load
     
    222224                   tend(k,j,i) = tend(k,j,i) -                                 &
    223225                                ( pt_d_t(k) / ( rho_surface * cp ) *           &
    224                                   ( df_p - df_m ) / dzw(k) )
     226                                  ( df_p - df_m ) / dzw(k) )                   &
     227                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     228                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    225229
    226230                ENDIF
     
    250254
    251255       USE indices,                                                            &
    252            ONLY:  nzb, nzb_s_inner, nzt
     256           ONLY:  nzb, nzt, wall_flags_0
    253257
    254258       USE kinds
     
    295299       blackbody_emission(nzb) = sigma * temperature**4
    296300
    297        DO  k = nzb_s_inner(j,i)+1, nzt
     301       DO  k = nzb+1, nzt
    298302          k_help = ( nzt+nzb+1 ) - k
    299303          lwp_ground(k)   = lwp_ground(k-1) + rho_surface * ql(k,j,i) * dzw(k)
     
    303307
    304308          temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
    305           blackbody_emission(k) = sigma * temperature**4
     309          blackbody_emission(k) = sigma * temperature**4                       &
     310                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     311                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    306312
    307313       ENDDO
     
    318324       impinging_flux_at_top = blackbody_emission(nzb) - 100.0_wp
    319325
    320        DO  k = nzb_s_inner(j,i)+1, nzt
     326       DO  k = nzb+1, nzt
    321327!
    322328!--       Store some computational time,
     
    356362!--          Compute tendency term         
    357363             tend(k,j,i) = tend(k,j,i) - ( pt_d_t(k) / ( rho_surface * cp ) *  &
    358                                          ( df_p - df_m ) / dzw(k) )
     364                                         ( df_p - df_m ) / dzw(k) )            &
     365                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     366                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    359367
    360368          ENDIF
Note: See TracChangeset for help on using the changeset viewer.