import visualization
from abaqus import *
from abaqusConstants import *
import xyPlot
import displayGroupOdbToolset as dgo
from odbAccess import openOdb
from viewerModules import *
from driverUtils import executeOnCaeStartup  
import sys
import os
import visualization
from abaqusConstants import *
#-----------------------------
l1 = [200.0,200.0]  #Bottom flange length
height = [1000.0] #Total height web+flange..5l1,5l1,6l1,7l1,8l1
l2 = [200.0,200.0]  #Top flange length
T1 = [20.0,25.0]   #Extruded Bottom flange thickness
delta_t= [20.0,25.0]# Difference of thickness
multiplied_list = [element * 4 for element in delta_t]
partition_datum_plane_XY = [750.0,1080.0,1250.0,1786.0]  #offsets in Vertical Direction
#------------------------------

# list of ODB-Files
odb_names = ['Para_Analysis_H_1000_W_200_t1_40_T1_20.odb','Para_Analysis_H_1000_W_200_t1_50_T1_30.odb',
'Para_Analysis_H_1000_W_200_t1_30_T1_15.odb']  # Dummy list with length of odb_names list. In this we can store the ODB objects 
path='./'                          #I don't know writing is right or not
odbs = [1]*len(odb_names); 
vips = [1]*len(odb_names); 
for s,odb_name in enumerate(odb_names):    
	print('\n\nlooping for %s ODB...'%odb_names)     
	vips[s] = session.Viewport(name='Girder Viewport_%d'%(s+1))     
	odbs[s] = session.openOdb(name=odb_name)

	path=session.Path(name='Path-1', type=POINT_LIST, expression=((l1[i], T1[i],partition_datum_plane_XY[0]), 
	(l1[i], T1[i],partition_datum_plane_XY[1]), (l1[i], t1[i],multiplied_list[i]+partition_datum_plane_XY[1]),
	(l1[i], t1[i],partition_datum_plane_XY[2])))    #giving error as it cannot find the coordinates
	
	pth=session.path['Path-1']
	vips[s].setValues(displayedObject=odbs[s])
	newXYData=session.XYDataFromPath(name=xyName,path=pth,
	includeIntersections=TRUE,
	shape=UNDEFORMED,
	labelType=TRUE_DISTANCE,
	step=1,
	frame=1,
	variable=(('S',INTEGRATION_POINT,((COMPONENT, 'S11' ),)),))
	session.writeXYReport(fileName='%s.txt'%name, xyData=(newXYData,))