Skip to content

Instantly share code, notes, and snippets.

View ridlo's full-sized avatar
🎯
Focusing

Ridlo W. Wibowo ridlo

🎯
Focusing
  • Atmospheric and Planetary Sciences, ITERA
  • Indonesia
View GitHub Profile
#!/bin/bash
# Post installing script for Ubuntu 14.04 Trusty Tahr Desktop
# First, check your internet connection (e.g. proxy)
#
# Check for updates
sudo apt-get -y update
# Upgrade Packages
sudo apt-get -y upgrade
# Major Upgrades
sudo apt-get -y dist-upgrade
import sys
import math
filename = 'lalapo.txt'
param = math.pi
while (len(sys.argv) > 1):
option = sys.argv[1]; del sys.argv[1]
if option == '-f':
filename = sys.argv[1]; del sys.argv[1]
@ridlo
ridlo / music.conf
Created November 30, 2016 14:28
MUSIC input file
[setup]
boxlength = 60 #Mpc
zstart = 50 #redshift at which simulation starts
levelmin = 7 #2^X
levelmin_TF = 7 #only zoom, otherwise same as others
levelmax = 7 #only zoom, otherwise same as others
padding = 7 #padding between different zones
overlap = 4 #number of padding grids when computing transfer function
ref_center = 0.0, 0.0, 0.0 #refinement
ref_extent = 0.0, 0.0, 0.0
InitCondFile /home/ridlo/Gadget-2.0.7/run4/ics_60_gadget_128.dat
OutputDir /home/ridlo/Gadget-2.0.7/run4/
EnergyFile energy.txt
InfoFile info.txt
TimingsFile timings.txt
CpuFile cpu.txt
RestartFile restart
SnapshotFileBase snapshot
@ridlo
ridlo / rockstar.cfg
Last active November 30, 2016 16:35
#Rockstar Halo Finder
#Parallel config file for multi-cpu, multi-snapshot halo finding
FILE_FORMAT = "GADGET2"
# For GADGET2, you may need to specify conversion parameters.
# Rockstar's internal units are Mpc/h (lengths) and Msun/h (masses)
GADGET_LENGTH_CONVERSION = 0.001
GADGET_MASS_CONVERSION = 1e+10
#!/bin/bash
#$ -N cosmo_gadget2
#$ -cwd
#$ -j y
#$ -m ea
#$ -M email@host.com
#$ -pe mpifill 16
#$ -V
#$ -o gadget-log
J1229+0203 / 3C273
J2253+1608 / 3C454.3
J2232+1143
J1924-2914
J1256-0547
J1058+0133
2017-04-13 20:23:30 INFO vishead ##########################################
2017-04-13 20:23:30 INFO vishead ##### Begin Task: vishead #####
2017-04-13 20:23:30 INFO vishead vishead(vis="uid___A002_X9998b8_X11eb.ms.split.cal",mode="list",listitems=['field', 'ptcs'],hdkey="",hdindex="",
2017-04-13 20:23:30 INFO vishead hdvalue="")
2017-04-13 20:23:30 INFO vishead field: ['J0423-0120', 'J0423-013', 'J0517-0520', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2', 'OMC-2']
2017-04-13 20:23:30 INFO vishead ptcs: ['UNKNOWN 04:23:15.80 -001.20.32.991', 'UNKNOWN 05:35:23.73 -005.01.57.682', 'UNKNOWN 05:35:25.44 -005.02.30.478', 'UNKNOWN 05:35:26.03 -005.03.10.493', 'UNKNOWN 05:35:26.34 -005.03.51.334', 'UNKNOWN 05:35:26.30 -005.05.29.722', 'UNKNOWN 05:35:26.82 -005.04.43.931', 'UNKNOWN 04:23:15.80 -001.20.32.991', 'UNKNOWN 05:17:28.11 -005.20.40.794', 'UNKNOWN 05:35:25.26 -005.02.34.397', 'UNKNOWN 05:35:14.52 -004.59.32.884', 'UNKNOWN 05:35:15.85 -005
def estimate_sensitivity(msName):
# find integration time
time_os = au.timeOnSource(msName)
print("-------------------------------------------------")
int_time = time_os[0]['minutes_on_source'] # always in field 0 (from split)
print "Integration time = ", int_time, " min"
vm = au.ValueMapping(msName)
nant = vm.numAntennas # find number of antenna
npol = vm.nPolarizations # number of polarizations
'2013.1.01099.S/science_goal.uid___A001_X13e_X168/group.uid___A001_X13e_X169/member.uid___A001_X13e_X16a', '2013.1.01099.S/science_goal.uid___A001_X13e_X16c/group.uid___A001_X13e_X16d/member.uid___A001_X13e_X16e', '2013.1.01099.S/science_goal.uid___A001_X13e_X164/group.uid___A001_X13e_X165/member.uid___A001_X13e_X166'
'2013.1.00020.S/science_goal.uid___A001_X122_X524/group.uid___A001_X122_X525/member.uid___A001_X122_X526', '2013.1.00020.S/science_goal.uid___A001_X122_X530/group.uid___A001_X122_X531/member.uid___A001_X122_X532', '2013.1.00020.S/science_goal.uid___A001_X122_X518/group.uid___A001_X122_X519/member.uid___A001_X122_X51a', '2013.1.00020.S/science_goal.uid___A001_X122_X520/group.uid___A001_X122_X521/member.uid___A001_X122_X522', '2013.1.00020.S/science_goal.uid___A001_X122_X52c/group.uid___A001_X122_X52d/member.uid___A001_X122_X52e', '2013.1.00020.S/science_goal.uid___A001_X122_X51c/group.uid___A001_X122_X51d/member.uid___A001_X122_X51e', '2013.1.00020.S/science_goal.uid___A001_X122_X528/group.uid__