爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 9796|回复: 12

[作图] ERA-20C和ERA-Interim资料拼接画图出问题,求指教

[复制链接]
发表于 2016-7-3 10:33:44 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 subtropical 于 2016-7-3 10:35 编辑

我需要1901-2000年的感热资料,然后好像只搜到ERA20C有这么长的。下载了ERA-20C的1901-2000年的资料但是发现只有1901-1978年的是有数据的,也不知道是为什么,然后就把ERA-Interim的1979年的资料接在后面,两个资料都是Instantaneous surface sensible heat flux精度都是1*1的,我想把我定义的一个指数回归到感热场上。但是结果却空白的
111.png
112.png
但是我把这两种资料分别单独看的时候都是有数据的,想请教各路大神,帮忙指教一下。不胜感激。谢谢
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-7-3 16:18:38 | 显示全部楼层
有人用过这种资料吗?实在是不知道哪出错了,因为单独拿ERA20C的1901-1978年的资料做结果是没问题的,不知道为什么加了后面的79-00年就成这样了。。。{:5_275:}{:5_275:}
密码修改失败请联系微信:mofangbao
 成长值: 19710
发表于 2016-7-4 07:04:34 | 显示全部楼层
只能说先检查拼接后的数据是否正常,也许你拼接的时候就弄错了。。。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-7-4 09:49:48 | 显示全部楼层
兰溪之水 发表于 2016-7-4 07:04
只能说先检查拼接后的数据是否正常,也许你拼接的时候就弄错了。。。

谢谢版主的回复,我检查了一下,的确拼接后的数组好像没有数据,可是为什么呢?
可以麻烦您帮忙看一下吗,谢谢版主
1.png 2.png
min_lat =  -30.0
  max_lat =90.0
  min_lon = 30.0
  max_lon = 180.0
  yrStrt1 = 1901
  yrLast1 = 1978
   yrStrt2 = 1979
  yrLast2 = 2000
season = "DJF"    ; choose chun seasonal
f1=addfile("./era20monofdaily.nc","r")
  f2=addfile("./era.interim.monofdaily1.1.nc","r")
  TIME1   = f1->time
  TIME2   = f2->time
YYYY1   = cd_calendar(TIME1,-1)/100                 ; entire file
  iYYYY1  = ind(YYYY1.ge.yrStrt1 .and. YYYY1.le.yrLast1)
   YYYY2   = cd_calendar(TIME2,-1)/100                 ; entire file
  iYYYY2 = ind(YYYY2.ge.yrStrt2 .and. YYYY2.le.yrLast2)
   ishf1 = short2flt(f1->ishf)
      ishf2 = short2flt(f2->ishf)
      gt1=ishf1(iYYYY1,:,:)
        gt2=ishf2(iYYYY2,:,:)
  printVarSummary(gt1)
  printVarSummary(gt2)
  lat=f1->latitude
lon=f1->longitude   
  
dsizes_x1=dimsizes(gt1)
ntim1=dsizes_x1(0)
nlat1=dsizes_x1(1)
nlon1=dsizes_x1(2)
dsizes_x2=dimsizes(gt2)
ntim2=dsizes_x2(0)
nlat2=dsizes_x2(1)
nlon2=dsizes_x2(2)
gt = new((/ntim1+ntim2,nlat1,nlon1/),"float")
time=fspan(0,ntim1+ntim2-1,ntim1+ntim2)
gt!0="time"  
       gt!1="lat"
      gt!2="lon"
      gt&time=time
      gt&lat=gt1&latitude
      gt&lon=gt2&longitude
      time!0         = "time"
time@long_name = "time"
time@units     = "yr"
  lat!0 = "lat"
        lat@long_name = "latitude"
lat@units     = "degrees_north"
      lon!0 = "lon"
  lon@long_name = "longitude"
  lon@units     = "degrees_east"
  printVarSummary(gt)
      
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-7-4 15:14:51 | 显示全部楼层
兰溪之水 发表于 2016-7-4 07:04
只能说先检查拼接后的数据是否正常,也许你拼接的时候就弄错了。。。

版主,我发现我的问题了,原。。。。来。。。。是。。。。没。。。赋。。。值。。。。被我自己打败了。麻烦您了,谢谢您{:eb511:}
密码修改失败请联系微信:mofangbao
发表于 2016-7-5 14:18:48 | 显示全部楼层
这两种资料可以接在一起用么?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-7-5 21:36:59 | 显示全部楼层
66kqkq 发表于 2016-7-5 14:18
这两种资料可以接在一起用么?

应该可以的,因为它们的资料类型是一样的,我在想ERA20C会不会是考虑到已经有1979年后的这些资料了所以就做到1978年,我也不知道,只是猜测
密码修改失败请联系微信:mofangbao
发表于 2016-7-7 21:39:57 | 显示全部楼层
subtropical 发表于 2016-7-5 21:36
应该可以的,因为它们的资料类型是一样的,我在想ERA20C会不会是考虑到已经有1979年后的这些资料了所以就 ...

不可以的吧   资料的同化方法什么的  应该不一样吧
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-7-8 09:50:37 | 显示全部楼层
66kqkq 发表于 2016-7-7 21:39
不可以的吧   资料的同化方法什么的  应该不一样吧

但是结果做出来没问题啊
密码修改失败请联系微信:mofangbao
发表于 2017-9-15 08:13:43 | 显示全部楼层
请问ERA-20C和ERA Interim的拼接效果可以吗? 我也需要用到数据是1949年-2016年不知道应该怎么选择
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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