Changeset 1560 for palm


Ignore:
Timestamp:
Mar 6, 2015 10:48:54 AM (9 years ago)
Author:
keck
Message:

implemented possibility of adding a y shift to the recycled inflow turbulence

Location:
palm/trunk/SOURCE
Files:
5 edited

Legend:

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

    r1558 r1560  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! output for recycling y shift
    2323!
    2424! Former revisions:
     
    11011101       WRITE ( io, 318 )  use_cmax, pt_damping_width, pt_damping_factor       
    11021102       IF ( turbulent_inflow )  THEN
    1103           WRITE ( io, 319 )  recycling_width, recycling_plane, &
    1104                              inflow_damping_height, inflow_damping_width
     1103          IF ( .NOT. recycling_yshift ) THEN
     1104             WRITE ( io, 319 )  recycling_width, recycling_plane, &
     1105                                inflow_damping_height, inflow_damping_width
     1106          ELSE
     1107             WRITE ( io, 322 )  recycling_width, recycling_plane, &
     1108                                inflow_damping_height, inflow_damping_width
     1109          END IF
    11051110       ENDIF
    11061111    ENDIF
     
    20472052321 FORMAT (//' Initial profiles:'/ &
    20482053              ' ----------------')
     2054322 FORMAT ('       turbulence recycling at inflow switched on'/ &
     2055            '       y shift of the recycled inflow turbulence switched on'/ &
     2056            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
     2057            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m'/q)
    20492058325 FORMAT (//' List output:'/ &
    20502059             ' -----------'//  &
  • palm/trunk/SOURCE/inflow_turbulence.f90

    r1354 r1560  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Option recycling_yshift added. If this option is switched on, the turbulence
     23! data, which is mapped from the recycling plane to the inflow, is shifted in
     24! y direction (by ny * dy / 2 )
    2325!
    2426! Former revisions:
     
    5961       
    6062    USE control_parameters,                                                    &
    61         ONLY:  recycling_plane
     63        ONLY:  recycling_plane, recycling_yshift
    6264       
    6365    USE cpulog,                                                                &
     
    8183    INTEGER(iwp) ::  k        !:
    8284    INTEGER(iwp) ::  l        !:
     85    INTEGER(iwp) ::  next     !:
    8386    INTEGER(iwp) ::  ngp_ifd  !:
    8487    INTEGER(iwp) ::  ngp_pr   !:
     88    INTEGER(iwp) ::  prev     !:
    8589
    8690    REAL(wp), DIMENSION(nzb:nzt+1,5,nbgp)           ::                         &
    8791       avpr, avpr_l  !:
    8892    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,5,nbgp) ::                         &
    89        inflow_dist   !:
     93       inflow_dist, local_inflow_dist  !:
    9094
    9195    CALL cpu_log( log_point(40), 'inflow_turbulence', 'start' )
     
    190194    IF ( myidx == id_recycling  .AND.  myidx /= id_inflow )  THEN
    191195
    192        CALL MPI_SEND( inflow_dist(nzb,nysg,1,1), ngp_ifd, MPI_REAL, &
     196       CALL MPI_SEND( inflow_dist(nzb,nysg,1,1), ngp_ifd, MPI_REAL,            &
    193197                      id_inflow, 1, comm1dx, ierr )
    194198
     
    196200
    197201       inflow_dist = 0.0_wp
    198        CALL MPI_RECV( inflow_dist(nzb,nysg,1,1), ngp_ifd, MPI_REAL, &
     202       CALL MPI_RECV( inflow_dist(nzb,nysg,1,1), ngp_ifd, MPI_REAL,            &
    199203                      id_recycling, 1, comm1dx, status, ierr )
    200204
    201205    ENDIF
     206
     207   
     208    IF ( recycling_yshift .AND. myidx == id_inflow ) THEN
     209
     210       IF ( pdims(2) >= 2 ) THEN
     211 
     212          IF ( myidy >= INT( pdims(2) / 2 ) ) THEN
     213             prev = myidy - INT( pdims(2) / 2 )
     214          ELSE
     215             prev = pdims(2) - ( INT( pdims(2) / 2 ) - myidy )
     216          ENDIF
     217       
     218          IF ( myidy < pdims(2) - INT( pdims(2) / 2 ) ) THEN
     219             next = myidy + INT( pdims(2) / 2 )
     220          ELSE
     221             next = INT( pdims(2) / 2 ) - ( pdims(2) - myidy )
     222          ENDIF
     223
     224       ENDIF
     225
     226       local_inflow_dist = 0.0_wp
     227   
     228       CALL MPI_SENDRECV( inflow_dist(nzb,nysg,1,1), ngp_ifd, MPI_REAL,        &
     229                          next, 1, local_inflow_dist(nzb,nysg,1,1), ngp_ifd,   &
     230                          MPI_REAL, prev, 1, comm1dy, status, ierr )
     231       
     232    ENDIF
     233
    202234#endif
    203235
     
    205237!-- Add the disturbance at the inflow
    206238    IF ( nxl == 0 )  THEN
    207 
    208        DO  j = nysg, nyng
    209           DO  k = nzb, nzt + 1
    210 
    211               u(k,j,-nbgp+1:0) = mean_inflow_profiles(k,1) + &
    212                            inflow_dist(k,j,1,1:nbgp) * inflow_damping_factor(k)
    213               v(k,j,-nbgp:-1)  = mean_inflow_profiles(k,2) + &
    214                            inflow_dist(k,j,2,1:nbgp) * inflow_damping_factor(k)
    215               w(k,j,-nbgp:-1)  =                             &
     239       
     240       IF ( recycling_yshift ) THEN       
     241
     242          DO  j = nysg, nyng
     243             DO  k = nzb, nzt + 1
     244
     245                u(k,j,-nbgp+1:0) = mean_inflow_profiles(k,1) +                 &
     246                                   local_inflow_dist(k,j,1,1:nbgp) *           &
     247                                   inflow_damping_factor(k)
     248                v(k,j,-nbgp:-1)  = mean_inflow_profiles(k,2) +                 &
     249                                   local_inflow_dist(k,j,2,1:nbgp) *           &
     250                                   inflow_damping_factor(k)
     251                w(k,j,-nbgp:-1)  =                                             &
     252                                   local_inflow_dist(k,j,3,1:nbgp) *           &
     253                                   inflow_damping_factor(k)
     254                pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +                 &
     255                                   local_inflow_dist(k,j,4,1:nbgp) *           &
     256                                   inflow_damping_factor(k)
     257                e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                 &
     258                                   local_inflow_dist(k,j,5,1:nbgp) *           &
     259                                   inflow_damping_factor(k)
     260                e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
     261
     262             ENDDO
     263          ENDDO
     264
     265       ELSE
     266
     267          DO  j = nysg, nyng
     268             DO  k = nzb, nzt + 1
     269 
     270                 u(k,j,-nbgp+1:0) = mean_inflow_profiles(k,1) +                   &
     271                              inflow_dist(k,j,1,1:nbgp) * inflow_damping_factor(k)
     272                 v(k,j,-nbgp:-1)  = mean_inflow_profiles(k,2) +                   &
     273                              inflow_dist(k,j,2,1:nbgp) * inflow_damping_factor(k)
     274                 w(k,j,-nbgp:-1)  =                                               &
    216275                           inflow_dist(k,j,3,1:nbgp) * inflow_damping_factor(k)
    217               pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) + &
     276                 pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +                  &
    218277                           inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
    219               e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) + &
     278                 e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                  &
    220279                           inflow_dist(k,j,5,1:nbgp) * inflow_damping_factor(k)
    221               e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
    222 
    223           ENDDO
    224        ENDDO
    225 
    226     ENDIF
     280                 e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
     281
     282             ENDDO
     283          ENDDO
     284
     285       ENDIF
     286   
     287    ENDIF
     288
    227289
    228290    CALL cpu_log( log_point(40), 'inflow_turbulence', 'stop' )
  • palm/trunk/SOURCE/modules.f90

    r1558 r1560  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! +recycling_yshift
    2323!
    2424! Former revisions:
     
    669669                prandtl_layer = .TRUE., &
    670670                precipitation = .FALSE., &
    671                 random_heatflux = .FALSE., &
     671                random_heatflux = .FALSE., recycling_yshift = .FALSE.,&
    672672                run_control_header = .FALSE., run_coupled = .TRUE., &
    673673                scalar_rayleigh_damping = .TRUE., sloping_surface = .FALSE., &
  • palm/trunk/SOURCE/parin.f90

    r1497 r1560  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! +recycling_yshift
    2323!
    2424! Former revisions:
     
    214214               q_vertical_gradient_level, random_generator,                    &
    215215               random_heatflux, rayleigh_damping_factor,                       &
    216                rayleigh_damping_height, recycling_width, reference_state,      &
    217                residual_limit, restart_time, return_addres, return_username,   &
     216               rayleigh_damping_height, recycling_width, recycling_yshift,     &
     217               reference_state, residual_limit,                                &
     218               restart_time, return_addres, return_username,                   &
    218219               revision, rif_max, rif_min, roughness_length, runnr,            &
    219220               run_identifier, sa_surface, sa_vertical_gradient,               &
     
    307308             q_vertical_gradient, q_vertical_gradient_level,                   &
    308309             random_generator, random_heatflux,                                &
    309              rayleigh_damping_factor, rayleigh_damping_height, recycling_width,&
     310             rayleigh_damping_factor, rayleigh_damping_height,                 &
     311             recycling_width, recycling_yshift,                                &
    310312             reference_state, residual_limit,                                  &
    311313             rif_max, rif_min, roughness_length, sa_surface,                   &
  • palm/trunk/SOURCE/read_var_list.f90

    r1551 r1560  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! +recycling_yshift
    2323!
    2424! Former revisions:
     
    528528          CASE ( 'recycling_width' )
    529529             READ ( 13 )  recycling_width
     530          CASE ( 'recycling_yshift' )
     531             READ ( 13 ) recycling_yshift
    530532          CASE ( 'reference_state' )
    531533             READ ( 13 )  reference_state
Note: See TracChangeset for help on using the changeset viewer.