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.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( np.sum(np.in1d(geo.IKLE2, np.sort(BOR - 1)).reshape(geo.NELEM2, geo.NDP2),
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 ini.NPOIN2 = geo.NPOIN2 ini.NPOIN3 = geo.NPOIN2*ini.NPLAN ini.NELEM2 = geo.NELEM2 ini.NELEM3 = ini.NELEM2*(ini.NPLAN-1) print ' +> setting connectivity' ini.IKLE3 = \ np.repeat(geo.NPOIN2*np.arange(ini.NPLAN-1),geo.NELEM2*ini.NDP3).reshape((geo.NELEM2*(ini.NPLAN-1),ini.NDP3)) + \ np.tile(np.add(np.tile(geo.IKLE2,2),np.repeat(geo.NPOIN2*np.arange(2),geo.NDP2)),(ini.NPLAN-1,1))
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) atm.IKLE2 = geo.IKLE2 atm.IKLE3 = \
#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] #print ' +> Set SELAFIN mesh'