Ignore:
Timestamp:
Aug 21, 2017 2:59:59 PM (4 years ago)
Author:
kanani
Message:

Vertical nesting implemented (SadiqHuq?)

File:
1 edited

Legend:

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

    r2298 r2365  
    2525! ------------------
    2626! $Id$
     27! Vertical nesting implemented (SadiqHuq)
     28!
     29! 2298 2017-06-29 09:28:18Z raasch
    2730! MPI2 coupling removed
    2831!
     
    6770    USE pegrid
    6871
     72    USE vertical_nesting_mod
     73
    6974    IMPLICIT NONE
    7075
     
    96101       IF ( TRIM( coupling_mode ) == 'coupled_run' )  THEN
    97102          i = 1
     103       ELSEIF ( TRIM( coupling_mode ) == 'vnested_twi' )  THEN
     104          i = 9
    98105       ELSE
    99106          i = 0
     
    111118    i = bc_data(0)
    112119
    113     IF ( i == 0 ) THEN
     120    IF ( i == 0 )  THEN
    114121       coupling_mode = 'uncoupled'
    115122!
     
    121128          CLOSE ( 90 )
    122129       ENDIF
     130    ELSEIF ( i == 9 )  THEN
     131
     132!
     133!--    Set a flag to identify runs with vertical nesting
     134       vnested = .TRUE.
     135       
     136       comm_inter = MPI_COMM_WORLD
     137       
     138!
     139!--    Split the total available PE's into two groups
     140!--    numprocs for Coarse and Fine Grid are specified via mrun argument -N
     141       IF ( myid < bc_data(1) )  THEN
     142          inter_color     = 0
     143          numprocs        = bc_data(1)
     144          coupling_mode   = 'vnested_crse'
     145       ELSE
     146          inter_color     = 1
     147          numprocs        = bc_data(2)
     148          coupling_mode   = 'vnested_fine'
     149       ENDIF
     150       
     151       CALL MPI_COMM_SPLIT( MPI_COMM_WORLD, inter_color, 0, comm_palm, ierr )
     152       comm2d = comm_palm
     153       
     154       OPEN( 90, FILE='VNESTING_PORT_OPENED', FORM='FORMATTED' )
     155       WRITE ( 90, '(''TRUE'')' )
     156       CLOSE ( 90 )
     157     
    123158    ELSE
    124159       comm_inter = MPI_COMM_WORLD
     
    154189    ENDIF
    155190
     191    IF (  TRIM( coupling_mode ) == 'vnested_fine' )  THEN
     192!
     193!-- Set file extension for vertical nesting
     194       coupling_char = '_N'
     195    ENDIF
     196
    156197 END SUBROUTINE init_coupling
Note: See TracChangeset for help on using the changeset viewer.