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
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
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
bnd.NPOIN2 = len(BOR)
bnd.NPOIN3 = bnd.NPOIN2 * slf.NPLAN
bnd.IPARAM = [0, 0, 0, 0, 0, 0, bnd.NPLAN, 0, 0, 0]
bnd.IPOB2 = BOR # /!\ Note that IPOBO keeps the original numbering
print ' +> masking and setting connectivity'
# Set the array that only includes elements of geo.IKLE2 with at least two nodes in BOR
MASK = geo.IKLE2[np.where(
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
atm.NELEM2 = geo.NELEM2
print ' +> setting connectivity'
if atm.NPLAN > 1:
atm.NELEM3 = geo.NELEM2 * (atm.NPLAN - 1)