?? basdy.f90
字號(hào):
#include <misc.h>#include <params.h>subroutine basdy(phi ,lbasdy )!----------------------------------------------------------------------- ! ! Purpose: ! Compute weights for the calculation of derivative estimates at the two! center points of the four point stencil for each interval in the! unequally spaced latitude grid. Estimates are from differentiating! a Lagrange cubic polynomial through the four point stencil.! ! Method: ! phi Latitude values in the extended grid.! lbasdy Weights for derivative estimates based on Lagrange cubic! polynomial on the unequally spaced latitude grid.! If grid interval j (in extended grid) is surrounded by! a 4 point stencil, then the derivative at the "bottom"! of the interval uses the weights lbasdy(1,1,j),! lbasdy(2,1,j), lbasdy(3,1,j), and lbasdy(4,1,j).! The derivative at the "top" of the interval! uses lbasdy(1,2,j), lbasdy(2,2,j), lbasdy(3,2,j),! and lbasdy(4,2,j).! ! Author: J. Olson! !-----------------------------------------------------------------------!! $Id: basdy.F90,v 1.1 2001/11/06 18:42:46 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) ! latitude coordinates of model grid real(r8), intent(out) :: lbasdy(4,2,platd) ! derivative estimate weights!-----------------------------------------------------------------------!---------------------------Local variables----------------------------- integer jj ! index!-----------------------------------------------------------------------! do jj = jfirst,jlast call lcdbas( phi(jj-1), lbasdy(1,1,jj), lbasdy(1,2,jj) ) end do! returnend subroutine basdy
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -