From 44c15d39c716fa0a8b02b9950f011ae2cd354d26 Mon Sep 17 00:00:00 2001
From: Douglas Guptill <douglas.guptill@dal.ca>
Date: Mon, 29 Jun 2009 18:03:19 +0000
Subject: [PATCH] split long lines
---
slices.f90 | 84 +++++++++++++++++++++++++++++++-----------------------
1 file changed, 48 insertions(+), 36 deletions(-)
diff --git a/slices.f90 b/slices.f90
index 50bbcbd2..5d899c58 100644
--- a/slices.f90
+++ b/slices.f90
@@ -454,36 +454,36 @@ do ie=1,osolve%nleaves ! loop on all leaves
!*******************************************************************
case (1) !******************************************************
!*******************************************************************
- valpha1=( (section%slice-osolve%x(node1)) * ov%vnode(node1) + (osolve%x(node2)-section%slice) * ov%vnode(node2) )/maxfield
- valpha2=( (section%slice-osolve%x(node3)) * ov%vnode(node3) + (osolve%x(node4)-section%slice) * ov%vnode(node4) )/maxfield
- valpha3=( (section%slice-osolve%x(node5)) * ov%vnode(node5) + (osolve%x(node6)-section%slice) * ov%vnode(node6) )/maxfield
- valpha4=( (section%slice-osolve%x(node7)) * ov%vnode(node7) + (osolve%x(node8)-section%slice) * ov%vnode(node7) )/maxfield
- vbeta1 =( (section%slice-osolve%x(node1)) * ov%wnode(node1) + (osolve%x(node2)-section%slice) * ov%wnode(node2) )/maxfield
- vbeta2 =( (section%slice-osolve%x(node3)) * ov%wnode(node3) + (osolve%x(node4)-section%slice) * ov%wnode(node4) )/maxfield
- vbeta3 =( (section%slice-osolve%x(node5)) * ov%wnode(node5) + (osolve%x(node6)-section%slice) * ov%wnode(node6) )/maxfield
- vbeta4 =( (section%slice-osolve%x(node7)) * ov%wnode(node7) + (osolve%x(node8)-section%slice) * ov%wnode(node8) )/maxfield
+ valpha1=( (section%slice-osolve%x(node1)) * ov%vnode(node1) + (osolve%x(node2)-section%slice) * ov%vnode(node2) )/maxfield
+ valpha2=( (section%slice-osolve%x(node3)) * ov%vnode(node3) + (osolve%x(node4)-section%slice) * ov%vnode(node4) )/maxfield
+ valpha3=( (section%slice-osolve%x(node5)) * ov%vnode(node5) + (osolve%x(node6)-section%slice) * ov%vnode(node6) )/maxfield
+ valpha4=( (section%slice-osolve%x(node7)) * ov%vnode(node7) + (osolve%x(node8)-section%slice) * ov%vnode(node7) )/maxfield
+ vbeta1 =( (section%slice-osolve%x(node1)) * ov%wnode(node1) + (osolve%x(node2)-section%slice) * ov%wnode(node2) )/maxfield
+ vbeta2 =( (section%slice-osolve%x(node3)) * ov%wnode(node3) + (osolve%x(node4)-section%slice) * ov%wnode(node4) )/maxfield
+ vbeta3 =( (section%slice-osolve%x(node5)) * ov%wnode(node5) + (osolve%x(node6)-section%slice) * ov%wnode(node6) )/maxfield
+ vbeta4 =( (section%slice-osolve%x(node7)) * ov%wnode(node7) + (osolve%x(node8)-section%slice) * ov%wnode(node8) )/maxfield
!*******************************************************************
case (2) !******************************************************
!*******************************************************************
- valpha1=( (section%slice-osolve%y(node1)) * ov%unode(node1) + (osolve%y(node3)-section%slice) * ov%unode(node3) )/maxfield
- valpha2=( (section%slice-osolve%y(node2)) * ov%unode(node2) + (osolve%y(node4)-section%slice) * ov%unode(node4) )/maxfield
- valpha3=( (section%slice-osolve%y(node5)) * ov%unode(node5) + (osolve%y(node7)-section%slice) * ov%unode(node7) )/maxfield
- valpha4=( (section%slice-osolve%y(node6)) * ov%unode(node6) + (osolve%y(node8)-section%slice) * ov%unode(node8) )/maxfield
- vbeta1 =( (section%slice-osolve%y(node1)) * ov%wnode(node1) + (osolve%y(node3)-section%slice) * ov%wnode(node3) )/maxfield
- vbeta2 =( (section%slice-osolve%y(node2)) * ov%wnode(node2) + (osolve%y(node4)-section%slice) * ov%wnode(node4) )/maxfield
- vbeta3 =( (section%slice-osolve%y(node5)) * ov%wnode(node5) + (osolve%y(node7)-section%slice) * ov%wnode(node7) )/maxfield
- vbeta4 =( (section%slice-osolve%y(node6)) * ov%wnode(node6) + (osolve%y(node8)-section%slice) * ov%wnode(node8) )/maxfield
+ valpha1=( (section%slice-osolve%y(node1)) * ov%unode(node1) + (osolve%y(node3)-section%slice) * ov%unode(node3) )/maxfield
+ valpha2=( (section%slice-osolve%y(node2)) * ov%unode(node2) + (osolve%y(node4)-section%slice) * ov%unode(node4) )/maxfield
+ valpha3=( (section%slice-osolve%y(node5)) * ov%unode(node5) + (osolve%y(node7)-section%slice) * ov%unode(node7) )/maxfield
+ valpha4=( (section%slice-osolve%y(node6)) * ov%unode(node6) + (osolve%y(node8)-section%slice) * ov%unode(node8) )/maxfield
+ vbeta1 =( (section%slice-osolve%y(node1)) * ov%wnode(node1) + (osolve%y(node3)-section%slice) * ov%wnode(node3) )/maxfield
+ vbeta2 =( (section%slice-osolve%y(node2)) * ov%wnode(node2) + (osolve%y(node4)-section%slice) * ov%wnode(node4) )/maxfield
+ vbeta3 =( (section%slice-osolve%y(node5)) * ov%wnode(node5) + (osolve%y(node7)-section%slice) * ov%wnode(node7) )/maxfield
+ vbeta4 =( (section%slice-osolve%y(node6)) * ov%wnode(node6) + (osolve%y(node8)-section%slice) * ov%wnode(node8) )/maxfield
!*******************************************************************
case (3) !******************************************************
!*******************************************************************
- valpha1=( (section%slice-osolve%z(node1)) * ov%unode(node1) + (osolve%z(node5)-section%slice) * ov%unode(node5) )/maxfield
- valpha2=( (section%slice-osolve%z(node2)) * ov%unode(node2) + (osolve%z(node6)-section%slice) * ov%unode(node6) )/maxfield
- valpha3=( (section%slice-osolve%z(node3)) * ov%unode(node3) + (osolve%z(node7)-section%slice) * ov%unode(node7) )/maxfield
- valpha4=( (section%slice-osolve%z(node4)) * ov%unode(node4) + (osolve%z(node8)-section%slice) * ov%unode(node8) )/maxfield
- vbeta1 =( (section%slice-osolve%z(node1)) * ov%vnode(node1) + (osolve%z(node5)-section%slice) * ov%vnode(node5) )/maxfield
- vbeta2 =( (section%slice-osolve%z(node2)) * ov%vnode(node2) + (osolve%z(node6)-section%slice) * ov%vnode(node6) )/maxfield
- vbeta3 =( (section%slice-osolve%z(node3)) * ov%vnode(node3) + (osolve%z(node7)-section%slice) * ov%vnode(node7) )/maxfield
- vbeta4 =( (section%slice-osolve%z(node4)) * ov%vnode(node4) + (osolve%z(node8)-section%slice) * ov%vnode(node8) )/maxfield
+ valpha1=( (section%slice-osolve%z(node1)) * ov%unode(node1) + (osolve%z(node5)-section%slice) * ov%unode(node5) )/maxfield
+ valpha2=( (section%slice-osolve%z(node2)) * ov%unode(node2) + (osolve%z(node6)-section%slice) * ov%unode(node6) )/maxfield
+ valpha3=( (section%slice-osolve%z(node3)) * ov%unode(node3) + (osolve%z(node7)-section%slice) * ov%unode(node7) )/maxfield
+ valpha4=( (section%slice-osolve%z(node4)) * ov%unode(node4) + (osolve%z(node8)-section%slice) * ov%unode(node8) )/maxfield
+ vbeta1 =( (section%slice-osolve%z(node1)) * ov%vnode(node1) + (osolve%z(node5)-section%slice) * ov%vnode(node5) )/maxfield
+ vbeta2 =( (section%slice-osolve%z(node2)) * ov%vnode(node2) + (osolve%z(node6)-section%slice) * ov%vnode(node6) )/maxfield
+ vbeta3 =( (section%slice-osolve%z(node3)) * ov%vnode(node3) + (osolve%z(node7)-section%slice) * ov%vnode(node7) )/maxfield
+ vbeta4 =( (section%slice-osolve%z(node4)) * ov%vnode(node4) + (osolve%z(node8)-section%slice) * ov%vnode(node8) )/maxfield
end select
CALL PGSCI (1)
CALL PGARRO (alpha_pos,beta_pos,alpha_pos+valpha1,beta_pos+vbeta1)
@@ -504,30 +504,42 @@ do ie=1,osolve%nleaves ! loop on all leaves
case (1) !***************************************************
!****************************************************************
do ilsf=1,nlsf
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,section%slice,y0 ,z0 ,f1)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,section%slice,y0+dxyz,z0 ,f2)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,section%slice,y0 ,z0+dxyz,f3)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,section%slice,y0+dxyz,z0+dxyz,f4)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ section%slice,y0 ,z0 ,f1)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ section%slice,y0+dxyz,z0 ,f2)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ section%slice,y0 ,z0+dxyz,f3)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ section%slice,y0+dxyz,z0+dxyz,f4)
!call find_lsf_points (section%xyz,osolve,alpha_lsf(:,ilsf),beta_lsf(:,ilsf),ilsf,f1,f2,f3,f4,section%slice,x0,y0,z0,dxyz,lilflag(ilsf))
end do
!****************************************************************
case (2) !***************************************************
!****************************************************************
do ilsf=1,nlsf
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0 ,section%slice,z0 ,f1)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0+dxyz,section%slice,z0 ,f2)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0 ,section%slice,z0+dxyz,f3)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0+dxyz,section%slice,z0+dxyz,f4)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0 ,section%slice,z0 ,f1)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0+dxyz,section%slice,z0 ,f2)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0 ,section%slice,z0+dxyz,f3)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0+dxyz,section%slice,z0+dxyz,f4)
!call find_lsf_points (section%xyz,osolve,alpha_lsf(:,ilsf),beta_lsf(:,ilsf),ilsf,f1,f2,f3,f4,section%slice,x0,y0,z0,dxyz,lilflag(ilsf))
end do
!****************************************************************
case (3) !***************************************************
!****************************************************************
do ilsf=1,nlsf
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0 ,y0 ,section%slice,f1)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0+dxyz,y0 ,section%slice,f2)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0 ,y0+dxyz,section%slice,f3)
- call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,x0+dxyz,y0+dxyz,section%slice,f4)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0 ,y0 ,section%slice,f1)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0+dxyz,y0 ,section%slice,f2)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0 ,y0+dxyz,section%slice,f3)
+ call octree_interpolate (osolve%octree,osolve%noctree,osolve%icon,osolve%nleaves,osolve%lsf(:,ilsf),osolve%nnode,&
+ x0+dxyz,y0+dxyz,section%slice,f4)
!call find_lsf_points (section%xyz,osolve,alpha_lsf(:,ilsf),beta_lsf(:,ilsf),ilsf,f1,f2,f3,f4,section%slice,x0,y0,z0,dxyz,lilflag(ilsf))
end do
end select
--
GitLab