i20xsave.F File Reference

#include "implicit_f.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i20xsave (itask, xa, nty, nsn, nmn, nsne, nmne, nln, nsv, msr, xsav, nsve, msre, xsave, xmin, ymin, zmin, xmax, ymax, zmax, c_max, curv_max, icurv, irect, nrtm_t, xmine, ymine, zmine, xmaxe, ymaxe, zmaxe)

Function/Subroutine Documentation

i20xsave()

subroutine i20xsave	(	integer	itask,
			xa,
		integer	nty,
		integer	nsn,
		integer	nmn,
		integer	nsne,
		integer	nmne,
		integer	nln,
		integer, dimension(*)	nsv,
		integer, dimension(*)	msr,
			xsav,
		integer, dimension(*)	nsve,
		integer, dimension(*)	msre,
			xsave,
			xmin,
			ymin,
			zmin,
			xmax,
			ymax,
			zmax,
			c_max,
			curv_max,
		integer	icurv,
		integer, dimension(4,*)	irect,
		integer	nrtm_t,
			xmine,
			ymine,
			zmine,
			xmaxe,
			ymaxe,
			zmaxe )

Definition at line 28 of file i20xsave.F.

C sauvegarde des XSAV
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NSN,NMN,ITASK,NTY, ICURV, NRTM_T,NLN,
     .        NSV(*),MSR(*)
      INTEGER NSNE,NMNE,NSVE(*),MSRE(*), IRECT(4,*)
      my_real
     .   xmin, ymin, zmin, xmax, ymax, zmax, c_max,
     .   xmine, ymine, zmine, xmaxe, ymaxe, zmaxe,
     .   xa(3,*),xsav(3,*),xsave(3,*),
     .   curv_max(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NSNF,NMNF,NSNL,NMNL,I, J, II
      INTEGER NSNEF,NMNEF,NSNEL,NMNEL
      my_real
     .        xxx, yyy, zzz 
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
 
      nsnf=1+itask*nsn/nthread
      nsnl=(itask+1)*nsn/nthread
      nmnf=1+itask*nmn/nthread
      nmnl=(itask+1)*nmn/nthread
 
      nsnef=1+itask*nsne/nthread
      nsnel=(itask+1)*nsne/nthread
      nmnef=1+itask*nmne/nthread
      nmnel=(itask+1)*nmne/nthread
 
      xmin=ep30
      xmax=-ep30
      ymin=ep30
      ymax=-ep30
      zmin=ep30
      zmax=-ep30
 
      xmine=ep30
      xmaxe=-ep30
      ymine=ep30
      ymaxe=-ep30
      zmine=ep30
      zmaxe=-ep30
 
 
      DO i=nsnf,nsnl
           j=abs(nsv(i))
           xsav(1,i)=xa(1,j)
           xsav(2,i)=xa(2,j)
           xsav(3,i)=xa(3,j)
      END DO
      DO i=nmnf,nmnl
           ii = i+nsn
           j=msr(i)
C msr < 0 <=> shooting nodes
           IF(j>0) THEN 
             xmin= min(xmin,xa(1,j))
             ymin= min(ymin,xa(2,j))
             zmin= min(zmin,xa(3,j))
             xmax= max(xmax,xa(1,j))
             ymax= max(ymax,xa(2,j))
             zmax= max(zmax,xa(3,j))
             xsav(1,ii)=xa(1,j)
             xsav(2,ii)=xa(2,j)
             xsav(3,ii)=xa(3,j)
           END IF
      END DO
 
      c_max = zero
      IF(icurv/=0)THEN
          DO i=1,nrtm_t
            xxx=max(xa(1,irect(1,i)),xa(1,irect(2,i)),
     .              xa(1,irect(3,i)),xa(1,irect(4,i)))
     .         -min(xa(1,irect(1,i)),xa(1,irect(2,i)),
     .              xa(1,irect(3,i)),xa(1,irect(4,i)))
            yyy=max(xa(2,irect(1,i)),xa(2,irect(2,i)),
     .              xa(2,irect(3,i)),xa(2,irect(4,i)))
     .         -min(xa(2,irect(1,i)),xa(2,irect(2,i)),
     .              xa(2,irect(3,i)),xa(2,irect(4,i)))
            zzz=max(xa(3,irect(1,i)),xa(3,irect(2,i)),
     .              xa(3,irect(3,i)),xa(3,irect(4,i)))
     .         -min(xa(3,irect(1,i)),xa(3,irect(2,i)),
     .              xa(3,irect(3,i)),xa(3,irect(4,i)))
            curv_max(i) = half * max(xxx,yyy,zzz)
            c_max = max(c_max,curv_max(i))
          ENDDO
        ELSE
          DO i=1,nrtm_t
            curv_max(i)=zero
          ENDDO
      ENDIF
 
      DO i=nsnef,nsnel
        j=abs(nsve(i))
        xsave(1,i)=xa(1,j)
        xsave(2,i)=xa(2,j)
        xsave(3,i)=xa(3,j)
      END DO
 
      DO i=nmnef,nmnel
        ii = i+nsne
        j=msre(i)
C msr < 0 <=> shooting nodes
        IF(j>0) THEN 
          xmine= min(xmine,xa(1,j))
          ymine= min(ymine,xa(2,j))
          zmine= min(zmine,xa(3,j))
          xmaxe= max(xmaxe,xa(1,j))
          ymaxe= max(ymaxe,xa(2,j))
          zmaxe= max(zmaxe,xa(3,j))
          xsave(1,ii)=xa(1,j)
          xsave(2,ii)=xa(2,j)
          xsave(3,ii)=xa(3,j)
        END IF
      END DO
 
 
      RETURN