10分钟入门NCL

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
现在不需要密码了，直接看

好的呢~谢谢楼主啦~

十月的天空

厉害！学习