예제 #1
0
파일: GB_lib.py 프로젝트: gbolzon/ogstm
    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]
예제 #2
0
파일: aveScan.py 프로젝트: inogs/bit.sea
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)