def dealias_berrimah_volume(filename, **kwargs):
pattern=kwargs.get('pattern', 'BerrimaVol')
deal_add=kwargs.get('deal_add', '_deal')
raw_path=kwargs.get('raw_path', '/data/uf_ber/')
deal_path=kwargs.get('deal_path', '/data/deal_ber/')
deal_files=os.listdir(deal_path)
raw_files=os.listdir(raw_path)
#deal_files=[]
#raw_files=[]
#for file in uf_files:
# if 'deal' in file:
# deal_files.append(file)
# else:
# raw_files.append(file)
#
if not(filename in raw_files):
raise IOError, 'File not there'
deal_files.sort()
raw_files.sort()
this_date_str=filename[len(pattern):-3].replace('_', ' ')
this_date_str=this_date_str[0:this_date_str.find(' ')+5]
sonde_name=read_sounding.make_deal_sonde(this_date_str)
deal_date=sonde_name[0:sonde_name.find('_')]
date_num_list=[]
for file in deal_files:
deal_str=file[len(pattern):-8].replace('_', ' ')
deal_str=deal_str[0:deal_str.find(' ')+5]
date_num_list.append(datestr2num(deal_str))
offset=array(date_num_list)-datestr2num(this_date_str)
idec_where=nwhere(offset < 0.0)[0]
offset_least=100.0
if len(idec_where)!=0: offset_least=offset[idec_where[-1]]
print "**********************"
if abs(offset_least) < 30.0/(60.0*24.0):
deal_fname=deal_files[idec_where[-1]]
print "Found a previous de-aliased file ", deal_fname
else:
print "No de-aliased file found within 30 minutes, only using sounding"
deal_fname='dummy'
print "********************"
outfile=filename[0:-3]+deal_add+".uf"
cwd=os.getcwd()
os.chdir('/home/scollis/bom_mds/dealias/')
execbl='./FourDD_berrimah '
command=execbl+deal_path+deal_fname+' '+raw_path+filename+' '+deal_path+outfile+' '+ sonde_name+' '+deal_date+' '+'0 1 1 1'
print command
os.system(command)
os.chdir(cwd)
return outfile
def dealias_single_volume(filename, **kwargs):
pattern=kwargs.get('pattern', 'Gunn_pt_')
deal_add=kwargs.get('deal_add', '_deal')
lassen_path=kwargs.get('lassen_path', '/bm/gscratch/scollis/lassen_cpol/')
uf_path=kwargs.get('lassen_path', '/bm/gscratch/scollis/uf_cpol/')
deal_files=os.listdir(uf_path)
raw_files=os.listdir(lassen_path)
#deal_files=[]
#raw_files=[]
#for file in uf_files:
# if 'deal' in file:
# deal_files.append(file)
# else:
# raw_files.append(file)
#
if not(filename in raw_files):
raise IOError, 'File not there: '+filename
deal_files.sort()
raw_files.sort()
this_date_str=filename[len(pattern):-11]
#this_date_str=this_date_str[0:this_date_str.find(' ')+5]
sonde_name=read_sounding.make_deal_sonde(this_date_str)
deal_date=sonde_name[0:sonde_name.find('_')]
date_num_list=[]
for file in deal_files:
deal_str=file[len(pattern):-8]
#print deal_str
#deal_str=deal_str[0:deal_str.find(' ')+5]
#print deal_str
date_num_list.append(datestr2num(deal_str))
offset=array(date_num_list)-datestr2num(this_date_str)
idec_where=nwhere(offset < 0.0)[0]
offset_least=100.0
if len(idec_where)!=0: offset_least=offset[idec_where[-1]]
print "**********************"
if abs(offset_least) < 30.0/(60.0*24.0):
deal_fname=deal_files[idec_where[-1]]
print "Found a previous de-aliased file ", deal_fname
else:
print "No de-aliased file found within 30 minutes, only using sounding"
deal_fname='dummy'
print "********************"
outfile=filename[0:-11]+deal_add+".uf"
cwd=os.getcwd()
os.chdir('/flurry/home/scollis/bom_mds/dealias/')
execbl='./FourDD_lassen '
command=execbl+uf_path+deal_fname+' '+lassen_path+filename+' '+uf_path+outfile+' '+ sonde_name+' '+deal_date+' '+'1 1 1 1'
print command
os.system(command)
os.chdir(cwd)
return outfile
