爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
楼主: 小小人

[源代码] 双线性插值用于Micaps

[复制链接]
 楼主| 发表于 2014-1-20 16:45:44 | 显示全部楼层

还有就是您说的已知四个站点插值一个站点,这个我还没有考虑过的。不过,我会努力的。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-1-21 11:04:33 | 显示全部楼层
嗯嗯,谢谢大家的厚爱,我下去后有改进了下这个程序,是他可以直接在终端输入时显示坐标。源代码我把它贴出来。
program main
implicit none
integer::i,j
character*3 judge
real::x=101.87,y=25.73,x1(2),y1(2)
real::w(2,2),Vq(2,2),Vr(2),Vp,g,p
DATA(x1(i),i=1,2)/99.98,102.49/
DATA(y1(j),j=1,2)/25.02,27.52/
333  write(*,*)"Please enter the Value of Q"
do i=1,2
do j=1,2
write(*,"('Vq(',I1,',',I1,')=')")i,j
        read(*,*) Vq(i,j)
        end do
        end do
g=((x-x1(1))**2)+((x-x1(2))**2)
p=((y-y1(1))**2)+((y-y1(2))**2)
do i=1,2
w(1,i)=1-(((x-x1(i))**2)/g)
end do
do j=1,2
w(2,j)=1-(((y-y1(j))**2)/p)
end do

do i=1,2
Vr(i)=Vq(1,i)*w(1,1)+Vq(2,i)*w(1,2)
end do
Vp=Vr(1)*w(2,1)+Vr(2)*w(2,2)
write(*,*)  vp
!write(*,*)'是否需要继续?,是请按Y,其他键退出'
!  read(*,*) judge
  ! if(judge=="Y".or.judge=="y") then
   ! goto 333
   !endif

end
其中,被注释的那几句,我仿照了下论坛里面的,我想把它改为循环输入,
密码修改失败请联系微信:mofangbao
发表于 2014-1-21 16:54:46 | 显示全部楼层
本帖最后由 kc121 于 2014-1-21 16:56 编辑

我不是很懂fortran,但我想
DATA(x1(i),i=1,2)/99.98,102.49/
DATA(y1(j),j=1,2)/25.02,27.52/

这里改一下,改成四个值如:
DATA(x1(i),i=1,4)/99,102,100,103/
DATA(y1(j),j=1,4)/25,27,24,26/
其他不变,这样是否能实现之前我说的任意四个站插值到一个站的问题呢?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-1-21 20:17:08 | 显示全部楼层
本帖最后由 小小人 于 2014-1-21 20:23 编辑

嗯嗯,谢谢您的关注。这样的四个站点在我的这个程序里是不行,我在micaps里的插值,插的是由两条经纬线之间那个点,用图来表示就,它是在兰伯特投影下的四个点,所以在X方向上是在同一条经度,在y方向上是在同一条维度上。所以虽然有四个点,但是,他的值分别只有两个。
插值原理图.png
插值原理图.png
密码修改失败请联系微信:mofangbao
发表于 2014-2-18 21:19:25 | 显示全部楼层
LZ有木有不用输入这四个点,直接调用生成的几个网格点数据的方法?
密码修改失败请联系微信:mofangbao
发表于 2014-4-9 16:52:50 | 显示全部楼层
贡献不够,先收藏了啊!
密码修改失败请联系微信:mofangbao
发表于 2014-4-9 16:53:08 | 显示全部楼层
贡献不够,先收藏了啊!
密码修改失败请联系微信:mofangbao
发表于 2014-4-17 06:15:50 | 显示全部楼层
呵呵,等待改进
密码修改失败请联系微信:mofangbao
发表于 2014-5-5 15:21:08 | 显示全部楼层
谢谢楼主分享
密码修改失败请联系微信:mofangbao
发表于 2014-5-21 15:32:21 | 显示全部楼层
好东东,值得分享啊!!!!!!!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表