python读取grd数据_一个关于grads显示grd格式数据的问题。

是合成差值场的t检验程序。

program main

implicit none

integer i,j,t,irec

integer,parameter :: x=145

integer,parameter :: y=73

integer,parameter :: n=58 !总体样本数

real,parameter :: n1=10.0 !偏早年份

real,parameter :: n2=10.0 !偏晚年份

real :: dat(x,y,n),earlyave(x,y),lateave(x,y),earlyvar(x,y),latevar(x,y),ts(x,y),dif(x,y)

!资料读取

open(100,file='d:\48-05u200june.grd',form='unformatted',access='direct',recl=x*y*4)

irec=1

do t=1,n

read(100,rec=irec) ((dat(i,j,t),i=1,x),j=1,y)

irec=irec+1

enddo

close(100)

!资料读取完毕

!偏早年和偏晚年的均值

do j=1,y

do i=1,x

earlyave(i,j)=(dat(i,j,9)+dat(i,j,14)+dat(i,j,24)+dat(i,j,27)+dat(i,j,33)+dat(i,j,35)+dat(i,j,37)+dat(i,j,42)+dat(i,j,44)+dat(i,j,52))/n1

lateave(i,j)=(dat(i,j,11)+dat(i,j,12)+dat(i,j,17)+dat(i,j,18)+dat(i,j,19)+dat(i,j,20)+dat(i,j,22)+dat(i,j,45)+dat(i,j,51)+dat(i,j,58))/n2

enddo

!!!!!!!!!!!!!!!!!!!!!

!偏早年和偏晚年的差值

do j=1,y

do i=1,x

dif(i,j)=earlyave(i,j)-lateave(i,j)

enddo

!!!!!!!!!!!!!!!!!!!!!

!偏早年和偏晚年的方值

do j=1,y

do i=1,x

earlyvar(i,j)=((dat(i,j,9)-earlyave(i,j))**2+(dat(i,j,14)-earlyave(i,j))**2+(dat(i,j,24)-earlyave(i,j))**2+(dat(i,j,27)-earlyave(i,j))**2+(dat(i,j,33)-earlyave(i,j))**2+(dat(i,j,35)-earlyave(i,j))**2+(dat(i,j,37)-earlyave(i,j))**2+(dat(i,j,42)-earlyave(i,j))**2+(dat(i,j,44)-earlyave(i,j))**2+(dat(i,j,52)-earlyave(i,j))**2)/n1

latevar(i,j)=((dat(i,j,11)-lateave(i,j))**2+(dat(i,j,12)-lateave(i,j))**2+(dat(i,j,17)-lateave(i,j))**2+(dat(i,j,18)-lateave(i,j))**2+(dat(i,j,19)-lateave(i,j))**2+(dat(i,j,20)-lateave(i,j))**2+(dat(i,j,22)-lateave(i,j))**2+(dat(i,j,45)-lateave(i,j))**2+(dat(i,j,51)-lateave(i,j))**2+(dat(i,j,58)-lateave(i,j))**2)/n2

enddo

!!!!!!!!!!!!!!!!!!!!

!计算t统计量

do j=1,y

do i=1,x

ts(i,j)=(earlyave(i,j)-lateave(i,j))/(sqrt(((n1-1)*earlyvar(i,j)+(n2-1)*latevar(i,j))/(n1+n2-2))*sqrt(1/n1+1/n2))

enddo

!!!!!!!!!!!!!!!!!!!!

!资料输出

open(200,file='d:\u200junedif.grd',form='binary')

do j=1,y

do i=1,x

write(200) dif(i,j)

enddo

close(200)

open(300,file='d:\u200junet-test.grd',form='binary')

do j=1,y

do i=1,x

write(300) ts(i,j)

enddo

close(300)

!资料输出完毕

stop

end program main

python读取grd数据_一个关于grads显示grd格式数据的问题。相关推荐

Python语言编程学习：numpy中的array格式数据切片与pandas中的dataframe格式数据切片、相互转换
Python语言编程学习:numpy中的array格式数据切片与pandas中的dataframe格式数据切片.相互转换目录 numpy中的array格式数据切片与pandas中的dataframe ...
Excel合并数据查找函数VLOOKUP()一直显示最后一行数据或者一直报错的解决方法
Excel合并数据查找函数VLOOKUP()一直显示最后一行数据或者一直报错的解决方法大家好,我叫亓官劼(qí guān jié ),在CSDN中记录学习的点滴历程,时光荏苒,未来可期,加油~博客地 ...
Burp Post、Get数据包转为上传multipart/form-data格式数据包
方法一: 新建一个网页进行上传,代码代码如下: <html> <head></head> <body> <form method="po ...
python读取第二行_使用Python操作Excel（二）：读取数据表
上一节我们提到,使用openpyxl可以方便的对数据表进行操作,例如:抽象Excel数据并存入数据库将数据库数据导出到Excel 给一个已存在的数据表追加信息我们还介绍了一些Excel的基本术语, ...
python读取pdf文件_深入学习python解析并读取PDF文件内容的方法
这篇文章主要学习了python解析并读取PDF文件内容的方法,包括对学习库的应用,python2.7和python3.6中python解析PDF文件内容库的更新,包括对pdfminer库的详细解释和应 ...
Python数据挖掘与机器学习_通信信用风险评估实战(1)——读数据
案例及数据说明本选题用到的数据包括公共数据集.用户通话数据.违约行为数据. 训练集数据:用户通话数据(201703-201706),违约行为数据(201707),公共数据-基本信息(201703-2 ...
python 读取excel图片_如何用Python读取Excel中图片？
公众号: 早起Python 作者:刘早起大家好,在使用Python进行办公自动化操作时,一定少不了与Excel表格的交互,我们通常是用pandas处理表格数据,但大多数情况下,都是读取表格中的数值进 ...
怎么用python读取excel图_如何用Python读取Excel中图片？
公众号: 早起Python 作者:刘早起大家好,在使用Python进行办公自动化操作时,一定少不了与Excel表格的交互,我们通常是用pandas处理表格数据,但大多数情况下,都是读取表格中的数值进 ...
python发明小故事_一个与python有关的故事
背景近期系统需求上来迟缓,也腾出了一些时间来梳理沉淀.对系统的代码进行审视,sonar静态检测,消除检测问题.也是还债的过程.回顾过往的坑,有些东西还是要沉淀下来.接下来节听听故事的来龙去脉. 一个 ...
python入侵数据库数据库_一个简单的Python访问Mysql数据库例子
2020/11/3操作记录搭建好Python的数据环境之后,接下来就是在Python代码中访问数据库我先在Navicat图形化界面创建一个数据库命名为pythontest,再在数据库中创建了一个表s ...

python读取grd数据_一个关于grads显示grd格式数据的问题。

python读取grd数据_一个关于grads显示grd格式数据的问题。相关推荐

最新文章

热门文章