v1.2.0 (Compatibilidad con IRIS 2020)

Uso directo de archivos SAC y SACPZ y cambios menores.

1.0_rdseed2sac_v1.0.csh

@@ -29,6 +29,7 @@ rm *.50.BH*.SAC
 rm *.60.BH*.SAC
 rm *.70.BH*.SAC
 rm *.80.BH*.SAC
+rm *.XX.BH*.SAC
 rm RESP*.01.BH*
 rm RESP*.02.BH*
 rm RESP*.03.BH*
@@ -42,6 +43,7 @@ rm RESP*.50.BH*
 rm RESP*.60.BH*
 rm RESP*.70.BH*
 rm RESP*.80.BH*
+rm RESP*.XX.BH*
 
 #$rdseed -f pisco2007.seed -p -d -o 1
 #rm *.10.BH*.SAC

1.0_rm_sac_inv_v1.0.csh

@@ -0,0 +1,63 @@
+#!/bin/csh -f
+# Script para cambiar de formato SEED a SAC
+# Nestor Luna - 20 Jan 2019
+
+echo '**************Script para eliminar archivos SAC no usados*************'
+echo '************Creado por: Nestor Luna Diaz, 20 de enero de 2019*********'
+echo '**************Ultima modificacion, 11 de julio de 2020****************'
+echo '**********************************************************************'
+echo ''
+
+# Ejecutable de rdseed. Modificarlo de ser necesario.
+#set rdseed='/opt/rdseed/rdseed.rh6.linux_64'
+
+#set archivo_seed = '2007-08-15_mww8.0_Pisco.seed'
+
+#$rdseed -f $1 -R -d -o 1
+#echo '->Archivos a eliminar (registros de aceleración y redundantes):'
+rm *.01.BH*.SAC
+rm *.02.BH*.SAC
+rm *.03.BH*.SAC
+rm *.04.BH*.SAC
+rm *.05.BH*.SAC
+rm *.10.BH*.SAC
+rm *.20.BH*.SAC
+rm *.30.BH*.SAC
+rm *.40.BH*.SAC
+rm *.50.BH*.SAC
+rm *.60.BH*.SAC
+rm *.70.BH*.SAC
+rm *.80.BH*.SAC
+rm *.XX.BH*.SAC
+rm SACPZ*.01.BH*
+rm SACPZ*.02.BH*
+rm SACPZ*.03.BH*
+rm SACPZ*.04.BH*
+rm SACPZ*.05.BH*
+rm SACPZ*.10.BH*
+rm SACPZ*.20.BH*
+rm SACPZ*.30.BH*
+rm SACPZ*.40.BH*
+rm SACPZ*.50.BH*
+rm SACPZ*.60.BH*
+rm SACPZ*.70.BH*
+rm SACPZ*.80.BH*
+rm SACPZ*.XX.BH*
+
+#$rdseed -f pisco2007.seed -p -d -o 1
+#rm *.10.BH*.SAC
+#rm SAC_PZs*_10_*99999
+
+echo ''
+echo '*******************************Resumen********************************'
+echo ''
+echo '->Archivos ".SAC" (formato SAC binario) a utilizar:'
+ls *.SAC | wc -l
+ls *.SAC
+echo ''
+echo '->Archivos "SACPZ*.BH*" (Polos y zeros) a utilizar:'
+ls SACPZ*.BH* | wc -l
+ls SACPZ*.BH*
+
+echo ''
+echo '**************************Fin del programa****************************'

1.1_renombrar_BHZ_v1.0.sh

@@ -5,38 +5,86 @@
 
 echo '**************Renombrar archivos SAC a NET.STN.BH[Z].SAC**************'
 echo '*****************Creado por: Zhigang Peng, 2014***********************'
-echo '******Ultima modificacion: Nestor Luna Diaz, 23 de enero de 2019******'
+echo '******Ultima modificacion: Nestor Luna Diaz, 30 de junio de 2020******'
 echo '**********************************************************************'
 echo ''
 
 backup_renombrar='backup_renombrar'
 STN_ZC='stn_zc.id'
+PZ_ZC='pz_zc.id'
 
-#echo '->Archivos a eliminar (otras componentes):'
-rm *.BH1*.SAC *.BH2*.SAC
-rm *.BHE*.SAC *.BHN*.SAC
-rm RESP.*BH1 RESP.*BH2
-rm RESP.*BHE RESP.*BHN
-echo ''
+if [ $1 -eq 1 ];
+then
+    #echo '->Archivos a eliminar (otras componentes):'
+    rm *.BH1*.SAC *.BH2*.SAC
+    rm *.BHE*.SAC *.BHN*.SAC
+    rm RESP.*BH1 RESP.*BH2
+    rm RESP.*BHE RESP.*BHN
+    echo ''
 
-echo '->Lista de archivos antes de renombrar:'
-ls *.SAC | wc -l
-ls *.SAC
+    echo '->Lista de archivos antes de renombrar:'
+    ls *.SAC | wc -l
+    ls *.SAC
 
-rm -r $backup_renombrar
-mkdir $backup_renombrar
-cp *.SAC $backup_renombrar
+    rm -r $backup_renombrar
+    mkdir $backup_renombrar
+    cp *.SAC $backup_renombrar
 
-ls *.SAC | awk -F"." '{print $7,$8}' | sort | uniq -c | awk '{if ($1 == 1) print $2,$3}' | sort | uniq > $STN_ZC
-# find the station name which has three seismograms (three components)
-ls `awk '{print "*"$1"."$2".*SAC"}' $STN_ZC` | awk -F"." '{print "mv "$0,$7"."$8"."$10"."$12}' | sh
-# rename the SAC filename to NET.STN.BH[ZNE].SAC
-mv $STN_ZC $backup_renombrar
+    ls *.SAC | awk -F"." '{print $7,$8}' | sort | uniq -c | awk '{if ($1 == 1) print $2,$3}' | sort | uniq > $STN_ZC
+    # find the station name which has three seismograms (three components)
+    ls `awk '{print "*"$1"."$2".*SAC"}' $STN_ZC` | awk -F"." '{print "mv "$0,$7"."$8"."$10"."$12}' | sh
+    # rename the SAC filename to NET.STN.BH[ZNE].SAC
+
+    mv $STN_ZC $backup_renombrar
+
+    echo ''
+    echo '->Lista de archivos despues de renombrar:'
+    ls *.SAC | wc -l
+    ls *.SAC
+fi
+
+if [ $1 -eq 2 ];
+then
+    #echo '->Archivos a eliminar (otras componentes):'
+    rm *.BH1*.SAC *.BH2*.SAC
+    rm *.BHE*.SAC *.BHN*.SAC
+    rm SACPZ.*BH1 SACPZ.*BH2
+    rm SACPZ.*BHE SACPZ.*BHN
+    echo ''
+
+    echo '->Lista de archivos antes de renombrar:'
+    ls *.SAC | wc -l
+    ls *.SAC
+    ls SACPZ.* | wc -l
+    ls SACPZ.*
+
+    rm -r $backup_renombrar
+    mkdir $backup_renombrar
+    cp *.SAC $backup_renombrar
+    cp SACPZ.* $backup_renombrar
+
+    ls *.SAC | awk -F"." '{print $1,$2}' | sort | uniq -c | awk '{if ($1 == 1) print $2,$3}' | sort | uniq > $STN_ZC
+    # find the station name which has three seismograms (three components)
+    ls `awk '{print "*"$1"."$2".*SAC"}' $STN_ZC` | awk -F"." '{print "mv "$0,$1"."$2"."$4"."$9}' | sh
+    # rename the SAC filename to NET.STN.BH[ZNE].SAC
+
+    ls SACPZ.* | awk -F"." '{print $2,$3}' | sort | uniq -c | awk '{if ($1 == 1) print $2,$3}' | sort | uniq > $PZ_ZC
+    # find the station name which has three seismograms (three components)
+    ls `awk '{print "*"$1"."$2".*BHZ"}' $PZ_ZC` | awk -F"." '{print "mv "$0,$2"."$3"."$5".SAC.PZ"}' | sh
+    # rename the SACPZ filename to NET.STN.BH[ZNE].SAC.PZ
+
+    mv $STN_ZC $backup_renombrar
+    mv $PZ_ZC $backup_renombrar
+
+    echo ''
+    echo '->Lista de archivos despues de renombrar:'
+    ls *.SAC | wc -l
+    ls *.SAC
+    echo ''
+    ls *.SAC.PZ | wc -l
+    ls *.SAC.PZ
+fi
 
-echo ''
-echo '->Lista de archivos despues de renombrar:'
-ls *.SAC | wc -l
-ls *.SAC
 
 echo ''
 echo '**************************Fin del programa****************************'

4.0_deconvolucion_integrando_v1.5.csh

@@ -3,16 +3,17 @@
 echo '************Deconvolucion - Remover respuesta instrumental************'
 echo '************rdseed 5.3.1 & SAC 101.6a utilizados**********************'
 echo '************Creado por: Nestor Luna Diaz, 5 de octubre de 2016********'
-echo '***************Ultima modificacion, 01 de febrero de 2019*************'
+echo '***************Ultima modificacion, 01 de julio de 2020***************'
 echo '**********************************************************************'
 echo ''
-echo "SELECCION: f1 = $1, f2 = $2, Np = $3"
+echo "SELECCION: f1 = $2, f2 = $3, Np = $4"
 echo ''
 
 set N = `ls *.SAC | wc -l`
 set i = 1
 
 setenv SAC_DISPLAY_COPYRIGHT 0
+if ($1 == 1) then
 foreach file ( *.SAC )
 echo "  -> ($i/$N). Archivo SAC: $file"
 @ i++
@@ -48,7 +49,7 @@ echo "  -> ($i/$N). Archivo SAC: $file"
 #   *mul 1.0e-9
 #  **********************************************************************
 #  ********* Filtro Pasa-Banda, fc=[f1 f2], Butterworth n polos: ********
-   bp n $3 corner $1 $2
+   bp n $4 corner $2 $3
 #  **********************************************************************
 #  ********* Correcciones generales de señales 3: ***********************
    rmean
@@ -60,6 +61,58 @@ echo "  -> ($i/$N). Archivo SAC: $file"
 #  **********************************************************************
    q
 EOF
+endif
+
+if ($1 == 2) then
+foreach file ( *.SAC )
+echo "  -> ($i/$N). Archivo SAC: $file"
+@ i++
+   sac <<EOF
+#  **********************************************************************
+#  ********* [conf.] Desactivar "Quick and Dirty Plot" mode: ************
+   qdp off
+#  **********************************************************************
+#  ***** [conf.] Frecuencias de corte fc para el filtro Pasa-Banda: *****
+#   *SETBB f1 0.002
+#   *SETBB f2 0.1
+#   *SETBB np 5
+#  **********************************************************************
+#  ***** Leer archivo SAC: **********************************************
+   r $file
+#  **********************************************************************
+#  ********* Correcciones generales de señales 1: ***********************
+   rmean
+   rtrend
+   taper
+#  **********************************************************************
+#  ********* Remover respuesta instrumental - deconvolucion, ************
+#  ********* convertir a Desp: 'NONE', Vel: 'VEL', Acel: 'ACC' **********
+   trans from polezero subtype ${file}.PZ to vel
+#  **********************************************************************
+#  ********* Correcciones generales de señales 2: ***********************
+   rmean
+   rtrend
+   taper
+#  **********************************************************************
+#  *********** Integrando a desplazamiento en nm: ***********************
+   int
+   mul 1.0e9
+#  **********************************************************************
+#  ********* Filtro Pasa-Banda, fc=[f1 f2], Butterworth n polos: ********
+   bp n $4 corner $2 $3
+#  **********************************************************************
+#  ********* Correcciones generales de señales 3: ***********************
+   rmean
+   rtrend
+   taper
+#  **********************************************************************
+#  ********* Escribe archivo SAC filtrado: ******************************
+   w ${file}.deconv
+#  **********************************************************************
+   q
+EOF
+endif
+
 echo ''
 
 end