Changeset 1276 for palm/trunk/SOURCE/ls_forcing.f90
 Timestamp:
 Jan 15, 2014 1:40:41 PM (11 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/ls_forcing.f90
r1250 r1276 238 238 fac = ( time time_surf(t) ) / ( time_surf(t+1)  time_surf(t) ) 239 239 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 243 267 surface_pressure = p_surf(t) + fac * ( p_surf(t+1)  p_surf(t) ) 244 268
Note: See TracChangeset
for help on using the changeset viewer.