vertices.F90

Go to the documentation of this file.

module vertices

  use gengrid_kinds, only : dbl_kind
  use grdvars,       only : ni,nj,nv

  implicit none

contains

  subroutine fill_vertices(jbeg,jend,iVert,jVert,lat,lon,latvert,lonvert)

    integer, intent( in) :: jbeg,jend
    integer, intent( in) :: iVert(nv), jVert(nv)
    real(dbl_kind), dimension(ni,nj),    intent( in) ::  lat, lon

    real(dbl_kind), dimension(ni,nj,nv), intent(out) :: latvert,lonvert

    ! local variables
    integer :: i,j,n,ii,jj

    do j = jbeg,jend
       do i = 1,ni
          do n = 1,nv
             ii = i + ivert(n); jj = j + jvert(n)
             if(ii .eq.    0)ii = ni
             if(ii .eq. ni+1)ii = 1
             latvert(i,j,n)   = lat(ii,jj)
             lonvert(i,j,n)   = lon(ii,jj)
          enddo
       enddo
    enddo
  end subroutine fill_vertices


  subroutine fill_bottom(iVert,jVert,lat,lon,latvert,lonvert,dlat)

    integer, intent( in) :: iVert(nv), jVert(nv)
    real(dbl_kind), dimension(ni,nj),    intent( in) ::  lat, lon
    real(dbl_kind), dimension(ni),       intent( in) ::  dlat

    real(dbl_kind), dimension(ni,nj,nv), intent(out) :: latvert,lonvert

    ! local variables
    integer :: i,j,n,ii,jj

    ! fill in grid bottom (j=1)
    ! vertices 1,2 are available
    ! vertices 3,4 must be set manually
    j = 1
    do i = 1,ni
       do n = 1,2
          ii = i + ivert(n); jj = j + jvert(n)
          if(ii .eq.    0)ii = ni
          if(ii .eq. ni+1)ii = 1
          latvert(i,j,n)   = lat(ii,jj)
          lonvert(i,j,n)   = lon(ii,jj)
       enddo
       do n = 3,4
          ii = i + ivert(n)
          if(ii .eq.    0)ii = ni
          if(ii .eq. ni+1)ii = 1
          latvert(i,j, n) =  dlat(ii)
       enddo
       lonvert(i,j, 3) = lonvert(i,j,2)
       lonvert(i,j, 4) = lonvert(i,j,1)
    enddo
  end subroutine fill_bottom


  subroutine fill_top(iVert,jVert,lat,lon,latvert,lonvert,xlat,xlon)

    integer, intent( in) :: iVert(nv), jVert(nv)
    real(dbl_kind), dimension(ni,nj),    intent( in) ::  lat,  lon
    real(dbl_kind), dimension(ni),       intent( in) :: xlat, xlon

    real(dbl_kind), dimension(ni,nj,nv), intent(out) :: latvert,lonvert

    ! local variables
    integer :: i,j,n,ii,jj

    ! fill in grid top (j=nj)
    ! vertices 3,4 are available
    ! vertices 1,2 must be set manually using 'across seam' values
    j = nj
    do i = 1,ni
       do n = 3,4
          ii = i + ivert(n); jj = j + jvert(n)
          if(ii .eq.    0)ii = ni
          if(ii .eq. ni+1)ii = 1
          latvert(i,j,n)  = lat(ii,jj)
          lonvert(i,j,n)  = lon(ii,jj)
       enddo
       do n = 1,2
          ii = i + ivert(n)
          if(ii .eq.    0)ii = ni
          if(ii .eq. ni+1)ii = 1
          latvert(i,j,n)  = xlat(ii)
          lonvert(i,j,n)  = xlon(ii)
       enddo
    enddo
    !latCv_vert(i,j, 1) = latCv_vert(i,j,4)
    !latCv_vert(i,j, 2) = latCv_vert(i,j,3)
    !lonCv_vert(i,j, 1) = lonCv_vert(i,j,4)+240.d0
    !lonCv_vert(i,j, 2) = lonCv_vert(i,j,3)+240.d0
  end subroutine fill_top
end module vertices