tadd more simulations with different viscosity - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit 98d189ae93f391c456434bf8030619bc8d30d15e
 (DIR) parent 1df3a06a8dd92aef65a43ed30b7d7527c4750445
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 13 Apr 2015 14:11:28 +0200
       
       add more simulations with different viscosity
       
       Diffstat:
         M python/halfshear-darcy-strength-di… |       8 ++++----
         A python/halfshear-darcy-strength-di… |      53 ++++++++++++++++++++++++++++++
       
       2 files changed, 57 insertions(+), 4 deletions(-)
       ---
 (DIR) diff --git a/python/halfshear-darcy-strength-dilation-rate.py b/python/halfshear-darcy-strength-dilation-rate.py
       t@@ -25,8 +25,8 @@ k_c = 3.5e-15
        
        # 5.0e-8 results present
        #mu_f_vals = [1.797e-06, 1.204e-06, 5.0e-8, 1.797e-08]
       -#mu_f_vals = [1.797e-06, 1.204e-06, 3.594e-07, 1.797e-08]
       -mu_f_vals = [1.797e-06, 1.204e-06, 1.797e-08]
       +mu_f_vals = [1.797e-06, 1.204e-06, 3.594e-07, 1.797e-08]
       +#mu_f_vals = [1.797e-06, 1.204e-06, 1.797e-08]
        velfac = 1.0
        
        
       t@@ -209,8 +209,8 @@ ax1.legend(loc='upper right', prop={'size':18}, fancybox=True,
        #ax4.legend(loc='best', prop={'size':18}, fancybox=True,
                #framealpha=legend_alpha)
        
       -ax1.set_xlim([0.0, 0.09])
       -ax2.set_xlim([0.0, 0.09])
       +#ax1.set_xlim([0.0, 0.09])
       +#ax2.set_xlim([0.0, 0.09])
        #ax2.set_xlim([0.0, 0.2])
        
        ax1.set_ylim([-7, 45])
 (DIR) diff --git a/python/halfshear-darcy-strength-dilation-rate.sh b/python/halfshear-darcy-strength-dilation-rate.sh
       t@@ -0,0 +1,53 @@
       +#!/bin/sh
       +#PBS -N hs-d-r
       +#PBS -l nodes=1:ppn=3
       +#PBS -l walltime=48:00:00
       +#PBS -q qfermi
       +#PBS -M adc@geo.au.dk
       +#PBS -m abe
       +
       +# Grendel CUDA
       +source /com/gcc/4.6.4/load.sh
       +CUDAPATH=/com/cuda/5.5.22
       +export PATH=$HOME/bin:$PATH
       +export PATH=$CUDAPATH/bin:$PATH
       +export LD_LIBRARY_PATH=$CUDAPATH/lib64:$CUDAPATH/lib:$LD_LIBRARY_PATH
       +
       +# Manually installed Python modules
       +export PYTHONPATH=$HOME/.local/lib/python:$PYTHONPATH
       +export PYTHONPATH=$HOME/.local/lib64/python:$PYTHONPATH
       +
       +# Manually installed Python
       +#export PATH=/home/adc/.local/bin:$PATH
       +
       +# Shared Python2.7
       +PYTHON=/com/python/2.7.6
       +export PYTHONPATH=$PYTHON/lib:$PYTHONPATH
       +export PATH=$PYTHON/bin:$PATH
       +
       +echo "`whoami`@`hostname`"
       +echo "Start at `date`"
       +
       +ORIGDIR=/home/adc/code/sphere
       +#WORKDIR=/scratch/$PBS_JOBID
       +WORKDIR=$ORIGDIR
       +
       +#cp -r $ORIGDIR/* $WORKDIR
       +
       +cd $WORKDIR
       +nvidia-smi
       +rm CMakeCache.txt
       +cmake . && make
       +cd python
       +# $ python shear-starter.py <DEVICE> <FLUID> <C_PHI> <C_GRAD_P> <SIGMA_0>
       +#python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 20000.0 1.797e-6 1.0 & from hs-d-v1
       +#python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 20000.0 1.204e-6 1.0 & from iddqd
       +#python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 20000.0 3.594e-7 1.0 & from iddqd
       +python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 20000.0 6.110e-7 1.0 &
       +python halfshear-darcy-starter.py 1 1 1.0 3.5e-15 20000.0 1.797e-7 1.0 &
       +python halfshear-darcy-starter.py 2 1 1.0 3.5e-15 20000.0 1.797e-8 1.0 &
       +wait
       +
       +#cp $WORKDIR/output/* $ORIGDIR/output/
       +
       +echo "End at `date`"