CALPOT.f

Go to the documentation of this file.

      SUBROUTINE calpot(P1D,T1D,THETA)

!     
      use ctlblk_mod, only: jsta, jend, spval, im,  ista, iend
!- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
      implicit none
!     
!     SET REQUIRED CONSTANTS.
      real,PARAMETER :: capa=0.28589641,p1000=1000.e2
!
!     DECLARE VARIABLES.
!     
      real,dimension(ista:iend,jsta:jend),intent(in)    :: p1d,t1d
      real,dimension(ista:iend,jsta:jend),intent(inout) :: theta

      integer i,j
!     
!**********************************************************************
!     START CALPOT HERE.
!     
!     COMPUTE THETA
!     
!$omp parallel do private(i,j)
      DO j=jsta,jend
        DO i=ista,iend
          IF(t1d(i,j) < spval) THEN
!           IF(ABS(P1D(I,J)) > 1.0) THEN
            IF(p1d(i,j) > 1.0) THEN
              theta(i,j) = t1d(i,j) * (p1000/p1d(i,j))**capa
            ELSE
              theta(i,j) = 0.0
            ENDIF
          ELSE
            theta(i,j) = spval
          ENDIF
        ENDDO
      ENDDO
!     do j = 180, 185
!        print *, ' me, j, p1d,t1d,theta = ',
!    *   me, j, p1d(10,j),t1d(10,j),theta (10,j)
!     end do
!       stop
!     
!     END OF ROUTINE.
!     
      RETURN
      END