TRS command file

In the following example, the ``EAB'' neutron configuration is specified (codes: E=1, A=13, B=9). Lines to edit are indicated by ``<<<<<<<<''. The output file hg193pt.tab lists many variables, while the results are shown pictorially as contour plots in the hg193pt.ps (postscript) file. The output table is wide and should be printed with the command ``mpage -PQMS860 -b -l -2'' (set up an alias!).

#!/bin/csh
setenv TRSDAT  /net/nsb/spec5/esp/trs/meshes/
setenv TRSBIN  /net/nsb/spec5/esp/trs/trsource/
setenv TRSPLOT /net/nsb/spec5/esp/trs/ok/
#
rm -f hg193pt.dat                 <<<<<<<< globally edit ``hg193pt'' filename
rm -f hg193pt.tab
rm -f hg193pt.log
#
echo 'Checking TIME'
date
#
$TRSBIN'gamlatzn' << temp0
 &INDATA NX=9, NY=11, N4=4, ILOM=18,  NNMIN=94, NNMAX= 118,
  NZMIN=72, NZMAX=82, IZ0=76, IN0=104, IXORIG=0, IYORIG=5 &

 &RUN IMESH=8, IZ=80, IN=113 &             <<<<<<< Z, N
3                                          <<<<<<< number of neutrons
1                                          <<<<<<< code for orbital
13                                         <<<<<<< code for orbital
9                                          <<<<<<< code for orbital
0                                          <<<<<<< number of protons
hg193pt.dat
QUIT
#
echo 'End of gamlatzn !!!!!!!!'
skip:
#
rm -f for007.dat
mv -f for014.dat hg193pt.tab
rm -f for015.dat
#
date
#
#two sheets will be created by this commandfile
$TRSPLOT'fixdat'
hg193pt.dat
temp0
#
mv -f temp.dat hg193pt.dat
rm -f for007.dat
mv -f for014.dat hg193pt.tab
rm -f for015.dat
date
$TRSPLOT'fixdat' << temp2
hg193pt.dat
temp2
mv -f temp.dat hg193pt.dat
#
plot:
ln -s $TRSPLOT'init_laser.ps' fort.28
#
$TRSPLOT'newcsm4_ptss.com' << end2
hg193pt.dat
18                              !number of frequencies produced
1                               !first frequency to be plotted
4                               !last frequency to be plotted
1                               !start with the first
R                               !plot routhians
1                               !step in one

end2
mv -f fort.29 f1hg193pt.ps
#
$TRSPLOT'newcsm4_ptss.com' << end3
hg193pt.dat
18     
5
8
1
R
1

end3
mv -f fort.29 f2hg193pt.ps
rm -f hg193pt.ps
cat f1hg193pt.ps > hg193pt.ps
cat f2hg193pt.ps >> hg193pt.ps
rm -f f1hg193pt.ps
rm -f f2hg193pt.ps
rm -f fort.28
rm -f PLO.TMP
echo 'Checking TIME'
date
done:
exit
$exit

Document Last Modified