Commit 7da6ee8a authored by Elia Merzari's avatar Elia Merzari
Browse files

Adding time reset

parent cfe3cb17
......@@ -145,6 +145,8 @@ c endif
t(i,1,1,1,1)=min(t(i,1,1,1,1),773.0)
enddo
if (istep.eq.0) time=0.0
return
end
c-----------------------------------------------------------------------
......@@ -191,7 +193,7 @@ c implicit none
cbc(ifc,iel,2) = 't '
elseif (id_face.eq.2) then ! SideSet 2 outlet
cbc(ifc,iel,1) = 'O '
cbc(ifc,iel,2) = 'O '
cbc(ifc,iel,2) = 'I '
else
cbc(ifc,iel,1) = 'E '
cbc(ifc,iel,2) = 'E '
......
......@@ -145,6 +145,8 @@ c endif
t(i,1,1,1,1)=min(t(i,1,1,1,1),773.0)
enddo
if (istep.eq.0) time=0.0
return
end
c-----------------------------------------------------------------------
......
......@@ -139,13 +139,15 @@ c endif
ntot=lx1*ly1*lz1*nelt
tmax=glmax(t(1,1,1,1,1),ntot)
tmin=glmin(t(1,1,1,1,1),ntot)
if (nid.eq.0) write(6,*)"T bounds > ",tmin,tmax
if (nid.eq.0) write(6,*)"T bounds > ",time,tmin,tmax
do i=1,ntot
t(i,1,1,1,1)=max(t(i,1,1,1,1),573.0)
t(i,1,1,1,1)=min(t(i,1,1,1,1),773.0)
enddo
if (istep.eq.0) time=0.0
return
end
c-----------------------------------------------------------------------
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment