Subject: Integrated OpemModelica Changes into develop branch

Description: Integrated OpemModelica Changes into develop branch

2 files changed, 640 insertions, 0 deletions

diff --git a/src/ngspicetoModelica/NgspicetoModelica.py b/src/ngspicetoModelica/NgspicetoModelica.py
new file mode 100644
index 00000000..0e326f17
--- /dev/null
+++ b/src/ngspicetoModelica/NgspicetoModelica.py
@@ -0,0 +1,640 @@
+
+
+
+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()
+
+
diff --git a/src/ngspicetoModelica/__init__.py b/src/ngspicetoModelica/__init__.py
new file mode 100644
index 00000000..e69de29b
--- /dev/null
+++ b/src/ngspicetoModelica/__init__.py