爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 7926|回复: 0

[其他] HDFEXP文件不支持读取

[复制链接]
发表于 2019-4-17 20:40:06 | 显示全部楼层 |阅读模式

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

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

x
我生成文件,用HDFEXP不能读取是怎样?希望大家给建议。代码如下
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

    outfile = "/output/pr/pr_regrid.nc"
    ;diri = "/input/data/pr/"
    fil = "/input/data/pr/pr_GPCP-SG_L3_v2.3_197901-201710.nc"
    f = addfile(fil, "r")
    tim = f->time
    lat = f->lat
    lon = f->lon
    n1 = dimsizes(tim)
    n2 = dimsizes(lat)
    n3 = dimsizes(lon)



;printVarSummary(pr1)
    pr = new((/n1,n2,n3/), float)

    do i = 0, n1-1
        pr(i,:,:) = f->pr(i,:,:)
        ; write loop content
    end do


    Opt                = True

    Opt@SrcFileName    = "src_SCRIP.nc"
    Opt@DstFileName    = "dst_SCRIP.nc"
    Opt@WgtFileName    = "NCEP_2_Rect.nc"
    Opt@ForceOverwrite = True

    Opt@SrcGridMask    = where(.not.ismissing(pr),1,0)   ; use SrcMask2D in NCL V6.2.1 and earlier

    Opt@DstGridType    = "1x1"

    ;Opt@DstLLCorner    = (/-89.5.d, 89.d/)
    ;Opt@DstURCorner    = (/-180.d, 180.d/)

    pr_regrid = new((/n1,181,361/), float)

     do i = 0, n1-1
       pr_regrid(i,:,:) =ESMF_regrid(pr(i,:,:),Opt)
        ; write loop content
    end do

    copy_VarAtts(pr, pr_regrid)



    printVarSummary(pr_regrid)
    ;printVarSummary(src_SCRIP.nc)

    ;............................................
    ;write .nc
    ;............................................
   fout = addfile(outfile,"c")
   filedimdef(fout,"time",-1,True) ;-- make time and UNLIMITED dimension
fout->pr = pr_regrid





密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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