try:
os.chdir(outdir)
except OSError:
print('Could not change directory to {}'.format(outdir))
else:
if cloud in ['A439', 'GMC1', 'GMC104', 'N59C']:
f12 = 115.2712
f13 = 110.2013
else:
f12 = 230.538
f13 = 220.399
for line in ['12','13']:
label = cloud+'_'+line
if line == '12':
pltprops(label, fghz=f12, dvkms=0.2, beam=3.5,
xplot=['rad_pc', 'vrms_k', 'area_pc2','area_pc2','siglum', '8um_avg'],
yplot=['vrms_k', 'mlumco', 'mlumco', 'mvir','sigvir', 'sigvir'],
xlims=[[-1,1.5], [-2,2], [-1,3], [-1,3], [-1,4], [-1,1]],
ylims=[[-2,1.5], [-1,5], [-1,6], [-1,6], [-1,4], [-1,4]],
pltname=[ 'rdv', 'dvflux', 'areaflux', 'areamv', 'bnd', 'sfnorm'])
else:
pltprops(label, fghz=f13, dvkms=0.5, beam=3.5,
xplot=['rad_pc', 'vrms_k', 'area_pc2' ],
yplot=['vrms_k', 'mlumco', 'mlumco' ],
xlims=[[-1,1.5], [-2,2], [-1,3] ],
ylims=[[-2,1.5], [-1.5,4.5], [-1,4] ],
pltname=[ 'rdv', 'dvflux', 'areaflux' ])
# Exit statement
os.chdir(old_dir)
#!/usr/bin/env python
import os
from pltprops import pltprops
workdir = '../props'
old_dir = os.getcwd() # returns absolute path
if not os.path.isdir(workdir):
os.makedirs(workdir)
try:
os.chdir(workdir)
for line in ['12', '13']:
label = 'pcc_'+line
if line == '12':
pltprops(label, dvkms=0.2, beam=1.5,
xplot=['rad_pc', 'vrms_k', 'area_pc2', 'siglum', 'mlumco', 'mlte'],
yplot=['vrms_k', 'mlumco', 'mlumco', 'sigvir', 'mvir', 'mvir'],
xlims=[[-1,1.5], [-2,2], [-1,3], [-1,4], [-0.5,5.5], [-0.5,5.5]],
ylims=[[-2,1.5], [-1,5], [-1,5], [-1,4], [-0.5,5.5], [-0.5,5.5]],
pltname=[ 'rdv', 'dvflux', 'areaflux', 'bnd', 'virial', 'ltevir'])
else:
pltprops(label, fghz=220.399, dvkms=0.2, beam=1.5,
xplot=['rad_pc', 'vrms_k', 'area_pc2', 'mlte'],
yplot=['vrms_k', 'mlumco', 'mlumco', 'mvir'],
xlims=[[-1,1.5], [-2,2], [-1,3], [-0.5,4]],
ylims=[[-2,1.5], [-1.5,4.5], [-1,4], [-0.5,4]],
pltname=[ 'rdv', 'dvflux', 'areaflux', 'ltevir'])
finally:
os.chdir(old_dir)