Changeset 1276 for palm/trunk/SOURCE/ls_forcing.f90

Timestamp:

Jan 15, 2014 1:40:41 PM (10 years ago)

Author:

heinze

Message:

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

File:

• palm/trunk/SOURCE/ls_forcing.f90 (modified) (1 diff)

palm/trunk/SOURCE/ls_forcing.f90

@@ -238 +238 @@
        fac = ( time -time_surf(t) ) / ( time_surf(t+1) - time_surf(t) )
 
-       shf              = shf_surf(t) + fac * ( shf_surf(t+1) - shf_surf(t) )
-       qsws             = qsws_surf(t) + fac * ( qsws_surf(t+1) - qsws_surf(t) )
-       pt_surface       = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
+       IF ( ibc_pt_b == 0 )  THEN
+!
+!--       In case of Dirichlet boundary condition shf must not 
+!--       be set - it is calculated via MOST in prandtl_fluxes
+          pt_surface = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
+
+       ELSEIF ( ibc_pt_b == 1 )  THEN
+!
+!--       In case of Neumann boundary condition pt_surface is needed for 
+!--       calculation of reference density
+          shf        = shf_surf(t) + fac * ( shf_surf(t+1) - shf_surf(t) )
+          pt_surface = pt_surf(t) + fac * ( pt_surf(t+1) - pt_surf(t) )
+
+       ENDIF
+
+       IF ( ibc_q_b == 0 )  THEN
+!
+!--       In case of Dirichlet boundary condition qsws must not 
+!--       be set - it is calculated via MOST in prandtl_fluxes
+          q_surface = q_surf(t) + fac * ( q_surf(t+1) - q_surf(t) )
+
+       ELSEIF ( ibc_pt_b == 1 )  THEN
+
+          qsws = qsws_surf(t) + fac * ( qsws_surf(t+1) - qsws_surf(t) )
+
+       ENDIF
+
        surface_pressure = p_surf(t) + fac * ( p_surf(t+1) - p_surf(t) )