hourly daily pentad monthly

在知道calculate_monthly_values前我是这么做的:

3 hourly to daily

var_name=get_script_prefix_name
year=1958
fils=systemfunc("ls /data1/DuanGroup/zy/jra55/55/4times/"+var_name+"/*"+year+"*.nc")
a_in=addfiles(fils,"r")
names = getfilevarnames(a_in[0])       ; get CCM file's variable names and print them
print (names)

that_var=a_in[:]->$names(3)$
printVarSummary(that_var)
dims=dimsizes(that_var)
ntime=dims(0)
nlat=dims(1)
nlon=dims(2)
that_var_onedim=ndtooned(that_var)
that_var_ndim=onedtond(that_var_onedim,(/ntime/8,8,nlat,nlon/))
that_var_daily=dim_avg_n_Wrap(that_var_ndim,1)
copy_VarCoords(that_var(::8,:,:),that_var_daily)
printVarSummary(that_var_daily)

pentad to daily

target_var_num=7

file_list=systemfunc("ls /public/damtest/zy/data/godas_pentad/"+year+"/*.grb")
f_in=addfile(file_list(0),"r")
var_names=getfilevarnames(f_in)
print(var_names)
x5_0=f_in->$var_names(target_var_num)$
dim_x5_0 = dimsizes( x5_0 ) 
x5=new((/73,dim_x5_0(0),dim_x5_0(1)/),typeof(x5_0),x5_0@_FillValue)
do i=0,72
    f_in=addfile(file_list(i),"r")
    var_names=getfilevarnames(f_in)
    x5(i,:,:)=f_in->$var_names(target_var_num)$
end do
x5!0="ncl_join"

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;process
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
dim_x5 = dimsizes( x5 ) 
ntim5 = dim_x5(0) 
nlat = dim_x5(1) 
nlon = dim_x5(2) 
 
ntim = 5*ntim5 ; 5x the number of pentads = # days 
 
time5 = ispan(3,ntim,5) 
time5!0 = "time5" 
time5_at_units = "pentad day" 
time5&time5 = time5 
 
time = ispan(1,ntim,1) ; create time/day coordinate variable 
time!0 = "time" 
time_at_units = "day" 
time&time = time 
 
xTemp = linint1_Wrap (time5, x5(g0_lat_1|:,g0_lon_2|:,ncl_join|:), False, time, 0) 
xDay = xTemp(time|:,g0_lat_1|:,g0_lon_2|:) 
delete(xTemp) 
copy_VarAtts(x5, xDay) 
 
xDay!0 = "time" 
xDay!1 = "lat" 
xDay!2 = "lon" 
xDay&time= time 
xDay&lat = x5&g0_lat_1 
xDay&lon = x5&g0_lon_2 
 
printVarSummary( xDay ) 
printMinMax( xDay,False ) 

pentad to monthly

file_list=systemfunc("ls /data1/DuanGroup/zy/godas/pentad/equ/thermo_equ_*heat_budget_detail_glo_tro_.nc")
a=addfiles(file_list,"r")
sst1=a[:]->sw
 utc_date = cd_calendar(a[:]->time, 0)
 year_   = tointeger(utc_date(:,0))    ; Convert to integer for
 month  = tointeger(utc_date(:,1))    ; use sprinti 
 day    = tointeger(utc_date(:,2))
 hour   = tointeger(utc_date(:,3))
 minute = tointeger(utc_date(:,4))
 second = utc_date(:,5)

sst_6=new((/year_num,181,361/),"double")
do i=year_start,year_end
    month_id:=ind(month.eq.target_month.and.year_.eq.i)
    tmp=dim_avg_n_Wrap(sst1(month_id,:,:),0)
    sst_6(i-year_start,:,:)=tmp
end do
printVarSummary(sst_6)
printMinMax(sst_6, False)
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 220,492评论 6 513
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,048评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 166,927评论 0 358
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,293评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,309评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,024评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,638评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,546评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,073评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,188评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,321评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,998评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,678评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,186评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,303评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,663评论 3 375
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,330评论 2 358

推荐阅读更多精彩内容

  • 原文:How CALCULATE works in DAX链接:http://sqlblog.com/blogs/...
    PowerBI战友联盟阅读 19,288评论 5 21
  • 背景 一年多以前我在知乎上答了有关LeetCode的问题, 分享了一些自己做题目的经验。 张土汪:刷leetcod...
    土汪阅读 12,747评论 0 33
  • 爱对了是爱情,爱错了是青春。 ——左耳
    3350a63f8487阅读 131评论 0 0
  • 卖电子产品的声音 是黑夜的口舌 霓虹打算放大瞳孔 让属于无限的神秘 铺满脂粉与不安 远方不是地理名词 远方传递一种...
    孟小繁阅读 180评论 0 0
  • 荠菜,每地的叫法不同,枕头菜,菱角菜,护生草,都是它的名字。荠菜翠绿,叶片光滑而呈齿状,可拌、可炒、可烩、可汤。 ...
    糖糖婉儿阅读 324评论 1 1