10分钟入门NCL

风子

xuebiz 发表于 2016-3-22 15:32
#5.1
台风例子中第三行
    f = asciiread("megi.1.txt", -1, "string") ; 已字符串格式读取文件所有行

已更正，希望你最后学完能给个各方面的评估，如难度、学下去的兴趣！

sunnygirlbo

哈哈，在群里面多次观摩了大神，，学习学习{:5_231:}

xuebiz

#第五部分
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>

lucky_girl

{:eb511:}{:eb511:}{:eb511:}

婧观气变

风子

王磊 发表于 2016-3-22 15:08
原文转自家园@MSN强

看版权声明
禁止转载
哈哈

还有是 MCS ， MCS哦， 不是MSN， 不是MSN， 不是MSN

王磊

原文转自家园@MSN强

南宫海语

支持支持写的不错学习学习

南宫海语

王磊 发表于 2016-3-22 15:02
谢强哥

哟！这不是蕾蕾嘛

王磊

风子 发表于 2016-3-22 15:01
恩，最后肯定整一份发给你

谢强哥{:eb511:}{:eb511:}