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
