Ignore:
Timestamp:
Jan 15, 2014 1:40:41 PM (8 years ago)
Author:
heinze
Message:

Usage of Dirichlet bottom boundary condition for scalars in conjunction with large scale forcing enabled

File:
1 edited

Legend:

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

    r1250 r1276  
    238238       fac = ( time -time_surf(t) ) / ( time_surf(t+1) - time_surf(t) )
    239239
    240        shf              = shf_surf(t) + fac * ( shf_surf(t+1) - shf_surf(t) )
    241        qsws             = qsws_surf(t) + fac * ( qsws_surf(t+1) - qsws_surf(t) )
    242        pt_surface       = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
     240       IF ( ibc_pt_b == 0 )  THEN
     241!
     242!--       In case of Dirichlet boundary condition shf must not
     243!--       be set - it is calculated via MOST in prandtl_fluxes
     244          pt_surface = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
     245
     246       ELSEIF ( ibc_pt_b == 1 )  THEN
     247!
     248!--       In case of Neumann boundary condition pt_surface is needed for
     249!--       calculation of reference density
     250          shf        = shf_surf(t) + fac * ( shf_surf(t+1) - shf_surf(t) )
     251          pt_surface = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
     252
     253       ENDIF
     254
     255       IF ( ibc_q_b == 0 )  THEN
     256!
     257!--       In case of Dirichlet boundary condition qsws must not
     258!--       be set - it is calculated via MOST in prandtl_fluxes
     259          q_surface = q_surf(t) + fac * ( q_surf(t+1) - q_surf(t) )
     260
     261       ELSEIF ( ibc_pt_b == 1 )  THEN
     262
     263          qsws = qsws_surf(t) + fac * ( qsws_surf(t+1) - qsws_surf(t) )
     264
     265       ENDIF
     266
    243267       surface_pressure = p_surf(t) + fac * ( p_surf(t+1) - p_surf(t) )
    244268
Note: See TracChangeset for help on using the changeset viewer.