-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathProcessKITP.tar
More file actions
141 lines (125 loc) · 20 KB
/
ProcessKITP.tar
File metadata and controls
141 lines (125 loc) · 20 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
ProcessKITP/ 0042700 1567261 0127643 00000000000 11316765254 012131 5 ustar ux454321 uic221 ProcessKITP/process 0100700 1567261 0127643 00000002330 11316765237 013527 0 ustar ux454321 uic221 #!/bin/tcsh
if( $1 == "" )then
echo "process n"
echo " where n is the zero-padded dump number (i.e. 0001)"
echo " setups up all the Aake analysis tools."
exit
endif
set num = $1
# Datafile setup (.dim)
# Powerspectrum file setup.
# power spectrum computation
# structure function file setup
# structure function computation.
set dir = DD$num.products
if( !( -e $dir ) ) then
print "No directory $dir"
exit
endif
cd $dir
echo Analize: `pwd`
set power = 0
set sf = 0
set pdf = 1
set size = 256
@ snap = $num / 20
if( $power == 1)then
cp ../AakeFiles/base.dim ./data$num.dim
sed -e 's:input:'data$num.dat':' ../AakeFiles/g2p.in > ./g2p.in
date > start_power
g2p.x g2p.in
date > end_power
endif
if( $sf == 1) then
echo "structure function"
sed -e 's:input:'data$num.dat':' ../AakeFiles/g2sf.in > ./g2sf.in
date > start_sf
g2sf.x g2sf.in
date > end_sf
endif
cd -
cd DD$num
if( $pdf == 1) then
pwd
unpoo
set file = data$num
poe pdfgen.exe $file -nodes 1 -tasks_per_node 8
mv $file.PDF.log_density ../final_products/mhd"$size"snapshot"$snap"densityPDF.txt
mv $file.PDF.energy ../final_products/mhd"$size"snapshot"$snap"energyPDF.txt
endif
cd ..
#end