starting test on geodyncomp

src/Makefile.inc.geodyncomp

+# Define our compilers.................................
+#
+# on geodyncomp, the defaults are OpenMPI (1.4.3), 
+# built with the GNU compilers
+F90 = mpif90
+F77 = mpif90
+cc  = mpicc
+
+# we take the standard ar
+AR=ar
+
+# There are no BITS
+BITS=
+
+# no MKL
+# IMKL = /scinet/gpc/intel/Compiler/11.1/056/mkl/include
+# IMKL = /scinet/gpc/intel/Compiler/11.1/072/mkl/include
+# INCLUDE=-I$(IMKL)
+
+# Define FORTRAN compiler flags.......................
+#
+F90FLAGS=-ffree-line-length-0
+
+# Define C compiler flags.............................
+#
+CFLAGS=
+
+# no Intel MKL library
+#
+# LMKL = /scinet/gpc/intel/Compiler/11.1/056/mkl/lib
+# LMKL = /scinet/gpc/intel/Compiler/11.1/072/mkl/lib
+
+# MPI libraries
+# (use the default)
+
+# WSMP Library
+#
+# WSMP = /home/dwhipp/software/wsmp/wsmp-Linux64/lib/Intel
+# WSMP = /home/dwhipp/software/wsmp/wsmp-Linux64-MPI2
+# WSMP = /home/dwhipp/software/wsmp/wsmp-Linux64-IMPI
+WSMP = -L/opt/wsmp/wsmp-Linux64/lib/GNU/openmpi -lpwsmp64
+
+# Not using Google's TCMalloc
+# TCML   = /home/dwhipp/software/tcml
+
+LIBS = \
+ -LOCTREE -lOctree$(BITS) \
+ -LNN -lnn_f$(BITS) \
+ -LNN -lnn_c$(BITS) \
+ -LCASCADE -lcascade$(BITS) \
+ -LRESAMPLE -lresample$(BITS) \
+ $(WSMP) -lpthread
+
+# The Makefiles look at compile output with this
+PAGER=more
+
+# compile rules.......................................
+#
+COMPILE_OUT=$(NAME)$(BITS).compile
+.SUFFIXES:
+.SUFFIXES: .o .c .f .f90
+.f90.o:
+	@echo "--"                                   >>$(COMPILE_OUT) 2>&1
+	@echo "$(F90) $(F90FLAGS) $(INCLUDES) -c $<" >>$(COMPILE_OUT) 2>&1
+	$(F90) $(F90FLAGS) $(INCLUDES) -c $<         >>$(COMPILE_OUT) 2>&1
+
+.f.o:
+	@echo "--"                                   >>$(COMPILE_OUT) 2>&1
+	@echo "$(F77) $(F77FLAGS) $(INCLUDES) -c $<" >>$(COMPILE_OUT) 2>&1
+	$(F77) $(F77FLAGS) $(INCLUDES) -c $<         >>$(COMPILE_OUT) 2>&1
+
+.c.o:
+	@echo "--"                                 >>$(COMPILE_OUT) 2>&1
+	@echo "$(CC) $(CFLAGS) -c $<"              >>$(COMPILE_OUT) 2>&1
+	$(CC) $(CFLAGS) -c $<                      >>$(COMPILE_OUT) 2>&1

src/define_bc_2011-01-21.f90

+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+!                                                                              |
+!              ||===\\                                                         | 
+!              ||    \\                                                        |
+!              ||     ||   //==\\   ||  ||   //==||  ||/==\\                   |
+!              ||     ||  ||    ||  ||  ||  ||   ||  ||    ||                  |
+!              ||    //   ||    ||  ||  ||  ||   ||  ||                        |
+!              ||===//     \\==//    \\==\\  \\==\\  ||                        |
+!                                                                              |
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+!                                                                              |
+!              DEFINE_BC   Feb. 2009                                           |
+!                                                                              |
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+
+subroutine define_bc (params,osolve,vo,zi) 
+
+!------------------------------------------------------------------------------|
+!((((((((((((((((  Purpose of the routine  )))))))))))))))))))))))))))))))))))))
+!------------------------------------------------------------------------------|
+! this routine assigns the boundary condition for the Tibet experiment
+
+!------------------------------------------------------------------------------|
+!((((((((((((((((  declaration of the subroutine arguments  ))))))))))))))))))))
+!------------------------------------------------------------------------------|
+
+use definitions
+
+implicit none
+
+type (parameters) params
+type (octreesolve) osolve
+type (void) vo
+type (ziso) zi
+
+!------------------------------------------------------------------------------|
+!(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
+!------------------------------------------------------------------------------|
+
+integer i,iproc,nproc,ierr
+double precision eps,lsf0,pi,lorig,h,x1,x2,phi,spu,spv,yend,cper,cscl,xstart
+double precision ystart,theta,l,vin,vzfluxscl,cntvel,dxy,xend,xsym,ymax,xwidth
+double precision ywidth
+double precision,dimension(:),allocatable :: x0,ldisp
+integer ie,ij,j,jp,nelemx,nelemz
+logical firstcall
+
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|
+
+INCLUDE 'mpif.h'
+
+call mpi_comm_size (mpi_comm_world,nproc,ierr)
+call mpi_comm_rank (mpi_comm_world,iproc,ierr)
+
+eps=1.d-10
+osolve%kfix=0
+osolve%kfixt=0
+pi=atan(1.d0)*4.d0
+
+select case(trim(params%infile))
+
+case('input.txt','input.small.txt')
+   do i=1,osolve%nnode
+     if (osolve%x(i).lt.eps) then
+       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+     endif
+
+     if (osolve%x(i).gt.1.d0-eps) then
+       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+     endif
+
+     if (osolve%y(i).lt.eps) then
+!       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+!       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+     endif
+
+     if (osolve%y(i).gt.1.d0-eps) then
+!       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+!       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+     endif
+
+     if (osolve%z(i).lt.eps) then
+       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+       osolve%kfixt(i)=1
+       osolve%temp(i)=1.d0
+     endif
+
+     if (osolve%z(i).gt.1.d0-eps) then
+       osolve%kfix((i-1)*3+1)=1 ; osolve%u(i)=0.d0
+       osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+       osolve%kfix((i-1)*3+3)=1 ; osolve%w(i)=0.d0
+       osolve%kfixt(i)=1
+       osolve%temp(i)=0.d0
+     endif
+
+     if (.not.vo%influid(i)) then
+       osolve%kfixt(i)=1
+       osolve%temp(i)=0.d0
+     endif
+
+     ! Set velocity B/C everywhere to disallow flow in the y-direction
+     osolve%kfix((i-1)*3+2)=1 ; osolve%v(i)=0.d0
+   enddo
+
+case ('input.jgr')
+   call define_bc_jgr (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.sphere')
+   call define_bc_sphere (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.3Dpunch')
+   call define_bc_3Dpunch (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.2Dpunch')
+   call define_bc_2Dpunch (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.2Dpunch_vert')
+   call define_bc_2Dpunch_vert (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.folding')
+   call define_bc_folding (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.pipo')
+   call define_bc_pipo (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.ritske')
+   call define_bc_ritske (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.ritske_isurf')
+   call define_bc_ritske_isurf (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.subduction')
+   call define_bc_subduction (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.riedel')
+   call define_bc_riedel (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case ('input.parallipipede')
+   call define_bc_parallipipede (osolve%nnode,osolve%kfix,osolve%kfixt,osolve%x,osolve%y,osolve%z,osolve%u,osolve%v,osolve%w,osolve%temp,vo) 
+
+case default
+   if (iproc.eq.0) print *,params%infile
+   call stop_run (' pb with input file in define_bc.f90')
+
+end select
+
+! here we need to impose that the bc satisfy the bad faces constrains
+
+do ie=1,osolve%nface
+  do j=1,4
+  jp=1+mod(j,4)
+  i=osolve%iface(j+4,ie)
+  if (osolve%kfix((i-1)*3+1).eq.1) osolve%u(i)=(osolve%u(osolve%iface(j,ie))+osolve%u(osolve%iface(jp,ie)))/2.d0
+  if (osolve%kfix((i-1)*3+2).eq.1) osolve%v(i)=(osolve%v(osolve%iface(j,ie))+osolve%v(osolve%iface(jp,ie)))/2.d0
+  if (osolve%kfix((i-1)*3+3).eq.1) osolve%w(i)=(osolve%w(osolve%iface(j,ie))+osolve%w(osolve%iface(jp,ie)))/2.d0
+  if (osolve%kfixt(i).eq.1) osolve%temp(i)=(osolve%temp(osolve%iface(j,ie))+osolve%temp(osolve%iface(jp,ie)))/2.d0
+  enddo
+i=osolve%iface(9,ie)
+if (osolve%kfix((i-1)*3+1).eq.1) osolve%u(i)=sum(osolve%u(osolve%iface(1:4,ie)))/4.d0
+if (osolve%kfix((i-1)*3+2).eq.1) osolve%v(i)=sum(osolve%v(osolve%iface(1:4,ie)))/4.d0
+if (osolve%kfix((i-1)*3+3).eq.1) osolve%w(i)=sum(osolve%w(osolve%iface(1:4,ie)))/4.d0
+if (osolve%kfixt(i).eq.1) osolve%temp(i)=sum(osolve%temp(osolve%iface(1:4,ie)))/4.d0
+enddo
+
+if (params%debug==2) call output_bc (osolve)
+
+return
+end
+!------------------------------------------------------------------------------|
+!------------------------------------------------------------------------------|