Ignore:
Timestamp:
Mar 13, 2017 8:18:57 AM (5 years ago)
Author:
maronga
Message:

nesting extended for cloud physics quantities, bugfixes in nesting, change of default most_method

File:
1 edited

Legend:

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

    r2119 r2174  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Added support for nesting with cloud microphysics
    2323!
    2424! Former revisions:
     
    250250
    251251    USE arrays_3d,                                                             &
    252         ONLY:  diss, dzu, e, e_p, nr_p, prho, pt, pt_p, pt_init, q_init, q,    &
    253                ql, ql_c, ql_v, ql_vp, qr_p, q_p, ref_state, rho_ocean, s, s_p, sa_p, &
    254                tend, u, u_p, v, vpt, v_p, w, w_p
     252        ONLY:  diss, dzu, e, e_p, nr, nr_p, prho, pt, pt_p, pt_init, q_init, q,&
     253               ql, ql_c, ql_v, ql_vp, qr, qr_p, q_p, ref_state, rho_ocean,    &
     254               s, s_p, sa_p, tend, u, u_p, v, vpt, v_p, w, w_p
    255255
    256256    USE calc_mean_profile_mod,                                                 &
     
    603603             CALL pmci_datatrans( nesting_mode )
    604604
    605              IF ( TRIM( nesting_mode ) == 'two-way' .OR.                               &
     605             IF ( TRIM( nesting_mode ) == 'two-way' .OR.                       &
    606606                  nesting_mode == 'vertical' )  THEN
    607607!
     
    611611                CALL exchange_horiz( v, nbgp )
    612612                CALL exchange_horiz( w, nbgp )
    613                 IF ( .NOT. neutral            )  CALL exchange_horiz( pt, nbgp )
    614                 IF ( humidity                 )  CALL exchange_horiz( q, nbgp  )
    615                 IF ( passive_scalar           )  CALL exchange_horiz( s, nbgp  )
    616                 IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp  )
     613                IF ( .NOT. neutral )  CALL exchange_horiz( pt, nbgp )
     614
     615                IF ( humidity )  THEN
     616
     617                   CALL exchange_horiz( q, nbgp )
     618
     619                   IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
     620!                        CALL exchange_horiz( qc, nbgp )
     621                       CALL exchange_horiz( qr, nbgp )
     622!                        CALL exchange_horiz( nc, nbgp )
     623                       CALL exchange_horiz( nr, nbgp )
     624                   ENDIF
     625
     626                ENDIF
     627
     628                IF ( passive_scalar )  CALL exchange_horiz( s, nbgp )
     629                IF ( .NOT. constant_diffusion )  CALL exchange_horiz( e, nbgp )
    617630             ENDIF
    618631!
     
    686699             CALL calc_liquid_water_content
    687700          ENDIF
    688 ! 
     701!
    689702!--       If required, compute virtual potential temperature
    690703          IF ( humidity )  THEN
Note: See TracChangeset for help on using the changeset viewer.