DIMS = AVE.dimensions
CHL.createDimension('time', None)
CHL.createDimension('lon', DIMS['lon'])
CHL.createDimension('lat', DIMS['lat'])
ncvar = CHL.createVariable("lchlm", 'f', ('lat', 'lon'))
ncvar[:] = a
setattr(ncvar, "missing_value", fillValue)
AVE.close()
CHL.close()
if __name__ == "__main__":
import read_descriptor
RD = read_descriptor.read_descriptor('VarDescriptor_2.xml')
var = 'P_i'
INPUT_AVEDIR = "MODEL/AVE_FREQ_1/"
AGGREGATE_AVEDIR = "POSTPROC/TMP/"
avelist = [
"ave.20130101-12:00:00.N1p.nc", "ave.20130102-12:00:00.N1p.nc",
"ave.20130103-12:00:00.N1p.nc"
]
avelist = [
"ave.20130101-12:00:00.nc", "ave.20130102-12:00:00.nc",
"ave.20130103-12:00:00.nc"
]
avelist = [INPUT_AVEDIR + f for f in avelist]
nsub = len(SUBlist)
SubDescr =""
CoastDescr =""
DepthDescr =""
for i in SUBlist : SubDescr +=str(i) + ", "
for i in COASTNESS_LIST : CoastDescr +=str(i) + ", "
for i in DEPTHlist : DepthDescr +=str(i) + ", "
RD = read_descriptor.read_descriptor(args.descriptor)
NATIVE_VARS = RD.NATIVE_VARS
AGGREGATE_DICT = RD.AGGREGATE_DICT
SOME_VARS = RD.SOME_VARS
OUTPUT_DIR_INT = BASEDIR + 'INTEGRALS/'
OUTPUT_DIR_MPR = BASEDIR + 'STAT_PROFILES/'
OUTPUT_DIR_PRO = BASEDIR + 'PROFILES/'
OUTPUT_DIR_PPN = BASEDIR + 'INTEGRALS/PPN/'
if rank==0 :
RD.printStatus()
print "INPUT_AVEDIR", INPUT_AVEDIR
for DIR in [OUTPUT_DIR_INT, OUTPUT_DIR_MPR, OUTPUT_DIR_PPN, TMPDIR]:
if doStatistics: os.system("mkdir -p " + DIR)
