Ignore:
Timestamp:
Feb 25, 2016 12:31:13 PM (5 years ago)
Author:
hellstea
Message:

Introduction of nested domain system

File:
1 edited

Legend:

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

    r1748 r1762  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Introduction of nested domain feature
    2222!
    2323! Former revisions:
     
    113113        ONLY:  coupling_char, coupling_mode, do2d_at_begin, do3d_at_begin,     &
    114114               io_blocks, io_group, large_scale_forcing, message_string,       &
    115                nudging, simulated_time, simulated_time_chr,                    &
     115               nest_domain, nudging, simulated_time, simulated_time_chr,       &
    116116               user_interface_current_revision,                                &
    117117               user_interface_required_revision, version, wall_heatflux,       &
     
    144144
    145145    USE pegrid
     146
     147#if defined( PMC_ACTIVE )
     148    USE pmc_interface,                                                         &
     149        ONLY:  cpl_id, pmci_init, pmci_modelconfiguration
     150#endif
    146151
    147152    USE statistics,                                                            &
     
    176181!-- it will be defined in init_pegrid but is used before in cpu_log.
    177182    CALL MPI_INIT( ierr )
     183
     184#if defined( PMC_ACTIVE )
     185!
     186!-- Initialize the coupling for nested-domain runs
     187    CALL pmci_init( comm_palm )
     188    comm2d = comm_palm
     189
     190    IF ( cpl_id >= 2 )  THEN
     191       nest_domain = .TRUE.
     192       WRITE( coupling_char, '(A1,I1.1)') '_', cpl_id
     193    ENDIF
     194
     195    CALL MPI_COMM_SIZE( comm_palm, numprocs, ierr )
     196    CALL MPI_COMM_RANK( comm_palm, myid, ierr )
     197#else
    178198    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, numprocs, ierr )
    179199    CALL MPI_COMM_RANK( MPI_COMM_WORLD, myid, ierr )
    180200    comm_palm = MPI_COMM_WORLD
    181201    comm2d    = MPI_COMM_WORLD
    182 
    183 !
    184 !-- Initialize PE topology in case of coupled runs
     202!
     203!-- Initialize PE topology in case of coupled atmosphere-ocean runs (comm_palm
     204!-- will be splitted in init_coupling)
    185205    CALL init_coupling
     206#endif
    186207#endif
    187208
     
    291312!-- Initialize all necessary variables
    292313    CALL init_3d_model
     314
     315#if defined( PMC_ACTIVE )
     316!
     317!-- Coupling protocol setup for nested-domain runs
     318    CALL pmci_modelconfiguration
     319#endif
    293320
    294321!
Note: See TracChangeset for help on using the changeset viewer.