convert2nc转tigge资料出现问题

小虹

大家好，我刚刚接触ncl，用convert2nc这个命令将grib格式的tigge资料转为nc格式，然后画图，觉得数据南北反向了
我用的是2008年1月27日的资料，ncl脚本如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
begin
;********************************************************
; open file

  ff_con = addfile ("20080127-con.nc", "r")
  con_tem  = ff_con->t_P1_L100_GLL0(0,7,:,:)  
   wks= gsn_open_wks("png","tigge_con")
  gsn_define_colormap(wks,"BlAqGrYeOrRe")

  res=True
  res@cnFillOn              = True            ; turn on color
  res@cnLinesOn             = False           ; turn contour lines
  res@gsnSpreadColors       = True            ; use entire color map
  res@lbLabelAutoStride     = True            ; let NCL determine label spacing
  res@tfDoNDCOverlay       = True   
   plot = gsn_csm_contour_map(wks,con_tem,res)

end

画出来底层的温度资料如下图

1月份应该北半球冷，南半球暖啊，可是这个图却相反啊，求大神们解答到底是哪里出问题了？

kongfeng0824

你说说绘图整个反了了还是解释结果绘出来的反了呢？

小虹

kongfeng0824 发表于 2013-10-4 15:34

                               
登录/注册后可看大图

你说说绘图整个反了了还是解释结果绘出来的反了呢？

就是画出来的图不对啊，应该不是这个样子的，应该是北半球冷一些，南半球暖一些

小虹

找到原因了，把这句去掉就可以了res@tfDoNDCOverlay       = True。
数据对应的纬度坐标起始是90N，终止是90S，用上述语句的时候要求数据与图的坐标要匹配，但是我数据的坐标与图的坐标南北向正好相反了。
不知道NCL画图时默认的坐标原点是什么，是（90S,180W）吗？还请大神们前来解释下。

homezhaow

小虹 发表于 2013-10-4 16:21

                               
登录/注册后可看大图

找到原因了，把这句去掉就可以了res@tfDoNDCOverlay       = True。
数据对应的纬度坐标起始是90N，终止是 ...

默认是 [0,0] 大西洋视角
要用太平洋视角的话
请用  res@mpCenterLonF  = 180

homezhaow

小虹 发表于 2013-10-4 16:21

                               
登录/注册后可看大图

找到原因了，把这句去掉就可以了res@tfDoNDCOverlay       = True。
数据对应的纬度坐标起始是90N，终止是 ...

res@tfDoNDCOverlay       = True
另外本来就是用来调整纬度的，因为有的人会自然把纬度和经度一样按照从小到大的顺序排列。但实际图上格点从上到下，对应的纬度是从大到小排的。所以发现图反了，就调一下，反正是二选一的事儿，ture不对就false

小虹

homezhaow 发表于 2013-10-6 20:07

                               
登录/注册后可看大图

res@tfDoNDCOverlay       = True
另外本来就是用来调整纬度的，因为有的人会自然把纬度和经度一样按照从 ...

嗯，好的，多谢指点

lisixiuhe

{:eb502:}{:eb502:}{:eb502:}搂住真棒

chenhongfengye

homezhaow 发表于 2013-10-6 19:52
默认是 [0,0] 大西洋视角
要用太平洋视角的话
请用  res@mpCenterLonF  = 180

请问我为什么EC的Total Precipitation都是0 ？这个是降水吗？

liulian

~~~~~~~~~~~