求助求助！！站点降水 报错问题

syl136

请问下，我想画站点降水的图，脚本是根据例子改的，可是一直提示这个错误，我是在大型机上面用的ncl，之前运行其他脚本好像也有这个错误，求大神解答，file:///C:/Users/gua/Desktop/捕获.PNG

syl136

图片没上传成功。。这个错误fatal:CvtStringGenArrayToColorIndexGenArray: Unable to convert string "LightGray" to requested type
warning:Error retrieving resource mpFillColors from args - Ignoring Arg
fatal:CvtStringToColorIndex: Unable to convert string "black" to requested type
warning:Error retrieving resource mpGeophysicalLineColor from args - Ignoring Arg
fatal:CvtStringToColorIndex: Unable to convert string "black" to requested type
warning:Error retrieving resource mpNationalLineColor from args - Ignoring Arg

syl136

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"                                                                                                                                                                                                                                                                                                                                            
   
5. 
   
6. begin
   
7. 
   
8. ;diri = "/home/shenyanling/program/ncl/14062608"
   
9. ;fils = systemfunc ("ls "+ diri + "140*.000")
   
10. ;f   = addfiles("14062608.asc","r")  
    
11.   ;  data = asciiread("fils",(/488,20/),"float")
    
12. c := asciiread("/home/shenyanling/program/ncl/14062608.000",-1,"string")                              
    
13. delim =" "
    
14. lon = stringtofloat(str_get_field(c(1:),2,delim))
    
15. lat = stringtofloat(str_get_field(c(1:),3,delim))
    
16. rain = stringtofloat(str_get_field(c(1:),5,delim))
    
17. olon=new(66,"float");
    
18.   olat=new(66,"float");
    
19.   data1=new((/66,66/),"float");
    
20.   do i=0,65
    
21.    olon(i)=72+i
    
22.   end do
    
23.   do j=0,65
    
24.    olat(j)=17+j
    
25.   end do
    
26.   
    
27. ;;;;;设置属性，符合netcdf规定的数据格式;;;;;;;;
    
28.   olon!0        ="lon"
    
29.   olon@long_name="lon"
    
30.   olon@units    ="degrees-east"
    
31.   olon&lon      =olon
    
32.   olat!0        ="lat"
    
33.   olat@long_name="lat"
    
34.   olat@units    ="degrees-north"
    
35.   olat&lat      =olat
    
36.   
    
37.   rain@_FillValue = 999999.000000  
    
38.   rscan = (/10,5,3/)   ;连续的有效半径大小，最大为10，依次递减
    
39.   data1 = obj_anal_ic_deprecated(lon,lat,rain,olon,olat,rscan, False)  ;Creanm插值
    
40. ;data1=linint2_points(lon,lat,pre10,True,olon,olat,0)
    
41.   printVarSummary(data1)
    
42.   ;exit
    
43. 
    
44. 
    
45. 
    
46. wks = gsn_open_wks("ps","station")
    
47. gsn_define_colormap(wks, "BlRe")  
    
48. 
    
49. res=True
    
50. ;res@gsnSpreadColors      = True         ; spread out color table
    
51. res@cnFillOn = True            ; turn on color
    
52. res@cnLinesOn = False   
    
53. res@lbLabelBarOn        = False
    
54. res@cnInfoLabelOn       = False
    
55. res@cnLineLabelsOn      = False
    
56. 
    
57. res@mpPerimOn   = True
    
58. res@gsnMaximize=True
    
59. res@pmTickMarkDisplayMode   = "Always"
    
60. res@gsnDraw   = False        ; don't draw yet
    
61. res@gsnFrame  = False
    
62. res@mpGeophysicalLineColor = "black"
    
63. res@mpNationalLineColor    = "black"
    
64. res@mpDataSetName         = "Earth..4"
    
65. res@mpDataBaseVersion     = "MediumRes"
    
66. res@mpOutlineSpecifiers = ("China:states")
    
67. res@mpNationalLineThicknessF= 0.1
    
68. res@mpOutlineDrawOrder    = "PostDraw"
    
69. res@mpLimitMode = "LatLon"    ; Limit the map view.
    
70. res@gsnAddCyclic = False
    
71. ;res@mpMinLonF   = min(lon2)
    
72. ;res@mpMaxLonF   = max(lon2)
    
73. ;res@mpMinLatF   = min(lat2)
    
74. ;res@mpMaxLatF   = max(lat2)
    
75. 
    
76. 
    
77. ;>============================================================<
    
78. ;                      add China map
    
79. ;>------------------------------------------------------------<
    
80. ;cnres           = True
    
81. ;cnres@china     = True       ;draw china map or not
    
82. ;cnres@river     = True       ;draw changjiang&huanghe or not
    
83. ;cnres@province  = True       ;draw province boundary or not
    
84. ;cnres@nanhai    = True       ;draw nanhai or not
    
85. ;cnres@diqu      = False       draw diqujie or not
    
86. plot=gsn_csm_contour_map(wks,data1,res)
    
87. end  
    

复制代码

糯米团子爱芥末

想问+1.。。

sarshar

gsn_define_colormap(wks, "BlRe") 这句里定义的颜色只用蓝色到红色，没有黑色。。。
请参考http://www.ncl.ucar.edu/Document/Graphics/named_colors.shtml
以及http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml比较简单的方法就是你换一个colortable的名字，或者把black换成blue

syl136

sarshar 发表于 2015-5-13 16:20
gsn_define_colormap(wks, "BlRe") 这句里定义的颜色只用蓝色到红色，没有黑色。。。
请参考http://www.nc ...

谢谢谢~~我试试！！纠结了好几天了

syl136

sarshar 发表于 2015-5-13 16:20
gsn_define_colormap(wks, "BlRe") 这句里定义的颜色只用蓝色到红色，没有黑色。。。
请参考http://www.nc ...

想问下为什么不管我换什么色板  gsn_define_colormap(wks,"rainbow") 都提示fatal:CvtStringGenArrayToColorIndexGenArray: Unable to convert string "LightGray" to requested type

sarshar

syl136 发表于 2015-5-14 20:25
想问下为什么不管我换什么色板  gsn_define_colormap(wks,"rainbow") 都提示fatal:CvtStringGenArrayTo ...

不好意思，上次只看到了“black”的错误提示，没看到“lightgray”的。。。
NCL的色名中应该是没有lightgray这个名字的，可以改用“gray”+数字的方式，数字越大颜色越浅。
另外关于颜色这方面我比较建议你自定义颜色而不使用色板，比如下面这个：
colors = (/"white","black","White","RoyalBlue","LightSkyBlue",\
             "PowderBlue","lightseagreen","PaleGreen","Wheat","Brown",\
             "Pink"/)
gsn_define_colormap(wks, colors)  
把自己要用到的颜色定义出来，比较不容易出错 :)

syl136

sarshar 发表于 2015-5-20 11:03
不好意思，上次只看到了“black”的错误提示，没看到“lightgray”的。。。
NCL的色名中应该是没有light ...

好的！！谢谢你！！

杨雨轩

syl136 发表于 2015-5-13 15:17

请问楼主，第14行c（1:）是什么意思