Ignore:
Timestamp:
Mar 30, 2011 9:31:40 AM (10 years ago)
Author:
raasch
Message:

formatting adjustments

File:
1 edited

Legend:

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

    r668 r709  
    44! Current revisions:
    55! -----------------
     6! formatting adjustments
    67!
    78! Former revisions:
     
    5455!
    5556!-- First, local averaging within the recycling domain
    56 
    5757    i = recycling_plane
    5858
     
    7979!-- Now, averaging over all PEs
    8080    IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    81     CALL MPI_ALLREDUCE( avpr_l(nzb,1,1), avpr(nzb,1,1), ngp_pr,  &
    82                     MPI_REAL, MPI_SUM, comm2d, ierr )
     81    CALL MPI_ALLREDUCE( avpr_l(nzb,1,1), avpr(nzb,1,1), ngp_pr, MPI_REAL, &
     82                        MPI_SUM, comm2d, ierr )
    8383
    8484#else
     
    165165          DO  k = nzb, nzt + 1
    166166
    167               u(k,j,-nbgp+1:0)   = mean_inflow_profiles(k,1) + &
     167              u(k,j,-nbgp+1:0) = mean_inflow_profiles(k,1) + &
    168168                           inflow_dist(k,j,1,1:nbgp) * inflow_damping_factor(k)
    169169              v(k,j,-nbgp:-1)  = mean_inflow_profiles(k,2) + &
    170170                           inflow_dist(k,j,2,1:nbgp) * inflow_damping_factor(k)
    171               w(k,j,-nbgp:-1)  = inflow_dist(k,j,3,1:nbgp) * inflow_damping_factor(k)
     171              w(k,j,-nbgp:-1)  =                             &
     172                           inflow_dist(k,j,3,1:nbgp) * inflow_damping_factor(k)
    172173              pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) + &
    173174                           inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
     
    181182    ENDIF
    182183
    183 !
    184 !-- Conserve the volume flow at the inflow in order to avoid generation of
    185 !-- waves in the stable layer
    186 !    IF ( conserve_volume_flow  .AND.  inflow_l )  THEN
    187 
    188 !       volume_flow(1)   = 0.0
    189 !       volume_flow_l(1) = 0.0
    190 
    191 !       i = 0
    192 
    193 !       DO  j = nys, nyn
    194 !
    195 !--       Sum up the volume flow through the south/north boundary
    196 !          DO  k = nzb_2d(j,i) + 1, nzt
    197 !             volume_flow_l(1) = volume_flow_l(1) + u(k,j,i) * dzu(k)
    198 !          ENDDO
    199 !       ENDDO
    200 
    201 #if defined( __parallel )   
    202 !       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    203 !       CALL MPI_ALLREDUCE( volume_flow_l(1), volume_flow(1), 1, MPI_REAL, &
    204 !                           MPI_SUM, comm1dy, ierr )   
    205 #else
    206 !       volume_flow = volume_flow_l 
    207 #endif
    208 !       volume_flow_offset(1) = ( volume_flow_initial(1) - volume_flow(1) )    &
    209 !                               / volume_flow_area(1)
    210 
    211 !       DO  j = nys-1, nyn+1
    212 !          DO  k = nzb_v_inner(j,i) + 1, nzt
    213 !             u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
    214 !          ENDDO
    215 !       ENDDO
    216 
    217 !    ENDIF
    218 
    219184    CALL cpu_log( log_point(40), 'inflow_turbulence', 'stop' )
    220185
Note: See TracChangeset for help on using the changeset viewer.