Ignore:
Timestamp:
Apr 13, 2018 3:09:23 PM (3 years ago)
Author:
suehring
Message:

Bugfix in large-scale forcing

File:
1 edited

Legend:

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

    r2938 r2970  
    2525! -----------------
    2626! $Id$
     27! Bugfix in old large-scale forcing mode
     28!
     29! 2938 2018-03-27 15:52:42Z suehring
    2730! Further improvements for nesting in larger-scale model
    2831!
     
    6265    USE control_parameters,                                                    &
    6366        ONLY:  bc_lr, bc_ns, bc_pt_b, bc_q_b, constant_diffusion,              &
     67               constant_heatflux, constant_waterflux,                          &
    6468               data_output_pr, dt_3d, end_time, forcing,                       &
    6569               force_bound_l, force_bound_n, force_bound_r, force_bound_s,     & 
    66                humidity, intermediate_timestep_count, ibc_pt_b, ibc_q_b,       &
     70               humidity, initializing_actions, intermediate_timestep_count,    &
     71               ibc_pt_b, ibc_q_b,                                              &
    6772               large_scale_forcing, large_scale_subsidence, lsf_surf, lsf_vert,&
    6873               lsf_exception, message_string, neutral, nudging, passive_scalar,&
    69                pt_surface, ocean, q_surface, surface_pressure, topography,     &
     74               pt_surface, ocean, q_surface, surface_heatflux,                 &
     75               surface_pressure, surface_waterflux, topography,                &
    7076               use_subsidence_tendencies
    7177
     
    729735
    730736       IF ( large_scale_forcing  .AND.  ( bc_lr /= 'cyclic'  .OR.              &
    731                                        bc_ns /= 'cyclic' ) )  THEN
     737                                          bc_ns /= 'cyclic' ) )  THEN
    732738          message_string = 'Non-cyclic lateral boundaries do not allow for &'//&
    733739                        'the usage of large scale forcing from external file.'
     
    14721478
    14731479       ENDIF
     1480!
     1481!--    Surface heat- and waterflux will be written later onto surface elements
     1482       IF ( .NOT.  neutral  .AND.  constant_heatflux  .AND.                    &
     1483            TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
     1484             surface_heatflux = shf_surf(1)
     1485       ENDIF
     1486       IF ( humidity  .AND.  constant_waterflux  .AND.                         &
     1487            TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
     1488             surface_waterflux = qsws_surf(1)
     1489       ENDIF
    14741490
    14751491       surface_pressure = p_surf(nt) + fac * ( p_surf(nt+1) - p_surf(nt) )
Note: See TracChangeset for help on using the changeset viewer.