modified sponge for airfoil case and added surfID driven forces for f…

…orces.dat

phSolver/compressible/e3source.f

@@ -75,12 +75,7 @@ subroutine e3source (ri, rmi, rlyi,
  src(:,5) = u1*src(:,2)
  endwhere
 
- else if(matflg(5,1).ge.4) then ! cool case (sponge outside of a
- ! revolved box defined from zinSponge to
- ! zoutSponge in axial extent and 0
- ! to radSponge in radial extent for
- ! all theta)
-
+ else if(matflg(5,1).ge.4) then 
 c determine coordinates of quadrature pt
  xx=zero
  do n = 1,nenl
@@ -97,8 +92,9 @@ subroutine e3source (ri, rmi, rlyi,
  enddo
  if (1.eq.1) then ! bringing in x BL sponge
  do id=1,npro
- if((xx(id,1).gt.zoutSponge)) then
- bcool(id)=grthOSponge*(xx(id,1)-zoutSponge)**2
+ rsq=xx(id,1)*xx(id,1) + xx(id,2)*xx(id,2)
+ if(rsq.gt.zoutSponge) then
+ bcool(id)=grthOSponge*(rsq-zoutSponge)**2
  bcool(id)=min(bcool(id),betamax)
 c Determine the resulting density and energies
  den = ytargeti(id,1) / (Rgas * ytargeti(id,5))

phSolver/compressible/rstat.f

@@ -82,10 +82,11 @@ subroutine rstat (res, ilwork, b)
  HFlux = Forout(4)
  if (myrank .eq. master) then
  write (iforce,1000) lstep+1, (Force(i), i=1,nsd), HFlux, 
- & spmasss
+ & spmasss,(ftot(i),i=1,nsd)
  call flush(iforce)
  endif
  endif
+
 c
 c.... -----------------------> Convergence <-------------------------
 c
@@ -143,7 +144,7 @@ subroutine rstat (res, ilwork, b)
 c
  return
 c
-1000 format(1p,i6,5e13.5)
+1000 format(1p,i6,8e13.5)
 2000 format(1p,i6,e10.3,e10.3,1x,'(',i4,')',
  & 1x,e10.3,1x,'(',i4,')',
  & ' [',i3,'-',i3,']',i10)
@@ -327,7 +328,7 @@ subroutine rstatp (resNrm,resNrmP)
  HFlux = Forout(4)
  if (myrank .eq. master) then
  write (iforce,1000) lstep+1, (Force(i), i=1,nsd), HFlux, 
- & spmasss
+ & spmasss,(ftot(i),i=1,nsd)
  call flush(iforce)
  endif
  endif