def ret_pr(prtype,yeart,mont,dayt,hourt): #--- if prtype =="GSMaP": a2pr = gsmap_func.timeave_gsmap_backward_saone(yeart,mont,dayt,hourt,1) #-- test -- #a2pr = ones([120,360],float32) if prtype =="GSMaP.dec": a2pr = gsmap_func.timeave_gsmap_backward_org(yeart,mont,dayt,hourt,1) #--- return ma.masked_equal(a2pr, miss) * calcoef
tempdir ="/home/utsumi/mnt/iis.data2/GSMaP/standard/v5/hourly/%04d/01/01"%(year+1) tempname = tempdir + "/gsmap_mvk.%04d0101.0000.v5.222.1.dat.gz"%(year+1) #tempdir ="/media/disk2/data/GSMaP/sa.one/1hr/ptot/%04d%02d/"%(year+1,1) #tempname = tempdir + "/gsmap_mvk.1rh.%04d0101.0000.v5.222.1.sa.one"%(year+1) if not os.access(tempname, os.F_OK): continue #------------------- print yfirst, year, mon, day, "nhour=",nhour, "percent=",percent #------------------- for hour in lhour: itimes = itimes + 1 itimes_seg = itimes_seg + 1 #-- load --------------- a2in = gsmap_func.timeave_gsmap_backward_org(year,mon,day,hour, nhour) #-- count -- a2count = a2count + ma.masked_where(a2in == miss, a2one).filled(0.0) #-- init --------------- a2in_seg = a2in[yfirst:ylast, :].copy() if itimes == 1: a2stck = a2in_seg else: a2stck = r_[a2stck, a2in_seg] #----------------------- if itimes_seg > (ntop * numrat): print "sort!! on ",year,mon,day,hour a2stck = a2stck.reshape(-1,ylast-yfirst,nx_org) for iiy in range( (ylast -yfirst)): iy = yfirst + iiy