basiy.f90

来自「CCSM Research Tools: Community Atmospher」· F90 代码 · 共 47 行

#include <misc.h>
#include <params.h>

subroutine basiy(phi     ,lbasiy  )

!----------------------------------------------------------------------- 
! 
! Purpose: 
! Compute weights used in Lagrange cubic polynomial interpolation in 
! the central interval of a four point stencil.  Done for each interval
! in the unequally spaced latitude grid.
! 
! Method: 
! 
! Author: J. Olson
! 
!-----------------------------------------------------------------------
!
! $Id: basiy.F90,v 1.1 2001/11/06 18:42:47 erik Exp $
! $Author: erik $
!
  use precision
  use pmgrid
  implicit none

!------------------------------Parameters-------------------------------
  integer, parameter ::  jfirst = nxpt + 1          ! first index to be computed
  integer, parameter ::  jlast  = platd - nxpt - 1  ! last  index to be computed
!-----------------------------------------------------------------------

!------------------------------Arguments--------------------------------
  real(r8), intent(in)  :: phi(platd)           ! grid values in extended grid
  real(r8), intent(out) :: lbasiy(4,2,platd)    ! Weights for Lagrange cubic interp
!-----------------------------------------------------------------------

!---------------------------Local variables-----------------------------
  integer jj                ! index
!-----------------------------------------------------------------------
!
  do jj = jfirst,jlast
     call lcbas( phi(jj-1),lbasiy(1,1,jj),lbasiy(1,2,jj) )
  end do
!
  return
end subroutine basiy