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

Timestamp:

Mar 6, 2014 1:15:21 PM (11 years ago)

Author:

heinze

Message:

enable usage of large_scale subsidence in combination with large_scale_forcing

File:

• palm/trunk/SOURCE/ls_forcing.f90 (modified) (10 diffs)

palm/trunk/SOURCE/ls_forcing.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! Ensure a zero large scale vertical velocity at the surface
+! Bugfix: typo in case of boundary condition in if-clause
 !
 ! Former revisions:
@@ -39 +40 @@
 ! ------------
 ! Calculates large scale forcings (geostrophic wind and subsidence velocity) as 
-! well as surfaces fluxes dependend on time given in an external file (LSF_DATA).
+! well as surfaces fluxes dependent on time given in an external file (LSF_DATA).
 ! Code is based in parts on DALES and UCLA-LES.
 !--------------------------------------------------------------------------------!
@@ -110 +111 @@
 
           IF ( ierrn < 0 )  THEN
-            WRITE ( message_string, * ) 'No time dependend surface variables ',&
+            WRITE ( message_string, * ) 'No time dependent surface variables ',&
                               'in&LSF_DATA for end of run found'
 
@@ -119 +120 @@
 
        IF ( time_surf(1) > end_time )  THEN
-          WRITE ( message_string, * ) 'No time dependend surface variables in ',&
+          WRITE ( message_string, * ) 'No time dependent surface variables in ',&
                                      '&LSF_DATA for end of run found - ',  &
                                      'lsf_surf is set to FALSE'
@@ -146 +147 @@
              READ ( finput, *, IOSTAT=ierrn ) hash, time_vert(t)
              IF ( ierrn < 0 )  THEN 
-                WRITE( message_string, * ) 'No time dependend vertical profiles',&
+                WRITE( message_string, * ) 'No time dependent vertical profiles',&
                                  ' in&LSF_DATA for end of run found'
                 CALL message( 'ls_forcing', 'PA0372', 1, 2, 0, 6, 0 )
@@ -158 +159 @@
           IF ( ierrn /= 0 )  THEN
              message_string = 'errors in file LSF_DATA'
-             CALL message( 'nudging', 'PA0369', 1, 2, 0, 6, 0 )
+             CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
           ENDIF
 
@@ -166 +167 @@
           IF ( ierrn /= 0 )  THEN
              message_string = 'errors in file LSF_DATA'
-             CALL message( 'nudging', 'PA0369', 1, 2, 0, 6, 0 )
+             CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
           ENDIF
 
@@ -183 +184 @@
                 IF ( ierrn /= 0 )  THEN
                    message_string = 'errors in file LSF_DATA'
-                   CALL message( 'nudging', 'PA0369', 1, 2, 0, 6, 0 )
+                   CALL message( 'ls_forcing', 'PA0369', 1, 2, 0, 6, 0 )
                 ENDIF
 
@@ -200 +201 @@
        ENDDO 
 
+!
+!--    Large scale vertical velocity has to be zero at the surface
+       wsubs_vert(nzb,:) = 0.0
+ 
        IF ( time_vert(1) > end_time )  THEN
           WRITE ( message_string, * ) 'Time dependent large scale profile ',&
@@ -262 +267 @@
           q_surface = q_surf(t) + fac * ( q_surf(t+1) - q_surf(t) )
 
-       ELSEIF ( ibc_pt_b == 1 )  THEN
+       ELSEIF ( ibc_q_b == 1 )  THEN
 
           qsws = qsws_surf(t) + fac * ( qsws_surf(t+1) - qsws_surf(t) )