#!/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"