10分钟入门NCL

xuebiz

#第二部分

ncl 6> arr1 = ispan(1,10,1)
ncl 7> print(arr1(4))


Variable: arr1 (subsection)
Type: integer
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes:        [1]
Coordinates:
(0)        5
ncl 8> print(hgt(0,4,1,10))


Variable: hgt (subsection)
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes:        [1]
Coordinates:
Number Of Attributes: 24
  lon :        25
  lat :        87.5
  level :        600
  time :        1832064
  missing_value :        -32767
  valid_range :        ( -1500, 35799 )
  missing_value_original :        32766
  _FillValue :        -32767
  _FillValue_original :        -32767
  long_name :        Monthly Geopotential Heights on Pressure Levels
  unpacked_valid_range :        ( -1500, 35800 )
  actual_range :        ( -209, 32301 )
  units :        m
  precision :        0
  least_significant_digit :        0
  GRIB_id :        7
  GRIB_name :        HGT
  var_desc :        Geopotential height
  dataset :        NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Monthly Averages
  level_desc :        Pressure Levels
  statistic :        Individual Obs
  parent_stat :        Other
  standard_name :        geopotential_height
  cell_methods :        time: mean (montly from 6-hourly values)
(0)        3751
ncl 9> print(arr1(2:4))


Variable: arr1 (subsection)
Type: integer
Total Size: 12 bytes
            3 values
Number of Dimensions: 1
Dimensions and sizes:        [3]
Coordinates:
(0)        3
(1)        4
(2)        5
ncl 10> print(arr1(1::2))


Variable: arr1 (subsection)
Type: integer
Total Size: 20 bytes
            5 values
Number of Dimensions: 1
Dimensions and sizes:        [5]
Coordinates:
(0)        2
(1)        4
(2)        6
(3)        8
(4)        10
ncl 11> print(hgt(:15:2,0,15,29))


Variable: hgt (subsection)
Type: float
Total Size: 32 bytes
            8 values
Number of Dimensions: 1
Dimensions and sizes:        [time | 8]
Coordinates:
            time: [1832064..1842240]
Number Of Attributes: 23
  lon :        72.5
  lat :        52.5
  level :        1000
  missing_value :        -32767
  valid_range :        ( -1500, 35799 )
  missing_value_original :        32766
  _FillValue :        -32767
  _FillValue_original :        -32767
  long_name :        Monthly Geopotential Heights on Pressure Levels
  unpacked_valid_range :        ( -1500, 35800 )
  actual_range :        ( -209, 32301 )
  units :        m
  precision :        0
  least_significant_digit :        0
  GRIB_id :        7
  GRIB_name :        HGT
  var_desc :        Geopotential height
  dataset :        NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Monthly Averages
  level_desc :        Pressure Levels
  statistic :        Individual Obs
  parent_stat :        Other
  standard_name :        geopotential_height
  cell_methods :        time: mean (montly from 6-hourly values)
(0)        220
(1)        189
(2)        120
(3)        46
(4)        151
(5)        228
(6)        253
(7)        177
-----------太激动，中间手误敲错2次，强迫症患者们，行号不连续请无视啊。。。
ncl 14> hgt_500 = hgt(0, {500}, :, :)
ncl 15> hgt_500_0 = hgt(time| 0, {level| 500:700}, {lon| 90:120}, lat| :)
ncl 16> idx = (/1,3,4,5,7,8/)
ncl 17> print(arr1(idx))


Variable: arr1 (subsection)
Type: integer
Total Size: 24 bytes
            6 values
Number of Dimensions: 1
Dimensions and sizes:        [6]
Coordinates:
(0)        2
(1)        4
(2)        5
(3)        6
(4)        8
(5)        9
ncl 18> a = (/1, 2, 3, 4, 5, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5/)
ncl 19> a@_FillValue = 5
ncl 20> valid = ind(ismissing(a))
ncl 21> a(valid) = 0
ncl 22> a(ind(a.ne.0)) = -a(ind(a.ne.0)) * 2 + 1
ncl 23> print(a)


Variable: a
Type: integer
Total Size: 60 bytes
            15 values
Number of Dimensions: 1
Dimensions and sizes:        [15]
Coordinates:
Number Of Attributes: 1
  _FillValue :        5
(0)        -1
(1)        -3
(2)        -5
(3)        -7
(4)        0
(5)        0
(6)        -7
(7)        -5
(8)        -3
(9)        -1
(10)        -1
(11)        -3
(12)        -5
(13)        -7
(14)        0

xuebiz

#第一部分

$ ncl
Copyright (C) 1995-2014 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.2.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
ncl 0> f = addfile("hgt_2009_2012_monthly.nc", "r")
ncl 1> print(f)

Variable: f
Type: file
filename:        hgt_2009_2012_monthly
path:        hgt_2009_2012_monthly.nc
   file global attributes:
   dimensions:
      time = 48
      level = 12
      lat = 37
      lon = 144
   variables:
      short hgt ( time, level, lat, lon )
         cell_methods :        time: mean (montly from 6-hourly values)
         standard_name :        geopotential_height
         parent_stat :        Other
         statistic :        Individual Obs
         level_desc :        Pressure Levels
         dataset :        NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Monthly Avera
ges
         var_desc :        Geopotential height
         GRIB_name :        HGT
         GRIB_id :        7
         least_significant_digit :        0
         precision :        0
         _FillValue :        -32767
         missing_value :        32766
         scale_factor :         1
         add_offset :        31265
         units :        m
         actual_range :        ( -209, 32301 )
         unpacked_valid_range :        ( -1500, 35800 )
         valid_range :        ( -32765, 4534 )
         long_name :        Monthly Geopotential Heights on Pressure Levels

      double time ( time )
         bounds :        time_bnds
         coordinate_defines :        start
         axis :        T
         standard_name :        time
         prev_avg_period :        0000-00-01 00:00:00
         avg_period :        0000-01-00 00:00:00
         delta_t :        0000-01-00 00:00:00
         actual_range :        ( 1569072, 1866384 )
         long_name :        Time
         units :        hours since 1800-1-1 00:00:00

      float level ( level )
         coordinate_defines :        point
         axis :        Z
         GRIB_name :        hPa
         GRIB_id :        100
         positive :        down
         long_name :        Level
         actual_range :        ( 1000, 10 )
         units :        millibar

      float lat ( lat )
         coordinate_defines :        point
         axis :        Y
         standard_name :        latitude
         long_name :        Latitude
         actual_range :        ( 90,  0 )
         units :        degrees_north

      float lon ( lon )
         coordinate_defines :        point
         axis :        X
         standard_name :        longitude
         actual_range :        (  0, 357.5 )
         long_name :        Longitude
         units :        degrees_east

ncl 2> list_filevars(f)

float        lon [ lon | 144 ]
        coordinate_defines
        axis
        standard_name
        actual_range
        long_name
        units

float        lat [ lat | 37 ]
        coordinate_defines
        axis
        standard_name
        long_name
        actual_range
        units

float        level [ level | 12 ]
        coordinate_defines
        axis
        GRIB_name
        GRIB_id
        positive
        long_name
        actual_range
        units

double        time [ time | 48 ]
        bounds
        coordinate_defines
        axis
        standard_name
        prev_avg_period
        avg_period
        delta_t
        actual_range
        long_name
        units

short        hgt [ time | 48 ] x [ level | 12 ] x [ lat | 37 ] x [ lon | 144 ]
        cell_methods
        standard_name
        parent_stat
        statistic
        level_desc
        dataset
        var_desc
        GRIB_name
        GRIB_id
        least_significant_digit
        precision
        _FillValue
        missing_value
        scale_factor
        add_offset
        units
        actual_range
        unpacked_valid_range
        valid_range
        long_name
ncl 3> hgt = short2flt(f->hgt)
ncl 4> printvarsummary(hgt)
fatal:syntax error: line 4 before or near \n
printvarsummary(hgt)
--------------------^

fatal:syntax error: possibly an undefined procedure
------------这里没有区分大小写，出错了！！！
ncl 5> printVarSummary(hgt)

Variable: hgt
Type: float
Total Size: 12275712 bytes
            3068928 values
Number of Dimensions: 4
Dimensions and sizes:        [time | 48] x [level | 12] x [lat | 37] x [lon | 144]
Coordinates:
            time: [1832064..1866384]
            level: [1000..100]
            lat: [90.. 0]
            lon: [ 0..357.5]
Number Of Attributes: 20
  cell_methods :        time: mean (montly from 6-hourly values)
  standard_name :        geopotential_height
  parent_stat :        Other
  statistic :        Individual Obs
  level_desc :        Pressure Levels
  dataset :        NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Monthly Averages
  var_desc :        Geopotential height
  GRIB_name :        HGT
  GRIB_id :        7
  least_significant_digit :        0
  precision :        0
  units :        m
  actual_range :        ( -209, 32301 )
  unpacked_valid_range :        ( -1500, 35800 )
  long_name :        Monthly Geopotential Heights on Pressure Levels
  _FillValue_original :        -32767
  _FillValue :        -32767
  missing_value_original :        32766
  valid_range :        ( -1500, 35799 )
  missing_value :        -32767

xuebiz

楼主辛苦啊，好好写啊，正在学，初学，

低调高姿态

好棒好棒。。。。。。

风子

网络小丁丁 发表于 2016-3-21 09:11
有没有PDF或者WORD档收藏

有一种东西，叫做 印象笔记·剪藏
跟着本文过一遍才是王道，收藏的东西很少有人再看

喜欢本文，请关注+收藏+分享+赞+评分

我会不断完善，毕竟要写出一个好东西不是一朝一夕的事

网络小丁丁

有没有PDF或者WORD档收藏

婧观气变

风子

李舒文 发表于 2016-3-20 22:47
楼主还没有将整个教程完成。。果断时间才行

时间太宝贵，帖子慢慢写

李舒文

Hvoice 发表于 2016-3-20 10:18
舒文快给我发一份啦

楼主还没有将整个教程完成。。果断时间才行

Hvoice

风子 发表于 2016-3-20 10:22
现在不需要密码了，直接看

好的呢~谢谢楼主啦~