NCL 问题求助..??? 面板图 地图设置 出现错误？

guizi

我在读取SODA 海流 和 温度数据时候，进行月平均处理

运行脚本错误如下？
Copyright (C) 1995-2013 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.2
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.

Variable: temp_s
Type: float
Total Size: 155520000 bytes
            38880000 values
Number of Dimensions: 4
Dimensions and sizes:    [time | 36] x [depth | 40] x [lat | 150] x [lon | 180]
Coordinates:
            time: [552.5..587.5]
            depth: [5.01..5375]
            lat: [-44.75..29.75]
            lon: [30.25..119.75]
Number Of Attributes: 7
  pointwidth :     1
  units :    Celsius_scale
  standard_name :    sea_water_temperature
  long_name :    seasonal means: TEMPERATURE
  missing_value :    -9.99e+33
  _FillValue :    -9.99e+33
  season :    <ARRAY of 12 elements>
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
warning:["SetValues.c":474]:SetValuesChild never occurred on vector_vector.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on vector_vector.Title: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on vector_vector.LabelBar: Error in plotManagerClass
fatal:VectorPlotSetValues: error setting overlay object view
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerSetValues: error setting transformation
fatal:PlotManagerSetValues: error setting overlay plot member view
warning:["SetValues.c":474]:SetValuesChild never occurred on map.TickMark: Error in plotManagerClass
warning:["SetValues.c":474]:SetValuesChild never occurred on map.Title: Error in plotManagerClass
fatal:MapPlotSetValues: error setting overlay object view
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA
fatal:PlotManagerPreDraw: error setting transformation
fatal:_NhlPlotManagerDraw: PreDraw error
....... balabala  还很多 一样的错误






NCL脚本 如下
;================================================;
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/csm/shea_util.ncl"
; ================================================;
begin
;************************************************
;*************  open files and  read the data **********                                    
  filsu = systemfunc("ls /home/Model/Test/Iod/test_data/U/*")
  fu    = addfiles(filsu, "r")
  ListSetType (fu,"cat")

  filsv = systemfunc("ls /home/Model/Test/Iod/test_data/V/*")
  fv    = addfiles(filsv, "r")
  ListSetType (fv,"cat")

  filst = systemfunc("ls /home/Model/Test/Iod/test_data/T/*")
  ft    = addfiles(filst, "r")
  ListSetType (ft,"cat")


  u    = fu[:]->u(:,:,{-45:30},{30:120})    ; read in example data [2D only here]
  v    = fv[:]->v(:,:,{-45:30},{30:120})
  temp = ft[:]->temp(:,:,{-45:30},{30:120})
  dep  = ft[:]->depth(:)


;*************   deal with the data *************
    dsizes_temp = dimsizes(temp)
    TIME        = dsizes_temp(0)
    LEV         = dsizes_temp(1)
    LAT         = dsizes_temp(2)
    LON         = dsizes_temp(3)

;**********  calculate the season mean
   temp_s = month_to_season12(temp)
   u_s    = month_to_season12(u)
   v_s    = month_to_season12(v)
   printVarSummary(temp_s)

   
   temp_s12 = new((/12,LEV,LAT,LON/),"float")
   u_s12    = new((/12,LEV,LAT,LON/),"float")
   v_s12    = new((/12,LEV,LAT,LON/),"float")

   do it=0,11
      u_s12(it,:,:,:)    = dim_avg_n_Wrap(u_s(it:TIME-1:12,:,:,:), 0)
      v_s12(it,:,:,:)    = dim_avg_n_Wrap(v_s(it:TIME-1:12,:,:,:), 0)
      temp_s12(it,:,:,:) = dim_avg_n_Wrap(temp_s(it:TIME-1:12,:,:,:), 0)
      
;************************************  No 3 month mean
;      u_s12(it,:,:,:)    = dim_avg_n_Wrap(u(it:TIME-1:12,:,:,:), 0)
;      v_s12(it,:,:,:)    = dim_avg_n_Wrap(v(it:TIME-1:12,:,:,:), 0)
;      temp_s12(it,:,:,:) = dim_avg_n_Wrap(temp(it:TIME-1:12,:,:,:), 0)  
  end do
;   printVarSummary(temp_s12)



;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","vector")             ; open a ps file
  gsn_define_colormap(wks,"BlAqGrYeOrRevi200")  ; choose color map
  plot = new(12,graphic)                          ; create a plot array
  
  res                      = True               ; plot mods desired
  res@gsnDraw               = False    ; don't draw
  res@gsnFrame              = False    ; don't advance frame
  res@cnFillOn             = True               ; turn on color for contours
  res@cnLinesOn            = False              ; turn off contour lines
  res@cnLineLabelsOn       = False              ; turn off contour line labels
  res@gsnScalarContour     = True               ; contours desired
  res@gsnSpreadColors      = True               ; use full color map
  res@gsnSpreadColorStart  = 17                 ; start at color 17
  res@gsnSpreadColorEnd    = 200                ; end at color 200
  res@gsnAddCyclic         = False
  res@gsnMaximize          = True               ; make ps, pdf, eps large
;
; This is not necessary in V6.1.0 and later. Named colors can
; be used without having to first add them to the color map.
;
  i = NhlNewColor(wks,0.7,0.7,0.7)              ; add gray to colormap
  res@mpLandFillColor       = "gray"            ; set land to be gray

  res@lbOrientation            = "Vertical"     ; vertical label bar
  res@pmLabelBarOrthogonalPosF = -0.01          ; move label bar closer
  res@lbLabelStride            = 4

; note, when doing a subregion, NCL determines the range of the data from
; the full domain. If you wish to just consider the domain you are plotting,
; you must manually set those levels.

  res@cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
  res@cnMinLevelValF       = 0.0               ; set min contour level
  res@cnMaxLevelValF       = 30.                 ; set max contour level
  res@cnLevelSpacingF      = 0.5               ; set contour spacing

  res@gsnLeftString   = " "              ; add the gsn titles
  res@gsnCenterString = " "
  res@gsnRightString  = " "

  res@vcRefMagnitudeF           = 0.3             ; define vector ref mag
  res@vcRefLengthF              = 0.05           ; define length of vec ref
  res@vcRefAnnoOrthogonalPosF   = -1.0            ; move ref vector
  res@vcRefAnnoArrowLineColor   = "black"         ; change ref vector color
  res@vcRefAnnoArrowUseVecColor = False           ; don't use vec color for ref

  res@vcGlyphStyle            = "CurlyVector"     ; turn on curley vectors
  res@vcLineArrowColor        = "white"           ; change vector color
  res@vcLineArrowThicknessF   = 2.0               ; change vector thickness
  res@vcVectorDrawOrder       = "PostDraw"        ; draw vectors last


; note that the gsn_csm_*map_ce templates automatically set
; res@mpLimitMode="LatLon" for you. If you are plotting a different projection,
; you may have to set this resource.

   res@mpMinLatF            = -45      ; range to zoom in on
   res@mpMaxLatF            =  30.
   res@mpMinLonF            =  30.
   res@mpMaxLonF            = 120.


   do ik=0,LEV-1   
      do it = 0, 11
         plot=gsn_csm_vector_scalar_map_ce(wks,u_s12(it,ik,::7,::7),v_s12(it,ik,::7,::7),\
                                    temp_s12(it,ik,:,:),res)

      end do
;************************************************
; create panel
;************************************************
    resP            = True                         ; modify the panel plot
    resP@gsnFrame         = False                  ; don't advance panel plot
;    resP@gsnMaximize           = True              ; make as large as possible
    resP@txString   = " "
    resP@gsnPanelLabelBar    = True                ; add common colorbar
    resP@lbLabelFontHeightF  = 0.007               ; make labels smaller
    resP@gsnPanelBottom   = 0.05                   ; add space at bottom
    resP@gsnPanelFigureStrings= (/"DJF)","JFM)","FMA)","MAM)","AMJ)","MJJ)","JJA)","JAS)","ASO)","SON)","OND)","NDJ)"/)      ; add strings to panel
    resP@amJust   = "TopLeft"   
   
    gsn_panel(wks,plot,(/4,3/),resP)               ; now draw as one plot
   

;***********************************************   
    txres               = True
    txres@txFontHeightF = 0.015
    gsn_text_ndc(wks,"Figure   "+(ik+1)+" : Tempreture with depth  "+dep(ik)+" m ",0.5,0.02,txres)

   frame(wks)

   end do

    delete(wks)   ; Make sure PS file is closed
    cmd = "convert  -rotate -90 -delay 25 vector.ps  vector.gif"
    system(cmd)

end



求大家 给看看撒 小弟感激不尽 o(╯□╰)o

guizi

已解决 低级失误  

stopraining

想请教lz错误
“warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: error initializing map: MAPINT/MDPINT - MAP HAS ZERO AREA”
究竟是什么问题造成的，画面板图遇到类似问题
谢谢lz！

guizi

stopraining 发表于 2013-9-8 23:46

                               
登录/注册后可看大图

想请教lz错误
“warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: err ...

do it = 0, 11
         plot=gsn_csm_vector_scalar_map_ce(wks,u_s12(it,ik,::7,::7),v_s12(it,ik,::7,::7),\
                                    temp_s12(it,ik,:,:),res)

      end do


我的 错误好像是这块  修改如下 ：当时
do it = 0, 11
         plot（it）=gsn_csm_vector_scalar_map_ce(wks,u_s12(it,ik,::7,::7),v_s12(it,ik,::7,::7),\
                                    temp_s12(it,ik,:,:),res)

      end do



脚本贴出来 呀

ling19880131

请问楼主，您的soda海流资料从哪里下载的？我现在怎么无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）了？！可否发个链接？

sth.special.

stopraining 发表于 2013-9-8 23:46
想请教lz错误
“warning:MapSetTrans: map limits invalid - using maximal area
fatal:MapSetTrans: err ...

刚学习ncl也遇到同样的问题 请问是怎么解决的呢

打破砂锅纹到底

感谢楼主，完全一样的低级错误。。害我找了大半天- - -