绘制wrfchemi_00z_d02文件变量E_PM_10图出错

悠然一梦123

绘制E_PM_10，出现下列错误提示
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 13 in file E_PM_10。
脚本如下：
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"

begin
a = addfile("/WRF/WRFV3/test/em_real/wrfchemi_00z_d02","r")

wks = gsn_open_wks(pdf,"E1")

ter = a->E_PM_10(0,0,:,:)

ter!0="lat"
ter&lat=E_PM_10&lat（错误提示行）
ter!1="lon"
ter&lon=E_PM_10&lon
ter&lat@units = "degrees_north"
ter&lon@units = "degrees_east"

res                   = True                  
res@cnFillOn          = True                                
res@cnLinesOn         = False                  
res@cnLevelSelectionMode = "ManualLevels"      
res@cnMinLevelValF       = 0                  
res@cnMaxLevelValF       =  0.5                  
res@cnLevelSpacingF      =  0.05                 
res@mpDataSetName              = "Earth..4"
res@mpDataBaseVersion          = "MediumRes"
res@mpAreaMaskingOn            = True
res@mpMaskAreaSpecifiers       = (/"China"/)
res@mpOutlineSpecifiers        = (/"China","China:Provinces"/)

plot = gsn_csm_contour_map(wks,ter,res)
end

E_PM_10变量如下：
float E_PM_10(Time, emissions_zdim, south_north, west_east) ;
                E_PM_10:FieldType = 104 ;
                E_PM_10:MemoryOrder = "XYZ" ;
                E_PM_10:description = "EMISSIONS" ;
                E_PM_10:units = "mol km^-2 hr^-1" ;
                E_PM_10:stagger = "" ;
                E_PM_10:coordinates = "XLONG XLAT" ;

qiu79879

PM_10=wrf_user_getvar（a，“E_PM_10”）
ter=PM_10(0,0,:,:)

悠然一梦123

qiu79879 发表于 2016-12-8 20:46
PM_10=wrf_user_getvar（a，“E_PM_10”）
ter=PM_10(0,0,:,:)

修改添加你所说的
E_PM_10=wrf_user_getvar(a，"E_PM_10",0)
ter=E_PM_10(0,:,:)，出现如下错误提示

(0)        check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.
(0)        A valid latitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)            'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)        check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.
(0)        A valid longitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)            'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

qiu79879

ter=wrf_user_getvar(a，"E_PM_10",0)
就行了，如果你不做循环

悠然一梦123

qiu79879 发表于 2016-12-8 20:59
ter=wrf_user_getvar(a，"E_PM_10",0)
就行了，如果你不做循环

数据就一个时次的，没做循环

sun92

wrfchem用ncl作图变量在哪里查啊？