处理海温陆地格点数据的一些小技巧

言深深

前一段时间一直在忙海温，温度场数据，在处理的过程中，发现noaa下载的海温资料中缺测部分（陆地数据）用-999代替的，比较讨厌，因为要做和温度场的相关（或者一些别的什么需要分别海洋和陆地部分的），这样就使得一一对应起来比较困难了，为此将海洋和陆地写了一个txt文件，陆地部分用“1”代替，海洋部分用“0”代替，这样就可以得到一张71*144的格点分布图啦
这里给出印度洋的txt文档截图，为了让大家看清楚，用黑色线将“1”圈起来了
全部数据放在 istandard.txt (40.08 KB, 下载次数: 77)

2011-10-25 10:07 上传

点击文件名下载附件 standard.txt (59.91 KB, 下载次数: 59)

2011-10-25 10:07 上传

点击文件名下载附件 注意istand中存放的格式是（71,144），即144列（经度），71行（维度）。stand中将所有数据都放在一列了。可以用下面的fortran小程序调用
integer stand(71*144),istand(71,144)
call get_istand(istand)
end
subroutine get_stand(stand)
integer stand(71*144)
        open(1,file="standard.txt")
        do i=1,71*144
        read (1,*)stand(i)
        enddo
        close(1)
end

subroutine get_istand(istand)
integer istand(71,144)
        open(1,file="istandard.txt")
        do i=1,71
        read (1,*)(istand(i,j),j=1,144)
        enddo
        close(1)
end
主程序中得到stand（71*144）（istand（71,144））数组后和需要处理的数据sst（71,144）/air（71,144）比较，添加一个判别式if（istand（i，j）.eq.1) then处理陆地数据，否则处理海洋数据。我用surfer将istand数据绘制出来如下：

尽头的尽头

支持深深~~~

言深深

尽头的尽头 发表于 2011-10-25 10:20

                               
登录/注册后可看大图

支持深深~~~

你也很给力啊，刚刚发出来就给你抢沙发了

mofangbao

来抢深深的沙发了 这个是要加分加精华滴

言深深

mofangbao 发表于 2011-10-25 10:25

                               
登录/注册后可看大图

来抢深深的沙发了 这个是要加分加精华滴

咣···只剩地板了，哈哈

言深深

mofangbao 发表于 2011-10-25 10:25

                               
登录/注册后可看大图

来抢深深的沙发了 这个是要加分加精华滴

咣···只剩地板了，哈哈

尽头的尽头

言深深 发表于 2011-10-25 10:24

                               
登录/注册后可看大图

你也很给力啊，刚刚发出来就给你抢沙发了

哈哈，抢沙发专业户。。

言深深

尽头的尽头 发表于 2011-10-25 11:00

                               
登录/注册后可看大图

哈哈，抢沙发专业户。。

好孩子啊···

天目神眉

深深给力哇，顺便来一张财神卡吧

天目神眉

木有