Multi-Processor age Calculation

Self-Scheduling Model for multiple processors.

src/Makefile.inc.esd.argand

@@ -22,10 +22,10 @@ INCLUDE = $(MPI_INC)
 # Define compiler flags:
 
 #F90FLAGS = -O3 -Wall -Werror -Warray-bounds -Wunderflow -fbacktrace -fdump-core -ffree-line-length-none -fcheck=all
-F90FLAGS = -O3 -Wall -Warray-bounds -Wunderflow -fbacktrace -fdump-core -ffree-line-length-none -fcheck=all
+F90FLAGS = -O2 -Wall -Warray-bounds -Wunderflow -fbacktrace -fdump-core -ffree-line-length-none -fcheck=all
 F77FLAGS = $(F90FLAGS)
 # CFLAGS = -O3 -Wall -Werror
-CFLAGS = -O3 -Wall
+CFLAGS = -O2 -Wall
 WSMP = -L ../../wsmp/wsmp-Linux64-GNU/lib/openmpi/ -lpwsmp64
 
 

src/TRPOST/input_tracking_outputs.txt

-'/wrk/schmiddu/DOUAR_trackT5/TROUT/' !directory for tracking output files, ending with a /
-'/homeappl/home/schmiddu/runs/VTK/DOUAR_trackT5/TROUT/' !directory for tracking output files, ending with a /
-3 !Output mode: 1=only ages, 2=ages and surfaced tracks, 3=all tracks
+'/data/data3/mschmiddunser/model_runs/DOUAR/DOUAR_trackT5/' !Simulation Directory
+'/data/data3/mschmiddunser/data_large/DOUAR/DOUAR_trackT5/' !directory for output files
+2 !Output mode: 1=only ages, 2=ages and surfaced tracks, 3=all tracks
 0 !Scale Output? 0=Douar Units, 1=Natural Units 
-10 !First timestep for which tracking output will be calculated
-10 !Last timestep for which tracking output will be calculated
-20 !Step width for tracking output
-0 !acuumulate surfaced particles from n timesteps
-3 !only work on surfaced tracks that are at least n timesteps long
-2 !First timestep that has tracking input files
-2 !Step width of tracking input files to be read
+150 !First timestep for which tracking output will be calculated
+150 !Last timestep for which tracking output will be calculated
+30 !Step width for tracking output
+29 !acuumulate surfaced particles from n timesteps
+1 !only work on surfaced tracks that are at least n timesteps long
+10 !First timestep that has tracking input files
+20 !Step width of tracking input files to be read
 !Age prediction options
+15.0 !pad temperatures back to this age (in Ma)
+1 !Write AHe Ages
+1 !Write ZHe Ages
+1 !Write AFT Ages        
+1 !Write ZFT Ages        
+1 !Write Muscovite Ages  
 !Track options
 1 !write velocity
 1 !write z velocity

src/TRPOST/tridag.f90

+module m_tridag
+contains
+  subroutine tridag(a,b,c,r,u,n)
+    implicit none
+
+! inpout + output parameter
+    integer(4), intent(in) :: n
+
+    real(8), intent(out), dimension(n) :: u
+
+    real(8), intent(in), dimension(n) :: a, b, c, r
+
+
+! local vars
+    real(8) :: bet
+    real(8), dimension(:), allocatable :: gam
+
+    integer(4) :: j
+
+    allocate(gam(n))
+
+    if (b(1) == 0) then
+       print *, "tridag.f90: b(1) == 0"
+       stop
+    end if
+
+    bet = b(1)
+    u(1)=r(1)/bet
+
+    do j=2,n
+       gam(j)=c(j-1)/bet
+       bet=b(j)-a(j)*gam(j)
+       if (bet == 0) then
+          print *, "tridag.f90: bet == 0"
+          stop
+       end if
+       u(j)=(r(j)-a(j)*u(j-1))/bet
+    end do
+
+    do j=n-1,1,-1
+       u(j)=u(j)-gam(j+1)*u(j+1)
+    end do
+
+    deallocate(gam)
+
+  end subroutine tridag
+end module m_tridag

test/input_schmiddu.txt

@@ -1356,7 +1356,7 @@ approximately 0.25*ntrpx*ntrpy*ntrpz.
       ntrpy = 5
       ntrpz = 3
 
-If gridtrack is false [int]ntrp defines the total number of particles to be
+If gridtrack is false, [int]ntrp defines the total number of particles to be
 injected.
 
       ntrp = 25