数组有值，却无法出图

longlivehj

追风的阳光 发表于 2014-12-5 21:54
Variable: dv
Type: float
Total Size: 469504 bytes

max=1.28247e+38
这个正常？

追风的阳光

longlivehj 发表于 2014-12-5 22:00
max=1.28247e+38
这个正常？

不好意思，回晚了，我查找了下值，那个异常值是在边界上
(0)        1.282467e+38
(0)        dv（0 ，0， 209）
其他dv都是在200以下，异常值在边界上也可以理解啊。。。

longlivehj

异常值不是缺测值，前者参与绘图，后者忽略。你把uv2dv_cdf的边界选项换一下看看，比如用0，直接设为缺测。

追风的阳光

longlivehj 发表于 2014-12-6 12:12
异常值不是缺测值，前者参与绘图，后者忽略。你把uv2dv_cdf的边界选项换一下看看，比如用0，直接设为缺测。

不知道什么情况了。。改为0依旧没用，wrf_contour输出图像大小依旧为2kb，gsn_csm_contour依旧输出不了图像（运行不终止），脚本之前都没有问题，都是输出语句这里出现了问题。。

longlivehj

先看看数据的范围printMinMax，再说改变边界选项有没有用。画图是下一步，数据问题还没有解决呢！

追风的阳光

longlivehj 发表于 2014-12-8 10:19
先看看数据的范围printMinMax，再说改变边界选项有没有用。画图是下一步，数据问题还没有解决呢！

除去 dv(0,0,209)=1.282467e+38，其他值都比较正常(由Max到Min）：

dv(0,14,218)=139.1
dv(0,20 ,219)=111.4231
dv(0,207,143)=104.9907
dv(0,163 ,150)=103.6495
.............
.............
dv(0,4 214)=-122.042
dv(0,153 40)=-122.4412
dv(0,221 76)=-137.763

追风的阳光

longlivehj 发表于 2014-12-8 10:19
先看看数据的范围printMinMax，再说改变边界选项有没有用。画图是下一步，数据问题还没有解决呢！

自己写中央差后可以绘出图像：
mapfactorm=wrffile->MAPFAC_M(0,:,:)

  dims=dimsizes(u)
  nyarps=dims(0)
  print(nyarps)
  nxarps=dims(1)
  print(nxarps)
;**********************divergence************
  d=4000.0
  div=new((/nyarps,nxarps/),"float")
  do i=0,nxarps-1
   do j=0,nyarps-1
     if(i.eq.0.and.j.eq.0) then
       div(j,i)=(u(j,i+1)-u(j,i))/(d*mapfactorm(j,i))+(v(j+1,i)-v(j,i))/(d*mapfactorm(j,i))
     end if
     if(i.eq.0.and.j.eq.nyarps-1) then
       div(j,i)=(u(j,i+1)-u(j,i))/(d*mapfactorm(j,i))+(v(j,i)-v(j-1,i))/(d*mapfactorm(j,i))
     end if
。。。。。。
但是绘出图像后，想要缩小区域就出现问题了，
我先后尝试了：
  res@mpLimitMode         = "Corners"
  res@mpLeftCornerLatF    = leftbottom_lat_c
  res@mpLeftCornerLonF    = leftbottom_lon_c
  res@mpRightCornerLatF   = rightup_lat_c
  res@mpRightCornerLonF   = rightup_lon_c
  ************************************
  res@mpLimitMode = "LatLon"
  res@mpMaxLatF = 38
  res@mpMinLatF = 32
  res@mpMinLonF = 116
  res@mpMaxLonF = 120
**************************************
及在输出函数中增加数值范围
plotdiv = gsn_csm_contour_map(fstwks,plotdiv(lat1:lat2,lon1:lon2),wres)
都没有作用，请问有什么好的办法呢。。或者哪里出现了问题。。