Attachment 'NIKA2_Time_Estimator4.py'

Download

   1 #!/usr/bin/env python
   2 ##############################################################################################
   3 ### Created by: P. Garcia                                                                  ###
   4 ### Revised by: A. Ritacco                                                                 ###
   5 ### Revised by: B. Ladjelate                                                               ###
   6 ### Title: NIKA2_Time_Estimator.py                                                         ###
   7 ### Purpose: Calculation of the total integration time for observing proposals             ###
   8 ### Creation Date: 2016.JULY.28                                                            ###
   9 ### Date of last revision: 2017.JAN.29                                                     ###
  10 ### Revision history:                                                                      ###
  11 ###   - 2016.JULY.29: comments from C. Kramer & A. Sievers implemented                     ###
  12 ###   - 2016.AUG.03 : changed formula to N. Billot expression                              ###
  13 ###   - 2017.JAN.29 : - fpix set to 0.75 (commisi. results)                                ###
  14 ###                   - NEFD0 @ 2 mm set to 15 mJy/sqrt(hz) (commisi. results)             ###
  15 ###                   - overhead set as used defined parameter between 1.5 - 2.0           ###
  16 ###   - 2017.JULY.28: - According to the NIKA2 technical paper submitted on July, 3rd 2017 ###
  17 ###                   - fpix @1mm, 2mm set to 0.84 and 0.90, respectively                  ###
  18 ###                   - NEFD0 @1mm, 2mm set to 20 and 6 mJy/sqrt(hz), respectively         ###
  19 ###                   - FWHM set to 11.2 @1mm and 17.7 @2mm                                ###
  20 ###   - 2017.SEPT.28  - An error in the data software led sensitivities too optimistics    ###
  21 ###                   - NEFD0 @1mm, 2mm set to 40 and 10 mJy/sqrt(hz), respectively        ###
  22 ###   - 2018.JAN.20:  - NEFD0 @1mm, 2mm set to 33 and 8 mJy/sqrt(hz), respectively         ###
  23 ##############################################################################################
  24 import numpy as np, sys
  25 import math
  26 import sys
  27 import os
  28 os.system("rm -f *.py~")
  29 
  30 #version = 'v 2017.JAN.29'
  31 #version = 'v 2017.JULY.28'
  32 version = 'v 2017.SEPT.28'
  33 
  34 
  35 class bcolors:
  36     HEADER = '\033[95m'
  37     OKBLUE = '\033[94m'
  38     OKGREEN = '\033[92m'
  39     WARNING = '\033[93m'
  40     FAIL = '\033[91m'
  41     ENDC = '\033[0m'
  42     BOLD = '\033[1m'
  43     UNDERLINE = '\033[4m'
  44 
  45 print ''
  46 print bcolors.UNDERLINE + bcolors.HEADER + bcolors.BOLD + 'Time-estimator-NIKA2 '+version+' '+bcolors.ENDC
  47 print ''
  48 print bcolors.HEADER + bcolors.BOLD + ' See the "Guidelines for observing time estimates with the NIKA2 continuum camera' + bcolors.ENDC
  49 print bcolors.HEADER + bcolors.BOLD + ' at the IRAM-30m Telescope" for details on used parameters and calculations.\n' + bcolors.ENDC
  50 
  51 ################################
  52 ##### CONVERSION FACTORS #######
  53 ################################
  54 
  55 ##########  K_cmb ############
  56 def b_v_cmb(freq,tcmb):
  57     h    = 6.626070040*10**(-34) # (J.s) Planck constant
  58     k    = 1.38064852*10**(-23)  # (J/K) Boltzmann constant
  59     c    = 299792458             # (m/s^2) light speed
  60     freq = freq*10**9            # Hz    
  61     bv = 1.0/(((2*h*freq**3)/((c**2)*(np.e**((h*freq)/(k*tcmb)) - 1)))*(np.e**((h*freq)/(k*tcmb))/(np.e**((h*freq)/(k*tcmb)) - 1))*(h*freq/(k*tcmb**2))) / 10**20
  62     return bv
  63 
  64 ##########  Ysz ##############
  65 def y_sz(freq,tcmb):
  66     h    = 6.626070040*10**(-34) # (J.s) Planck constant
  67     k    = 1.38064852*10**(-23)  # (J/K) Boltzmann constant
  68     c    = 299792458             # (m/s^2) light speed
  69     freq = freq*10**9            # Hz    
  70     Ysz  = 1.0
  71     bv   = 1.0/(((2*h*freq**3)/((c**2)*(np.e**((h*freq)/(k*tcmb)) - 1)))*(np.e**((h*freq)/(k*tcmb))/(np.e**((h*freq)/(k*tcmb)) - 1))*(h*freq/(k*tcmb**2)))
  72     bv   = 1.0/bv
  73     y_sz = (bv*tcmb)*(((h*freq/(k*tcmb))*((np.e**((h*freq)/(k*tcmb)) + 1)/(np.e**((h*freq)/(k*tcmb)) -1))) -4)*Ysz /10**(-20)
  74     y_sz = 1.0/y_sz
  75     return y_sz
  76 
  77 #%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  78 #%%%%%%%%%%%%%%%%%%%%% FIXED GENERAL PARAMETERS %%%%%%%%%%%%%%%%%%%%%%  
  79 #%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  80 
  81 ################################################
  82 ##### BEAM SIZES AND CONVERSION TO RADIANS #####
  83 ################################################
  84 res_band1_spec        = 11.2  #From NIKA2 1st paper #MEMO # 12.0  # [arcsec]
  85 res_band2_spec        = 17.7  #From NIKA2 1st paper #MEMO # 20.0  # [arcsec]
  86 teta_band1_spec       = res_band1_spec*math.pi/180.0/3600.0 # in sr units 
  87 teta_band2_spec       = res_band2_spec*math.pi/180.0/3600.0 # in sr units 
  88 area_teta_band1_spec  = ((teta_band1_spec/2.0)**2)*math.pi  # in sr units 
  89 area_teta_band2_spec  = ((teta_band2_spec/2.0)**2)*math.pi  # in sr units 
  90 ############################
  91 ##### OTHER PARAMETERS #####
  92 ############################
  93 #fpix_spec             = 0.75
  94 fpix_spec_1mm         = 0.84
  95 fpix_spec_2mm         = 0.90
  96 bv_1                  = 0.075   # band1, do not change
  97 cv_1                  = 0.001   # band1, do not change
  98 bv_2                  = 0.025   # band2, do not change
  99 cv_2                  = 0.001   # band2, do not change
 100 tiempo                = 1.0     # [hours]
 101 FoV                   = 6.5     # [arcmin] diameter
 102 FoVArea               = ((FoV/2.0)**2)*math.pi # [arcmin^2]
 103 narrays1              = 2
 104 narrays2              = 1
 105 verbose_flag          = 0
 106 #foverhead             = 2.0 # value fixed for all proposals in 2016
 107 ##########################################
 108 ##### NEFDo VALUES FROM OBSERVATIONS #####
 109 ##########################################
 110 NEFD_band1_spec       = 33 #Changed in Jan. 20 #From NIKA2 1st paper arXiv:1707.00908  #MEMO 
 111 NEFD_band2_spec       = 8  #Changed in Jan. 20 #From NIKA2 1st paper arXiv:1707.00908  #comissioning results 2017
 112 
 113 #######################################
 114 ##### SAMPLING ANGULAR VELOCITIES #####
 115 #######################################
 116 sampl_vel_fast = 48.0    # [arcsec/second] 
 117 sampl_vel_mid  = 24.0    # [arcsec/second] 
 118 sampl_vel_slow = 12.0    # [arcsec/second] 
 119 sampl_rate     = 23.84   # [hz]
 120 sampl_rate_pol = 2*23.84 # [hz]
 121 slewing_loss   = 1.0     # to account for slewing of single OTF line
 122 ##########################################
 123 ##########################################
 124 ##########################################
 125 
 126 rango = len(sys.argv)
 127 for i in range(rango):
 128 
 129     if sys.argv[i] == ("--help") or rango == 1:
 130 
 131         print ' USAGE:\n'
 132         print '   python NIKA2_Time_Estimator.py --help '
 133         print '   python NIKA2_Time_Estimator.py --band 1 --rms 2.00 --pwv 2 --elevation 40 --Xsize 6.5 --Ysize 6.5 --filter 1.0 --overhead 1.5'
 134         print '   python NIKA2_Time_Estimator.py --band 2 --rms 1.00 --pwv 4 --elevation 50 --Xsize 15  --Ysize 15  --filter 2.0 --overhead 2.0 --verbose '
 135         print '\n'
 136         print ' OPTIONS:\n'
 137         print '          help =>  This help.'
 138         print '          band =>  Set 1 or 2 for the 1 mm or the 2 mm bands, respectively.'
 139         print '           rms =>  Wanted flux density per beam. Any value above the confusion limit in [mJy/beam].'
 140         print '           pwv =>  Precipitable water vapor in [mm].'
 141         print '     elevation =>  Values from 15 to < 83 [deg].'
 142         print '   Xsize Ysize =>  Map lengths. Xsize and Ysize are in [arcmin]. Minimum map size is 6.5x6.5 [arcmin^2] for homogeneous RMS noise distribution.'
 143         print '        filter =>  Factor for post-processing noise filtering. Values from 1.0 (point-like source) to 2.0 (extended bright emission).'
 144         print '      overhead =>  Factor for telescope overheads between 1.5 and 2.0 (slewing, pointing, focusing, calibration), i.e. all telescope time which is not spend integrating on-source.'
 145         print '       verbose =>  Set to get list of parameters used in the calculations, RMS noise unit conversion, and allowed OTF scan speeds.\n'
 146         sys.exit()
 147 
 148     if sys.argv[i] == ("--band"):
 149         band  = int(sys.argv[i+1]) 
 150         if band == 1:
 151             narrays = narrays1
 152             mili    = 1.2
 153         if band == 2:
 154             narrays = narrays2
 155             mili    = 2.0
 156         if band < 1 or band > 2:
 157             print '    %%%%%%%%%%%%%%%%%%%%%%%'
 158             print bcolors.BOLD +  bcolors.FAIL +'    Band %2i is not defined.' % (band) + bcolors.ENDC
 159             print '    %%%%%%%%%%%%%%%%%%%%%%%'
 160             print ''
 161             sys.exit()
 162 
 163     if sys.argv[i] == ("--rms"):
 164         rms = float(sys.argv[i+1])
 165 
 166     if sys.argv[i] == ("--pwv"):
 167         pwv = float(sys.argv[i+1])
 168 
 169     if sys.argv[i] == ("--elevation"):
 170         elevation = float(sys.argv[i+1])
 171         if elevation < 15 or elevation > 83:
 172             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 173             print bcolors.BOLD +  bcolors.FAIL +'    Elevation %3i [degrees] is outside the telescope`s limits: 15 - 83 [degrees].' % (elevation) + bcolors.ENDC
 174             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 175             print ''
 176             sys.exit()
 177 
 178     if sys.argv[i] == ("--Xsize"):
 179         dx = float(sys.argv[i+1])
 180 
 181     if sys.argv[i] == ("--Ysize"):
 182         dy = float(sys.argv[i+1])
 183 
 184     if sys.argv[i] == ("--filter"):
 185         filtering = float(sys.argv[i+1])
 186         if filtering < 1 or filtering > 2:
 187             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 188             print bcolors.BOLD +  bcolors.FAIL +'    Filter value %2.1f is outside the standard limits: 1.0 - 2.0' % (filtering) + bcolors.ENDC
 189             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 190             print ''
 191             sys.exit()
 192 
 193     if sys.argv[i] == ("--overhead"):
 194         foverhead = float(sys.argv[i+1])
 195         if foverhead < 1.5 or foverhead > 2:
 196             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 197             print bcolors.BOLD +  bcolors.FAIL +'    Overhead value %2.1f is outside the standard limits: 1.5 - 2.0' % (foverhead) + bcolors.ENDC
 198             print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 199             print ''
 200             sys.exit()
 201 
 202     if sys.argv[i] == ("--verbose"):
 203         verbose_flag = 1
 204 
 205 area_obs        = dx*dy
 206 
 207 if area_obs < 4.0:
 208     print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 209     print bcolors.BOLD +  bcolors.FAIL +'    Introduced map size is %3.1f [arcmin^2]. Minimum value is 4.0 [arcmin^2]' % (area_obs) + bcolors.ENDC
 210     print '    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 211     print ''
 212     sys.exit()
 213 
 214 ###############################################################################################
 215 ### SIGMA ESTIMATION CALCULATION, EFFECTIVE NEFD FROM PWV GIVEN VALUES AND UNITS CONVERSION ###
 216 ###############################################################################################
 217         
 218 if band == 1: 
 219     tau = bv_1*pwv + cv_1
 220     NEFD_spec_eff      = NEFD_band1_spec*(np.e**(tau/np.sin(elevation*math.pi/180.)))
 221     rms_point_MJy_spec = rms/area_teta_band1_spec/10**9
 222     rms_point_cmb_spec = rms_point_MJy_spec*b_v_cmb(260,2.726)*10**6
 223     rms_point_ys_spec  = rms_point_MJy_spec*y_sz(260,2.726)*10**6
 224     rms_point_ys_spec  = np.abs(rms_point_ys_spec)
 225     #########################################
 226     ### ratio of areas with fraction of valid pixel in the 2mm band ###       
 227     #########################################
 228     factor_area_spec   = (1+(area_obs/(FoVArea*fpix_spec_1mm))) # from N. Billot Doc.
 229     ##########################################################
 230     ### Integration Time calculations for given conditions ###
 231     ##########################################################
 232     t_spec     = ((NEFD_spec_eff*filtering/rms)**2)*factor_area_spec*(foverhead)/3600.0 #hours
 233     s_map_spec = ((dx*dy)/(((NEFD_spec_eff*filtering)**2)*foverhead))*3600.0            # arcmin^2 / hour^-1/ mJy^-2 
 234 if band == 2:
 235     tau = bv_2*pwv + cv_2
 236     NEFD_spec_eff      = NEFD_band2_spec*(np.e**(tau/np.sin(elevation*math.pi/180.)))
 237     rms_point_MJy_spec = rms/area_teta_band2_spec/10**9
 238     rms_point_cmb_spec = rms_point_MJy_spec*b_v_cmb(150,2.726)*10**6   
 239     rms_point_ys_spec  = rms_point_MJy_spec*y_sz(150,2.726)*10**6     
 240     rms_point_ys_spec  = np.abs(rms_point_ys_spec)
 241     #########################################
 242     ### ratio of areas with fraction of valid pixel in the 2mm band ###       
 243     #########################################
 244     factor_area_spec   = (1+(area_obs/(FoVArea*fpix_spec_2mm))) # from N. Billot Doc.
 245     ##########################################################
 246     ### Integration Time calculations for given conditions ###
 247     ##########################################################
 248     t_spec     = ((NEFD_spec_eff*filtering/rms)**2)*factor_area_spec*(foverhead)/3600.0 #hours
 249     s_map_spec = ((dx*dy)/(((NEFD_spec_eff*filtering)**2)*foverhead))*3600.0            # arcmin^2 / hour^-1/ mJy^-2 
 250 
 251 #########################################
 252 #### TIME PER OTF LINE TO REACH TOTAL ###
 253 #########################################
 254 
 255 T_SLOW = slewing_loss*((dx*60.0)/sampl_vel_slow)/60.0  # [minutes] 
 256 T_MID  = slewing_loss*((dx*60.0)/sampl_vel_mid)/60.0   # [minutes] 
 257 T_FAST = slewing_loss*((dx*60.0)/sampl_vel_fast)/60.0  # [minutes] 
 258 
 259 if verbose_flag == 1:
 260 
 261     print ''
 262     print '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 263     print '%%%%%%%%%% General Parameters Used for the Calculations %%%%%%%%%%%'
 264     print '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 265     print ''
 266     print '------------------------------------------------------------------------'
 267     print '|      Parameters      |    Band 1 (1.2 mm)    |    Band 2 (2.0 mm)    |'
 268     print '------------------------------------------------------------------------'
 269     print '|      opacity         |  %6.3fxpwv + %6.3f  |  %6.3fxpwv + %6.3f  |' % (bv_1,cv_1,bv_2,cv_2)
 270     print '------------------------------------------------------------------------'
 271     print '| HPBW [arcsec]        |       %6.1f          |       %6.1f          |' % (res_band1_spec,res_band2_spec)
 272     print '------------------------------------------------------------------------'
 273     print '| NEFDo [mJy.s^0.5]    |       %6.1f          |       %6.1f          |' % (NEFD_band1_spec,NEFD_band2_spec) 
 274     print '------------------------------------------------------------------------'
 275     print '| fpix                 |       %6.2f          |       %6.2f          |' % (fpix_spec_1mm, fpix_spec_2mm)
 276     print '------------------------------------------------------------------------'
 277     if band  == 1:
 278         print '| rms [mJy/beam]       |       %6.2f          |                       |' % (rms)
 279     if band  == 2:
 280         print '| rms [mJy/beam]       |                       |       %6.2f          |' % (rms)
 281     print '------------------------------------------------------------------------'
 282     print '| FoV [arcmin]         |                   %6.1f                      |' % (FoV) 
 283     print '------------------------------------------------------------------------'
 284     print '| h-filtering          |                   %6.2f                      |' % (filtering) 
 285     print '------------------------------------------------------------------------'
 286     print '| h-overhead           |                   %6.2f                      |' % (foverhead) 
 287     print '------------------------------------------------------------------------'
 288     print '| OTF slow [arcsec/s]  |                    %4i                       |' % (sampl_vel_slow)
 289     print '------------------------------------------------------------------------'
 290     print '| OTF mid  [arcsec/s]  |                    %4i                       |' % (sampl_vel_mid)
 291     print '------------------------------------------------------------------------'
 292     print '| OTF fast [arcsec/s]  |                    %4i                       |' % (sampl_vel_fast)
 293     print '------------------------------------------------------------------------'
 294     print '| Dump [hz]            |                   %6.2f                      |' % (sampl_rate)
 295     print '------------------------------------------------------------------------'
 296     print '| Dump POL [hz]        |                   %6.2f                      |' % (sampl_rate_pol)
 297     print '------------------------------------------------------------------------'
 298     print ''
 299 
 300     print '     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 301     print '     %%%%%%%%%%%%%%%%%   Units Conversion   %%%%%%%%%%%%%%%%%'
 302     print '     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'
 303     print ''    
 304     print '     --------------------------------------------------------'
 305     print '     |         (User Defined)             |    Band %1i       |' % (band)
 306     print '     --------------------------------------------------------'
 307     print '     |rms [mJy/beam]                      |  %8.2f       |' % (rms) 
 308     print '     --------------------------------------------------------'
 309     print '     |rms point-like source [MJy/sr]      |     %6.2f      |'  % (rms_point_MJy_spec) 
 310     print '     |rms point-like source [K_CMB]x10^-6 |     %5i       |'  % (rms_point_cmb_spec) 
 311     print '     |*rms point.like source [Ysz]x10^-6  |     %5i       |'  % (rms_point_ys_spec) 
 312     print '     --------------------------------------------------------'
 313     print '     |* for Ysz = 1.0                                       |'
 314     print '     --------------------------------------------------------'
 315     print ''
 316 
 317 
 318 #############################
 319 #### TABLE FINAL RESULTS ####
 320 #############################
 321         
 322 print '     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'    
 323 print '     %%%%%%%%%%%%%%%%%%%%%%   Results  %%%%%%%%%%%%%%%%%%%%%%'    
 324 print '     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'    
 325 print ''    
 326 print '     for: El = %2i (deg), PWV = %1i (mm), tau = %4.2f, rms = %6.2f'  %  (elevation,pwv,tau,rms) 
 327 print '     Area = %4.1f (arcmin^2), Filter = %4.1f, Overhead = %4.1f'  %  (area_obs,filtering,foverhead) 
 328 print '  ---------------------------------------------------------------'
 329 print '  |    (User Defined)      |                          |         |'
 330 print '  ---------------------------------------------------------------'
 331 print '  | Band                   |                          |    %1i    |'  % (band) 
 332 print '  ---------------------------------------------------------------'
 333 if band  == 1:
 334     print '  | Total Integration Time at 1mm |         [hours]          | %6.1f  |'  % (t_spec)
 335     print '  ---------------------------------------------------------------'
 336     print '  | Mapping speed          | [arcmin^-2/hour/mJy^-2 ] | %6.1f  |'  % (s_map_spec) 
 337 if band  == 2:
 338     print '  | Total Integration Time at 2mm |         [hours]          | %6.1f  |'  % (t_spec)
 339     print '  ---------------------------------------------------------------'
 340     print '  | Mapping speed          | [arcmin^-2/hour/mJy^-2 ] | %6.1f  |'  % (s_map_spec) 
 341 print '  ---------------------------------------------------------------'
 342 print ' '
 343 
 344 if verbose_flag == 1:
 345 
 346     print bcolors.WARNING + '  WARNING1: For proposed maps taking longer than 40 minutes, the total' 
 347     print bcolors.WARNING + '  integration time should be split into several maps.\n' + bcolors.ENDC
 348     print ' '
 349     print bcolors.WARNING + '  WARNING2: Note that the confusion noise is not included in this Time' 
 350     print bcolors.WARNING + '  Estimator and it should be discussed in the proposal whether appro-'
 351     print bcolors.WARNING + '  priate. The RMS noise does not decrease indefinitely with integra- '
 352     print bcolors.WARNING + '  tion time but stops at the confusion limit caused by unresolved sour-'
 353     print bcolors.WARNING + '  ces within the beam. The exact threshold at which the RMS noise    '
 354     print bcolors.WARNING + '  reaches the confusion limit will vary with wavelength, beam size,  '
 355     print bcolors.WARNING + '  and the type of astronomical source (Galactic or Extra-Galactic). For' 
 356     print bcolors.WARNING + '  instance, for the GOODS-N field (part of the Deep Field GT Proposal),' 
 357     print bcolors.WARNING + '  a preliminary 1 sigma confusion limit around ~ 0.090 mJy and ~ 0.056' 
 358     print bcolors.WARNING + '  mJy at 1.2 mm and 2.0 mm, respectively, has been estimated (A. Beelen,'
 359     print bcolors.WARNING + '  private communication).\n' + bcolors.ENDC
 360 
 361     print ' '
 362     print '             Time per OTF line for allowed scanning speeds   '
 363     print '        -----------------------------------------------------'
 364     print '        |          | OTF Scan Velocity  | TIME PER OTF LINE |'
 365     print '        -----------------------------------------------------'
 366     print '        |          |     [arcsec/s]     |    [minutes]      |'
 367     print '        -----------------------------------------------------'
 368     print '        | OTF-SLOW |       %6.1f       |     %6.2f        |' % (sampl_vel_slow,T_SLOW)
 369     print '        | OTF-MID  |       %6.1f       |     %6.2f        |' % (sampl_vel_mid,T_MID)
 370     print '        | OTF-FAST |       %6.1f       |     %6.2f        |' % (sampl_vel_fast,T_FAST)
 371     print '        -----------------------------------------------------'
 372     print ' '
 373 
 374 ###############################
 375 #### CONVERSION TO MINUTES ####
 376 ###############################
 377 
 378 if t_spec >= 1.:
 379    time_str_spec=str(round(t_spec,1))+' hours'
 380 
 381 if t_spec < 1. and t_spec >= 1./60.:
 382    time_str_spec=str(round(t_spec*60,1))+' minutes'
 383 
 384 if t_spec < 1./60. and t_spec >= 0.1/3600.:
 385    time_str_spec=str(int(t_spec*3600))+' seconds'
 386 
 387 if t_spec < 0.1/3600.:
 388    time_str_spec='0.1 (seconds)'
 389 
 390 slew_overhe = str(int((slewing_loss-1.0)*100.0))
 391 
 392 ##############################
 393 #### OUTPUT FINAL RESULTS ####
 394 ##############################
 395     
 396 print bcolors.HEADER + '**********************************************'
 397 print '***  Total Integration Time => ' + time_str_spec + ' ***'
 398 print '**********************************************\n'    
 399 print ''
 400 
 401 print bcolors.HEADER +  bcolors.UNDERLINE + 'Please include the following text into your proposal:\n' + bcolors.ENDC
 402 
 403 print 'According to the specifications of the NIKA2 instrument, the total observing time using the NIKA2 '+str(band)+' mm band to map a region of '+str(round(area_obs,1))+' [arcmin^2] to reach an rms of '+str(rms)+' [mJy/beam], assuming '+str(pwv)+' [mm] pwv, '+str(elevation)+' [deg] elevation, Filter = '+str(filtering)+', Overhead = '+str(foverhead)+', was estimated to be *'+time_str_spec+'*, using the time estimator '+version+'.\n'

Attached Files

To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.
  • [get | view] (2017-07-28 12:35:36, 21.8 KB) [[attachment:NIKA2_Time_Estimator.py]]
  • [get | view] (2017-07-28 12:39:29, 21.8 KB) [[attachment:NIKA2_Time_Estimator1.py]]
  • [get | view] (2017-09-28 09:12:07, 22.1 KB) [[attachment:NIKA2_Time_Estimator2.py]]
  • [get | view] (2017-09-28 15:02:49, 22.1 KB) [[attachment:NIKA2_Time_Estimator3.py]]
  • [get | view] (2018-01-20 15:02:28, 22.3 KB) [[attachment:NIKA2_Time_Estimator4.py]]
  • [get | view] (2017-07-28 11:54:53, 276.4 KB) [[attachment:NIKA2_Time_Estimator_Guideline.pdf]]
 All files | Selected Files: delete move to page copy to page

You are not allowed to attach a file to this page.