爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5424|回复: 1

求助:ncl脚本出错

[复制链接]

新浪微博达人勋

发表于 2020-7-4 15:11:00 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-7-5 12:59:19 | 显示全部楼层
检查两个相乘的数组大小是否一样
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表