[ncl-talk] How to get the regional data of SMAP L4 *(HDF5 format)

[易路]

Hi, all
Sorry to disturb you for technical help. 
 
I am trying to extract the regional data from the soil product SMAP L4 which is stored in the format of HDF5. I wrote the ncl script taking referecnce to the example listed on the web of http://www.ncl.ucar.edu/Applications/HDF.shtml. But the extracted data were not in the demond region. So would you please help me to find the faults in my ncl script?  The main part of this script is as follow, and the detail information can be downloaded from the accessory. 
********main part of the script *********
 nFiles  = dimsizes(iFiles)  
  f       = addfile(iFiles(0), "r")
  lat2d     = f->cell_lat

  lon2d     = f->cell_lon
 
 nSmapLat = dimsizes(lat2d(:,0))
 nSmapLon = dimsizes(lon2d(0,:))
 minSmapLat = min(lat2d)
 maxSmapLat = max(lat2d)
 minSmapLon = min(lon2d)
 maxSmapLon = max(lon2d)
 
 lat=fspan(minSmapLat, maxSmapLat, nSmapLat)   
    lat!0="lat"
    lat&lat=lat
    lat at units="degrees_north"
    lat at long_name="Latitude"
 lon=fspan(minSmapLon,maxSmapLon, nSmapLon)
    lon!0="lon"
    lon&lon=lon
    lon at units="degrees_east"
    lon at long_name="Longitude"
   
  lat2d!0="lat"
  lat2d!1="lon"
  lat2d&lon=lon
  lat2d&lat=lat
 
  lon2d!0="lat"
  lon2d!1="lon"
  lon2d&lon=lon
  lon2d&lat=lat
 
  minLatBoundary  = 25
  maxLatBoundary  = 35
  minLonBoundary  = 110
  maxLonBoundary  = 120

  sm   = f->sm_rootzone
  printVarSummary(sm)
  sm!0="lat"
  sm!1="lon"
  sm&lon=lon
  sm&lat=lat
  printVarSummary(sm)
  
  sm_regional = sm({minLatBoundary:maxLatBoundary},{minLonBoundary:maxLonBoundary})
***********************************
 
Thanks a lot for your precious attentions and time!
 
Yi Lu
 






易路

南大邮件系统/学生/博士生/12级博士生

南京市汉口路22号 




 


从腾讯企业邮箱发来的超大附件

SMAP_L4_SM_gph_20150815T013000_Vv3030_001.h5 (136.33M, 2018年01月12日 11:27 到期)进入下载页面：http://exmail.qq.com/cgi-bin/ftnExs_download?k=523465623d8c4a92182f6ffe1038075719575c5a030c065706190304015b1d51000c504f0e0e55561b075257010c0452050d070051381e65657924326874043a65793a0547506f57060550520f0905310605565207086f334007555107670055181a4b4c5f0d3058&t=exs_ftn_download&code=64eb780e&s=email
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171214/120a2045/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smap_regional_results.pdf
Type: application/octet-stream
Size: 136196 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171214/120a2045/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smap_regional.ncl
Type: application/octet-stream
Size: 4055 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171214/120a2045/attachment-0001.obj>