cressman差值，程序、说明及例子
程序代码如下：
subroutine interp_cressman(lati,loni,data,nxi,nyi, !此处给定了nxi,nyi,nxo,nyo目的是给据网格上的i、j在主程序是方便计算经纬度的值，如经纬
* lato,lono,datao,nxo,nyo) !给定了 是否就不用再给出nxi,nyi,nxo,nyo呢？而只要给出待插值网格的循环次数呢？
!datao 是插值格点数据 !自己改写主程序和子程序的时候，可以使fvcom的经纬变成一维的L，就不必给出nxo,nyo了
c**** interpolates fields to another grid as specified
c**** by lat., lon. values.
c**** uses a Cressman interpolation with a specified
c**** search radius
c**** uses output from lookup.f to determine search points
parameter(pi=3.14159,num=500)
real lati(nxi,nyi),loni(nxi,nyi)
real lato(nxo,nyo),lono(nxo,nyo)
real data(nxi,nyi)
integer ilon(num),ilat(num)
real datao(nxo,nyo)
open(30,file='lookup.tab1',status='unknown')
open(31,file='lookup.tab2',status='unknown')
guess = 100. ! search radius
xmiss = -999.
scale = 1.
c************************************************************
c reinterpolate to NMC Octagonal grid using Cressman (应用cressman方法，逐步订正到NMC网格上)
c weights with a specified search radius
do 100 jj=1,nyo
do 100 ii=1,nxo
sum1=0.
sum2=0.
ilevs=0
c read in indices from lookup table
read(30,*) ilevs
c process indices
if(ilevs.ne.0.)then
read(31,*) (ilon(i),i=1,ilevs),(ilat(i),i=1,ilevs)
!ilon,ilat没有定义？ 有定义
do 90 kk=1,ilevs
i = ilon(kk)
j = ilat(kk)
c find distance between NMC point to be interpolated to
c and original point.
dim=(sin(pi*lati(i,j)/180.)*sin(pi*lato(ii,jj)/180.))
dam=(cos(pi*lati(i,j)/180.)*cos(pi*lato(ii,jj)/180.))
dam1=lono(ii,jj)-loni(i,j)
dam1=cos(pi*dam1/180.)
dist=dim+dam*dam1
dist=acos(dist)*110.949*(180./pi)
c dim=(sin(lati(i,j))*sin(lato(ii,jj)))
c dam=(cos(l
2021-11-23 19:37:51
1.78MB
1
本文针对现有四线法测量接触电阻过程中接入了被测电阻两侧部分引线电阻的问题，对四线法进行了改进，设计了一种能够直接读取接触电阻值的等阻差值四线法测量电路。利用等阻差值四线法测量电路对导线电阻进行了测量，通过与安捷伦34401A台式万用表测量结果的对比，证明了该电路的可行性。
2021-11-18 21:45:37
460KB
1
两条离散曲线求差，采用拟合的方式得到两条曲线的差值。程序中画了离散程度不同的曲线求差；两条离散曲线求差，采用拟合的方式得到两条曲线的差值。程序中画了离散程度不同的曲线求差
2021-11-08 10:59:59
325B
1
Python实现列表划分求子列表和之差最小值，从长度为n的列表中随机取m个元素，将取出的m个元素重新赋值给一个list，返回列表list，'将',list,'划分为',[l for l in mi if sum(l)==maxx],'中的任意一个子列表时,与列表剩余元素之和的差为最小值，为',sum(list)-2*sum(le[0])。
2021-10-28 15:51:54
2KB
1
使用拉格朗日差值多项式和牛顿差值多项式计算f(x)在x=1.682和1.813处的近似值。（计算结果取5位小数）
2021-10-21 10:38:27
56KB
1
MATLAB证明一阶RC电路中，任一点切线的横坐标差值A（t2-t1）等于时间常数τ
2021-10-18 21:03:38
2KB
1
最大间隙问题
最大间隙问题：给定n 个实数x1, x2,... , xn，求这n 个数在实轴上相邻2 个数之间的最大差值。假设对任何实数的下取整函数耗时O(1)，设计解最大间隙问题的线性时间算法。
编程任务：对于给定的n 个实数x1, x2,...,xn,编程计算它们的最大间隙。
Input
输入数据的第1 行有1 个正整数n。接下来的1 行中有n个实数x1, x2,... , xn。
Output
程序运行结束时，将找到的最大间隙输出。
Sample Input
5
2.3 3.1 7.5 1.5 6.3
Sample Output
3.2
2021-09-24 12:51:03
840B
1