请大家进来看看帮帮忙，谢谢大家

kanyu

longlivehj 发表于 2014-8-11 22:30
没有就自己添加吧！
首先，要把你的坐标顺序调整一下，把bottom_top换到最左边，即第一维；然后，要为第一 ...

你好，想请问你，我想讲实测数据插值到wrfout中，然后输出茶之后的便粮食局，编写程序如下，但总是出现这样的错误：
fatal:Could not create (rain.nc)
fatal:Coordinate variables must be the same dimension as their dimension
fatal:No coordinate variable exists for dimension (lat_0) in variable (en_tem2)
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 35 in file real7

35行就是en_tem2&lat_0=ff->lat这个地方，不知道是什么问题，想请你解决一下

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
;********************************************************
  f = addfile("/home/kanyu/shanghai/shanghai/20080825/mp6/wrfout_d01_2008-08-25_06.nc","r")

  xlon   = f->XLONG(0,:,:)              ;wrfout的经纬度，由于是兰伯特投影，因此不是均匀网格，经纬度是2维的
  xlat   = f->XLAT(0,:,:)
  
  ff = addfile("/home/kanyu/shanghai/shanghai/real_/082500.nc","r")

  lat_0=ff->lat  
  lon_0=ff->lon  

  en_tem1  = ff->crain (0,:,:)   
  ;printVarSummary(lat_0)
  en_tem1!0="lat_0"      ;赋予数据属性
  en_tem1!1="lon_0"
  en_tem1&lat_0=lat_0(:)
  en_tem1&lon_0=lon_0(:)
  
  en_tem2=rgrid2rcm(en_tem1&lat_0,en_tem1&lon_0,en_tem1,xlat,xlon,0)
  out=addfile("rain.nc","c")  
  en_tem2!0="lat_0"      ;赋予数据属性
  en_tem2!1="lon_0"
  en_tem2&lat_0=ff->lat  
  en_tem2&lon_0=ff->lon
  ;print(lat_0(:))

end

腐朽的我Zero

官网上copy的time_hight，网址有吗，同学

Maziy

腐朽的我Zero 发表于 2016-4-1 22:30
官网上copy的time_hight，网址有吗，同学

现在才看到帖子，不好意思，你现在还需要吗