-
Notifications
You must be signed in to change notification settings - Fork 1
/
run_make_mms.py
43 lines (34 loc) · 1.36 KB
/
run_make_mms.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
import inspect, os, sys, json
filename = inspect.getframeinfo(inspect.currentframe()).filename
sys.path.append(os.path.dirname(os.path.realpath(filename)))
from VLBI_pipe_functions import *
try:
# CASA 6
import casatools
from casatasks import *
casalog.showconsole(True)
casa6=True
except:
# CASA 5
from casac import casac as casatools
from taskinit import casalog
casa6=False
casalog.origin('vp_make_mms')
inputs = load_json('vp_inputs.json')
params = load_json(inputs['parameter_file_path'])
steps_run = load_json('vp_steps_run.json',Odict=True,casa6=casa6)
gt_r = load_json('vp_gaintables.last.json', Odict=True, casa6=casa6)
gt_r['make_mms'] = {'gaintable':[],'gainfield':[],'spwmap':[],'interp':[]}
rmdirs(['%s/%s.ms.flagversions'%(params['global']['cwd'],params['global']['project_code'])])
mmsfile='%s/%s.ms'%(params['global']['cwd'],params['global']['project_code'])
msfile='%s_2.ms'%(mmsfile.split('.ms')[0])
## Make mms data-set
os.system('mv %s %s'%(mmsfile,msfile))
partition(vis=msfile,\
outputvis=mmsfile,
separationaxis=params['make_mms']['separationaxis'],
numsubms= params['make_mms']['numsubms'])
rmdirs([msfile])
save_json(filename='%s/vp_gaintables.last.json'%(params['global']['cwd']), array=gt_r, append=False)
steps_run['make_mms'] = 1
save_json(filename='%s/vp_steps_run.json'%(params['global']['cwd']), array=steps_run, append=False)