diff options
Diffstat (limited to 'noao/onedspec/smw/smwswattrs.x')
| -rw-r--r-- | noao/onedspec/smw/smwswattrs.x | 162 |
1 files changed, 162 insertions, 0 deletions
diff --git a/noao/onedspec/smw/smwswattrs.x b/noao/onedspec/smw/smwswattrs.x new file mode 100644 index 00000000..ff859cfc --- /dev/null +++ b/noao/onedspec/smw/smwswattrs.x @@ -0,0 +1,162 @@ +include <error.h> +include <smw.h> + + +# SMW_SWATTRS -- Set spectrum attribute parameters +# This routine has the feature that if the coordinate system of a single +# spectrum in an EQUISPEC WCS is changed then the image WCS is changed +# to a MULTISPEC WCS. + +procedure smw_swattrs (smw, index1, index2, ap, beam, dtype, w1, dw, nw, z, + aplow, aphigh, coeff) + +pointer smw # SMW pointer +int index1 # Spectrum index +int index2 # Spectrum index +int ap # Aperture number +int beam # Beam number +int dtype # Dispersion type +double w1 # Starting coordinate +double dw # Coordinate interval +int nw # Number of valid pixels +double z # Redshift factor +real aplow[2], aphigh[2] # Aperture limits +char coeff[ARB] # Nonlinear coeff string + +bool fp_equald() +int i, j, sz_val, strlen() +double a, b +pointer sp, str, val, mw +errchk smw_mw + +define start_ 10 + +begin + + call smark (sp) + call salloc (str, SZ_LINE, TY_CHAR) + +start_ + switch (SMW_FORMAT(smw)) { + case SMW_ND: + if (!IS_INDEFI(SMW_DTYPE(smw)) && (!fp_equald(w1,SMW_W1(smw)) || + !fp_equald(dw,SMW_DW(smw)) || !fp_equald(z,SMW_Z(smw)))) { + call malloc (val, 15, TY_DOUBLE) + mw = SMW_MW(smw,0) + i = SMW_PDIM(smw) + j = SMW_PAXIS(smw,1) + call mw_gwtermd (mw, Memd[val], Memd[val+i], Memd[val+2*i], i) + Memd[val+j-1] = 1. + switch (dtype) { + case DCNO, DCLINEAR: + a = w1 / (1 + z) + b = dw / (1 + z) + case DCLOG: + a = log10 (w1 / (1 + z)) + b = log10 ((w1 + (nw - 1) * dw) / w1) / (nw - 1) + case DCFUNC: + call error (1, + "Nonlinear functions not allowed for NSPEC format") + } + Memd[val+i+j-1] = a + Memd[val+2*i+(i+1)*(j-1)] = b + call mw_swtermd (mw, Memd[val], Memd[val+i], Memd[val+2*i], i) + call mfree (val, TY_DOUBLE) + } + SMW_DTYPE(smw) = dtype + SMW_NW(smw) = nw + SMW_W1(smw) = w1 + SMW_DW(smw) = dw + SMW_Z(smw) = z + + case SMW_ES: + # Check for any changes to the dispersion system. + if (dtype == DCFUNC) { + call smw_esms(smw) + goto start_ + } + if (!IS_INDEFI(SMW_DTYPE(smw)) && (dtype != SMW_DTYPE(smw) || + nw != SMW_NW(smw) || !fp_equald(w1,SMW_W1(smw)) || + !fp_equald(dw,SMW_DW(smw)) || !fp_equald(z,SMW_Z(smw)))) { + if (SMW_NSPEC(smw) > 1 && index1 > 0) { + call smw_esms(smw) + goto start_ + } + if (!fp_equald(w1,SMW_W1(smw)) || !fp_equald(dw,SMW_DW(smw)) || + !fp_equald(z,SMW_Z(smw))) { + call malloc (val, 15, TY_DOUBLE) + mw = SMW_MW(smw,0) + i = SMW_PDIM(smw) + j = SMW_PAXIS(smw,1) + call mw_gwtermd (mw, Memd[val], Memd[val+i], + Memd[val+2*i], i) + Memd[val+j-1] = 1. + switch (dtype) { + case DCNO, DCLINEAR: + a = w1 / (1 + z) + b = dw / (1 + z) + case DCLOG: + a = log10 (w1 / (1 + z)) + b = log10 ((w1 + (nw - 1) * dw) / w1) / (nw - 1) + } + Memd[val+i+j-1] = a + Memd[val+2*i+(i+1)*(j-1)] = b + call mw_swtermd (mw, Memd[val], Memd[val+i], + Memd[val+2*i], i) + call mfree (val, TY_DOUBLE) + } + } + + SMW_DTYPE(smw) = dtype + SMW_NW(smw) = nw + SMW_W1(smw) = w1 + SMW_DW(smw) = dw + SMW_Z(smw) = z + + if (index1 > 0) { + Memi[SMW_APS(smw)+index1-1] = ap + Memi[SMW_BEAMS(smw)+index1-1] = beam + Memr[SMW_APLOW(smw)+2*index1-2] = aplow[1] + Memr[SMW_APHIGH(smw)+2*index1-2] = aphigh[1] + Memr[SMW_APLOW(smw)+2*index1-1] = aplow[2] + Memr[SMW_APHIGH(smw)+2*index1-1] = aphigh[2] + } + + case SMW_MS: + # We can't use SPRINTF for the whole string because it can only + # handle a limited length and trucates long coefficient strings. + # Use STRCAT instead. + + call smw_mw (smw, index1, index2, mw, i, j) + sz_val = strlen (coeff) + SZ_LINE + call salloc (val, sz_val, TY_CHAR) + call sprintf (Memc[str], SZ_LINE, "spec%d") + call pargi (i) + call sprintf (Memc[val], sz_val, + "%d %d %d %.14g %.14g %d %.14g %.2f %.2f") + call pargi (ap) + call pargi (beam) + call pargi (dtype) + if (dtype == DCLOG) { + call pargd (log10 (w1)) + call pargd (log10 ((w1+(nw-1)*dw)/w1)/(nw-1)) + } else { + call pargd (w1) + call pargd (dw) + } + call pargi (nw) + call pargd (z) + call pargr (aplow[1]) + call pargr (aphigh[1]) + if (coeff[1] != EOS) { + call strcat (" ", Memc[val], sz_val) + call strcat (coeff, Memc[val], sz_val) + } + call mw_swattrs (mw, 2, Memc[str], Memc[val]) + + if (SMW_APS(smw) != NULL) + Memi[SMW_APS(smw)+index1-1] = ap + } + + call sfree (sp) +end |