! Procedure onoff.nic calculates on-off flux and RA/Dec position of each ! bolometer in Counts/beam ! The high frequency skynoise is suppressed, using "Toon Snijders method" ! The called procedures are : flux.nic radec.nic wmean.nic and fitline.nic ! ! @ p_signal &1 define real flux_onf[nchan] flux_rms[nchan] /global define real Delta_RA[nchan] Delta_DEC[nchan] /global define real ra_s dec_s define character unity*3 dummy*3 let dummy "r&2" ! if (dummy.eq."r") then despike /iter 3 /interp say "! Despiking requested..." else say "! Despiking omitted..." endif set gain_channel set gain_amplifier say "! Calculation of the skynoise corrected flux of each channel [Counts/Jy]" say "! This will take about 1 minute..." var * if scan_type.eq.1 then let unity "Counts/beam" ! let unity "Jy" @ flux flux_onf flux_rms @ radec Delta_RA Delta_DEC ! if Flux_onf[1].le.2 then ! let flux_onf Flux_onf*1000 ! let flux_rms Flux_rms*1000 ! let unity "mJy" ! endif sic precision single ! let flux_onf INT(Flux_onf*1000)|1000 ! let flux_rms INT(Flux_rms*1000)|1000 let Delta_RA INT(Delta_RA*100)|100 let Delta_DEC INT(Delta_DEC*100)|100 let ra_s (LAMBDA*1000)|1000 let dec_s (BETA*1000)|1000 SAY "!Date="'DATE_OBSERV'" LST(sec)="'LST[1]|3600'" Chop_throw(arcsec)="'wobbler_throw' say "!SCAN:"'scan'" SOURCE:"'SOURCE'" RA(Deg):"'ra_s'" DEC(Deg):"'dec_s'" Bol_gain="'ref_gain' ! say "!Opacity="'scan_opacity'" Bol_gain="'ref_gain'" Counts/Jy="'counts_per_jy'" Gain_elv=" 'gaine_flag'" say "!Scan Channel Off(RA) Off(DEC) Flux["'unity'"] Error["'unity'"]" say " "'scan'" "1" "'Delta_RA[1]'" "'Delta_DEC[1]'" "'Flux_onf[1]'" "'Flux_rms[1]' for k 2 to nchan say " "'scan'" "'k'" "'Delta_RA[k]'" "'Delta_DEC[k]'" "'Flux_onf[k]'" "'Flux_rms[k]' next ! endif else say "! No meaning to process, this scan is NOT an on-off" endif del /var flux_onf flux_rms Delta_RA Delta_DEC ra_s dec_s