path: root/src/ngspicetoModelica/NgspicetoModelica.py

import sys
import os.path
from string import maketrans

def readNetlist(filename):
    """Read Ngspice Netlist"""
    if os.path.exists(filename):
        try:
            f = open(filename)
        except:
            print("Error in opening file")
            sys.exit()
    else:
        print filename + " does not exist"
        sys.exit()

    data = f.read()
#    data = data.translate(maketrans('\n+', '  '))
    f.close()
    return data.splitlines()


def separateNetlistInfo(data):
  """ Separate schematic data and option data"""
  optionInfo=[]
  schematicInfo=[]
 
  for eachline in data:
    if len(eachline) > 1:
##     if eachline[0] == '+':
##       eachline=eachline.translate(maketrans('\n+','  '))
     if eachline[0]=='*':
       continue
     elif eachline[0]=='.':
       optionInfo.append(eachline.lower())
     else:
       schematicInfo.append(eachline.lower())
  return optionInfo,schematicInfo


def addModel(optionInfo):
   """ Add model parameters in the modelica file and create dictionary of model parameters"""
   modelName = []
   modelInfo = {}
   subcktName = []
   paramInfo = []
   #modelInfo['paramInfo'] = {}
   for eachline in optionInfo:
    words = eachline.split()
    if words[0] == '.include':
      name = words[1].split('.')
      if name[1] == 'lib':
        modelName.append(name[0])
      if name[1] == 'sub':
  	subcktName.append(name[0])
    elif words[0] == '.param':
      paramInfo.append(eachline) 
   for eachmodel in modelName:
     filename = eachmodel + '.lib'
     if os.path.exists(filename):
        try:
            f = open(filename)
        except:
            print("Error in opening file")
            sys.exit()
     else:
        print filename + " does not exist"
        sys.exit()

     data = f.read()
     data = data.lower()
     newdata = data.split('(')
     newdata = newdata[1].split()
     modelInfo[eachmodel] = {}
     for eachline in newdata:
       if len(eachline) > 1:
        info = eachline.split('=')
       # modelInfo[eachmodel][info[0]] = {}
        for eachitem in info:
         modelInfo[eachmodel][info[0]] = info[1] #dic within a dic
     #modelInfo[eachmodel] = modelInfo[eachmodel].split()
    # modelInfo[eachmodel] = modelInfo[eachmodel].lower()
     f.close()
          
   return modelName, modelInfo, subcktName, paramInfo     
     

def processParam(paramInfo):
    """ Process parameter info and update in Modelica syntax"""
    modelicaParam = []
    for eachline in paramInfo:
      eachline = eachline.split('.param')
      stat = 'parameter Real ' + eachline[1] + ';'
      stat = stat.translate(maketrans('{}', '  '))
      modelicaParam.append(stat)
    return modelicaParam


def separatePlot(schematicInfo):
   """ separate print plot and component statements"""
   compInfo = []
   plotInfo = []
   
   for eachline in schematicInfo:
    words = eachline.split()
    if words[0] == 'run':
      continue
    elif words[0] == 'plot' or words[0] == 'print':
      plotInfo.append(eachline)
    else:
      compInfo.append(eachline)
   return compInfo, plotInfo

def separateSource(compInfo):
   """Find if dependent sources are present in the schematic and if so make a dictionary with source details"""
   sourceInfo = {}
   source = []
   for eachline in compInfo:
     if eachline[0] in ['f', 'h']:
       source.append(words[3])
   if len(source) > 0:
     for eachline in compInfo:
       words_s = eachline.split()
       if words_s[0] in source:
         sourceInfo[words_s[0]] = words_s[1:3]
   return sourceInfo


def splitIntoVal(val):
   """ Split the number k,u,p,t,g etc into powers e3,e-6 etc"""
   for i in range(0,len(val),1):
       if val[i] in ['k','u','p','t','g','m','n','f']:
         newval = val.split(val[i])
         if val[i] == 'k':
           value = newval[0] + 'e3'
         if val[i] == 'u':
           value = newval[0] + 'e-6'
         if val[i] == 'p':
           value = newval[0] + 'e-12'
         if val[i] == 't':  
           value = newval[0] + 'e12'
         if val[i] == 'g':
           value = newval[0] + 'e9'
         if val[i] == 'm':
           if i != len(val)-1:
             if val[i+1] == 'e':
              value = newval[0] + 'e6'
           else:
            value = newval[0] +'e-3'
         if val[i] == 'n':
           value = newval[0] + 'e-9'       
         if val[i] == 'f':
           value = newval[0] +'e-15'
       else:
         value = val
   return value

      
def compInit(compInfo, node, modelInfo, subcktName):
   """For each component in the netlist initialise it acc to Modelica format"""
#### initial processign to check if MOs is present. If so, library to be used is BondLib
   modelicaCompInit = []
   numNodesSub = {} 
   IfMOS = '0'
   for eachline in compInfo:
#     words = eachline.split()
     if eachline[0] == 'm':
       IfMOS = '1'
       break
   if len(subcktName) > 0:
     subOptionInfo = []
     subSchemInfo = []
     for eachsub in subcktName:
       filename_tem = eachsub + '.sub'
       data = readNetlist(filename_tem)
       subOptionInfo, subSchemInfo = separateNetlistInfo(data)
       for eachline in subSchemInfo:
#        words = eachline.split()
        if eachline[0] == 'm':
          IfMOS = '1'
          break
   for eachline in compInfo:
     words = eachline.split()
     val = words[3]
     value = splitIntoVal(val)
     if eachline[0] == 'r':
       stat = 'Analog.Basic.Resistor ' + words[0] + '(R = ' + value + ');'
       modelicaCompInit.append(stat)
     elif eachline[0] == 'c':
       stat = 'Analog.Basic.Capacitor ' + words[0] + '(C = ' + value + ');'
       modelicaCompInit.append(stat)
     elif eachline[0] == 'l':
       stat = 'Analog.Basic.Inductor ' + words[0] + '(L = ' + value + ');'
       modelicaCompInit.append(stat) 
     elif eachline[0] == 'e':
       stat = 'Analog.Basic.VCV ' + words[0] + '(gain = ' + splitIntoVal(words[5]) + ');'
       modelicaCompInit.append(stat) 
     elif eachline[0] == 'g':
       stat = 'Analog.Basic.VCC ' + words[0] + '(transConductance = ' + splitIntoVal(words[5]) + ');'
       modelicaCompInit.append(stat) 
     elif eachline[0] == 'f':
       stat = 'Analog.Basic.CCC ' + words[0] + '(gain = ' + splitIntoVal(words[4]) + ');'
       modelicaCompInit.append(stat) 
     elif eachline[0] == 'h':
       stat = 'Analog.Basic.CCV ' + words[0] + '(transResistance = ' + splitIntoVal(words[4]) + ');'
       modelicaCompInit.append(stat) 
     elif eachline[0] == 'd':
       if len(words) > 3:
        n = float(modelInfo[words[3]]['n'])
        vt_temp = 0.025*n
        vt = str(vt_temp)
        stat = 'Analog.Semiconductors.Diode ' + words[0] + '(Ids = ' + modelInfo[words[3]]['is'] + ', Vt = ' + vt + ', R = 1e12' +');'
       else:
        stat = 'Analog.Semiconductors.Diode ' + words[0] +';'
       modelicaCompInit.append(stat)
     elif eachline[0] == 'm':
       line_l = words[7].split('=')
       line_w = words[8].split('=')
       line_pd = words[9].split('=')
       line_ps = words[10].split('=')
       line_ad = words[11].split('=')
       line_as = words[12].split('=')
       if words[5] == "mos_n" or words[5] == "mosfet_n":
          start = 'BondLib.Electrical.Analog.Spice.Mn '
       if words[5] == "mos_p" or words[5] == "mosfet_p":
          start = 'BondLib.Electrical.Analog.Spice.Mp '
       stat = start + words[0] + '(Tnom = 300, VT0 = ' + modelInfo[words[5]]['vto'] + ', GAMMA = ' + modelInfo[words[5]]['gamma'] + ', PHI = ' + modelInfo[words[5]]['phi'] + ', LD = ' + splitIntoVal(modelInfo[words[5]]['ld']) + ', U0 = ' + str(float(splitIntoVal(modelInfo[words[5]]['uo']))*0.0001) + ', LAMBDA = ' + modelInfo[words[5]]['lambda'] + ', TOX = ' + splitIntoVal(modelInfo[words[5]]['tox']) + ', PB = ' + modelInfo[words[5]]['pb'] + ', CJ = ' + splitIntoVal(modelInfo[words[5]]['cj']) + ', CJSW = ' + splitIntoVal(modelInfo[words[5]]['cjsw']) + ', MJ = ' + modelInfo[words[5]]['mj'] + ', MJSW = ' + modelInfo[words[5]]['mjsw'] + ', CGD0 = ' + splitIntoVal(modelInfo[words[5]]['cgdo']) + ', JS = ' + splitIntoVal(modelInfo[words[5]]['js']) + ', CGB0 = ' + splitIntoVal(modelInfo[words[5]]['cgbo']) + ', CGS0 = ' + splitIntoVal(modelInfo[words[5]]['cgso']) + ', L = ' + splitIntoVal(line_l[1]) + ', W = ' + line_w[1] + ', Level = 1' + ', AD = ' + line_ad[1] + ', AS = ' + line_as[1] + ', PD = ' + line_pd[1] + ', PS = ' + line_pd[1] + ');'
       stat = stat.translate(maketrans('{}', '  '))
       modelicaCompInit.append(stat)
     elif eachline[0] == 'v':
       typ = words[3].split('(')
       if typ[0] == "pulse":
          per = words[9].split(')')
 #         if IfMOS == '0':
 #          stat = 'Spice3.Sources.V_pulse '+words[0]+'(TR = '+words[6]+', V2 = '+words[4]+', PW = '+words[8]+', PER = '+per[0]+', V1 = '+typ[1]+', TD = '+words[5]+', TF = '+words[7]+');'
 #         elif IfMOS == '1': 
          stat = 'Analog.Sources.TrapezoidVoltage '+words[0]+'(rising = '+words[6]+', V = '+words[4]+', width = '+words[8]+', period = '+per[0]+', offset = '+typ[1]+', startTime = '+words[5]+', falling = '+words[7]+');'
          modelicaCompInit.append(stat)
       if typ[0] == "sine":
          theta = words[7].split(')')
#          if IfMOS == '0':
#           stat = 'Spice3.Sources.V_sin '+words[0]+'(VO = '+typ[1]+', VA = '+words[4]+', FREQ = '+words[5]+', TD = '+words[6]+', THETA = '+theta[0]+');'
#          elif IfMOS == '1':
          stat = 'Analog.Sources.SineVoltage '+words[0]+'(offset = '+typ[1]+', V = '+words[4]+', freqHz = '+words[5]+', startTime = '+words[6]+', phase = '+theta[0]+');'
          modelicaCompInit.append(stat)
       if typ[0] == "pwl":
#          if IfMOS == '0':
#           keyw = 'Spice3.Sources.V_pwl '
#          elif IfMOS == '1':
          keyw = 'Analog.Sources.TableVoltage '
          stat = keyw + words[0] + '(table = [' + typ[1] + ',' + words[4] + ';' 
          length = len(words);
          for i in range(6,length,2):
             if i == length-2:
               w = words[i].split(')')
               stat = stat + words[i-1] + ',' + w[0] 
             else:
               stat = stat + words[i-1] + ',' + words[i] + ';'
          stat = stat + ']);'
          modelicaCompInit.append(stat) 
       if typ[0] == words[3] and typ[0] != "dc":
          val_temp = typ[0].split('v')
#          if IfMOS  == '0':
          stat = 'Analog.Sources.ConstantVoltage ' + words[0] + '(V = ' + val_temp[0] + ');'
#          elif IfMOS == '1':
#           stat = 'Analog.Sources.ConstantVoltage ' + words[0] + '(V = ' + val_temp[0] + ');'
          modelicaCompInit.append(stat)
       elif typ[0] == words[3] and typ[0] == "dc":
#          if IfMOS  == '0':
#           stat = 'Spice3.Sources.V_constant ' + words[0] + '(V = ' + words[4] + ');'    ### check this
#          elif IfMOS == '1':
          stat = 'Analog.Sources.ConstantVoltage ' + words[0] + '(V = ' + words[4] + ');'    ### check this
          modelicaCompInit.append(stat)
     elif eachline[0] == 'x':
       temp_line = eachline.split()
       temp = temp_line[0].split('x')
       index = temp[1]
       for i in range(0,len(temp_line),1):
         if temp_line[i] in subcktName:
           subname = temp_line[i]
           numNodesSub[subname] = i - 1
           point = i
       if len(temp_line) > point + 1:
         rem = temp_line[point+1:len(temp_line)]
         rem_new = ','.join(rem)
         stat = subname + ' ' + subname +'_instance' + index + '(' +  rem_new + ');'
       else:
         stat = subname + ' ' + subname +'_instance' + index + ';'
       modelicaCompInit.append(stat)
     else:
       continue
   if '0' in node:
      modelicaCompInit.append('Analog.Basic.Ground g;')
   return modelicaCompInit, numNodesSub

def getSubInterface(subname, numNodesSub):
   """ Get the list of nodes for subcircuit in .subckt line"""
   subOptionInfo_p = []
   subSchemInfo_p = []
   filename_t = subname + '.sub'
   data_p = readNetlist(filename_t)
   subOptionInfo_p, subSchemInfo_p = separateNetlistInfo(data_p)
   if len(subOptionInfo_p) > 0:
     newline = subOptionInfo_p[0]
     newline = newline.split('.subckt '+ subname)       
     intLine = newline[1].split()
     newindex = numNodesSub[subname]
     nodesInfoLine = intLine[0:newindex]
   return nodesInfoLine 


def getSubParamLine(subname, numNodesSub, subParamInfo):
   """ Take subcircuit name and give the info related to parameters in the first line and initislise it in """
#   nodeSubInterface = []
   subOptionInfo_p = []
   subSchemInfo_p = []
   filename_t = subname + '.sub'
   data_p = readNetlist(filename_t)
   subOptionInfo_p, subSchemInfo_p = separateNetlistInfo(data_p)
   if len(subOptionInfo_p) > 0:
     newline = subOptionInfo_p[0]
     newline = newline.split('.subckt '+ subname)       
     intLine = newline[1].split()
     newindex = numNodesSub[subname]
     appen_line = intLine[newindex:len(intLine)]
     appen_param = ','.join(appen_line)
     paramLine = 'parameter Real ' + appen_param + ';'
     paramLine = paramLine.translate(maketrans('{}', '  '))
     subParamInfo.append(paramLine)
   return subParamInfo

def nodeSeparate(compInfo, ifSub, subname, subcktName):
   """ separate the node numbers and create nodes in modelica file; the nodes in the subckt line should not be inside protected keyword. pinInit is the one that goes under protected keyword."""
   node = []
   nodeTemp = []
   nodeDic = {}
   pinInit = 'Modelica.Electrical.Analog.Interfaces.Pin '
   pinProtectedInit = 'Modelica.Electrical.Analog.Interfaces.Pin '
   protectedNode = []
   for eachline in compInfo:
     words = eachline.split()
     if eachline[0] in ['m', 'e', 'g', 't']:
      nodeTemp.append(words[1])
      nodeTemp.append(words[2])
      nodeTemp.append(words[3])
      nodeTemp.append(words[4])
     elif eachline[0] in ['q', 'j']:
      nodeTemp.append(words[1])
      nodeTemp.append(words[2])
      nodeTemp.append(words[3])
     elif eachline[0] == 'x':
      templine = eachline.split()
      for i in range(0,len(templine),1):
        if templine[i] in subcktName:
          point = i   
      nodeTemp.extend(words[1:point])
     else:
      nodeTemp.append(words[1])
      nodeTemp.append(words[2])
   for i in nodeTemp:
    if i not in node:
      node.append(i)
   for i in range(0, len(node),1):
     nodeDic[node[i]] = 'n' + node[i]
     if ifSub == '0':
       if i != len(node)-1:
        pinInit = pinInit + nodeDic[node[i]] + ', '
       else:
        pinInit = pinInit + nodeDic[node[i]] 
     else:
       nonprotectedNode = getSubInterface(subname, numNodesSub)
       if node[i] in nonprotectedNode:
        continue
       else:
        protectedNode.append(node[i])
   if ifSub == '1': 
     if len(nonprotectedNode) > 0:    
      for i in range(0, len(nonprotectedNode),1):
        if i != len(nonprotectedNode)-1:
         pinProtectedInit = pinProtectedInit + nodeDic[nonprotectedNode[i]] + ','
        else:
         pinProtectedInit = pinProtectedInit + nodeDic[nonprotectedNode[i]]
     if len(protectedNode) > 0:
      for i in range(0, len(protectedNode),1):
        if i != len(protectedNode)-1: 
         pinInit = pinInit + nodeDic[protectedNode[i]] + ','
        else:
         pinInit = pinInit + nodeDic[protectedNode[i]] 
   pinInit = pinInit + ';'
   pinProtectedInit = pinProtectedInit + ';'
   return node, nodeDic, pinInit, pinProtectedInit
  
def connectInfo(compInfo, node, nodeDic, numNodesSub):
   """Make node connections in the modelica netlist"""
   connInfo = []
   sourcesInfo = separateSource(compInfo)
   for eachline in compInfo:
     words = eachline.split()
     if eachline[0] == 'r' or eachline[0] == 'c' or eachline[0] == 'd' or eachline[0] == 'l' or eachline[0] == 'v':
       conn = 'connect(' + words[0] + '.p,' + nodeDic[words[1]] + ');'
       connInfo.append(conn)
       conn = 'connect(' + words[0] + '.n,' + nodeDic[words[2]] + ');'
       connInfo.append(conn)
     elif eachline[0] == 'm':
       conn = 'connect(' + words[0] + '.D,' + nodeDic[words[1]] + ');'
       connInfo.append(conn)
       conn = 'connect(' + words[0] + '.G,' + nodeDic[words[2]] + ');'
       connInfo.append(conn)
       conn = 'connect(' + words[0] + '.S,' + nodeDic[words[3]] + ');'
       connInfo.append(conn)
       conn = 'connect(' + words[0] + '.B,' + nodeDic[words[4]] + ');'
       connInfo.append(conn)
     elif eachline[0] in ['f','h']:
      vsource = words[3]
      sourceNodes = sourcesInfo[vsource]
      sourceNodes = sourceNodes.split()
      conn = 'connect(' + words[0] + '.p1,'+ nodeDic[sourceNodes[0]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.n1,'+ nodeDic[sourceNodes[1]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.p2,'+ nodeDic[words[1]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.n2,'+ nodeDic[words[2]] + ');'
      connInfo.append(conn)
     elif eachline[0] in ['g','e']:
      conn = 'connect(' + words[0] + '.p1,'+ nodeDic[words[3]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.n1,'+ nodeDic[words[4]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.p2,'+ nodeDic[words[1]] + ');'
      connInfo.append(conn)
      conn = 'connect(' + words[0] + '.n2,'+ nodeDic[words[2]] + ');'
      connInfo.append(conn)
     elif eachline[0] == 'x':
      templine = eachline.split()
      temp = templine[0].split('x')
      index = temp[1]
      for i in range(0,len(templine),1):
        if templine[i] in subcktName:
          subname = templine[i]
      nodeNumInfo = getSubInterface(subname, numNodesSub)
      for i in range(0, numNodesSub[subname], 1):
#        conn = 'connect(' + subname + '_instance' + index + '.' + nodeDic[nodeNumInfo[i]] + ',' + nodeDic[words[i+1]] + ');'
        conn = 'connect(' + subname + '_instance' + index + '.' + 'n'+ nodeNumInfo[i] + ',' + nodeDic[words[i+1]] + ');'
        connInfo.append(conn)              
     else:
     #elif eachline[0] == 'q':
     #elif eachline[0] == 'j':
       continue
   if '0' in node:
     conn = 'connect(g.p,n0);'
     connInfo.append(conn)
   return connInfo
## For testing
     

if len(sys.argv) < 2:
  filename=raw_input('Enter file name: ')
else:
  filename=sys.argv[1]


# get all the required info
lines=readNetlist(filename)
optionInfo, schematicInfo=separateNetlistInfo(lines)
modelName, modelInfo, subcktName, paramInfo = addModel(optionInfo)
modelicaParamInit = processParam(paramInfo)
compInfo, plotInfo = separatePlot(schematicInfo)
IfMOS = '0'
for eachline in compInfo:
 words = eachline.split()
 if eachline[0] == 'm':
  IfMOS = '1'
  break
if len(subcktName) > 0:
 subOptionInfo = []
 subSchemInfo = []
 for eachsub in subcktName:
  filename_temp = eachsub + '.sub'
  data = readNetlist(filename_temp)
  subOptionInfo, subSchemInfo = separateNetlistInfo(data)
  for eachline in subSchemInfo:
   words = eachline.split()
   if eachline[0] == 'm':
    IfMOS = '1'
    break
node, nodeDic, pinInit, pinProtectedInit = nodeSeparate(compInfo, '0', [], subcktName)
modelicaCompInit, numNodesSub  = compInit(compInfo,node, modelInfo, subcktName)
connInfo = connectInfo(compInfo, node, nodeDic, numNodesSub)

####Extract subckt data
def procesSubckt(subcktName):
   """ Process the subcircuit file .sub in the project folder"""
#   subcktDic = {}
   subOptionInfo = []
   subSchemInfo = []
   subModel = []
   subModelInfo = {}
   subsubName = [] 
   subParamInfo = []
   nodeSubInterface = []
   nodeSub = []
   nodeDicSub = {}
   pinInitsub = []
   connSubInfo = []
   if len(subcktName) > 0:
    for eachsub in subcktName:
     filename = eachsub + '.sub'
     data = readNetlist(filename)
     subOptionInfo, subSchemInfo = separateNetlistInfo(data)
     if len(subOptionInfo) > 0:
       newline = subOptionInfo[0]
       subInitLine = newline
       newline = newline.split('.subckt')       
       intLine = newline[1].split()
       for i in range(0,len(intLine),1):
         nodeSubInterface.append(intLine[i])
     subModel, subModelInfo, subsubName, subParamInfo = addModel(subOptionInfo)
     IfMOSsub = '0'
     for eachline in subSchemInfo:
#      words = eachline.split()
      if eachline[0] == 'm':
        IfMOSsub = '1'
        break
     if len(subsubName) > 0:
      subsubOptionInfo = []
      subsubSchemInfo = []
      for eachsub in subsubName:
       filename_stemp = eachsub + '.sub'
       data = readNetlist(filename_stemp)
       subsubOptionInfo, subsubSchemInfo = separateNetlistInfo(data)
       for eachline in subsubSchemInfo:
#        words = eachline.split()
        if eachline[0] == 'm':
         IfMOSsub = '1'
         break
     modelicaSubParam =  processParam(subParamInfo)
     nodeSub, nodeDicSub, pinInitSub, pinProtectedInitSub = nodeSeparate(subSchemInfo, '1', eachsub, subsubName)
     modelicaSubCompInit, numNodesSubsub = compInit(subSchemInfo, nodeSub, subModelInfo, subsubName)
     modelicaSubParamNew = getSubParamLine(eachsub, numNodesSub, modelicaSubParam)
     connSubInfo = connectInfo(subSchemInfo, nodeSub, nodeDicSub, numNodesSubsub)
     newname = filename.split('.')
     newfilename = newname[0]
     outfilename = newfilename+ ".mo"
     out = open(outfilename,"w")
     out.writelines('model ' + os.path.basename(newfilename))
     out.writelines('\n')
     if IfMOSsub == '0':
      out.writelines('import Modelica.Electrical.*;')
     elif IfMOSsub == '1':
      out.writelines('import BondLib.Electrical.*;')
     out.writelines('\n') 
     for eachline in modelicaSubParamNew:
       if len(subParamInfo) == 0:
         continue
       else:
        out.writelines(eachline) 
        out.writelines('\n') 
     for eachline in modelicaSubCompInit:
      if len(subSchemInfo) == 0:
       continue
      else:
       out.writelines(eachline)
       out.writelines('\n')
     out.writelines(pinProtectedInitSub)
     out.writelines('\n')
     if pinInitSub != 'Modelica.Electrical.Analog.Interfaces.Pin ;':
      out.writelines('protected')
      out.writelines('\n')
      out.writelines(pinInitSub)
      out.writelines('\n')
     out.writelines('equation')
     out.writelines('\n')
     for eachline in connSubInfo:
       if len(connSubInfo) == 0:
        continue
       else:
        out.writelines(eachline)
        out.writelines('\n')
     out.writelines('end '+ os.path.basename(newfilename) + ';')
     out.writelines('\n')
     out.close()
   
   return data, subOptionInfo, subSchemInfo, subModel, subModelInfo, subsubName, subParamInfo, modelicaSubCompInit, modelicaSubParam, nodeSubInterface, nodeSub, nodeDicSub, pinInitSub, connSubInfo 

if len(subcktName) > 0:
 data, subOptionInfo, subSchemInfo, subModel, subModelInfo, subsubName, subParamInfo, modelicaSubCompInit, modelicaSubParam,  nodeSubInterface, nodeSub, nodeDicSub, pinInitSub, connSubInfo = procesSubckt(subcktName)

# creating final output

newfile = filename.split('.')
newfilename = newfile[0]
outfile = newfilename + ".mo"
out = open(outfile,"w")
out.writelines('model ' + os.path.basename(newfilename))
out.writelines('\n')
if IfMOS == '0':
 out.writelines('import Modelica.Electrical.*;')
elif IfMOS == '1':
 out.writelines('import BondLib.Electrical.*;')
#out.writelines('import Modelica.Electrical.*;')
out.writelines('\n')

for eachline in modelicaParamInit:
  if len(paramInfo) == 0:
    continue
  else:
    out.writelines(eachline)
    out.writelines('\n')
for eachline in modelicaCompInit:
  if len(compInfo) == 0:
    continue
  else:
    out.writelines(eachline)
    out.writelines('\n')

out.writelines('protected')
out.writelines('\n')
out.writelines(pinInit)
out.writelines('\n')
out.writelines('equation')
out.writelines('\n')

for eachline in connInfo:
  if len(connInfo) == 0:
    continue
  else:
    out.writelines(eachline)
    out.writelines('\n')

out.writelines('end '+ os.path.basename(newfilename) + ';')
out.writelines('\n')


out.close()