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