-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path4.1_graficar_SAC_single_v1.0.csh
executable file
·121 lines (108 loc) · 4.06 KB
/
4.1_graficar_SAC_single_v1.0.csh
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
#!/bin/csh -f
echo '************Graficar - señal original y deconvolucionada*************'
echo '************rdseed 5.3.1 & SAC 101.6a utilizados**********************'
echo '************Creado por: Nestor Luna Diaz, 5 de octubre de 2016********'
echo '*****************Ultima modificacion, 23 de marzo de 2019*************'
echo '**********************************************************************'
echo ''
set CARPETA_GRAF = 'Graficas_SAC_single_01'
rm -r $CARPETA_GRAF
mkdir $CARPETA_GRAF
set N = `ls *.SAC | wc -l`
set i = 1
setenv SAC_DISPLAY_COPYRIGHT 0
foreach file ( *.SAC )
echo " -> ($i/$N). Archivo SAC: $file. (Bucle 1/2)"
@ i++
sac <<EOF
# **********************************************************************
# ********* [conf.] Desactivar "Quick and Dirty Plot" mode: ************
qdp off
# **********************************************************************
# ********* [conf.] Remover archivos sgf: ******************************
sc rm -f *.sgf
# **********************************************************************
# ********* [conf.] Relación de aspecto y/x: ***************************
vspace 0.25
# **********************************************************************
# ***** Leer archivo SAC: **********************************************
r $file
# **********************************************************************
# ********* Correcciones generales de señales 1: ***********************
*rmean
*rtrend
*taper
# **********************************************************************
# ********* Etiquetas de ejes, grillas y color: ************************
xlabel "Tiempo (s)"
ylabel "Amplitud (Cuentas) Vel."
grid on
color blue
# **********************************************************************
# *************** Convertir sgf a ps: **********************************
bd sgf
p
saveimg ${file}.original.ps
# **********************************************************************
# ********* Remover archivos sgf: **************************************
sc rm -f *.sgf
# **********************************************************************
q
EOF
ps2eps ${file}.original.ps
echo ''
end
echo '-------------------------------------------------------------------'
echo ''
set N = `ls *.deconv | wc -l`
set i = 1
foreach file ( *.deconv )
echo " -> ($i/$N). Archivo SAC: $file. (Bucle 2/2)"
@ i++
sac <<EOF
# **********************************************************************
# ********* [conf.] Desactivar "Quick and Dirty Plot" mode: ************
qdp off
# **********************************************************************
# ********* [conf.] Remover archivos sgf: ******************************
sc rm -f *.sgf
# **********************************************************************
# ********* [conf.] Relación de aspecto y/x: ***************************
vspace 0.25
# **********************************************************************
# ***** Leer archivo SAC: **********************************************
r $file
# **********************************************************************
# ********* Etiquetas de ejes grillas y color: *************************
xlabel "Tiempo (s)"
ylabel "Amplitud (nm) Desp."
grid on
color red
# **********************************************************************
# *************** Convertir sgf a ps: **********************************
bd sgf
p
saveimg ${file}.ps
# **********************************************************************
# ********* Remover archivos sgf: **************************************
sc rm -f *.sgf
# **********************************************************************
q
EOF
ps2eps ${file}.ps
echo ''
end
rm *.ps
echo ''
echo '*******************************Resumen********************************'
echo ''
echo '->Archivos "original.eps" creados:'
ls *.original.eps | wc -l
ls *.original.eps
echo ''
echo '->Archivos "deconv.eps" creados:'
ls *.deconv.eps | wc -l
ls *.deconv.eps
mv *.eps $CARPETA_GRAF
echo ''
echo '**************************Fin del programa****************************'