vertvel.f File Reference

Go to the source code of this file.

Functions/Subroutines
real function	dubyou (UT, VT, I, J, K)

Function/Subroutine Documentation

real function dubyou	(		UT,
			VT,
			I,
			J,
			K
	)

Definition at line 3 of file vertvel.f.

Referenced by bal5().

C********************************************************************
C
C  GETTING VERTICAL MOTION W FROM HORIZONTAL MOTION & FLOW SFC SLOPE. 
C       F LUDWIG 11/89
C
C  REVISED AS A FUNCTION 3/2000
C
          include 'NGRIDS.PAR'
C
          include 'ANCHOR.CMM'
          include 'FLOWER.CMM'
          include 'LIMITS.CMM'
C
    SAVE
C
C  GETTING FLOW SURFACE SLOPE
C
    IF (k.GT.1) THEN
       IF (rhs(i,j,k).GT.0.0) THEN
C
C  CHECK TO SEE IF WE ARE AT EAST OR WEST EDGE, THEN GET FLOW SURFACE 
C  HEIGHTS TO EAST & WEST OF THIS PT.  HTS ABOVE SFC (RHS)
C  ARE ADDED TO SFCHT TO GET MSL HT.  SURFACE HEIGHT IS USED FOR
C  1ST LEVEL (K=1) WHICH REPRESENTS ANEMOMETER HEIGHT WINDS.
C
C
          IF (i.GT.1 .AND. i.LT.ncol) THEN
             dgridx=2000.0*dscrs
             zp10=max(rhs(i+1,j,k),0.0)
             zm10=max(rhs(i-1,j,k),0.0)
             zp10=0.0
             zm10=0.0
             hsige=sfcht(i+1,j)+zp10
             hsigw=sfcht(i-1,j)+zm10
C
C  USE ADJACENT PTS IF AT EDGE
C
          ELSE IF (i.EQ.1) THEN
             dgridx=1000.0*dscrs
         zp10=max(rhs(i+1,j,k),0.0)
         zm10=max(rhs(i,j,k),0.0)
         hsige=sfcht(i+1,j)+zp10
         hsigw=sfcht(i,j)+zm10
          ELSE IF (i.EQ.ncol) THEN
             dgridx=1000.0*dscrs
         zp10=max(rhs(i,j,k),0.0)
         zm10=max(rhs(i-1,j,k),0.0)
         hsige=sfcht(i,j)+zp10
         hsigw=sfcht(i-1,j)+zm10
          END IF
C
C  NOW DO THE SAME FOR  NORTH SOUTH
C
          IF (j.GT.1 .AND. j.LT.nrow) THEN
             dgridy=2000.0*dscrs
             z0p1=max(rhs(i,j+1,k),0.0)
             z0m1=max(rhs(i,j-1,k),0.0)
             hsign=sfcht(i,j+1)+z0p1
             hsigs=sfcht(i,j-1)+z0m1
C
C  USE ADJACENT PTS IF AT EDGE
C
          ELSE IF (j.EQ.1) THEN
             dgridy=1000.0*dscrs
         z0p1=max(rhs(i,j+1,k),0.0)
         z0m1=max(rhs(i,j,k),0.0)
         hsign=sfcht(i,j+1)+z0p1
         hsigs=sfcht(i,j)+z0m1
          ELSE IF (j.EQ.nrow) THEN
             dgridy=1000.0*dscrs
             z0p1=max(rhs(i,j,k),0.0)
             z0m1=max(rhs(i,j-1,k),0.0)
             hsign=sfcht(i,j)+z0p1
             hsigs=sfcht(i,j-1)+z0m1
          END IF
C
C  NOW GET W-E & S-N SLOPES
C
          dhdx=(hsige-hsigw)/dgridx
          dhdy=(hsign-hsigs)/dgridy
C
C  VERTICAL MOTION IS INNER PRODUCT OF HORIZ WIND & HT GRADIENT
C  AT ABOVE GROUND PTS.
C
          dubyou=ut*dhdx+vt*dhdy
C
       ELSE
C
C  VERTICAL MOTION IS ZERO FOR BELOW GROUND PTS
C
          dubyou=0.0
C
       END IF
C
C  USE SFC HTS TO GET VERTICAL MOTION IN 1ST (ANEMOMETER HT) LAYER
C
    ELSE
       IF (i.GT.1 .AND. i.LT.ncol) THEN
          dgridx=2000.0*dscrs
          hsige=sfcht(i+1,j)
          hsigw=sfcht(i-1,j)
C
C  USE ADJACENT PTS IF AT EDGE
C
       ELSE IF (i.EQ.1) THEN
          dgridx=1000.0*dscrs
          hsige=sfcht(i+1,j)
          hsigw=sfcht(i,j)
       ELSE IF (i.EQ.ncol) THEN
          dgridx=1000.0*dscrs
          hsige=sfcht(i,j)
          hsigw=sfcht(i-1,j)
       END IF
C
C  NOW DO THE SAME FOR  NORTH SOUTH
C
       IF (j.GT.1 .AND. j.LT.nrow) THEN
          dgridy=2000.0*dscrs
          hsign=sfcht(i,j+1)
          hsigs=sfcht(i,j-1)
C
C  USE ADJACENT PTS IF AT EDGE
C
       ELSE IF (j.EQ.1) THEN
          dgridy=1000.0*dscrs
          hsign=sfcht(i,j+1)+z0p1
          hsigs=sfcht(i,j)+z0m1
       ELSE IF (j.EQ.nrow) THEN
          dgridy=1000.0*dscrs
          hsign=sfcht(i,j)+z0p1
          hsigs=sfcht(i,j-1)+z0m1
       END IF
C
C  NOW GET W-E & S-N SLOPES
C
       dhdx=(hsige-hsigw)/dgridx
       dhdy=(hsign-hsigs)/dgridy
C
C  VERTICAL MOTION IS INNER PRODUCT OF HORIZ WIND & HT GRADIENT
C  OF SFC HTS.
C
       dubyou=ut*dhdx+vt*dhdy
C
    END IF
C
    RETURN
C

Here is the caller graph for this function: