利用 ncl 自带海陆数据分别绘制海陆图

qdzwxz

如果你有一个全球分布的数据，比如地表温度 T（：，：），想分别画出海洋和陆地的等值线图，用 ncl 该怎么做呢？

ncl 自带了一个1X1 度的海陆数据，可以实现以上功能。总共分三步：

1，需要在文件开头加上：

     load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

2， 打开自带海陆数据，并按照自己的地表温度数据产生用来区分海陆的新数据：

     b = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")

     lsm_T= landsea_mask(b->LSMASK,T&lat,T&lon)

     其中 landsea_mask 是shea_util.ncl中的一个函数。

3，用新数据mask out 海洋或者陆地，得到自己想要的那部分：

   T_land = mask(T,lsm_T.eq.0, False)      ; mask out ocean points
          T_ocean = mask(T,lsm_T.eq.1, False)   ; mask out ocean points

          copy_VarMeta(T,T_land)        ；建议把这两句加上，这样新的land 和 ocean都有
          copy_VarMeta(T,T_ocean)     ；了和原始数据 T一样的meta信息，并且海洋(或陆地)设为了却测值-999

      大功告成。现在可以用T_land 和T_ocean去画自己想要的图了。

artmunich

用NCL画的图好简洁啊，美观~楼主能否给一下软件，我是win7，在这边没注册成功

mofangbao

NCL版区有希望啦！

qdzwxz

@artmunich没问题啊~ 能否给个邮箱？这个压缩包有50m左右，传不上来。
我是在Linux下用的，Windows下安装貌似有点麻烦，你可以参考附件上的（不知适合Win7不，成功的话别忘分享下经验哦）

artmunich

qdzwxz 发表于 2011-9-16 19:46

                               
登录/注册后可看大图

@artmunich没问题啊~ 能否给个邮箱？这个压缩包有50m左右，传不上来。
我是在Linux下用的，Windows下安装貌 ...

好的，我邮箱 artmunich@qq.com， 麻烦你了，试着安装一下~有问题再共同探讨~

artmunich

mofangbao 发表于 2011-9-16 19:03

                               
登录/注册后可看大图

NCL版区有希望啦！

代表管理员清风问下这位兄弟@qdzwxz愿意做下NCL版的版主吗？
不介意吧清风？

qdzwxz

artmunich 发表于 2011-9-16 20:08

                               
登录/注册后可看大图

好的，我邮箱 ， 麻烦你了，试着安装一下~有问题再共同探讨~

已发送，注意查收哦

mofangbao

artmunich 发表于 2011-9-16 20:11

                               
登录/注册后可看大图

代表管理员清风问下这位兄弟@qdzwxz愿意做下NCL版的版主吗？
不介意吧清风？

嘿嘿 偶表示木有问题
@qdzwxz自己呢？
可以的话联系一下我吧，偶滴Q:937062711

qdzwxz

artmunich 发表于 2011-9-16 20:11

                               
登录/注册后可看大图

代表管理员清风问下这位兄弟@qdzwxz愿意做下NCL版的版主吗？
不介意吧清风？

版主？？？ 我怕是不太会打理事物，凌乱了本版责任可是重大啊~~ 谢谢厚爱了

qdzwxz

mofangbao 发表于 2011-9-16 20:48

                               
登录/注册后可看大图

嘿嘿 偶表示木有问题
@qdzwxz自己呢？
可以的话联系一下我吧，偶滴Q:937062711

真对不起啦，我没啥经验，不堪此重任啊~
还是写点经验分享给大家比较在行，呵呵。谢谢厚爱