ncl怎么给数组增加坐标属性

浅水叶

pfile=addfile("/data4/TRMM/1999/3B42.19990625.12.7.HDF","r")
prec=pfile->precipitation

lon=fspan(-180,180,1440)
lat=fspan(50,-50,400)
lon@units="degrees_east"
lat@units="degrees_north"
prec!0="lon"
prec!1="lat"
prec&lon=lon
prec&lat=lat
res =True
res@mpMinLatF = -10.0
res@mpMaxLatF = 55.0
res@mpMinLonF = 60.
res@mpMaxLonF = 180.0
wks=gsn_open_wks("X11","pre")
plot=gsn_csm_contour_map(wks,prec,res)


我这代码想画一下TRMM的降水，但是提示没有有效的坐标属性。。。不知道怎么弄，初次使用，求教

若兰

prec应该自带经纬度吧，你这样重置不了，如果需要插值重设经纬度需要用linint2函数

longlivehj

能把错误提示信息、相关变量信息发上来么？

浅水叶

longlivehj 发表于 2015-11-10 20:50
能把错误提示信息、相关变量信息发上来么？

原来是prec的经纬顺序搞错了，gsn_contour(prec（lon|：，lat|：）)这样会报错，好像画等值线需要先lat，再lon的顺序。。

gyq1991

浅水叶 发表于 2015-11-11 16:15
原来是prec的经纬顺序搞错了，gsn_contour(prec（lon|：，lat|：）)这样会报错，好像画等值线需要先lat， ...

能把正确脚本贴上来吗？

gyq1991

浅水叶 发表于 2015-11-11 16:15
原来是prec的经纬顺序搞错了，gsn_contour(prec（lon|：，lat|：）)这样会报错，好像画等值线需要先lat， ...

能把正确脚本贴上来吗？

gyq1991

浅水叶 发表于 2015-11-11 16:15
原来是prec的经纬顺序搞错了，gsn_contour(prec（lon|：，lat|：）)这样会报错，好像画等值线需要先lat， ...

能把正确脚本贴上来吗？

巧克力瑞纳冰

我也遇到了这个问题，求解
f = addfile("/home/twhh/data/landcovertype.hdf","r")

LCT = f->Majority_Land_Cover_Type_1(:,:)

LCT!0 = "lat"
LCT!1 = "lon"

lat = fspan(0,90,1801)
lon = fspan(0,180,3601)

lat@units = "degrees_north"
lon@units = "degrees_east"

LCT&lon = lon
LCT&lat = lat

printVarSummary(lat)
printVarSummary(lon)
printVarSummary(LCT)

之后提示错误：
fatal:Coordinate variables must be the same dimension as their dimension
fatal:No coordinate variable exists for dimension (lon) in variable (LCT)
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 14 in file landcovertype.ncl

fatal:Coordinate variables must be the same dimension as their dimension
fatal:No coordinate variable exists for dimension (lat) in variable (LCT)
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 15 in file landcovertype.ncl

变量如下：
Variable: lat
Type: float
Total Size: 7204 bytes
            1801 values
Number of Dimensions: 1
Dimensions and sizes:   [1801]
Coordinates:
Number Of Attributes: 1
  units :       degrees_north

Variable: lon
Type: float
Total Size: 14404 bytes
            3601 values
Number of Dimensions: 1
Dimensions and sizes:   [3601]
Coordinates:
Number Of Attributes: 1
  units :       degrees_east

Variable: LCT
Type: ubyte
Total Size: 25920000 bytes
            25920000 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 3600] x [lon | 7200]

95niande77

巧克力瑞纳冰 发表于 2020-8-15 14:54
我也遇到了这个问题，求解
f = addfile("/home/twhh/data/landcovertype.hdf","r")

后来咋解决的呀