GET_PNT_POS_PGM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_pnt_pos_pgm (nb_el, aa11, aa12, aa13, aa21, aa22, aa23, aa31, aa32, aa33, bb11, bb12, bb13, bb21, bb22, bb23, bb31, bb32, bb33, gg1, gg2, gg3, dd1, dd2, dd3, x_true, y_true, z_true, el, xref, yref, zref,
	Findes node coordinates in the reference element.

Function/Subroutine Documentation

get_pnt_pos_pgm()

subroutine get_pnt_pos_pgm	(	integer*4	nb_el,
		real*8, dimension(nb_el)	aa11,
		real*8, dimension(nb_el)	aa12,
		real*8, dimension(nb_el)	aa13,
		real*8, dimension(nb_el)	aa21,
		real*8, dimension(nb_el)	aa22,
		real*8, dimension(nb_el)	aa23,
		real*8, dimension(nb_el)	aa31,
		real*8, dimension(nb_el)	aa32,
		real*8, dimension(nb_el)	aa33,
		real*8, dimension(nb_el)	bb11,
		real*8, dimension(nb_el)	bb12,
		real*8, dimension(nb_el)	bb13,
		real*8, dimension(nb_el)	bb21,
		real*8, dimension(nb_el)	bb22,
		real*8, dimension(nb_el)	bb23,
		real*8, dimension(nb_el)	bb31,
		real*8, dimension(nb_el)	bb32,
		real*8, dimension(nb_el)	bb33,
		real*8, dimension(nb_el)	gg1,
		real*8, dimension(nb_el)	gg2,
		real*8, dimension(nb_el)	gg3,
		real*8, dimension(nb_el)	dd1,
		real*8, dimension(nb_el)	dd2,
		real*8, dimension(nb_el)	dd3,
		real*8	x_true,
		real*8	y_true,
		real*8	z_true,
		integer*4	el,
		real*8	xref,
		real*8	yref,
		real*8	zref
	)

Findes node coordinates in the reference element.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_el	number of elements
[in]	A11	costant values for the bilinear map
[in]	A12	costant values for the bilinear map
[in]	A13	costant values for the bilinear map
[in]	A21	costant values for the bilinear map
[in]	A22	costant values for the bilinear map
[in]	A23	costant values for the bilinear map
[in]	A31	costant values for the bilinear map
[in]	A32	costant values for the bilinear map
[in]	A33	costant values for the bilinear map
[in]	B11	costant values for the bilinear map
[in]	B12	costant values for the bilinear map
[in]	B13	costant values for the bilinear map
[in]	B21	costant values for the bilinear map
[in]	B22	costant values for the bilinear map
[in]	B23	costant values for the bilinear map
[in]	B31	costant values for the bilinear map
[in]	B32	costant values for the bilinear map
[in]	B33	costant values for the bilinear map
[in]	GG1	costant values for the bilinear map
[in]	GG2	costant values for the bilinear map
[in]	GG3	costant values for the bilinear map
[in]	DD1	costant values for the bilinear map
[in]	DD2	costant values for the bilinear map
[in]	DD3	costant values for the bilinear map
[in]	x_true	x-coordinate of the node
[in]	y_true	y-coordinate of the node
[in]	z_true	z-coordinate of the node
[in]	highest	elevation of the elemnt
[in]	elevation	maximal elevation
[out]	el	element index
[out]	xref	x-coordinate in the refernce element
[out]	yref	y-coordinate in the refernce element
[out]	zref	z-coordinate in the refernce element

Definition at line 58 of file GET_PNT_POS_PGM.f90.

                                 highest,elevation)
      
      implicit none
      
      integer*4 :: nb_el,el
      integer*4 :: inr,ie
 
      real*8 :: x_true,y_true,z_true,xref,yref,zref
      real*8 :: csi,eta,zeta,delta_csi,delta_eta,delta_zeta,dist,dist_min
      real*8 :: termx,termy,termz,det_j
      real*8 :: dxdx,dxdy,dxdz
      real*8 :: dydx,dydy,dydz
      real*8 :: dzdx,dzdy,dzdz
      real*8 :: a11,a12,a13
      real*8 :: a21,a22,a23
      real*8 :: a31,a32,a33
      real*8 :: elevation
 
      real*8, dimension(nb_el) :: aa11,aa12,aa13
      real*8, dimension(nb_el) :: aa21,aa22,aa23
      real*8, dimension(nb_el) :: aa31,aa32,aa33
      real*8, dimension(nb_el) :: bb11,bb12,bb13
      real*8, dimension(nb_el) :: bb21,bb22,bb23
      real*8, dimension(nb_el) :: bb31,bb32,bb33
      real*8, dimension(nb_el) :: gg1,gg2,gg3
      real*8, dimension(nb_el) :: dd1,dd2,dd3
      real*8, dimension(nb_el) :: highest
      
      dist_min = 1.0e20
      xref = 0.0;      yref = 0.0;      zref = 0.0
      el = 1
      
      do ie = 1,nb_el
         if (highest(ie) .gt. elevation) then
            csi = 0.0;             eta = 0.0;     zeta = 0.0
            
            do inr = 1, 100
                dxdx = aa11(ie) + bb12(ie)*zeta + bb13(ie)*eta + gg1(ie)*eta*zeta
                dydx = aa21(ie) + bb22(ie)*zeta + bb23(ie)*eta + gg2(ie)*eta*zeta
                dzdx = aa31(ie) + bb32(ie)*zeta + bb33(ie)*eta + gg3(ie)*eta*zeta
          
                dxdy = aa12(ie) + bb11(ie)*zeta + bb13(ie)*csi + gg1(ie)*zeta*csi
                dydy = aa22(ie) + bb21(ie)*zeta + bb23(ie)*csi + gg2(ie)*zeta*csi
                dzdy = aa32(ie) + bb31(ie)*zeta + bb33(ie)*csi + gg3(ie)*zeta*csi
          
                dxdz = aa13(ie) + bb11(ie)*eta + bb12(ie)*csi  + gg1(ie)*csi*eta
                dydz = aa23(ie) + bb21(ie)*eta + bb22(ie)*csi  + gg2(ie)*csi*eta
                dzdz = aa33(ie) + bb31(ie)*eta + bb32(ie)*csi  + gg3(ie)*csi*eta
          
                det_j = dxdz * (dydx*dzdy - dzdx*dydy) - dydz * (dxdx*dzdy - dzdx*dxdy) &
                          + dzdz * (dxdx*dydy - dydx*dxdy)
          
          
                a11 = +dydy*dzdz -dydz*dzdy
                a12 = -dydx*dzdz +dydz*dzdx
                a13 = +dydx*dzdy -dydy*dzdx
          
                a21 = -dxdy*dzdz +dxdz*dzdy
                a22 = +dxdx*dzdz -dxdz*dzdx
                a23 = -dxdx*dzdy +dxdy*dzdx
          
                a31 = +dxdy*dydz -dxdz*dydy
                a32 = -dxdx*dydz +dxdz*dydx
                a33 = +dxdx*dydy -dxdy*dydx
          
                termx = aa11(ie)*csi + aa12(ie)*eta &
                   + aa13(ie)*zeta + bb11(ie)*eta*zeta &
                   + bb12(ie)*csi*zeta + bb13(ie)*csi*eta &
                   + gg1(ie)*csi*eta*zeta + dd1(ie) -x_true
          
                termy = aa21(ie)*csi + aa22(ie)*eta &
                   + aa23(ie)*zeta + bb21(ie)*eta*zeta &
                   + bb22(ie)*csi*zeta + bb23(ie)*csi*eta &
                   + gg2(ie)*csi*eta*zeta + dd2(ie) -y_true
            
                termz = aa31(ie)*csi + aa32(ie)*eta &
                   + aa33(ie)*zeta + bb31(ie)*eta*zeta &
                   + bb32(ie)*csi*zeta + bb33(ie)*csi*eta &
                   + gg3(ie)*csi*eta*zeta + dd3(ie) -z_true
            
                    delta_csi  = -(a11*termx +a21*termy +a31*termz) / det_j
                delta_eta  = -(a12*termx +a22*termy +a32*termz) / det_j
                delta_zeta = -(a13*termx +a23*termy +a33*termz) / det_j
            
                csi = csi + delta_csi
                eta = eta + delta_eta
                zeta = zeta + delta_zeta
                    
                    if (csi.lt.-1.0) csi = -1.0
                if (csi.gt. 1.0) csi =  1.0
                if (eta.lt.-1.0) eta = -1.0
                if (eta.gt. 1.0) eta =  1.0
                if (zeta.lt.-1.0) zeta = -1.0
                if (zeta.gt. 1.0) zeta =  1.0
            
            enddo
         
                 dist = termx*termx +termy*termy +termz*termz
         
                 if (dist.lt.dist_min) then
                    el = ie
                    xref = csi
                    yref = eta
                    zref = zeta
                    dist_min = dist
                 endif
 
        endif !if (highest(ie).gt.elevation) then
         
      enddo !do ie = 1,nb_el
      
      return
      

Referenced by find_monitor_position(), and read_system_position().

Here is the caller graph for this function: