コード例 #1
0
   def copyCommonData(self):

      SLFn = SELAFIN('')
      #   Meta data
      SLFn.TITLE = self.slf.TITLE
      SLFn.file = self.slf.file
      SLFn.IPARAM = self.slf.IPARAM
      #   Time
      SLFn.DATETIME = self.slf.DATETIME
      SLFn.tags = self.slf.tags
      #   Variables
      SLFn.NBV1 = self.slf.NBV1
      SLFn.VARNAMES = self.slf.VARNAMES
      SLFn.VARUNITS = self.slf.VARUNITS
      SLFn.NBV2 = self.slf.NBV2
      SLFn.CLDNAMES = self.slf.CLDNAMES
      SLFn.CLDUNITS = self.slf.CLDUNITS
      SLFn.NVAR = self.slf.NVAR
      SLFn.VARINDEX = range(self.slf.NVAR)
      #   Unchanged numbers
      SLFn.NPLAN = self.slf.NPLAN
      SLFn.NDP2 = self.slf.NDP2
      SLFn.NDP3 = self.slf.NDP3

      return SLFn
コード例 #2
0
ファイル: partitioning.py プロジェクト: nicogodet/PostTelemac
    def copyCommonData(self):

        SLFn = SELAFIN("")
        #   Meta data
        SLFn.TITLE = self.slf.TITLE
        SLFn.file = self.slf.file
        SLFn.IPARAM = self.slf.IPARAM
        #   Time
        SLFn.DATETIME = self.slf.DATETIME
        SLFn.tags = self.slf.tags
        #   Variables
        SLFn.NBV1 = self.slf.NBV1
        SLFn.VARNAMES = self.slf.VARNAMES
        SLFn.VARUNITS = self.slf.VARUNITS
        SLFn.NBV2 = self.slf.NBV2
        SLFn.CLDNAMES = self.slf.CLDNAMES
        SLFn.CLDUNITS = self.slf.CLDUNITS
        SLFn.NVAR = self.slf.NVAR
        SLFn.VARINDEX = range(self.slf.NVAR)
        #   Unchanged numbers
        SLFn.NPLAN = self.slf.NPLAN
        SLFn.NDP2 = self.slf.NDP2
        SLFn.NDP3 = self.slf.NDP3

        return SLFn
コード例 #3
0
    support3d = zip(support2d, len(xys) * [range(slf.NPLAN)])

    # <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
    # ~~~~ writes BND header ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    bndFile = options.args[3]
    bnd = SELAFIN('')
    bnd.fole = {}
    bnd.fole.update({'hook': open(bndFile, 'wb')})
    bnd.fole.update({'name': bndFile})
    bnd.fole.update({'endian': ">"})  # big endian
    bnd.fole.update({'float': ('f', 4)})  # single precision

    # Meta data and variable names
    bnd.TITLE = ''
    bnd.NBV1 = 5
    # /!\ ELEVATION has to be the first variable
    # (for possible vertical re-interpolation within TELEMAC)
    bnd.VARNAMES = ['ELEVATION Z     ', \
                    'VELOCITY U      ','VELOCITY V      ', \
                    'SALINITY        ','TEMPERATURE     ' ]
    bnd.VARUNITS = ['M               ', \
                    'M/S             ','M/S             ', \
                    '                ','                ' ]
    bnd.NVAR = bnd.NBV1
    bnd.VARINDEX = range(bnd.NVAR)

    # Sizes and mesh connectivity
    bnd.NPLAN = slf.NPLAN
    bnd.NDP2 = 2
    bnd.NDP3 = 4
コード例 #4
0
   support3d = zip(support2d,len(xys)*[range(slf.NPLAN)])

# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
# ~~~~ writes INI header ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

   iniFile = options.args[2]
   ini = SELAFIN('')
   ini.fole = {}
   ini.fole.update({ 'hook': open(iniFile,'wb') })
   ini.fole.update({ 'name': iniFile})
   ini.fole.update({ 'endian': ">" })     # big endian
   ini.fole.update({ 'float': ('f',4) })  # single precision

   # Meta data and variable names
   ini.TITLE = ''
   ini.NBV1 = 5
   # /!\ ELEVATION has to be the first variable
   # (for possible vertical re-interpolation within TELEMAC)
   ini.VARNAMES = ['ELEVATION Z     ', \
                   'VELOCITY U      ','VELOCITY V      ', \
                   'SALINITY        ','TEMPERATURE     ' ]
   ini.VARUNITS = ['M               ', \
                   'M/S             ','M/S             ', \
                   '                ','                ' ]
   ini.NVAR = ini.NBV1
   ini.VARINDEX = range(ini.NVAR)

   # Sizes and mesh connectivity
   ini.NPLAN = slf.NPLAN
   ini.NDP2 = 3
   ini.NDP3 = 6
コード例 #5
0
ファイル: convertToATM.py プロジェクト: ogoe/OpenTelemac
    support3d = zip(support2d, len(xys) * [range(slf.NPLAN)])

    # <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
    # ~~~~ writes ATM header ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    atmFile = options.args[2]
    atm = SELAFIN('')
    atm.fole = {}
    atm.fole.update({'hook': open(atmFile, 'wb')})
    atm.fole.update({'name': atmFile})
    atm.fole.update({'endian': ">"})  # big endian
    atm.fole.update({'float': ('f', 4)})  # single precision

    # Meta data and variable names
    atm.TITLE = ''
    atm.NBV1 = 4
    atm.VARNAMES = ['SURFACE PRESSURE', \
       'WIND VELOCITY U ','WIND VELOCITY V ', \
       'AIR TEMPERATURE ']
    atm.VARUNITS = ['UI              ', \
       'M/S             ','M/S             ', \
       'DEGREES         ']
    atm.NVAR = atm.NBV1
    atm.VARINDEX = range(atm.NVAR)

    # Sizes and mesh connectivity
    atm.NPLAN = slf.NPLAN  # it should be 2D but why the heack not ...
    atm.NDP2 = slf.NDP2
    atm.NDP3 = slf.NDP3
    atm.NPOIN2 = geo.NPOIN2
    atm.NPOIN3 = geo.NPOIN2 * atm.NPLAN
コード例 #6
0
#fout = open("junk.out","w")

#for i in range(len(b)):
#	fout.write(str(b[i]) + '\n')

# now we can delete the temp file
os.remove(temp_nodes_file)
os.remove(temp_elements_file)
os.remove(temp_boundaries_file)

# now to write the SELAFIN mesh file
slf2d = SELAFIN('')

#print '     +> Set SELAFIN variables'
slf2d.TITLE = 'Converted from gmsh'
slf2d.NBV1 = 1
slf2d.NVAR = 1
slf2d.VARINDEX = range(slf2d.NVAR)
slf2d.VARNAMES.append('BOTTOM          ')
slf2d.VARUNITS.append('M               ')

#print '     +> Set SELAFIN sizes'
slf2d.NPLAN = 1
slf2d.NDP2 = 3
slf2d.NDP3 = 3
slf2d.NPOIN2 = n
slf2d.NPOIN3 = n
slf2d.NELEM2 = e
slf2d.NELEM3 = e
slf2d.IPARAM = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0]
コード例 #7
0
ファイル: gmsh2slf.py プロジェクト: pprodano/pputils
#fout = open("junk.out","w")

#for i in range(len(b)):
#	fout.write(str(b[i]) + '\n')

# now we can delete the temp file
os.remove(temp_nodes_file)
os.remove(temp_elements_file)
os.remove(temp_boundaries_file)

# now to write the SELAFIN mesh file
slf2d = SELAFIN('')

#print '     +> Set SELAFIN variables'
slf2d.TITLE = 'Converted from gmsh'
slf2d.NBV1 = 1 
slf2d.NVAR = 1
slf2d.VARINDEX = range(slf2d.NVAR)
slf2d.VARNAMES.append('BOTTOM          ')
slf2d.VARUNITS.append('M               ')

#print '     +> Set SELAFIN sizes'
slf2d.NPLAN = 1
slf2d.NDP2 = 3
slf2d.NDP3 = 3
slf2d.NPOIN2 = n
slf2d.NPOIN3 = n
slf2d.NELEM2 = e
slf2d.NELEM3 = e
slf2d.IPARAM = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0]