Changeset 2104 for palm/trunk/SOURCE


Ignore:
Timestamp:
Jan 6, 2017 4:01:15 PM (7 years ago)
Author:
knoop
Message:

Bugfix for large scale forcing.

File:
1 edited

Legend:

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

    r2101 r2104  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Bugfix for approximation related flux input conversion
    2323!
    2424! Former revisions:
     
    115115           ONLY:  p_surf, pt_surf, q_surf, qsws_surf, shf_surf, td_lsa_lpt,    &
    116116                  td_lsa_q, td_sub_lpt, td_sub_q, time_surf, time_vert,        &
    117                   heatflux_input_conversion, waterflux_input_conversion,       &
    118117                  ug_vert, vg_vert, wsubs_vert, zu
    119118
     
    226225          CALL message( 'ls_forcing', 'PA0371', 0, 0, 0, 6, 0 )
    227226          lsf_surf = .FALSE.
    228        ELSE
    229           shf_surf  = shf_surf  * heatflux_input_conversion(nzb)
    230           qsws_surf = qsws_surf * waterflux_input_conversion(nzb)
    231227       ENDIF
    232228
     
    358354       USE arrays_3d,                                                          &
    359355           ONLY:  p_surf, pt_surf, q_surf, qsws, qsws_surf, shf, shf_surf,     &
     356                  heatflux_input_conversion, waterflux_input_conversion,       &
    360357                  time_surf, time_vert, ug, ug_vert, vg, vg_vert
    361358
     
    363360           ONLY:  bc_q_b, ibc_pt_b, ibc_q_b, pt_surface, q_surface,            &
    364361                  surface_pressure
     362
     363       USE indices,                                                            &
     364           ONLY:  nzb
    365365
    366366       USE kinds
     
    395395!--       In case of Neumann boundary condition pt_surface is needed for
    396396!--       calculation of reference density
    397           shf        = shf_surf(nt) + fac * ( shf_surf(nt+1) - shf_surf(nt) )
     397          shf        = ( shf_surf(nt) + fac * ( shf_surf(nt+1) - shf_surf(nt) )&
     398                       ) * heatflux_input_conversion(nzb)
    398399          pt_surface = pt_surf(nt) + fac * ( pt_surf(nt+1) - pt_surf(nt) )
    399400
     
    408409       ELSEIF ( ibc_q_b == 1 )  THEN
    409410
    410           qsws = qsws_surf(nt) + fac * ( qsws_surf(nt+1) - qsws_surf(nt) )
     411          qsws = ( qsws_surf(nt) + fac * ( qsws_surf(nt+1) - qsws_surf(nt) )   &
     412                 ) * waterflux_input_conversion(nzb)
    411413
    412414       ENDIF
Note: See TracChangeset for help on using the changeset viewer.