Changeset 1276 for palm/trunk/SOURCE/prandtl_fluxes.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/prandtl_fluxes.f90 (modified) (3 diffs)

palm/trunk/SOURCE/prandtl_fluxes.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! Use LSF_DATA also in case of Dirichlet bottom boundary condition for scalars
 !
 ! Former revisions:
@@ -120 +120 @@
 !--    For a given surface temperature:
 !--    (the Richardson number is still the one from the previous time step)
+   
+       IF ( large_scale_forcing .AND. lsf_surf )  THEN
+          pt(0,:,:) = pt_surface
+       ENDIF
+
        !$OMP PARALLEL DO PRIVATE( a, b, k, z_p )
        !$acc kernels loop
@@ -334 +339 @@
        ELSE
           coupled_run = ( coupling_mode == 'atmosphere_to_ocean' .AND. run_coupled )
+
+           IF ( large_scale_forcing .AND. lsf_surf )  THEN
+              q(0,:,:) = q_surface
+           ENDIF
+
           !$OMP PARALLEL DO PRIVATE( a, b, k, z_p )
           !$acc kernels loop independent