GrADS中批量描述的基础形式和基础模板

river

  前些天发了一个nc文件描述文件编写和批量描述nc文件的帖子。那个是针对nc文件的，但是其他的资料做起来也是一样的思路。今天这个就是想简单的介绍一下有关于GrADS批量描述的通用的部分。  需要了解批量描述中常用的替代模板，如下：

时间的替换模板

%x1   1 digit decade 1位数的十年

%x3   3 digit decade 3位数的十年

%y2   2 digit year 2位数的年

%y4   4 digit year 4位数的年

%m1   1 or 2 digit month 1或者2位数的月

%m2   2 digit month (leading zero if needed) 2位数的月（如果需要前面补0）

%mc   3 character month abbreviation 3个英文字母的月的缩写

%d1   1 or 2 digit day 1位或2位数的天

%d2   2 digit day (leading zero if needed) 2位数的月（如果需要前面补0）

%h1   1 or 2 digit hour  1位或2位数的小时

%h2   2 digit hour 2位数的小时%h3   3 digit hour (e.g., 120 or 012) 3位数的小时（例如：120 或者012）

%n2   2 digit minute; leading zero if needed  2位数的分钟（如果需要前面补0）

还有很多，不一一列举了。想看的去官方文档Using Templates或者我的附件 批量描述基本样式.docx (26.25 KB, 下载次数: 68)

2013-5-28 16:39 上传

点击文件名下载附件

再说一下批量描述文件编写的思路：

• 首先编写可以正确描述一个资料的ctl，并能正确出图。
• 修改该ctl里dest后的文件名，使用合适的替代格式替代
• 增加options template，这个是批量描述必须加上的一句
• 计算出所要批量描述的文件的总时次，修改ctl里的总时次
• 画图进行验证，如果不对，再根据具体情况做出相应的修改
• 大功告成，可以用了                                                                                             
  

再强调一下，所有的批量描述文件都是在可以正确描述一个文件的基础上得来的。


下面列举几个GrADS官方的例子：
1.  下面是一组二进制文件共一个月，每一天的每小时数据单独存在一个文件中:
  1may92.dat  2may92.dat
   ...
   31may92.dat
必须在数据描述文件中修改三个记录。请注意 TDEF 项要反映整个月的数据时间长度：      
   DSET ^%d1may92.dat
    OPTIONS template
   TDEF 744 linear 0z1may1992 1hr
2.  如果以上数据集继续扩展其他年份的每个月 ：
   1jun92.dat
   2jun92.dat
   ...
   1jan93.dat
那就需要在DSET项中增加年和月的替代模板，并且相应的延长TDEF 项中的时间长度：
   DSET^%d1%mc%y2.dat
   OPTIONS template
   TDEF 6624 linear 0z1may1992 1hr

  3.假设有一组七个 netcdf 文件，每个包含跨越十年的月度数据：   
   pr.1880_1889.nc
   pr.1890_1899.nc
   pr.1900_1909.nc
   pr.1910_1919.nc
   pr.1920_1929.nc
   pr.1930_1939.nc
   pr.1940_1949.nc
那就需要描述文件中包括下面的几项：   
   DSET ^pr.%x30_%x39.nc
   OPTIONStemplate
   DTYPE netcdf
   TDEF 840 linearjan1880 1mo

4.  有两个netcdf文件, 一个包含50年的月度数据(共600个月，也就是600 个时间步长),另一个包含100年的月度数据（共1200个月，也就是1200 个时间步长):
   pr.1851-1900.nc
   pr.1901-2000.nc
描述文件就需要包括下面几项：
   DSET^pr.%ch.nc
   CHSUB  1  600 1851-1900
   CHSUB  601 1800 1901-2000
   OPTIONS template
   DTYPE netcdf
   TDEF 1800 linearjan1851 1mo
注：其实红色部分很好理解，1 600表示第一个资料里数据的时间起点和结束点，1851-1900就是文件名中的%ch。


如果这两个文件在不同的磁盘路径下，那就需要把那几项用下面的写法代替   
   DSET %ch
   CHSUB 1  600/disk1/pr.1851-1900.nc
   CHSUB 601 1800/disk2/pr.1901-2000.nc

5.  预报模式输出结果的资料如下：  
  MMOUT_DOMAIN1_00
  MMOUT_DOMAIN1_01
  MMOUT_DOMAIN1_02
那 DSET 项就要这样写：
  DSET^MMOUT_DOMAIN1_%tm2

river

mofangbao 发表于 2013-5-28 17:26

                               
登录/注册后可看大图

我曾经也干过这种事儿

最近每天看着要到暴雨了，要到暴雨了的。实在等不及了用这办法赚点积分

天道酬勤

@river的帖子，必属精华呀！

river

天道酬勤 发表于 2013-5-28 16:58

                               
登录/注册后可看大图

@river的帖子，必属精华呀！

这个是从GrADS官网翻过来的，属于搬运啦

mofangbao

river 发表于 2013-5-28 17:16

                               
登录/注册后可看大图

这个是从GrADS官网翻过来的，属于搬运啦

我曾经也干过这种事儿

mofangbao

river 发表于 2013-5-28 19:28

                               
登录/注册后可看大图

最近每天看着要到暴雨了，要到暴雨了的。实在等不及了用这办法赚点积分

看来你成功了！

传说中的谁

你这刚发完，后面就有人提问这样的问题了。。。。

jeffggg

river 发表于 2013-5-28 19:28

                               
登录/注册后可看大图

最近每天看着要到暴雨了，要到暴雨了的。实在等不及了用这办法赚点积分

我顶你！！！！！

善人/jw

必须顶起啊~~

liutaoz21

顶river~成功暴雨~