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)