EOF分解（ncl）

安之若素

做了一个温度的EOF分解，不知道对不对。用的是每年冬季的月数据，共268个月。想问的是如果我想要年数据，可以用月数据做分解得出的时间序列在进行每个冬季月份平均得到每年冬季的时间序列么。

lat  =   dble2flt(temper&lat)
    rad     =   4.*atan(1.)/180.
    clat    =   cos(lat*rad)
    clat!0  =   "lat"
    clat&lat=   lat
    clat    =   where(clat .le. 0.0, 0.0, clat)
    print("min and max of clat: ")
    printMinMax(clat,0)
    plat    =   clat
    plat    =   where(ismissing(clat),clat,sqrt(clat))
    printVarSummary(plat)

    temper    = temper*conform(temper,plat,0)

          neof   = 3        ; number of EOFs
    optEOF = True
    optEOF@jopt = 0   ; This is the default; most commonly used; no need to specify.
    optETS = False
    eof_ref     =   eofunc_Wrap(temper, neof, optEOF)
    eof_ts_ref  =   eofunc_ts_Wrap(temper, eof_ref, optETS)

    eof_ts_ref=-1*dim_standardize_n( eof_ts_ref, 0, 1)

    eof_norm_ref = eof_ref                               ; create an array w meta data
    do ne = 0, 2
        eof_norm_ref(ne,:,:) = (/ regCoef(eof_ts_ref(ne,:), temper) /)
    end do

追风灬少年

一直都看不懂，实在太难了

dqhahaha

EOF分解做完也不知道是什么 稀里糊涂的

lgz

关注一下

owl113

   rad     =   4.*atan(1.)/180. 麻烦问一下这个什么意思