REDUCAO POLARIMETRIA ===================== Atualizacao: 03/Dezembro/2006 - CVR Este texto refere-se a reducao de dados obtidos com a calcita, que eh o elemento que divide o feixe de um dado objeto em dois. Assim, se a imagem que voce esta reduzindo tem duas imagens para cada objeto, tenha certeza que ela foi obtida com a calcita. Os dados polarimetricos podem ser obtidos com dois tipos de analisadores: lamina de meia-onda (l/2) ou lamina de quarto-de-onda (l/4). No primeiro caso, apenas a polarizacao linear pode ser obtida e o numero minimo de imagens para uma medida da polarizacao eh 4. Com l/4, podemos obter as polarizacoes linear e circular e precisamos de, pelo menos, 8 imagens. DICAS GERAIS ============ * estamos colocando no final deste arquivo, os ¨lpar¨ das tasks usadas * epar ccdred verbose=yes * para ler fita remotamente, use a dica que existe no FAQ do suporte da DAS * criar arquivo com lista de arquivos files *.fits > arq tirar .fits da lista acima * Acrescentando HJD no header de todas as imagem hedit @arq observat lna (aqui se corrigi o header das imagens) setjd @arq (aqui se calcula o hjd e atualiza-se o header) hselect @arq hjd > hjd.lis (cria-se o arquivo hjd.lis com o hjd de cada imagem) -> Confira: 1. se no no setjd as variaveis sao: time = ut exposure = exptime 2. se no cabecalhos das imagens existem os parametros UT e EXPTIME 3. se a data esta no formato: yyyy-mm-dd - se nao estiver, eh necessario usar o hedit para corrigir - a keyword da data eh date-obs * acrescentando extensao (fits, por exemplo) em arquivos rename extn -> ver epar do rename * TRIMSEC e outras informacoes sobre detetores do LNA /home1/claudiavr/iraf/lembretes * ccdlist mostra lista de imagens e flags de reducao REDUCAO BASICA ============== A reducao basica aplica-se para qualquer conjunto de dados. Passos: - criar um arquivo para corrigir do ruido de leitura do CCD: bias ou zero - criar um arquivo para corrigir diferenca de resposta dos pixels do CCD: flatfield - fazer as correcoes acima + ¨trim¨ - deslocar todas as imagens para a mesma referencia: registro Abaixo a execucao de cada um desses passos eh descrita. * Combinar bias Usar task ZEROCOMBINE * Parametros importantes: (combine = "average") Type of combine operation (reject = "avsigclip") Type of rejection (ccdtype = "") CCD image type to combine * Combinar flat flatcombine (vide lpar) * Antes de passar o ccdproc (abaixo) nas imagens, fazer backup * Reducao do trim, bias e flat ccdproc ATENCAO! Se imagens forem com l/4, NAO USAR CCDPROC PARA CORRECAO DE FLATFIELD! Comentarios: NAO usar o overscan! Testes feitos para imagens do LNA mostraram que utilizacao da area de overscan para o calculo/correcao do overscan em separado do bias so piora o resultado. TRIMSEC e outras informacoes sobre detetores do LNA /home1/claudiavr/iraf/lembretes * Correcao do flat-field para l/4 - ccdproc (bias e trim) no flatave; - calcular valor medio do plato (sem vignetting) do flat-field medio; - criar centro do flat (testar com tvmark) - fazer fotometria com aberturas de 20,30,...120; - escolher abertura com valor medio maximo - dividir flat medio pelo valor acima. Desse modo, o plato tera um valor proxima a 1; imarith flatave / [valor medio] flatave_norm - fazer a correcao do flat de cada imagem efetuando A DIVISAO DE CADA IMAGEM PELO FLAT CALCULADO DO MODO ACIMA. imarith @arq / flatave_norm @arq Foram feitos testes corrigindo-se o flat pelo modo acima e tambem acima o ccdproc (usando a normalizacao por toda area da imagem). No caso acima, os erros estao mais proximos do ruido de fotons esperado. * Fazer backup antes do xregister * Registrar imagens xregister Dicas: - parametros do xregister: Considerando deltax e deltay os deslocamentos maximos no conjunto de imagens, - xwindow e ywindow devem ser 2 vezes maior que os deltax e deltay; - regions: com dimensoes 2 vezes maiores que xwindow e ywindow. - os deltax e deltay de uma serie de imagens podem ser estimados escolhendo uma estrela de referencia cujas coordenadas serao verificadas em algumas imagens ao longo da sequencia (isso pode ser feito com a fotometria (opcao ¨a¨) do imexamine). As variacoes maximas nos eixos x e y sao os deltax e deltay. - visualizacao da correcao: A task disp_var (/home1/claudiavr/iraf/tasks) possibilita o display de varias imagens consecutivamente. Pode ser util para verificar se os deslocamentos aplicados pelo xregister foram os corretos. FOTOMETRIA ========== LAMINA DE MEIA-ONDA =================== * anotar FWHM e Sigma (raiz quadrado do ceu) calcular valores com imexam (a: fotometria de abertura - r: perfil radial - ?: help) tome cuidado com estrelas saturadas. Normalmente, um CCD tem contagem maxima em ~32000. Se a imagem possui muitos objetos saturados, pode ser conveniente configurar ¨datamax¨ no datapars. * editar parametros utilizados por diferentes tasks (tarefas) do apphot. Exemplos de tasks: phot, daofind. Exemplos de parametros: FWHM (largura meia altura), constante de magnitude fitskypars + fitskypars: parametros para ajuste de ceu - normalmente nao variamos os valores durante uma missao * annulus: deve ser tal que compreenda as duas imagens de uma mesma estrela (~50) * dannulus: configuramos como 10 * algoritmo: mode datapars * editar FWHM e Sigma de acordo com caracteristicas da imagem photpars * editar aberturas * achar estrelas ATENCAO! Usar sempre daofind, mesmo que seja estrela unica!! Esse cuidado eh necessario, pois o angulo da polarizacao eh dependente da ordem das duas imagens (de uma mesma estrela) no arquivo de coordenadas. Por exemplo, convencione que a imagem de baixo eh a primeira imagem do par a ser listada no arquivo de coordenadas. Siga esse procedimento EM TODO O TURNO! Caso isso nao seja feito, os angulos de polarizacao estarao ERRADOS!! Usando o daofind esse procedimento eh feito automaticamente - daofind -> acha estrelas na imagem Temos dois casos: 1. polarizacao de 1 unico objeto 2. polarizacao de todas as estrelas do campo 1. APENAS UMA ESTRELA DO CAMPO: Esse caso corresponde a estrelas-padroes e/ou interesse em 1 unico objeto usar como limite de deteccao no daofind as contagens do objeto 2. TODAS ESTRELAS DO CAMPO usar limite de deteccao em torno de 4 ou 5 sigmas - ordem (pccdpack) -> encontra os pares de imagens ordinaria e extraordinaria de um mesmo objeto e cria um novo arquivo de coordenada -> usar ordem2 (task da Claudia, nao do pccdpack) - tvmark -> plota centros no display * rodar phot - tasks auxiliares com parametros de entrada que devem ser editados: + photpars: aberturas + fitskypars: parametros para ajuste de ceu - normalmente nao variamos os valores. * annulus: deve ser tal que compreenda as duas imagens de uma mesma estrela (~50) * dannulus: configuramos como 10 * algoritmo: mode + centerpars: normalmente editado uma unica vez na vida * txdump ou cria_dat (este ultimo eh mais simples...) Nesta etapa, selecionados alguns resultados do phot (contidos no .mag.?) e jogamos para um arquivo de entrada do pccdgen. Parametros do txdump: textfiles = "*.mag.1" Input apphot/daophot text database(s) fields = "image,msky,nsky,rapert,sum,area" Fields to be extracted expr = "yes" Boolean expression for record selection (headers = no) Print the field headers ? (parameters = yes) Print the parameters if headers is yes ? (mode = "ql") Mode of task ** cria_dat - criar mag com ¨files *mag.1 > mag¨ criar serie de dats usando, no exemplo abaixo, conjunto de 8 imagens varim = "@mag" Input mag list (outdat = "dat") Raiz do arquivo de saida (interval = 8) Numero de imagens em 1 datfile (flistvar = "tmpvar9215co") (mode = "ql") * pccdgen (vide lpar) Esta task eh a que calcula a polarizacao. Cria um *.log de acordo com o definido nos parametros de entrada usando como executavel (fileexe = "/home1/claudiavr/iraf/pccdpack/pccd/pccd2000gen.exe") pccd execut => procurar a abertura que fornece o menor erro do ajuste. Eh o ¨SIGMA¨ da saida do pccd ou pccdgen * graf Essa tarefa permite visualizar/imprimir a razao (F_ord-F_ext)/(F_ord+F_ext) como funcao da posicao da lamina retardadora, alem do valor da polarizacao para uma dada abertura. * imprimindo o grafico do graf selecione no graf a criacao do metacode (.mc) impressora epar stdplot -> device = stdplot stdplot *.mc criar ps epar stdplot -> device = epsfl (landscape) ou epsf (portrait) stdplot *.mc * nao esqueca de deletar os diretorios bck! ==>> PARA UMA UNICA ESTRELA NO CAMPO NAO EH NECESSARIO REALIZAR AS DUAS ETAPAS ABAIXO! * macrol ou macrol_v (l/4) seleciona abertura de menor erro arquivo de entrada: .log (saida do pccdgen ou pccd) arquivo de saida: *.out. Neste arquivo cada objeto corresponde a uma linha * select (funciona apenas para l/2) visualiza resultado do campo * plotando vetores sobre uma imagem refer vecplot Comentarios/testes: para estrelas brilhantes escolha do annulus e dannulus nao eh fundamental para estrelas fracas escolha do ceu tambem nao eh crucial. Melhor usar annulus que inclua as duas imagens de cada objeto. LAMINA DE QUARTO-DE-ONDA ======================== com os mag criados: * crie um arquivo_lista que contenha os nomes de todos os mags: files *mag.1 > mag * cria_dat criar serie de dats usando, no exemplo abaixo, conjunto de 8 imagens varim = "@mag" Input mag list (outdat = "dat") Raiz do arquivo de saida (interval = 8) Numero de imagens em 1 datfile (flistvar = "tmpvar9215co") (mode = "ql") * acha_zero eu rodo o acha_zero para cada conjunto de dados em que a primeira imagem seja igual a (multiplo de 8)+1 NAO ESQUECA DE ATUALIZAR PARAMETROS DO PCCD!! em particular, use como executavel: /home1/claudiavr/pccd/vccd.exe indat = "bd32_r.dat" Arquivo/lista de entrada (outmac = "acha") Arquivo/lista de saida do macrol (platebeg = 0) Posicao da lamina - inicio (plateend = 90) Posicao da lamina - final (passo = 1) Passo para variar lamina (mode = "ql") O acha_zero roda tambem para uma lista de arquivos. Util para uma sequencia de imagens de candidatas a polar. um programa alternatico eh o zerofind do Antonio Pereyra. Porem, esse programa escolhe a abertura de minimo erro a partir do erro da polarizacao linear. indat = "dat.073" Input file (outmac = "zeroq") Macrol output file (emin = "first") macrol: minimum? (first|full) (platebeg = 0) Initial zero (degrees) (plateend = 90) Final zero (degrees) (passo = 1) step for zero (degrees) (mode = "ql") ==> NAO ESQUECA DE CONFIGURAR pccdgen * definir zero em cima dos resultados acima de toda a missao configure o pccd com esse valor no parametro ¨zero¨ * pccd_var em cada noite (ja cria o out) ==>> ATENCAO! Nao esqueca de editar o pccdgen corretamente! ==>> Leia secao sobre pccd_var no final deste arquivo! * plota_pol eh necessario criar arquivo hjd.txt Atencao: o hjd.txt deve ter o mesmo numero de linhas que datfiles para o plota_pol funcionar corretamente. o plota_pol permite criar um arquivo de saida com hjd,polarizacao que pode ser usado como entrada do diagfase2 * caso queira usar o SELECT com uma saida de l/4, crie um novo *.out com o comando fields [out.original] 3-11 > [out.select] SEQUENCIA TEMPORAL DE IMAGENS ============================= * com os mags criados (saida do phot) rodar cria_dat -> criar dats rodar pccd_var (leia secao sobre pccd_var no final deste arquivo) -> nao esquecer de editar pccdgen corretamente! Em particular, verifique se o numero de imagens (posicoes da lamina esta correto). O pccd_var usa o pccdgen, e nao o pccd, entao os parametros devem estar configurados corretamente no PCCDgen! -> caso a primeira imagem do primeiro dat nao for feita na posicao 1 (das 16 possiveis posicoes da lamina) - para l/2 -, temos que o parametro delta_theta no pccd_var deve ser configurado como: [(#pos. da lamina)-1]*(-45) Usar valores entre -360 e 0. * criar hjd.txt hselect @arq hjd > hjd.lis Atencao: o hjd.txt deve ter o mesmo numero de linhas que datfiles para o plota_pol funcionar corretamente. FOTOMETRIA - L/2 ou l/4 ======================= * depois de criar os *.mag.* com a rotina phot * criar um arquivo .pht com txdump - aqui estamos selecionando os parametros contidos nos *.mag.* para calculo da fotometria txdump > .pht Parametros do txdump: textfiles = "*.mag.1" Input apphot/daophot text database(s) fields = "image,msky,nsky,rapert,sum,area" Fields to be extracted expr = "yes" Boolean expression for record selection (headers = no) Print the field headers ? (parameters = yes) Print the parameters if headers is yes ? (mode = "ql") Mode of task * criar arquivo .luz com a task phot_pol - aqui estamos calculando as curvas de luz propriamente ditas usando como entrada .pht cada estrela corresponde a duas imagens. Assim, um arquivo com 64 imagens corresponde a 32 estrelas cl> lpar phot_pol file_in = "rej1007.pht" Arquivo com saida txdump file_out = "rej1007.luz" Nome archivo de saida (nstars = 32) Numero de estrelas (nhw = 112) Numero de posicoes da lamina (nap = 10) Numero de aberturas (star = 1) Numero da estrela a se fazer (fotometria (star_out = 0) Estrela que nao eh incluida na soma dos fluxos, (ganho = 5.5) Ganho - e/adu (mode = "ql") * criar hjd.txt (caso ele ainda nao tenha sido criado) hselect @arq hjd > hjd.lis Atencao: o hjd.txt deve ter o mesmo numero de linhas que o numero de laminas utilizados no phot_pol * visualizar as curvas de luz obtidas com o plota_luz cl> lpar plota_luz arqpht = "rej1007.luz" Input pht file (tempo = "hjd.lis") File with time input (comp = 1) Number of the comparison star (aper = 1) Ordinal number of the aperture (conecta = yes) Connect the points (pontos = yes) Plot points (title = "REJ 1007-20 - 18/marco/2004 - R") Title of the graphics (phase = no) Convert HJD to orbital phase (to = 0.) To das efemerides (per = 0.0605) Periodo (dias) (ffile = no) Create hjd,mag file? (mmagfile = "v1432aql_24set03.pht") Name of the hjd,mag file (metafile = no) Create mc file (eps = no) Create eps file (flist = "hjd.lis") (mode = "ql") o parametro aper corresponde a ordem da abertura de acordo com o que foi configurado como entrada do phot, nao ao seu valor real. o plota_luz permite criar um arquivo de saida com hjd,mad que pode ser usado como entrada do diagfase2 * pdm k -> faz o ajuste p -> plota binado no periodo h -> plota em hjd ****************************************************************** SOBRE AS TASKS ============== PCCD_VAR -------- O pccd_var usa o pccdgen para calcular a polarizacao de uma serie de datfiles. Assim, eh fundamental editar o pccdgen corretamente: numero de laminas, tipo de lamina, erro de leitura, etc. O pccd_var calcula a polarizacao de uma serie de datfiles (dat.001, dat.002, etc) e produz um log para cada dat (log.001, log.002,...). Ao fazer a fotometria, usamos diferentes aberturas. Assim, cada datfile contem o resultado de varias aberturas (tipicamente 10) e portanto cada log tem a polarizacao de cada objeto calculada usando 10 diferentes aberturas. Mas, qual abertura usaremos como a melhor estimativa do valor da polarizacao? Escolhemos a abertura que apresenta o menor erro do ajuste. Essa variavel eh representada por ¨SIGMA¨ nos arquivos de resultado (.log). O pccd_var usa uma rotina chamada ¨macrol¨ para escolher a abertura de menor erro dentro de cada log e cria um log final onde cada linha corresponde a melhor abertura de cada datfile. Esse ultimo eh o log.out. A rotina macrol tem uma opcao (minimum) que pode ser full ou first. No caso de estar configurada como first, o primeiro minimo eh escolhido. No caso full, o minimo absoluto eh escolhido. Veja o exemplo abaixo: | S | * * I | * * * G | M | * * * A | * | * | |-------------------------------- 2 3 4 5 6 7 8 9 10 11 Se o macrol estiver com minimum=first, a abertura escolhida sera a 3. Porem, se o minimum = full, a abertura escolhida sera a 10. Eu, Claudia, costumo usar o macrol configurado na opcao full. ****************************************************************** ****************************************************************** PARAMETROS DAS TASKS UTILIZADAS =============================== *** PACKAGE = ccdred TASK = zerocombine input = @arq List of zero level images to combine (output = biasave) Output zero level name (combine= average) Type of combine operation (reject = avsigclip) Type of rejection (ccdtype= ) CCD image type to combine (process= no) Process images before combining? (delete = no) Delete input images after combining? (clobber= no) Clobber existing output image? (scale = none) Image scaling (statsec= ) Image section for computing statistics (nlow = 0) minmax: Number of low pixels to reject (nhigh = 1) minmax: Number of high pixels to reject (nkeep = 1) Minimum to keep (pos) or maximum to reject (neg) (mclip = yes) Use median in sigma clipping algorithms? (lsigma = 3.) Lower sigma clipping factor (hsigma = 3.) Upper sigma clipping factor (rdnoise= 0.) ccdclip: CCD readout noise (electrons) (gain = 1.) ccdclip: CCD gain (electrons/DN) (snoise = 0.) ccdclip: Sensitivity noise (fraction) (pclip = -0.5) pclip: Percentile clipping parameter (blank = 0.) Value if there are no pixels (mode = ql) ***** PACKAGE = ccdred TASK = flatcombine input = @arq List of flat field images to combine (output = flatave) Output flat field root name (combine= average) Type of combine operation (reject = avsigclip) Type of rejection (ccdtype= ) CCD image type to combine (process= no) Process images before combining? (subsets= no) Combine images by subset parameter? (delete = no) Delete input images after combining? (clobber= no) Clobber existing output image? (scale = mode) Image scaling (statsec= ) Image section for computing statistics (nlow = 1) minmax: Number of low pixels to reject (nhigh = 1) minmax: Number of high pixels to reject (nkeep = 1) Minimum to keep (pos) or maximum to reject (neg) (mclip = yes) Use median in sigma clipping algorithms? (lsigma = 3.) Lower sigma clipping factor (hsigma = 3.) Upper sigma clipping factor (rdnoise= 0.) ccdclip: CCD readout noise (electrons) (gain = 1.) ccdclip: CCD gain (electrons/DN) (snoise = 0.) ccdclip: Sensitivity noise (fraction) (pclip = -0.5) pclip: Percentile clipping parameter (blank = 1.) Value if there are no pixels (mode = ql) ***** cl> lpar ccdproc images = "@arq" List of CCD images to correct (output = "") List of output CCD images (ccdtype = " ") CCD image type to correct (max_cache = 0) Maximum image caching memory (in Mbytes) (noproc = no) List processing steps only?\n (fixpix = no) Fix bad CCD lines and columns? (overscan = no) Apply overscan strip correction? (trim = yes) Trim the image? (zerocor = yes) Apply zero level correction? (darkcor = no) Apply dark count correction? (flatcor = yes) Apply flat field correction? (illumcor = no) Apply illumination correction? (fringecor = no) Apply fringe correction? (readcor = no) Convert zero level image to readout correction? (scancor = no) Convert flat field image to scan correction?\n (readaxis = "line") Read out axis (column|line) (fixfile = "") File describing the bad lines and columns (biassec = " ") Overscan strip image section (trimsec = "[10:521,1:512]") Trim data section (zero = "../../bias/biasave") Zero level calibration image (dark = "") Dark count calibration image (flat = "../../flat/v/flatvave") Flat field images (illum = "") Illumination correction images (fringe = "") Fringe correction images (minreplace = 1.) Minimum flat field value (scantype = "shortscan") Scan type (shortscan|longscan) (nscan = 1) Number of short scan lines\n (interactive = no) Fit overscan interactively? (function = "legendre") Fitting function (order = 1) Number of polynomial terms or spline pieces (sample = "*") Sample points to fit (naverage = 1) Number of sample points to combine (niterate = 1) Number of rejection iterations (low_reject = 3.) Low sigma rejection factor (high_reject = 3.) High sigma rejection factor (grow = 0.) Rejection growing radius (mode = "ql") cl> **************** PCCDGEN NO MODO l/2 cl> lpar pccdgen filename = "dat.001" input file (.dat) (nstars = 1) number of stars (max. 2000) (wavetype = "half") wave-plate used ? (half,quarter,other) (retar = 180.) retardance of waveplate (degrees) (nhw = 16) number of total wave-plate positions in input f (pospars = "") wave-plate positions used to calculus? :e (nap = 10) number of apertures (max. 10) (calc = "c") analyser: calcite (c) / polaroid (p) (readnoise = 0.82) CCD readnoise (adu) (ganho = 5.) CCD gain (e/adu) (deltatheta = 0.) correction in polarization angle (degrees) (zero = 50.) Zero of waveplate (fileout = "teste.log") output file (.log) (fileexe = "/home1/claudiavr/iraf/pccdpack/pccd/pccd2000gen.exe") pccd execut (norm = yes) include normalization? (flist = "pccdgen8522ij") (line1 = "") (line2 = "") (mode = "al") cl> ***** PCCDGEN NO MODO l/4 cl> lpar pccdgen filename = "dat.001" input file (.dat) (nstars = 1) number of stars (max. 2000) (wavetype = "quarter") wave-plate used ? (half,quarter,other) (retar = 90.) retardance of waveplate (degrees) (nhw = 16) number of total wave-plate positions in input f (pospars = "") wave-plate positions used to calculus? :e (nap = 10) number of apertures (max. 10) (calc = "c") analyser: calcite (c) / polaroid (p) (readnoise = 0.82) CCD readnoise (adu) (ganho = 5.) CCD gain (e/adu) (deltatheta = 0.) correction in polarization angle (degrees) (zero = 50.) Zero of waveplate (fileout = "teste.log") output file (.log) (fileexe = "/home1/claudiavr/iraf/pccdpack/pccd/pccd2000gen.exe") pccd execut (norm = yes) include normalization? (flist = "pccdgen8522ij") (line1 = "") (line2 = "") (mode = "al") cl> ***** PCCD para l/2 cl> lpar pccd filename = input file (.dat) (nstars = 1) number of stars (maximum 2000) !! numero de estrelas (nhw = 16) number of postions of wave-plate (maximum 16) !! numero de posicoes da lamina (nap = 10) number of apertures (maximum 10) !! numero de aberturas, normalmente 10 (calc = "c") analyser: calcite (c) / polaroid (p) !! deixe c de calcita (readnoise = 5.) CCD readnoise (adu) !! coloque o valor adequado para o CCD usado, veja /home1/claudiavr/iraf/lembretes (ganho = 0.82) CCD gain (e/adu) !! idem (deltatheta = 0.) correction in polarization angle (degrees) !! configurar zero (zero = 0.) Position zero of the l/4 plate !! configurar zero (fileout = "teste") output file (.log) !! nome do arquivo de saida, sera apendado ¨.log¨ (fileexe = "/home1/claudiavr/pccd/pccd.exe") pccd execute file (.exe) (mode = "ql") cl> ***** PCCD para l/4 cl> lpar pccd filename = input file (.dat) (nstars = 1) number of stars (maximum 2000) !! numero de estrelas (nhw = 16) number of postions of wave-plate (maximum 16) !! numero de posicoes da lamina (nap = 10) number of apertures (maximum 10) !! numero de aberturas, normalmente 10 (calc = "c") analyser: calcite (c) / polaroid (p) !! deixe c de calcita (readnoise = 5.) CCD readnoise (adu) !! coloque o valor adequado para o CCD usado, veja /home1/claudiavr/iraf/lembretes (ganho = 0.82) CCD gain (e/adu) !! idem (deltatheta = 0.) correction in polarization angle (degrees) !! configurar zero (zero = 0.) Position zero of the l/4 plate !! configurar zero (fileout = "teste") output file (.log) !! nome do arquivo de saida, sera apendado ¨.log¨ (fileexe = "/home1/claudiavr/pccd/vccd.exe") pccd execute file (.exe) (mode = "ql") cl> ***** TASK = fitskypars (salgori= mode) Sky fitting algorithm (annulus= 50.) Inner radius of sky annulus in scale units (dannulu= 10.) Width of sky annulus in scale units (skyvalu= 0.) User sky value (smaxite= 10) Maximum number of sky fitting iterations (sloclip= 0.) Lower clipping factor in percent (shiclip= 0.) Upper clipping factor in percent (snrejec= 50) Maximum number of sky fitting rejection iteratio (sloreje= 3.) Lower K-sigma rejection limit in sky sigma (shireje= 3.) Upper K-sigma rejection limit in sky sigma (khist = 3.) Half width of histogram in sky sigma (binsize= 0.1) Binsize of histogram in sky sigma (smooth = no) Boxcar smooth the histogram (rgrow = 0.) Region growing radius in scale units (mksky = no) Mark sky annuli on the display (mode = ql) ***** TASK = macrol (file_in = "imagem_4.log") pccd output file (.log) or list(@*) (file_out = "imagem_4.out") output file (.out) (minimun = "full") minimum? (first,full) (flist = "") (flistvar = "") (mode = "ql") *****