def main(*args):
if int(MON) < 10:
m = '0' + str(MON)
else:
m = str(MON)
if int(DAY) < 10:
d = '0' + str(DAY)
else:
d = str(DAY)
if int(TIME) < 10:
t = '0' + str(TIME)
else:
t = str(TIME)
### Read args ###
Data = fromfile(
DATADIR + 'Ens_' + En + '_' + str(Year) + str(m) + str(d) + str(t) +
'.bin', 'float32')
varName = 'Ens' + '_' + str(En)
itemName = 'PRCP'
gtUnit = 'kg/m2/s'
### Setting ###
#-- Resolution
nX = 500
nY = 440
nZ = 1
#-- Time Step and deltT
nT = 39
dT = 3600
#-- Coordination
GLON = 'GLON500M'
GLAT = 'GLAT440IM'
GATM = 'SFC1'
#missing_value = 1.E20
#missing_value = -999.
if os.path.isdir(outDir) != True:
raise IOError
### Pre-Process ###
os.environ['F_UFMTENDIAN'] = 'big' # don't need ??
outPath = os.path.join(outDir, '%s_%s%s%s%s.gt' % (varName, Year, m, d, t))
Data = Data.reshape(nT, nZ, nY, nX)
# DTIME = [ datetime(Year,1,1,0)+timedelta(seconds=dT)*i for i in range(nT) ]
DTIME = [
datetime.datetime(Year, MON, DAY, TIME) +
datetime.timedelta(seconds=dT) * i for i in range(0, nT)
]
### Writing ###========================
gtOut = gtool(outPath, 'ow') # ow: over write
gtOut[itemName] = Data, DTIME
gtOut[itemName].TITL1 = varName
gtOut[itemName].UNIT = gtUnit
gtOut[itemName].AEND1 = nX
gtOut[itemName].AEND2 = nY
gtOut[itemName].AEND3 = nZ
gtOut[itemName].AITM1 = GLON
gtOut[itemName].AITM2 = GLAT
gtOut[itemName].AITM3 = GATM
gtOut.save(10)
print 'outPath: ', outPath
print gtOut[itemName].header
def __getitem__(self,fName):
srcPath = os.path.join(self.srcDir,fName)
return gtool(srcPath)