10分钟入门NCL

风子

松下风 发表于 2016-3-22 14:59
强哥出品，比属精品

多谢支持！！
喜欢就点个赞、关注一下

风子

王磊 发表于 2016-3-22 14:59
强哥最后整一个pdf的版本就更好了。

恩，最后肯定整一份发给你

Alliris

新学NCL的福音，跟着楼主学

66kqkq

太详细了，赞！

王磊

强哥最后整一个pdf的版本就更好了。

松下风

强哥出品，比属精品

hnqx

好帖子，楼主有心了

xuebiz

风子 发表于 2016-3-22 10:08
哇，你这学习态度，我给一百分！

谢谢老师，板书有个别地方写错了，我标红了，你要改过来哦~~
还有，拜托你啦，早日写完哦，人家还要好好学习呢。。。

xuebiz

#第三部分

ncl 24> arr2 = (/(/1, 2, 3, 4/), (/5, 6, 7, 8/)/)
ncl 25> arr3 = transpose(arr2)
ncl 26> print(arr2)


Variable: arr2
Type: integer
Total Size: 32 bytes
            8 values
Number of Dimensions: 2
Dimensions and sizes:        [dim0 | 2] x [dim1 | 4]
Coordinates:
(0,0)        1
(0,1)        2
(0,2)        3
(0,3)        4
(1,0)        5
(1,1)        6
(1,2)        7
(1,3)        8
ncl 27> print(arr3)


Variable: arr3
Type: integer
Total Size: 32 bytes
            8 values
Number of Dimensions: 2
Dimensions and sizes:        [dim1 | 4] x [dim0 | 2]
Coordinates:
(0,0)        1
(0,1)        5
(1,0)        2
(1,1)        6
(2,0)        3
(2,1)        7
(3,0)        4
(3,1)        8
ncl 28> arr1 := array_append_record(arr1, (/11, 12/), 0)
ncl 29> print(arr1)


Variable: arr1
Type: integer
Total Size: 48 bytes
            12 values
Number of Dimensions: 1
Dimensions and sizes:        [dim0 | 12]
Coordinates:
(0)        1
(1)        2
(2)        3
(3)        4
(4)        5
(5)        6
(6)        7
(7)        8
(8)        9
(9)        10
(10)        11
(11)        12
ncl 30> x1 = (/(/-3.71, -3.70, -3.03/),\
ncl 30>        (/-1.72, -3.70, -3.64/),\
ncl 30>        (/-1.91, -4.92, -4.85/),\
ncl 30>        (/-4.17, -4.68, -4.41/)/)
ncl 31> x2 = (/(/14.31, 15.95, 19.75/),\
ncl 31>        (/14.00, 10.11, 12.65/)/)
ncl 32> a = table_attach_row (x1, x2, 0)
fatal:Undefined identifier: (table_attach_row) is undefined, can't continue
fatal:["Execute.c":8565]:Execute: Error occurred at or near line 32

#3.2将第二个数组与第一个数组按行连结，这里出错了，应该是
a = table_attach_rows (x1, x2, 0)
------------row后面少了个“s”

后面有输出的，俺就pass了

#3.2将包含12个元素的一维数组arr1扩展到多维（4行12列），这里出错了，应该是
    write_matrix (arr4, "12I3", False)
------------应该是arr4不是arr2哦~~


#3.5

• arr6 = where(arr.lt.0, arr+256, arr)  ; 将小于0的值加上256
  
• arr6_inv = 1. / where(arr6.ne.0, arr6, arr6@_FillValue)  ; 禁用0除
  
• 
  
• ; 将数组某块方形区域以外置为缺测
  
• hgt_mask = where((hgt@lat.ge.20 .and. hgt@lat.le.25  .and. \
  
•                   hgt@lon.ge.110 .and. hgt@lon.le.120), hgt, hgt@_FillValue)
  

其中第一行应为：arr6 = where(arr6.lt.0, arr6+256, arr6) ，即后面的3个arr改为arr6
第二行和第六行，arr6@_FillValue和hgt@_FillValue没有定义

风子

xuebiz 发表于 2016-3-22 09:42
#第一部分

$ ncl

哇，你这学习态度，我给一百分！