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

最新推荐文章于 2025-04-29 09:22:12 发布

阿水水水水水

最新推荐文章于 2025-04-29 09:22:12 发布

阅读量311

点赞数 1

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

本文链接：https://blog.youkuaiyun.com/qq_40256777/article/details/122244891

版权

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

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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[:]-