Ignore:
Timestamp:
Mar 22, 2007 9:54:05 AM (15 years ago)
Author:
raasch
Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:
1 edited

Legend:

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

    r73 r75  
    55! -----------------
    66! uxrp, vynp eliminated,
    7 ! dirichlet/neumann changed to dirichlet/radiation, etc.
     7! dirichlet/neumann changed to dirichlet/radiation, etc.,
     8! poisfft_init is only called if fft-solver is switched on
    89!
    910! Former revisions:
     
    493494    nzta = nz
    494495    nnz  = nz
    495 
    496 !
    497 !-- For non-cyclic boundaries extend array u (v) by one gridpoint
    498 !    IF ( bc_lr /= 'cyclic' )  uxrp = 1
    499 !    IF ( bc_ns /= 'cyclic' )  vynp = 1
    500496
    501497!
     
    757753          inflow_r  = .TRUE.
    758754       ENDIF
    759 !       uxrp      = 1
    760755    ENDIF
    761756
     
    774769          outflow_n = .TRUE.
    775770       ENDIF
    776 !       vynp      = 1
    777     ENDIF
    778 
    779 !
    780 !-- Additional MPI derived data type for the exchange of ghost points along x
    781 !-- needed in case of non-cyclic boundary conditions along y on the northmost
    782 !-- processors (for the exchange of the enlarged v array)
    783     IF ( bc_ns /= 'cyclic'  .AND.  pnorth == MPI_PROC_NULL )  THEN
    784        ngp_yz_p = ( nzt - nzb + 2 ) * ( nyn + vynp - nys + 3 )
    785        CALL MPI_TYPE_VECTOR( nxr-nxl+3, nzt-nzb+2, ngp_yz_p, &
    786                              MPI_REAL, type_xz_p, ierr )
    787        CALL MPI_TYPE_COMMIT( type_xz_p, ierr )
    788     ENDIF
     771    ENDIF
     772
    789773#else
    790774    IF ( bc_lr == 'dirichlet/radiation' )  THEN
    791775       inflow_l  = .TRUE.
    792776       outflow_r = .TRUE.
    793 !       uxrp      = 1
    794777    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
    795778       outflow_l = .TRUE.
     
    803786       outflow_n = .TRUE.
    804787       inflow_s  = .TRUE.
    805 !       vynp      = 1
    806788    ENDIF
    807789#endif
     
    809791    IF ( psolver == 'poisfft_hybrid' )  THEN
    810792       CALL poisfft_hybrid_ini
    811     ELSE
     793    ELSEIF ( psolver == 'poisfft' )  THEN
    812794       CALL poisfft_init
    813795    ENDIF
Note: See TracChangeset for help on using the changeset viewer.