- 积分
- 1901
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-6
- 最后登录
- 1970-1-1
|
发表于 2016-3-22 15:32:14
|
显示全部楼层
本帖最后由 xuebiz 于 2016-3-23 10:56 编辑
#第五部分
ncl 0> fg = addfile("hgt_2009_2012_monthly.nc", "r")
ncl 1> hgt = short2flt(fg->hgt)
##------------ 因为不是从头做的,需要用到hgt,就从前面抄了这两句。
##------------ f 后面用到了,会出错,所以改用 fg
ncl 2>
ncl 3> load "cd_inv_string.ncl"
ncl 4> begin
ncl 5> f = asciiread("megi.txt", -1, "string")
##------------ 原文: f = asciiread("megi.1.txt", -1, "string") ; 已字符串格式读取文件所有行
##------------ 作者链接下载的数据文件名字是"megi.txt",这里数据名写的是"megi.1.txt"
ncl 6> time_str = str_get_cols(f, 0, 13)
ncl 7> time = cd_inv_string(time_str, "%Y-%N-%D %H")
ncl 8>
ncl 9> lat = tofloat(str_get_field(f, 4, " "))*0.1
ncl 10> lon = tofloat(str_get_field(f, 5, " "))*0.1
ncl 11>
ncl 12> print(time+" "+lon+" "+lat)
ncl 13> printVarSummary(time)
ncl 14> end
(0) 1847664 141.4 11.8
(1) 1847670 141.1 11.9
(2) 1847676 140.9 12
(3) 1847682 140.6 12.1
(4) 1847688 140 12.2
(5) 1847694 139.2 12.7
(6) 1847700 138.5 13.1
(7) 1847706 137.8 13.6
(8) 1847712 137.2 14.3
(9) 1847718 136.6 14.9
(10) 1847724 135.6 15.7
(11) 1847730 134.3 16.6
(12) 1847736 133 17.5
(13) 1847742 131.6 18.1
(14) 1847748 130.2 18.4
(15) 1847754 128.8 18.7
(16) 1847760 127.5 18.7
(17) 1847766 126.2 18.5
(18) 1847772 125.1 18.1
(19) 1847778 124.2 17.6
(20) 1847784 123.3 17.5
(21) 1847790 121.8 17.2
(22) 1847796 120.6 16.8
(23) 1847802 119.5 16.6
(24) 1847808 118.9 16.5
(25) 1847814 118.4 16.6
(26) 1847820 117.9 16.7
(27) 1847826 117.6 16.9
(28) 1847832 117.3 17.2
(29) 1847838 117.2 17.8
(30) 1847844 117.2 18.4
(31) 1847850 117.4 18.8
(32) 1847856 117.5 19.4
(33) 1847862 117.6 19.9
(34) 1847868 117.8 20.2
(35) 1847874 118 20.6
(36) 1847880 118.2 21
(37) 1847886 118.3 21.6
(38) 1847892 118.3 22.3
(39) 1847898 118.2 22.8
(40) 1847904 118 23.4
(41) 1847910 117.7 24
(42) 1847916 117.8 24.4
(43) 1847922 118 24.8
(44) 1847928 118.4 25.3
(45) 1847934 118.6 26.2
Variable: time
Type: double
Total Size: 368 bytes
46 values
Number of Dimensions: 1
Dimensions and sizes: [46]
Coordinates:
Number Of Attributes: 3
_FillValue : 9.969209968386869e+36
calendar : standard
units : hours since 1800-01-01 00:00:00
ncl 15> load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
ncl 16> print(cd_string(hgt&time(:10), "%Y-%N-%D %H:%M:%S"))
Variable: unnamed (return)
Type: string
Total Size: 88 bytes
11 values
Number of Dimensions: 1
Dimensions and sizes: [11]
Coordinates:
Number Of Attributes: 1
_FillValue : missing
(0) 2009-01-01 00:00:00
(1) 2009-02-01 00:00:00
(2) 2009-03-01 00:00:00
(3) 2009-04-01 00:00:00
(4) 2009-05-01 00:00:00
(5) 2009-06-01 00:00:00
(6) 2009-07-01 00:00:00
(7) 2009-08-01 00:00:00
(8) 2009-09-01 00:00:00
(9) 2009-10-01 00:00:00
(10) 2009-11-01 00:00:00
ncl 17> ndims = dimsizes(hgt)
ncl 18> print(cd_string(hgt&time(ndims(0)-10:), "%c %Y"))
Variable: unnamed (return)
Type: string
Total Size: 80 bytes
10 values
Number of Dimensions: 1
Dimensions and sizes: [10]
Coordinates:
Number Of Attributes: 1
_FillValue : missing
(0) Mar 2012
(1) Apr 2012
(2) May 2012
(3) Jun 2012
(4) Jul 2012
(5) Aug 2012
(6) Sep 2012
(7) Oct 2012
(8) Nov 2012
(9) Dec 2012
ncl 19>
|
|