Ignore:
Timestamp:
Oct 10, 2007 12:03:15 AM (14 years ago)
Author:
raasch
Message:

preliminary updates for implementing buildings in poismg

File:
1 edited

Legend:

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

    r110 r114  
    44! Actual revisions:
    55! -----------------
     6! Allocation of wall flag arrays for multigrid solver
    67! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
    78!
     
    915916    ENDIF
    916917
     918!
     919!-- Allocate wall flag arrays used in the multigrid solver
     920    IF ( psolver == 'multigrid' )  THEN
     921
     922       DO  i = maximum_grid_level, 1, -1
     923
     924           SELECT CASE ( i )
     925
     926              CASE ( 1 )
     927                 ALLOCATE( wall_flags_1(nzb:nzt_mg(i)+1,         &
     928                                        nys_mg(i)-1:nyn_mg(i)+1, &
     929                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     930
     931              CASE ( 2 )
     932                 ALLOCATE( wall_flags_2(nzb:nzt_mg(i)+1,         &
     933                                        nys_mg(i)-1:nyn_mg(i)+1, &
     934                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     935
     936              CASE ( 3 )
     937                 ALLOCATE( wall_flags_3(nzb:nzt_mg(i)+1,         &
     938                                        nys_mg(i)-1:nyn_mg(i)+1, &
     939                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     940
     941              CASE ( 4 )
     942                 ALLOCATE( wall_flags_4(nzb:nzt_mg(i)+1,         &
     943                                        nys_mg(i)-1:nyn_mg(i)+1, &
     944                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     945
     946              CASE ( 5 )
     947                 ALLOCATE( wall_flags_5(nzb:nzt_mg(i)+1,         &
     948                                        nys_mg(i)-1:nyn_mg(i)+1, &
     949                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     950
     951              CASE ( 6 )
     952                 ALLOCATE( wall_flags_6(nzb:nzt_mg(i)+1,         &
     953                                        nys_mg(i)-1:nyn_mg(i)+1, &
     954                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     955
     956              CASE ( 7 )
     957                 ALLOCATE( wall_flags_7(nzb:nzt_mg(i)+1,         &
     958                                        nys_mg(i)-1:nyn_mg(i)+1, &
     959                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     960
     961              CASE ( 8 )
     962                 ALLOCATE( wall_flags_8(nzb:nzt_mg(i)+1,         &
     963                                        nys_mg(i)-1:nyn_mg(i)+1, &
     964                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     965
     966              CASE ( 9 )
     967                 ALLOCATE( wall_flags_9(nzb:nzt_mg(i)+1,         &
     968                                        nys_mg(i)-1:nyn_mg(i)+1, &
     969                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     970
     971              CASE ( 10 )
     972                 ALLOCATE( wall_flags_10(nzb:nzt_mg(i)+1,        &
     973                                        nys_mg(i)-1:nyn_mg(i)+1, &
     974                                        nxl_mg(i)-1:nxr_mg(i)+1) )
     975
     976              CASE DEFAULT
     977                 IF ( myid == 0 )  PRINT*, '+++ init_pegrid: more than 10 ', &
     978                                           ' multigrid levels'
     979                 CALL local_stop
     980
     981          END SELECT
     982
     983       ENDDO
     984
     985    ENDIF
     986
    917987 END SUBROUTINE init_pegrid
Note: See TracChangeset for help on using the changeset viewer.