爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6068|回复: 4

[作图] linint2_Wrap进行插值后画图,图是颠倒的。

[复制链接]
发表于 2016-3-13 15:40:12 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
本帖最后由 这个也有- 于 2016-3-13 15:48 编辑

以下是我的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"

begin

path      = "/cma/g3/ncl/tas_Amon_FGOALS2-s_rcp45_r1i1p1_200601-210012.nc"
nc_file   = addfile(path,"r")

;print(nc_file)

wks = gsn_open_wks("pdf","yuan")

temp = nc_file->tas(:,:,:)
lat1  = nc_file->lat(:)
lon1  = nc_file->lon(:)

res    = False
map = gsn_csm_contour_map_ce(wks,temp(5,:,:),res)

   latx = fspan(-90.0,90.0,362)
   lonx = fspan(0.0,360.0,720)

   tempx = new((/1140,362,720/),typeof(temp))

   tempx = linint2_Wrap(lon1,lat1,temp,True,lonx,latx,0)

;  tempx!0 = "lat"
;  tempx!1 = "lon"
;  tempx&lat = fspan(-90.0,90.0,362)
;  tempx&lon = fspan(0.0,360.0,720)
;  tempx&lat@units = "degrees_north"
;  tempx&lon@units = "degrees_east"

   map = gsn_csm_contour_map_ce(wks,tempx(5,:,:),res)

end

我是想将我的数据插值成为0.5*0.5的,但是我插值后发现画出来的图是有些颠倒的,第一张图是插值前的数据画的图,第二张是插值后画的。
请问,我插值后输出的数据和插值前数据的存放位置是否发生了变化?比如:原数据中x方向第一个格点表示西经180度,插值后x方向第一个格点表示的也是西经180度么?


插值前数据画的图

插值前数据画的图

插值后画的图

插值后画的图
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-13 15:46:52 | 显示全部楼层
还有一个问题,就是当我不注释这几行的时候,会出现错误。
  tempx!0 = "lat"
  tempx!1 = "lon"
  tempx&lat = fspan(-90.0,90.0,362)
  tempx&lon = fspan(0.0,360.0,720)
  tempx&lat@units = "degrees_north"
  tempx&lon@units = "degrees_east"
出现的问题如下图,第三十二行是tempx&lat = fspan(-90.0,90.0,362)这一行,但是这个的大小和上面定义的是一样的大小。

QQ图片20160313154843.png
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-15 16:17:44 | 显示全部楼层
有没有人知道怎么回事
密码修改失败请联系微信:mofangbao
发表于 2016-3-16 15:21:56 | 显示全部楼层
我知道,已回复
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-16 23:22:45 | 显示全部楼层
关于图是颠倒的问题解决了, lonx = fspan(0.0,359.5,720)改为这样就好了。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表