diff options
Diffstat (limited to 'Donjon/data/pulseTHM_0d.save')
| -rwxr-xr-x | Donjon/data/pulseTHM_0d.save | 54 |
1 files changed, 54 insertions, 0 deletions
diff --git a/Donjon/data/pulseTHM_0d.save b/Donjon/data/pulseTHM_0d.save new file mode 100755 index 0000000..a3cb9f6 --- /dev/null +++ b/Donjon/data/pulseTHM_0d.save @@ -0,0 +1,54 @@ +#!/bin/sh +if [ $# = 0 ] + then + echo "usage: pulseTHM_0d.save directory" 1>&2 + exit 1 +fi +echo access pulseTHM_0d.save +MACH=`uname -s` +Sysx="`echo $MACH | cut -b -6`" +if [ $Sysx = "CYGWIN" ]; then + MACH=`uname -o` +elif [ $Sysx = "AIX" ]; then + MACH=`uname -s` +else + MACH=`uname -sm | sed 's/[ ]/_/'` +fi +if [ $Sysx = "Linux" -o $Sysx = "Darwin" ]; then +mkdir Rundir +cd Rundir +grep ">|#PKINS#:" ../pulseTHM_0d.result > input_file +head input_file +cat >conv.f90 <<'eof2' +program conv + parameter(maxtim=3000) + real time(maxtim),power(maxtim) + character(len=8) :: text + it=0 + 10 it=it+1 + if(it.gt.maxtim) then + print *,'maxtim overflow' + stop + endif + open(7,file='input_file') + open(8,file='matlab_pulseTHM_0d.m') + read(7,*,end=20) text,itt,time(it),power(it) + go to 10 + 20 close(7) + write(8,200,advance='no') (time(i),i=1,it-1) + write(8,'(2h ])') + write(8,300,advance='no') (power(i),i=1,it-1) + write(8,'(2h ])') + close(8) + 200 format(7h time=[,1p,10e12.4,4h .../(7x,10e12.4,4h ...)) + 300 format(8h power=[,1p,10e12.4,4h .../(8x,10e12.4,4h ...)) +end program conv +eof2 +gfortran conv.f90 +./a.out <input_file +rm input_file conv.f90 a.out +mv matlab_pulseTHM_0d.m $1/"$MACH" +cd .. +rm -r Rundir +fi +echo "pulseTHM_0d.save completed" |