!&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& Subroutine Bocovs_oc(H,im,jm,lsd) ! ****************************************************************** ! * * ! * Boundary Conditions for Height Points on Octacon * ! * * ! ****************************************************************** integer,Parameter::nm=14 integer::im,jm,lsd !----------------------------------------------------------------------- Dimension & H(0:im+1,0:jm+1,nm,lsd) !----------------------------------------------------------------------- lm=im imp1=im+1 jmp1=jm+1 im1=im-1 jm1=jm-1 ! ! Domain 1 ! Do l=1,lm-1 H(l ,0 ,1,:) = H(im-l,1, 7,:) H(l ,jm,1,:) = H(l ,1, 3,:) H(0 ,l ,1,:) = H(l ,1, 8,:) H(im,l ,1,:) = H(im-l,1, 6,:) End Do H(0 ,0 , 1,:)=0. H(im,0 , 1,:)=0. H(0 ,jm, 1,:) = H(im1,1, 2,:) H(im ,jm, 1,:) = H(1 ,1, 4,:) ! ! Domain 2 ! Do l=1,lm-1 H(l ,0 , 2,:)=H(im1,jm-l, 8,:) H(l ,jm, 2,:)=H(im1,l ,14,:) H(0 ,l , 2,:)=H(im1,l ,11,:) H(im,l , 2,:)=H(1 ,l , 3,:) End Do H(0 ,0 , 2,:)=0. H(im,0 , 2,:)=H(1 ,jm1, 1,:) H(0 ,jm, 2,:)=0. H(im,jm, 2,:)=H(1 ,1 , 5,:) ! ! Domain 3 ! Do l=1,lm-1 H(l ,0 , 3,:)=H(l ,jm1, 1,:) H(l ,jm, 3,:)=H(l ,1 , 5,:) H(0 ,l , 3,:)=H(im1,l , 2,:) H(im,l , 3,:)=H(1 ,l , 4,:) End Do H(0 ,0 , 3,:)=H(im1,1 , 8,:) H(im,0 , 3,:)=H(1 ,1 , 6,:) H(0 ,jm, 3,:)=H(im1,jm1,14,:) H(im,jm, 3,:)=H(1 ,jm1,12,:) ! ! Domain 4 ! Do l=1,lm-1 H(l ,0 , 4,:)=H(1 ,l , 6,:) H(l ,jm, 4,:)=H(1 ,jm-l,12,:) H(0 ,l , 4,:)=H(im1,l , 3,:) H(im,l , 4,:)=H(1 ,l , 9,:) End Do H(0 ,0 , 4,:)=H(im1,jm1, 1,:) H(im,0 , 4,:)=0. H(0 ,jm, 4,:)=H(im1,1 , 5,:) H(im,jm, 4,:)=0. ! ! Domain 5 ! Do l=1,lm-1 H(l ,0 , 5,:)=H(l ,jm1, 3,:) H(l ,jm, 5,:)=H(im-l,jm1,13,:) H(0 ,l , 5,:)=H(im-l,jm1,14,:) H(im,l , 5,:)=H(l ,jm1,12,:) End Do H(0 ,0 , 5,:)=H(im1,jm1, 2,:) H(im,0 , 5,:)=H(1 ,jm1, 4,:) H(0 ,jm, 5,:)=0. H(im,jm, 5,:)=0. ! ! Domain 6 ! Do l=1,lm-1 H(l ,0 , 6,:)=H(im1,jm-l, 1,:) H(l ,jm, 6,:)=H(l ,1 , 9,:) H(0 ,l , 6,:)=H(l ,1 , 4,:) H(im,l , 6,:)=H(1 ,l , 7,:) End Do H(0 ,0 , 6,:)=H(im1,1, 3,:) H(im,0 , 6,:)=0. H(0 ,jm, 6,:)=0. H(im,jm, 6,:)=H(1 ,1,10,:) ! ! Domain 7 ! Do l=1,lm-1 H(l ,0 , 7,:)=H(im-l,2, 1,:) H(l ,jm, 7,:)=H(l ,1,10,:) H(0 ,l , 7,:)=H(im1 ,l, 6,:) H(im,l , 7,:)=H(1 ,l, 8,:) End Do H(0 ,0 , 7,:)=0. H(im,0 , 7,:)=0. H(0 ,jm, 7,:)=H(im1,1, 9,:) H(im,jm, 7,:)=H(1 ,1,11,:) ! ! Domain 8 ! Do l=1,lm-1 H(l ,0 , 8,:)=H(1 ,l, 1,:) H(l ,jm, 8,:)=H(l ,1,11,:) H(0 ,l , 8,:)=H(im1 ,l, 7,:) H(im,l , 8,:)=H(im-l,1, 2,:) End Do H(0 ,0 , 8,:)=0 H(im,0 , 8,:)=H(1 ,1, 3,:) H(0 ,jm, 8,:)=H(im1,1,10,:) H(im,jm, 8,:)=0. ! ! Domain 9 ! Do l=1,lm-1 H(l ,0 , 9,:)=H(l ,jm1, 6,:) H(l ,jm, 9,:)=H(l ,1 ,12,:) H(0 ,l , 9,:)=H(im1,l , 4,:) H(im,l , 9,:)=H(1 ,l ,10,:) End Do H(0 ,0 , 9,:)=0 H(im,0 , 9,:)=H(1,jm1, 7,:) H(0 ,jm, 9,:)=0. H(im,jm, 9,:)=H(1,1 ,13,:) ! ! Domain 10 ! Do l=1,lm-1 H(l ,0 ,10,:)=H(l ,jm1, 7,:) H(l ,jm,10,:)=H(l ,1 ,13,:) H(0 ,l ,10,:)=H(im1,l , 9,:) H(im,l ,10,:)=H(1 ,l ,11,:) End Do H(0 ,0 ,10,:)=H(im1,jm1, 6,:) H(im,0 ,10,:)=H(1 ,jm1, 8,:) H(0 ,jm,10,:)=H(im1,1 ,12,:) H(im,jm,10,:)=H(1 ,1 ,14,:) ! ! Domain 11 ! Do l=1,lm-1 H(l ,0 ,11,:)=H(l ,jm1, 8,:) H(l ,jm,11,:)=H(l ,1 ,14,:) H(0 ,l ,11,:)=H(im1,l ,10,:) H(im,l ,11,:)=H(1 ,l , 2,:) End Do H(0 ,0 ,11,:)=H(im1,jm1, 7,:) H(im,0 ,11,:)=0. H(0 ,jm,11,:)=H(im1,1 ,13,:) H(im,jm,11,:)=0. ! ! Domain 12 ! Do l=1,lm-1 H(l ,0 ,12,:)=H(l ,jm1, 9,:) H(l ,jm,12,:)=H(im1,l , 5,:) H(0 ,l ,12,:)=H(im-l,jm1, 4,:) H(im,l ,12,:)=H(1 ,l ,13,:) End Do H(0 ,0 ,12,:)=0. H(im,0 ,12,:)=H(1 ,jm1,10,:) H(0 ,jm,12,:)=H(im1,im1, 3,:) H(im,jm,12,:)=0. ! ! Domain 13 ! Do l=1,lm-1 H(l ,0 ,13,:)=H(l ,jm1,10,:) H(l ,jm,13,:)=H(im-l,jm1, 5,:) H(0 ,l ,13,:)=H(im1 ,l ,12,:) H(im,l ,13,:)=H(1 ,l ,14,:) End Do H(0 ,0 ,13,:)=H(im1,jm1, 9,:) H(im,0 ,13,:)=H(1 ,jm1,11,:) H(0 ,jm,13,:)=0. H(im,jm,13,:)=0. ! ! Domain 14 ! Do l=1,lm-1 H(l ,0 ,14,:)=H(l ,jm1,11,:) H(l ,jm,14,:)=H(1 ,jm-l, 5,:) H(0 ,l ,14,:)=H(im1 ,l ,13,:) H(im,l ,14,:)=H(l ,jm1, 2,:) End Do H(0 ,0 ,14,:)=H(im1,jm1,10,:) H(im,0 ,14,:)=0. H(0 ,jm,14,:)=0. H(im,jm,14,:)=H(1 ,jm1, 3,:) !----------------------------------------------------------------------- Return End