Ignore:
Timestamp:
Mar 17, 2017 11:07:39 AM (4 years ago)
Author:
hellstea
Message:

Nesting bugfixes

File:
1 edited

Legend:

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

    r2177 r2178  
    161161
    162162    USE control_parameters,                                                    &
    163         ONLY:  constant_diffusion, coupling_char, coupling_mode,               &
     163        ONLY:  cloud_physics, constant_diffusion, coupling_char, coupling_mode,&
    164164               do2d_at_begin, do3d_at_begin, humidity, initializing_actions,   &
    165165               io_blocks, io_group,                                            &
    166                large_scale_forcing, message_string, nest_domain, neutral,      &
     166               large_scale_forcing, message_string, microphysics_seifert,      &
     167               nest_domain, neutral,                                           &
    167168               nudging, passive_scalar, simulated_time, simulated_time_chr,    &
    168169               urban_surface,                                                  &
     
    203204    USE pmc_interface,                                                         &
    204205        ONLY:  cpl_id, nested_run, pmci_child_initialize, pmci_init,           &
    205                pmci_modelconfiguration, pmci_parent_initialize
     206               pmci_modelconfiguration, pmci_parent_initialize,                &
     207               pmci_ensure_nest_mass_conservation
    206208
    207209    USE radiation_model_mod,                                                   &
     
    352354!--       Send initial condition data from parent to children
    353355          CALL pmci_parent_initialize
    354        END IF
    355356!
    356357!--    Exchange_horiz is needed after the nest initialization
    357        IF ( nest_domain )  THEN
    358           CALL exchange_horiz( u, nbgp )
    359           CALL exchange_horiz( v, nbgp )
    360           CALL exchange_horiz( w, nbgp )
    361           IF ( .NOT. neutral )  THEN
    362              CALL exchange_horiz( pt, nbgp )
     358          IF ( nest_domain )  THEN
     359             CALL exchange_horiz( u, nbgp )
     360             CALL exchange_horiz( v, nbgp )
     361             CALL exchange_horiz( w, nbgp )
     362             IF ( .NOT. neutral )  THEN
     363                CALL exchange_horiz( pt, nbgp )
     364             ENDIF
     365             IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp )
     366             IF ( humidity )  THEN
     367                CALL exchange_horiz( q, nbgp )
     368                IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
     369!                   CALL exchange_horiz( qc, nbgp )
     370                   CALL exchange_horiz( qr, nbgp )
     371!                   CALL exchange_horiz( nc, nbgp )
     372                   CALL exchange_horiz( nr, nbgp )
     373                ENDIF
     374             ENDIF
     375             IF ( passive_scalar )  CALL exchange_horiz( s, nbgp )
     376
     377             CALL pmci_ensure_nest_mass_conservation
     378             CALL pres
    363379          ENDIF
    364           IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp )
    365           IF ( humidity       )  CALL exchange_horiz( q, nbgp )
    366           IF ( passive_scalar )  CALL exchange_horiz( s, nbgp )
    367380       ENDIF
    368381
Note: See TracChangeset for help on using the changeset viewer.