Ignore:
Timestamp:
Sep 14, 2012 2:35:53 PM (9 years ago)
Author:
raasch
Message:

subdomains must have identical size, i.e. grid_matching = "match" not allowed any more
parameter grid_matching removed
some obsolete variables removed

File:
1 edited

Legend:

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

    r779 r1003  
    44! Current revisions:
    55! -----------------
    6 !
     6! adjustment of array tend for cases with unequal subdomain sizes removed
    77!
    88! Former revisions:
     
    279279    ELSE
    280280       !$OMP PARALLEL DO SCHEDULE( STATIC )
    281        DO  i = nxl, nxra
    282           DO  j = nys, nyna
    283              DO  k = nzb+1, nzta
     281       DO  i = nxl, nxr
     282          DO  j = nys, nyn
     283             DO  k = nzb+1, nzt
    284284                d(k,j,i) = 0.0
    285285             ENDDO
     
    409409
    410410!
    411 !--    Enlarge the size of tend, used as a working array for the transpositions
    412        IF ( nxra > nxr  .OR.  nyna > nyn  .OR.  nza > nz )  THEN
    413           DEALLOCATE( tend )
    414           ALLOCATE( tend(1:nza,nys:nyna,nxl:nxra) )
    415        ENDIF
    416 
    417 !
    418411!--    Solve Poisson equation via FFT and solution of tridiagonal matrices
    419412       IF ( psolver == 'poisfft' )  THEN
     
    427420!--       are some optimization problems in poisfft
    428421          CALL poisfft_hybrid( d )
    429        ENDIF
    430 
    431 !
    432 !--    Resize tend to its normal size
    433        IF ( nxra > nxr  .OR.  nyna > nyn  .OR.  nza > nz )  THEN
    434           DEALLOCATE( tend )
    435           ALLOCATE( tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    436422       ENDIF
    437423
Note: See TracChangeset for help on using the changeset viewer.