绘制分布图（自由基、isop、O3）

阿水水水水水

于 2022-08-23 11:30:08 发布

阅读量381

点赞数

分类专栏：脚本积累文章标签： python 算法人工智能

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

版权

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"

wrffile=addfile("/data1/loucx/COPY/WRF/WRF/run/wrfout_d03_2020-08-01_00:00:00","r")
lat2d=wrffile->XLAT(0,3:176,3:170)
lon2d=wrffile->XLONG(0,3:176,3:170)
dims = dimsizes(lat2d)
geofile=addfile("/data1/loucx/COPY/CNAQ/data1/cctm/ZJ_CMAQ_appmap.nc","r")
flag=geofile->ZJ_1(0,0,:,:)
DataDir1 = "/data1/loucx/COPY/CNAQ/data/cctm/"
flist1 = systemfunc("ls "+DataDir1+"ACONC.*_d03")
Sfile1 = addfiles(flist1, "r")
ListSetType (Sfile1, "join")
isop=Sfile1[:]->OH ;HO2 ;[ncl_join | 31] x [TSTEP | 25] x [LAY | 15] x [ROW | 174] x [COL | 168]
Isop=dim_avg_n(isop(:,:,0,:,:),0)*1000  ;ppm->ppb Dimensions and sizes:[24] x [183] x [183]
ISop=dim_avg_n(Isop(:,:,:),0)*100000;*10^3;100000;*10^5

var_0=ISop(:,:)
num1=0
sumflag=0.0
w=new((/29232/),float)
do z=0,173
 do f=0,167
  if(flag(z,f) .ne. 0) then
  w(num1)=ISop