source: nappy/trunk/naFile2310.py @ 349

#   Copyright (C) 2004 CCLRC & NERC( Natural Environment Research Council ).
#   This software may be distributed under the terms of the
#   Q Public License, version 1.0 or later. http://ndg.nerc.ac.uk/public_docs/QPublic_license.txt

"""
naFile2310.py
=============

Container module for NAFile2310 class.

"""

# Imports from python standard library

# Imports from local package
from textParser import *
import naFile2110

class NAFile2310(naFile2110.NAFile2110):
    """
    Class to read, write and interact with NASA Ames files conforming to the
    File Format Index (FFI) 2310.
    """

    def readHeader(self):
        """
        Reads FFI-specifc header section.
        """        
        self._normalizedX="no"
        self._readCommonHeader()
        self.DX=readItemsFromLine(self.file.readline(), 1, float)
        self.XNAME=readItemsFromLines(self._readLines(self.NIV), self.NIV, str)
        self.XNAME.reverse()  # Reverse because C-type array is least-changing first          
        self._readVariablesHeaderSection()
        self._readAuxVariablesHeaderSection()
        self._readComments()
    
    def _readData1(self, datalines, ivar_count): 
        """
        Reads first line/section of current block of data.
        """        
        # Start with independent and Auxilliary vars
        (x_and_a, rtlines)=readItemsFromUnknownLines(datalines, self.NAUXV+1, float)
        (x, aux)=(x_and_a[0], x_and_a[1:])
        for a in range(self.NAUXV):
            self.A.append(aux[a])
        self.X.append([])
        self.X[ivar_count].append(x)
        # Set up list to take second changing independent variable
        self.X[ivar_count].append([aux[1]])  
        self.NX.append(int(aux[0]))
        self.DX.append(int(aux[2]))
        return rtlines
   
    def _readData2(self, datalines, ivar_count):
        """
        Reads second line/section (if used) of current block of data.
        """
        # Now get the dependent variables
        (v, rtlines)=readItemsFromUnknownLines(datalines, self.NV*self.NX[ivar_count], float)
        count=0
        for n in range(self.NV):
            self.V[n].append([])
            for i in range(self.NX[ivar_count]):
                self.V[n][ivar_count].append(v[count])
                count=count+1
        return rtlines