NCL读写.nc文件

批量解压,gz文件

(base) yfeng@������������:~$ cd /mnt/d/yff/Dataset/meteorological_forcing_TPDC/Temp/
(base) yfeng@������������:/mnt/d/yff/Dataset/meteorological_forcing_TPDC/Temp$ gunzip *.gz

查看.nc文件中存储的数据信息

 $ ncdump -h prec_ITPCAS-CMFD_V0106_B-01_03hr_010deg_200301.ncVariable: prec (subsection)
Type: short
Total Size: 2 bytes1 values
Number of Dimensions: 1
Dimensions and sizes:   [1]
Coordinates:
Number Of Attributes: 9lon : 70.15lat : 15.15time :        902880long_name :   Precipitation rateunits :       mm hr-1missing_value :       -32767_FillValue :  -32767add_offset :  50scale_factor :        0.0025
(0)     -32767

进入ncl状态

$ ncl
ncl 0> f = addfile("/文件路径/.nc","r")
ncl 1> AREA = f->area
ncl 2> printVarSummary(area)

输出：
Variable: AREA
Type: double
Total Size: 8 bytes
1 values
Number of Dimensions: 2
Dimensions and sizes: [nj | 1] x [ni | 1]
Coordinates:
Number Of Attributes: 3
long_name : area of grid cell in radians squared
coordinate : xc yc
units : radians2

print(AREA(0,0))

打印结果：
Variable: AREA (subsection)
Type: double
Total Size: 8 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 3
units : radians2
coordinate : xc yc
long_name : area of grid cell in radians squared
(0) 5.120350289773058e-05

cd $NCARG_ROOT/lib/ncarg/nclscripts/csm/
ll
total 2496
drwxr-xr-x 2 yfeng yfeng   4096 Dec  4 18:59 ./
drwxr-xr-x 6 yfeng yfeng   4096 Dec  4 18:59 ../
-rw-r--r-- 1 yfeng yfeng  35802 Jan  3  2022 bootstrap.ncl
-rw-r--r-- 1 yfeng yfeng  14623 Jan  3  2022 ccm_func.ncl
-rw-r--r-- 1 yfeng yfeng 729591 Jan  3  2022 contributed.ncl
-rw-r--r-- 1 yfeng yfeng  89921 Jan  3  2022 crop.ncl
-rw-r--r-- 1 yfeng yfeng 137924 Jan  3  2022 diagnostics_cam.ncl
-rw-r--r-- 1 yfeng yfeng  27956 Jan  3  2022 extval.ncl
-rw-r--r-- 1 yfeng yfeng 504377 Jan  3  2022 gsn_code.ncl
-rw-r--r-- 1 yfeng yfeng 555558 Jan  3  2022 gsn_csm.ncl
-rw-r--r-- 1 yfeng yfeng  29815 Jan  3  2022 heat_stress.ncl
-rw-r--r-- 1 yfeng yfeng  58262 Jan  3  2022 popRemap.ncl
-rw-r--r-- 1 yfeng yfeng 212631 Jan  3  2022 shea_util.ncl
-rw-r--r-- 1 yfeng yfeng  87018 Jan  3  2022 skewt_func.ncl
-rw-r--r-- 1 yfeng yfeng  39989 Jan  3  2022 wind_rose.ncl

指定经纬度
print(AREA({},{}))

画图

ncl 9> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
ncl 10> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

exercise 5 计算气象数据平均值例如某月prec

;计算某月气温均值并绘图
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;open multiple files
fns = systemfunc("ls /mnt/d/yff/Dataset/meteorological_forcing_TPDC/temp_CMFD_V0106_B-01_03hr_010deg_????01.nc")
print(fns)infile = addfiles(fns, "r")
ListSetType(infile, "cat");read dataT = short2flt(infile[:]->temp(:,:,:))printVarSummary(T);calculate long-term climatologyT_bar = dim_avg_n_Wrap(T, 0) - 273.15 ;K to °CprintVarSummary(T_bar);plotwks = gsn_open_wks("png", "T_bar") ;set workspaceres = True ;plot mode desiredres@cnFillOn            = True          ; turn on colorres@cnFillPalette       = "gui_default" ; set color mapres@vpWidthF            = 0.5           ; change aspect ratio of plotres@vpHeightF           = 0.3res@gsnMaximize         = True          ; maximize plot sizeres@gsnLeftString = "" ;clear up the long_name informationres@gsnCenterString = "temp Jun.climatology mean" ;setup plot titleplot = gsn_csm_contour(wks, T_bar, res) ;setup plot

NCL读写.nc文件相关推荐

Python 使用netCDF4读写nc文件以及截取指定经纬度范围内的数据生成新的nc文件
Python 使用netCDF4读写nc文件以及截取nc文件经纬度范围内的数据简单介绍nc文件的读写操作,以及实现输入nc文件和坐标范围,输出一个新的nc文件的功能环境 python3.8.13 ...
NCL学习笔记（1）——NCL读取nc文件一个变量存入一个新的nc文件
代码: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/nc ...
学习笔记_ncl_读取nc文件中的变量_制作nc文件的方法
由于生成的文件较大较多,在超算上不方便使用matlab,只好向ncl势力低头目的是在超算上对相关结果进行第一步简单处理后,下载到pc上再进行其他处理 ncl新手如果是从已有的nc文件中读取变量到新 ...
用NCL将GRIB/GRIB2文件转成nc文件（批量转），JRA-55再分析为例
背景: 答主最近需要用到JRA-55再分析数据,但在面对GRIB数据格式时犯难了,便想把GRIB格式转成nc文件再做分析.因为nc文件很容易用python或Matlab处理尝试了好些方法,发现先安装 ...
记录一次ncl从前端到后端出图过程：调用linux-shell执行ncl命令从nc文件出图
ncl出图大概长这样子: 数据文件: ncl脚本: ;************************************************ ; ; These files are load ...
NCL将Grib文件转换成nc文件 JRA55文件说明
NCL将GRIB文件转换为NC文件的函数具体说明 JRA下载手册下载好的JRA monthly 数据,如果是average数据就没有_var,如果是variance数据中有_var,在这里需要注意一 ...
ncl批量处理多个nc文件_利用MATLAB读取NC文件并绘图
最近需要对很多NC 文件进行绘图,通常情况下,我一直用GMT,不过不知道为什么,笔记本上的GMT 用不了了,可能是版本的问题,但是我下了很多版本都无法使用,画出来的图都有问题.由于笔记本安装硬盘也就2 ...
ncl批量处理多个nc文件_UG定制后处理入门必备基础
在UG中生成加工刀轨文件后,由于刀轨文件并不能直接驱动机床,所以必须针对机床来进行后处理.后处理的主要是把CAM软件生成的加工刀路转换成特定机床能识别的G代码,这样机床才能根据对应的指令来完成工件的加 ...
python读取nc文件转成img_使用python的netCDF4库读取.nc文件和创建.nc文件[转]
使用python netCDF4库读取.nc文件和创建.nc文件 1. 介绍 .nc(network Common Data Format)文件是气象上常用的数据格式,python上读取.nc使用 ...

NCL读写.nc文件

NCL读写.nc文件

NCL读写.nc文件相关推荐

最新文章

热门文章