cmorph降水融合产品作图出现错误

好久不见!

降雨落区和数据气象网给出的差距很大，官方给出的基本在中部和南部地区，我用出图都显示在新疆和东北地区呢？请问是哪里出现的问题呢还有想修改一下，出图的区域 ，想指定输出某经纬度范围内的数据请问应该修改哪里呢？谢谢大家！

脚本如下：res@mpMinLatF  =  18.
res@mpMaxLatF  =  56.
res@mpMinLonF  =  72.
res@mpMaxLonF  =  135.

res@tmXTOn  = False
res@tmYROn  = False
res@gsnMajorLonSpacing = 10

res@tmXBOn =  True
res@tmXBMajorThicknessF  = 2.0         ;default 2
res@tmXBMajorLengthF  = 0.01         ;default 0.02  

res@tmYLOn =  True
res@tmYLMajorThicknessF  = 2.0         ;default 2
res@tmYLMajorLengthF  = 0.01         ;default 0.02

res@tmYLMinorOn = True
res@tmYLMinorThicknessF  = 1.0
res@tmYLMinorLengthF   =  0.005     ;default 0.01
res@tmYLMinorPerMajor = 4

nlat = 440
nlon = 700
ncrain = 308000
ntim = 24
r = new((/ntim,nlat,nlon/),"double")

;time = new((/ntime/),"string")
;setfileoption("grd","ReadByteOrder","LittleEndian")

;reading files
        
     tim = (/ "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060400.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060401.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060402.grd",\
                      "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060403.grd",\
                   "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060404.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060405.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060406.grd",\
                   "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060407.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060408.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060409.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060410.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060411.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060412.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060413.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060414.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060415.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060416.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060417.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060418.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060419.grd",\
                   "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060420.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060421.grd",\
                    "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060422.grd",\
                     "SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-2017060423.grd"/)

lon = fspan(70.05,139.95,nlon)
lat = fspan(15.05,58.95,nlat)
lat!0 = "lat"
lat@long_name = "latitude"
lat@units = "degrees_north"
lon!0 = "lon"
lon@units = "degrees_east"
lon@long_name = "longitude"

do t = 0, ntim-1  
rain = fbindirread(tim(t),0,(/ncrain/),"double")
k=0
do i = 0, nlat-1
  do j = 0, nlon-1
   r(t,i,j) = rain(k)
   k = k+1
  end do
end do
end do

;r!0 = "time"
;r!1 = "lat"
;r!2 = "lon"
;r&time = tim
;r&lat = lat
;r&lon = lon
;r@units = "mm/hour"
;r@long_name = "china daily accumulated rain"
;r@_FillValue = -999.00
;n=num(ismissing(r))
;print("n="+n)
;printVarSummary(r)

res@gsnAddCyclic = False

wks = gsn_open_wks("pdf","rain")
cmap=(/(/1.,1.,1./),(/0.,0.,0./),(/1.,1.,1./),(/.26,.88,.95/),(/.153,.706,.733/),(/.08,.61,.26/),\
     (/.28,.90,.49/),(/.89,.88,.18/),(/.77,.75,.10/),(/.65,.45,.16/),(/.612,.125,.137/),\
     (/.8392,.2118,.2275/),(/.9059,.0745,.5098/),(/.65,.10,.86/),(/.24,.16,.45/)/)
gsn_define_colormap(wks,cmap)
;gsn_draw_colormap(wks)
;res@gsnSpreadColorStart = 0
res@cnLevelSelectionMode = "ExplicitLevels"   ; "ManualLevels"  ,"automaticlevels" , "explicitlevels"
res@cnLevels   = (/0.01,0.5,1,5,10,20,30,40,50,60,80,100/)

;res@cnMinLevelValF  =  .5
;res@cnMaxLevelValF  = 100.
;res@cnLevelSpacingF =  10.
res@gsnLeftString   = "rain"
res@gsnRightString  = "[mm~N~]"
res@tiMainString    = "2017-06-04-00:00 - 2017-06-04-23:00"
  
frain=r(0,:,:)
do t = 1,ntim-1
  srain = r(t,:,:)
  arain =srain + frain
  frain = arain
end do

print(frain)
filo="crain"
fbinrecwrite(filo,-1,frain)

frain!0 = "lat"
frain!1 = "lon"
frain&lat = lat
frain&lon = lon
frain@units = "mm/hour"
frain@long_name = "china daily accumulated rain"
frain@_FillValue = -999.00
n=num(ismissing(frain))
print("n="+n)
printVarSummary(frain)

plot = gsn_csm_contour_map(wks,frain,res)
   
   config      =  True
   config@china  = True
   config@river  = True
   config@province = True
   config@nanhai  = True
   config@diqu   = False
   chinamap = add_china_map(wks,plot,config)
   
   draw(plot)
   frame(wks)


end

;=============================================================================================
谢谢指导

煮茶的水

数据应该没错，看新疆那的海南岛形状。。。

哒哒哒

请问您问题解决了吗？我也是一样的问题，可以告知一下解决办法吗？谢谢！

wjn

同样的问题，不知道楼主还有没有后续了

wjn

最近我画这个资料也出现同样的问题，后来发现把读数据的格式double换成float就好了

whj112233

wjn 发表于 2020-4-5 15:30
最近我画这个资料也出现同样的问题，后来发现把读数据的格式double换成float就好了

warning:ContourPlotSetValues: Data values out of range of levels set by EXPLICITLEVELS       mode
那样改了，好像有这样的问题，怎么调大这个范围还是有warning，这个应该怎么解决鸭