Initial commit

1 files changed, 70 insertions, 0 deletions

diff --git a/noao/onedspec/smw/smwopen.x b/noao/onedspec/smw/smwopen.x
new file mode 100644
index 00000000..782c8749
--- /dev/null
+++ b/noao/onedspec/smw/smwopen.x
@@ -0,0 +1,70 @@
+include	<smw.h>
+
+
+# SMW_OPEN -- Open SMW structure.
+# The basic MWCS pointer and a template SMW pointer or image is input
+# and the SMW pointer is returned in its place.
+
+procedure smw_open (mw, smw1, im)
+
+pointer	mw		#U MWCS pointer input and SMW pointer output
+pointer	smw1		#I Template SMW pointer
+pointer	im		#I Template IMIO pointer
+
+int	i, nspec, nmw, format, strdic()
+pointer	sp, sys, smw, mw_sctran(), mw_newcopy()
+errchk	smw_daxis, smw_saxes, mw_sctran
+
+begin
+	call smark (sp)
+	call salloc (sys, SZ_FNAME, TY_CHAR)
+
+	call mw_gwattrs (mw, 0, "system", Memc[sys], SZ_FNAME)
+	format = strdic (Memc[sys], Memc[sys], SZ_FNAME, SMW_FORMATS)
+
+	call calloc (smw, SMW_LEN(1), TY_STRUCT)
+	call malloc (SMW_APID(smw), SZ_LINE, TY_CHAR)
+	SMW_FORMAT(smw) = format
+	SMW_DTYPE(smw) = INDEFI
+	SMW_NMW(smw) = 1
+	SMW_MW(smw,0) = mw
+
+	switch (format) {
+	case SMW_ND:
+	    call smw_daxis (smw, im, INDEFI, INDEFI, INDEFI)
+	    call smw_saxes (smw, smw1, im)
+
+	case SMW_ES:
+	    call smw_saxes (smw, smw1, im)
+
+	    nspec = SMW_NSPEC(smw)
+	    call calloc (SMW_APS(smw), nspec, TY_INT)
+	    call calloc (SMW_BEAMS(smw), nspec, TY_INT)
+	    call calloc (SMW_APLOW(smw), 2*nspec, TY_REAL)
+	    call calloc (SMW_APHIGH(smw), 2*nspec, TY_REAL)
+	    call calloc (SMW_APIDS(smw), nspec, TY_POINTER)
+	    if (SMW_PDIM(smw) > 1)
+		SMW_CTLP(smw) = mw_sctran (mw, "logical", "physical", 2)
+
+	case SMW_MS:
+	    call smw_saxes (smw, smw1, im)
+
+	    nspec = SMW_NSPEC(smw)
+	    call calloc (SMW_APIDS(smw), nspec, TY_POINTER)
+	    if (SMW_PDIM(smw) > 1)
+		SMW_CTLP(smw) = mw_sctran (mw, "logical", "physical", 2)
+
+	    nmw = 1 + (nspec - 1) / SMW_NSPLIT
+	    if (nmw > 1) {
+		call realloc (smw, SMW_LEN(nmw), TY_STRUCT)
+		call calloc (SMW_APS(smw), nspec, TY_INT)
+	    }
+	    do i = 1, nmw-1
+		SMW_MW(smw,i) = mw_newcopy (mw)
+	    SMW_NMW(smw) = nmw
+	}
+
+	mw = smw
+	    
+	call sfree (sp)
+end