summaryrefslogtreecommitdiff
path: root/Donjon/data/pulseTHM_0d.save
diff options
context:
space:
mode:
authorstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
committerstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
commit7dfcc480ba1e19bd3232349fc733caef94034292 (patch)
tree03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/data/pulseTHM_0d.save
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/data/pulseTHM_0d.save')
-rwxr-xr-xDonjon/data/pulseTHM_0d.save54
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"