求指教周变温的语句怎么写

zhangdapao

求指教：每周的变温怎么写？想写第二周的平均减去第一周的平均，然后通过差值的正负来看温度是升还是降。
如下：这是06年六月第一周最高温的平均。请问第二周的最高温平均减去第一周的最高温平均应该怎么写呢？
f = addfile("D:\Temdatabase\tmax.2006.nc", "r")
tmax = f->tmax(:,{minlat-1:maxlat+1},{minlon-1:maxlon+1})
   year = (/2006,2006/) month = (/6,6/) day = (/1,7/)   
doy = day_of_year(year, month, day)-1  
tmax_JJA_ave = dim_avg_n_Wrap(tmax(doy(0):doy(1),:,:), 0)
printVarSummary(tmax_JJA_ave)
因为是新手，只能在上文中去改。
doy = day_of_year(year, month, day1)-1  
tmax = f->tmax(:,{minlat-1:maxlat+1},{minlon-1:maxlon+1})
另外请问上面这两句是什么意思呢？

一大碗年糕

f = addfile("D:\Temdatabase\tmax.2006.nc", "r")
tmax = f->tmax(:,{minlat-1:maxlat+1},{minlon-1:maxlon+1})  ;读取2006年的逐日一定经纬度区域范围内的温度最大值数据
   year = (/2006,2006/) month = (/6,6/) day = (/1,7/)   
doy = day_of_year(year, month, day)-1  
;day_of_year是返回指定日期在当年所有天数中的序号（从1开始计数），用来定位6月第一周，-1是因为tmax这个数组第一维日期序列是0开始的
tmax_JJA_ave = dim_avg_n_Wrap(tmax(doy(0):doy(1),:,:), 0)；对第一周最高温作平均
printVarSummary(tmax_JJA_ave)

zhangdapao

我自己瞎捉摸的这样改并不行。doy这个可以改成doy1、doy2吗 ？
f = addfile("D:\Temdatabase\tmax.2006.nc", "r")
   
tmax = f->tmax(:,{minlat-1:maxlat+1},{minlon-1:maxlon+1})

   

year = (/2006,2006/)
    month = (/6,6/)
    day1 = (/1,7/)
   
doy1 = day_of_year(year, month, day1)-1
   
tmax_JJA_ave1 = dim_avg_n_Wrap(tmax(doy(0):doy(1),:,:), 0)

year = (/2006,2006/)
    month = (/6,6/)   day2 = (/8,14/)
   
doy2 = day_of_year(year, month, day2)-1
   
tmax_JJA_ave2 = dim_avg_n_Wrap(tmax(doy(0):doy(1),:,:), 0)
   
tmax_JJA_ave= (tmax_JJA_ave2-tmax_JJA_ave1)
printVarSummary(tmax_JJA_ave)

   

zhangdapao

报错是这样的
fatal:syntax error: line 33 in file /cygdrive/d/jiaoben/avetempv.ncl before or n                                                                                                                                  ear day2
    month = (/6,6/)   day2
-------------------------^

fatal:Syntax Error in block, block not executed
fatal:error at line 126 in file /cygdrive/d/jiaoben/avetempv.ncl

zhangdapao

顶顶 求指教

一大碗年糕

zhangdapao 发表于 2020-4-13 12:40
我自己瞎捉摸的这样改并不行。doy这个可以改成doy1、doy2吗 ？
f = addfile("D:\Temdatabase\tmax.2006.nc ...

你既然定义了doy1和doy2那么作平均的时候数组用的应该是doy1(0),doy1(1),doy2(0),doy2(1)把

zhangdapao

一大碗年糕 发表于 2020-4-13 22:15
你既然定义了doy1和doy2那么作平均的时候数组用的应该是doy1(0),doy1(1),doy2(0),doy2(1)把

了解了，谢谢您。因为这几天有事，好几天没看论坛了。非常感谢！

kermit

学习学习

DGGG

{:eb513:}