ncl计算最长连续无雨天数CDD

dh5002

请各位指教一下，求秋季最长连续无雨的天数，想知道哪里有问题

f = addfile("pre_daily_9_10_11.nc", "r")

pre = f->pre(:,{25:36},{100:111})
year=ispan(1961,2020,1)
nyr=dimsizes(year)

f1=addfile("pre_son.nc","r")
time = f1->time;此时间用于输出nc文件用

dim    = dimsizes(pre)   
ntim    = dim(0)   ;计算时间数
nlat    = dim(1)   ;计算纬度数
nlon    = dim(2)   ;计算经度数

npre=onedtond(ndtooned(pre),(/nyr,91,dimsizes(pre&lat),dimsizes(pre&lon)/)) ;onedtond将一维数组转为多维，ndtooned将多维转为一维
copy_VarMeta(pre(0,:,:),npre(0,0,:,:))
printVarSummary(npre)

;;---------------------------------------------------------------------------------------
CDD=new((/nyr,dimsizes(pre&lat),dimsizes(pre&lon)/),"float",-9999.0)
CDD=0.
n = new((/nyr,dimsizes(pre&lat),dimsizes(pre&lon)/),"float",-9999.0)
n = 0.
count = new((/nyr,dimsizes(pre&lat),dimsizes(pre&lon)/),"float",-9999.0)
count = 0.

do i=0,nlat-1
  do j=0,nlon-1
    do t = 0, nyr-1
      do d = 0, 90
        if (.not.ismissing(npre(t,d,i,j)).and.npre(t,d,i,j).lt.1) then
          n(t,i,j) = n(t,i,j)+1
        else
          count(t,i,j) = n(t,i,j)
          n(t,i,j) = 0
        end if
      end do
      CDD(t,i,j)=max(count(:,i,j))
    end do
  end do
end do

CDD!0="time"
CDD&time = time;输出nc时用f1的时间，其余用year,这样输出的文件才有年份
CDD!1="lat"
CDD&lat = pre&lat
CDD!2="lon"
CDD&lon = pre&lon
printVarSummary(CDD) ;min=2.61047   max=162.979

ncdf_out = addfile("CDD.nc" ,"c");先创建（输出这个数据，用cdo选出大于等于TH95的格点，然后再计算极端降水量）
ncdf_out->CDD= CDD

alextsui

你count的维度定义应该有问题，count应该一年秋季当中存在的个数是不定的啊，看你的定义是一年只有一个。
请问你现在问题解决了么？

dh5002

alextsui 发表于 2021-12-21 10:37
你count的维度定义应该有问题，count应该一年秋季当中存在的个数是不定的啊，看你的定义是一年只有一个。
...

嘿嘿，我没弄这个了，放弃了