(2011-11-08 04:02)
标签：
ferret
绘图
eof
太平洋
it
|
分类：
学术科研
|
下面使用ferret绘图分析前面eof分析的结果。使用的数据是cdo输出的结果eofs.nc、eigenvalues.nc，以及pc00000.nc-pc00010.nc。
我们使用循环，对每个变量分别绘制一张前4个eof的空间分布图，一张前4个pc的时间序列图，以及一张前10个eof的方差贡献图。
新建一个文件，命名为draweof.jnl。
1. 外层循环控制
let vps={'ul','ur','ll','lr'}
let vars={'tsw','slp','ustrw','vstrw'}
这是定义字符串数组。ferret的字符串要放在双引号中，而字符串数组则用大括号定义。这里定义的两个字符串数组分别存放子图的位置和变量的名字，下面会用到。
repeat/range=1:4/name=nvar \
(......... ;\
......... ;\
)
这是最外层循环，遍历所有4个要分析的变量。ferret的循环用repeat命令完成，其range属性定义循环变量的范围，name属性定义循环变量的名称，后面接一对括号，里面是
采用与前一篇文章相同的数据，这次改用ferret进行EOF分析。
1. 导入要分析的数据：
在ferret中输入命令
use tsw, slp, ustrw, vstrw
察看数据属性：
show data
ferret会显示当前打开的数据文件及其中所包含的数据变量的列表和信息。这些数据文件中只有一个处于当前激活状态，在列表中标记为（default），它所包含的数据变量可以直接访问。非激活状态的数据文件中的变量不能直接使用，而需要以var[d=file]的形式访问。其中var是变量名，file是数据文件名。要改变当前激活的数据文件，使用set
data file或set data n，其中file为数据文件名，n为数据文件的编号。
2. 设定时间和空间范围
set
reg/x=110s:70w/y=20s:60n/t=1-Jan-3500:31-Dec-3600
reg是region的简写。set命令用于设置属性，region是要设置的对象，/后面分隔的是不同的属性。ferret中的经纬度可以用110e,70w,20s,60n这样的写法表示，同时也支持用负数表示西经和南纬。时间的写法必须按照dd-mmm-yyyy的格式来写。这样设置了时空范
标签：
气候
数据处理
mpi
cdo
编译
安装
it
|
分类：
IT编程
|
CDO (Climate Data Operators) 是Max-Planck Institute for
Meteorology开发的一个气候数据处理软件。它提供非常简单明了的命令行操作方式，能够对气候（海洋、大气）领域科研中常用的数据格式（NetCDF、GRIB等）进行处理。支持超过400种运算，能够完成网格选择、转换、变换、比较、统计、回归、相关、插值、EOF分析等许多任务。
CDO以GPL开源协议发布，支持Unix/Linux、MacOS、Windows等多种平台。其项目首页为https://code.zmaw.de/projects/cdo，其中可以下载源码包。要想在CDO中操作NetCDF等数据格式，需要系统中已安装NetCDF等的支持库。如果没有安装，则可以使用开发小组提供的集成库libs4cdo（下载页面中提供），其中包含了CDO需要的各种外部支持库，安装后就可以启用CDO的全部功能。
下载libs4cdo-0.0.9后在终端中解压并进入源码目录：
tar -xzvf
libs4cdo-0.0.9.tar.gz
cd
libs4cdo-0.0.9
编译安装：
sudo make libs4cdo
PREFIX=/usr/local
其中P
加载中…
