- 积分
- 1339
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-5-29
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
想绘制污染物分布图,脚本中某个污染物的数值需要两个变量数值相乘,但是出错了。想请教一下前辈们怎么修改?
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
;grid file
grid_fn="/mnt/disk2/pxj/CMAQ/data/PRD12/met/mcip/GRIDCRO2D.nc"
;Fig title
fig_title="SO2 concentration (ppmv)"
g_font=0.02
;File IO
pdata_fn = "/mnt/disk2/pxj/CMAQ/data/output_CCTM_v52_gcc_PRD12/CCTM_CONC_v52_gcc_PRD12_20200324.nc"
ctrl_in = addfile(pdata_fn ,"r") ; open output netCDF file
var1 = ctrl_in->ASO4I(:,0,:,:)
var2 = ctrl_in->ANO3I(:,0,:,:)
var3 = ctrl_in->ANH4I(:,0,:,:)
var4 = ctrl_in->ANAI(:,0,:,:)
var5 = ctrl_in->ACLI(:,0,:,:)
var6 = ctrl_in->AECI(:,0,:,:)
var7 = ctrl_in->ALVPO1I(:,0,:,:)
var8 = ctrl_in->ASVPO1I(:,0,:,:)
var9 = ctrl_in->ASVPO2I(:,0,:,:)
var10 =var7/1.39+var8/1.32+var9/1.26
var11 = ctrl_in->ALVOO1I(:,0,:,:)
var12 = ctrl_in->ALVOO2I(:,0,:,:)
var13 = ctrl_in->ASVOO1I(:,0,:,:)
var14 = ctrl_in->ASVOO2I(:,0,:,:)
var15 = ctrl_in->AOTHRI(:,0,:,:)
var16 =var1+var2+var3+var4+var5+var6+var7+var8+var9+var10+var11+var12+var13+var14+var15
pdata_fn = "/mnt/disk2/pxj/CMAQ/data/output_CCTM_v52_gcc_PRD12/CCTM_PMDIAG_v52_gcc_PRD12_20200324.nc"
ctrl_in = addfile(pdata_fn ,"r") ; open output netCDF file
var17 = ctrl_in->PM10AT(:,0,:,:)
var1=var16*var17
latlon_in = addfile(grid_fn,"r")
lat = latlon_in->LAT(0,0,:,:)
lon = latlon_in->LON(0,0,:,:)
var1@lat2d=lat
var1@lon2d=lon
res = True ; use plot options
..........
fatal:Mul: Dimension size, for dimension number 0, of operands does not match, can't continue
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 42 in file cmax2.ncl
|
|