更新提取数据（单位统一）

阿水水水水水

于 2021-12-30 20:34:56 发布

阅读量277

点赞数 1

分类专栏：脚本积累文章标签：脚本语言

本文链接：https://blog.csdn.net/qq_40256777/article/details/122244891

版权

该博客介绍了如何使用NCL脚本读取CCTM CONC数据，并通过WRF文件获取经纬度信息，将大气污染物O3和ISOP浓度转换为站点数据。通过`rcm2points`函数将格点数据转换到特定站点，最后生成O3和ISOP站点对比CSV文件。

摘要由CSDN通过智能技术生成

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
DataDir1 = "/data1/loucx/COPY/CNAQ/data/cctm/"
flist1 = systemfunc("ls "+DataDir1+"CONC.*_d03")
Sfile1 = addfiles(flist1, "r")
ListSetType (Sfile1, "join")
isop=Sfile1[:]->ISOP
o3=Sfile1[:]-