ncl运行wrfout提取数据脚本出错

sushx

我用ncl提取wrfout文件脚本时出现如下错误：
[~/ssx/wrf3.8.1/WRFV3/run]$ncl zhandian.ncl
Copyright (C) 1995-2017 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.4.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
fatal:syntax error: line 1 in file zhandian.ncl before or near

^

fatal:error in statement

脚本内容如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
  a=addfile("/gpfshome/chenq/ssx/wrf3.8.1/WRFV3/run/wrfout_d02_2018-01-01_00_00_00","r")

  d1=wrf_user_getvar(a,"U10",-1)
  d2=wrf_user_getvar(a,"V10",-1)
  d=sqrt((d1)*(d1)+(d2)*(d2))
; RH=wrf_user_getvar(a,"RH",-1)
  T2=wrf_user_getvar(a,"T2",-1)
  ;d1=wrf_user_getvar(a,"DUST_1",-1)
  ;d2=wrf_user_getvar(a,"DUST_2",-1)
  ;d3=wrf_user_getvar(a,"DUST_3",-1)
  ;d4=wrf_user_getvar(a,"DUST_4",-1)
  ;d5=wrf_user_getvar(a,"DUST_5",-1)
  ;d=d1+d2+d3+d4+d5

; d=wrf_user_getvar(a,"O3",-1)
; d=d*1.96*1000  

  ;d=wrf_user_getvar(a,"PM2_5_DRY",-1)

  times = wrf_user_getvar(a,"times",-1)

  xlat=wrf_user_getvar(a,"XLAT",0)
  xlon=wrf_user_getvar(a,"XLONG",0)

  ;locationij=wrf_user_ll_to_ij(a,103.839462,36.039448,True);
  ;locationij=wrf_user_ll_to_ij(a,106.795478,39.668052,True) ;   
  ;locationij=wrf_user_ll_to_ij(a,106.822042,39.698685,True) ;   
  locationij=wrf_user_ll_to_ij(a,106.8,39.8,True) ;   

  ;print(locationij)
  lonlat=wrf_user_ij_to_ll(a,locationij(0),locationij(1),True);              
  ;print(lonlat)
  xx=xlat(locationij(1)-1,locationij(0)-1)                              
  ;print(xx)
  dust=d(:,locationij(1)-1,locationij(0)-1)
  T=T2(:,locationij(1)-1,locationij(0)-1)
;  R=RH(:,locationij(1)-1,locationij(0)-1)
;  print(dust)
asciiwrite("/gpfshome/chenq/ssx/result/Wuhai/V.xls", dust)
asciiwrite("/gpfshome/chenq/ssx/result/Wuhai/T.xls", T)
asciiwrite("/gpfshome/chenq/ssx/result/Wuhai/R.xls", R)
end

请问是哪里出了问题

小其其格

a=addfile("/gpfshome/chenq/ssx/wrf3.8.1/WRFV3/run/wrfout_d02_2018-01-01_00_00_00.nc","r")
试试这样呢

sushx

小其其格 发表于 2019-11-15 08:00
a=addfile("/gpfshome/chenq/ssx/wrf3.8.1/WRFV3/run/wrfout_d02_2018-01-01_00_00_00.nc","r")
试试这样 ...

也没用，最后重新输了一遍代码发现好了，不知道哪里的问题