Changeset 1380 for palm/trunk/SOURCE/nudging.f90

Timestamp:

Apr 28, 2014 12:40:45 PM (10 years ago)

Author:

heinze

Message:

Upper boundary conditions for pt and q in case of nudging adjusted

File:

• palm/trunk/SOURCE/nudging.f90 (modified) (3 diffs)

palm/trunk/SOURCE/nudging.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Subroutine nudge_ref added to account for proper upper scalar boundary 
+! conditions in case of nudging
 ! 
 ! Former revisions:
@@ -72 +73 @@
 
     PRIVATE
-    PUBLIC init_nudge, calc_tnudge, nudge 
+    PUBLIC init_nudge, calc_tnudge, nudge, nudge_ref
     SAVE
 
@@ -515 +516 @@
     END SUBROUTINE nudge_ij
 
+
+    SUBROUTINE nudge_ref ( time )
+
+       USE arrays_3d,                                                          &
+           ONLY:  time_vert, ptnudge, pt_init, qnudge, q_init
+
+       USE kinds
+
+
+       IMPLICIT NONE
+
+       INTEGER(iwp) ::  nt                    !:
+
+       REAL(wp)             ::  fac           !:
+       REAL(wp), INTENT(in) ::  time          !:
+
+!
+!--    Interpolation in time of NUDGING_DATA for pt_init and q_init. This is 
+!--    needed for correct upper boundary conditions for pt and q and in case that 
+!      large scale subsidence as well as scalar Rayleigh-damping are used
+       nt = 1
+       DO WHILE ( time > time_vert(nt) )
+          nt = nt + 1
+       ENDDO
+       IF ( time /= time_vert(nt) )  THEN
+        nt = nt - 1
+       ENDIF
+
+       fac = ( time-time_vert(nt) ) / ( time_vert(nt+1)-time_vert(nt) )
+
+       pt_init = ptnudge(:,nt) + fac * ( ptnudge(:,nt+1) - ptnudge(:,nt) )
+       q_init  = qnudge(:,nt) + fac * ( qnudge(:,nt+1) - qnudge(:,nt) )
+
+    END SUBROUTINE nudge_ref
+
+
  END MODULE nudge_mod