CALSTRM.f

Go to the documentation of this file.

      SUBROUTINE calstrm(Z1D,STRM)

!     
!
!     
!     
!     INCLUDE ETA GRID DIMENSIONS.  SET/DERIVE OTHER PARAMETERS.
!     
!      use vrbls2d, only:
      use params_mod, only: g
      use ctlblk_mod, only: jsta, jend, im, ista, iend
!- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      implicit none
!
      real,PARAMETER :: omega=7.292e-5,twomg=2*omega
!
!     DECLARE VARIABLES.
!     
!      LOGICAL FIRST,OLDRD,RESTRT,RUN,SIGMA,STRD
      REAL, dimension(ista:iend,jsta:jend), intent(in)    ::  z1d
      REAL, dimension(ista:iend,jsta:jend), intent(inout) ::  strm
!
      LOGICAL oldrd,strd
      integer imid,i,j
      real f0,gof0
!     
!***************************************************************************
!     START CALSTRM HERE.
!     
!     COMPUTE CORIOLIS PARAMETER AT 40N
!
      imid=im/2
      f0   = 1.454441e-4*sin(40.0*0.01745329)
      gof0 = g/f0
!     
!     COMPUTE GEOSTROPHIC STREAMFUNCTION.
!$omp  parallel do
      DO j=jsta,jend
        DO i=ista,iend
          strm(i,j) = gof0*z1d(i,j)
        ENDDO
      ENDDO
!     
!     END OF ROUTINE.
      RETURN
      END