Hi all,

I'm looking for some help with subroutines and multiple cpus. I'm using the disp subroutine to apply temperature boundary condition to a surface. When I try to run multiple cpus, abaqus hits an error when standard starts. I'm curious how I can fix this. Is there something in my subroutine that is preventing this? I do have tokens for the cpus I'm trying to run.

Please forgive me, I'm new to any sort of programming. Here is my subroutine for reference:

SUBROUTINE DISP(U,KSTEP,KINC,TIME,NODE,NOEL,JDOF,COORDS)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION U(3),TIME(2),COORDS(3),xi(10),zi(7)
real temperature(70,347),n,m,length,width
 
X=coords(1)
Y=coords(2)
Z=coords(3)
 
length=1.500
width=0.500
m=10.000
n=7.000
ti=time(2)/120
 
do i=1,m+1
xi(i)=(length/(m)*(i-1))
end do
xi(1)=xi(1)-.01
xi(m+1)=xi(m+1)+.01
do i=1,n+1
zi(i)=(width/(n)*(i-1))
end do
zi(1)=zi(1)-.01
zi(n+1)=zi(n+1)+.01
 
open(unit=2, file="C:\Makersat Thermal Model\Flight Model Temp BC\temperature.txt", status='old', action='read')
 
read(2,*)temperature
close(2)
 
do i=2,m+1
do j=2,n+1
if (x<=xi(i) .and. x>=xi(i-1) .and. z<=zi(j) .and. z>=zi(j-1)) then
xm=m*(j-2)+(i-1)
U(1)=temperature(xm,ti)
endif
enddo
enddo
 
C

RETURN
END