把中国以外区域设置为缺省值(问题已解决)

蓝馨蝶韵

新手编程序就是不细心，在朋友的帮助下搞定了，没什么大问题，现在的程序是修改正确的，跟大家分享下。
ccru-china-pre.nc中只有中国区域有值，其他都为缺省；现在把SPEI_01.nc写成与它一致的nc，只保留中国区域的值，把中国区域以外的值都写为缺省
; ==============================================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

           system ("rm spei-china.nc")

begin

     fi=addfile("ccru-china-pre.nc","r")
     pre1=fi->pre(:,:,:)

        time1  = fi->time
        lon1   = fi->lon  
        lat1   = fi->lat

     pre_china=new((/72,128/),float)   
     pre_china=pre1(0,:,:)
        copy_VarCoords(pre1(0,:,:),pre_china)
        printVarSummary(pre_china)

     in1 =addfile("/home/lenovo/SPEI_01.nc","r")   
     spei_new=in1->spei(:,216:287,504:631)

    mlon =128
    nlat =72
    nt =1332

     spei_china=new((/1332,72,128/),float)

do i=0,nlat-1
     do j=0,mlon-1
         if ismissing(pre_china(i,j))  then
            spei_china(:,i,j)= 1e+30
         else
            spei_china(:,i,j)= spei_new(:,i,j)
         end if
end do
end do
            spei_china@_FillValue    = 1e+30
            spei_china@missing_value = 1e+30

        copy_VarCoords(pre1(0:1331,:,:),spei_china)
        printVarSummary(spei_china)

     fo=addfile("spei-china.nc","c")
     fo->spei = spei_china
end

红尘_海

LZ 能否上传分享一下ccru-china-pre.nc文件？

风之牧语

学习了，能否上传分享一下ccru-china-pre.nc文件？

风之牧语

跪求ccru-china-pre.nc文件啊{:5_275:}

mxgang

楼主，如果前后两种数据的空间分辨率不一致，这种方法还可以生效吗？

zongqing123

能否上传分享一下ccru-china-pre.nc文件？

tongqing

看看

zongqing123

{:eb315:}{:eb315:}