ncl做时间平均图

不骗你哦

变量：
Variable: f
Type: file
filename:       1
path:   /cygdrive/e/datesurface/indicatedate/1.grib
   file global attributes:
   dimensions:
      initial_time0_hours = 6938
      g0_lat_1 = 361
      g0_lon_2 = 720
   variables:
      float LGWS_GDS0_SFC ( initial_time0_hours, g0_lat_1, g0_lon_2 )
         center :       European Center for Medium-Range Weather Forecasts (RSMC
)
         long_name :    East-West surface stress
         units :        N m**-2 s
         _FillValue :   1e+20
         level_indicator :      1
         gds_grid_type :        0
         parameter_table_version :      128
         parameter_number :     180
         forecast_time :        3
         forecast_time_units :  hours
脚本：
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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
; ================================================;
begin
;=================================================;
; open file and read in data
;=================================================;
  f  = addfile("/cygdrive/e/datesurface/indicatedate/1.grib","r")
  u  = f->LGWS_GDS0_SFC(6937,360,719)
aveu = dim_avg_n(u,0)

  wks  = gsn_open_wks("pdf","11111110")     ; open a ps file
gsn_define_colormap(wks,"gui_default")  ; choose a colormap
  res                     = True          ; plot mods desired
  res@cnFillOn            = True         ; turn on color
  
;---This resource not needed in NCL V6.1.0
  res@gsnSpreadColors     = True          ; use full range of colormap  
  plot = gsn_csm_contour_map_ce(wks,aveu,res)     
  end
不知道哪里出错了 求大神指点啊

river

发错版块了吧

尽头的尽头

楼主仔细看看求平均的函数的用法吧

freekiller

应该是这里的问题吧 u  = f->LGWS_GDS0_SFC(6937,360,719)
有6937个时次？

錯莂牸

u  = f->LGWS_GDS0_SFC(6937,360,719)改成u  = f->LGWS_GDS0_SFC
或者u  = f->LGWS_GDS0_SFC(：,：,：)
你只读了（360,719)这个格点的6937时次

不骗你哦

錯莂牸 发表于 2013-12-26 16:09

                               
登录/注册后可看大图

u  = f->LGWS_GDS0_SFC(6937,360,719)改成u  = f->LGWS_GDS0_SFC
或者u  = f->LGWS_GDS0_SFC(：,：,：)
你 ...

我改过你那个试过 直接来个错误 你的意思就是不定义 根据数据本身的数组定义 呗 错误好像是error[12]什么的

不骗你哦

freekiller 发表于 2013-12-26 09:29

                               
登录/注册后可看大图

应该是这里的问题吧 u  = f->LGWS_GDS0_SFC(6937,360,719)
有6937个时次？

对有这么多个啊 我前面不是有资料的说明么

freekiller

不骗你哦 发表于 2013-12-30 17:10

                               
登录/注册后可看大图

对有这么多个啊 我前面不是有资料的说明么

这态度我都不敢回复了，你哪句读取了一个时次的数据，还求个毛的时间平均

易小凯

只有一个时间一个点的数据啊