-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy path031_get_mseed-sac_parallel.sh
executable file
·141 lines (133 loc) · 5.78 KB
/
031_get_mseed-sac_parallel.sh
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
#!/bin/bash
#
#
# (c) 2024 Valentino Lauciani <valentino.lauciani@ingv.it>,
# Istituto Nazione di Geofisica e Vulcanologia.
#
#####################################################
# Import config file
. $(dirname $0)/config.sh
### START - Check parameters ###
TYPE=
INPUT_STRING=
while getopts :o:k:d:u:t:s:e: OPTION
do
case ${OPTION} in
o) FILE_OUTPUT_MSEED="${OPTARG}"
;;
d) FILE_OUTPUT_DLESS="${OPTARG}"
;;
k) INPUT_STRING="${OPTARG}"
;;
t) TYPE="${OPTARG}"
;;
u) DATASELECT_URL="${OPTARG}"
;;
s) STARTTIME="${OPTARG}"
;;
e) ENDTIME="${OPTARG}"
;;
\?) echo "Invalid option: -$OPTARG" >/dev/null
shift
;;
*) #echo $OPTARG >/dev/null
echo "Invalid OPTARG: -$OPTARG" >&2
;;
esac
done
### END - Check parameters ###
# Set var
BASENAME_MSEED=$( basename ${FILE_OUTPUT_MSEED} )
DIRNAME_MSEED=$( dirname ${FILE_OUTPUT_MSEED} )
DIR_NODE=$( cd ${DIRNAME_MSEED} && cd ../ && pwd )
DIR_MSEED_LOG=${DIRNAME_MSEED}/log
FILE_OUTPUT_MSEED_HTTPCODE_LOG=${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).httpcode
# create dir
if [ ! -d ${DIR_MSEED_LOG} ]; then
mkdir -p ${DIR_MSEED_LOG}
fi
echo "${INPUT_STRING} - retreving miniseed from \"${DATASELECT_URL}\""
# running process
COUNT=1
COUNT_LIMIT=5
HTTP_CODE=429
RET_CODE=-9
PREV=0
while ( (( ${HTTP_CODE} == 429 )) || (( ${HTTP_CODE} == 503 )) ) && (( ${COUNT} <= ${COUNT_LIMIT} )); do
curl --digest "${DATASELECT_URL}" -o "${FILE_OUTPUT_MSEED}" --write-out "%{http_code}\\n" > ${FILE_OUTPUT_MSEED_HTTPCODE_LOG} -s
RET_CODE=$?
HTTP_CODE=$( cat ${FILE_OUTPUT_MSEED_HTTPCODE_LOG} )
if (( ${HTTP_CODE} == 429 )); then
echo "TOO MANY REQUEST (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}. Tentative: ${COUNT}/${COUNT_LIMIT}"
sleep 5
PREV=1
elif (( ${HTTP_CODE} == 503 )); then
echo "SERVICE UNAVAILABLE (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}. Tentative: ${COUNT}/${COUNT_LIMIT}"
sleep 5
PREV=1
elif (( ${PREV} == 1 )); then
echo " DONE (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}. Tentative: ${COUNT}/${COUNT_LIMIT}"
fi
COUNT=$(( ${COUNT} + 1 ))
done
if (( ${RET_CODE} == 0 )); then
if (( ${HTTP_CODE} == 200 )); then
if [ -f ${FILE_OUTPUT_MSEED} ]; then
echo "OK - file \"${FILE_OUTPUT_MSEED}\" successfully downloaded." > /dev/null
# Use qmerge to cut file properly
if [[ ! -z ${STARTTIME} ]] && [[ ! -z ${ENDTIME} ]]; then
#echo " use qmerge to cut file properly"
qmerge -f ${STARTTIME} -t ${ENDTIME} ${FILE_OUTPUT_MSEED} > ${FILE_OUTPUT_MSEED}.new 2> ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).qmerge.log
mv ${FILE_OUTPUT_MSEED}.new ${FILE_OUTPUT_MSEED}
fi
if [[ "${TYPE}" == "sac" ]]; then
# create SAC dir
DIR_SAC_NODE=${DIR_NODE}/sac
if [ ! -d ${DIR_SAC_NODE} ]; then
mkdir -p ${DIR_SAC_NODE}
fi
# Convert to SAC.
python3 /opt/seed_handler.py --oud ${DIR_SAC_NODE} --fmtout SAC --filein ${FILE_OUTPUT_MSEED}
RET=$?
if [ $RET -ne 0 ]; then
echo " ERROR - converting \"${FILE_OUTPUT_MSEED}\" to SAC format."
else
echo " OK (for ${INPUT_STRING}) - converting \"${FILE_OUTPUT_MSEED}\" to SAC format."
fi
fi
else
echo " ERROR - skip conversion. File ${FILE_OUTPUT_DLESS} not found."
fi
elif (( ${HTTP_CODE} == 204 )); then
echo " NODATA (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
if [ -f ${FILE_OUTPUT_MSEED} ]; then
mv ${FILE_OUTPUT_MSEED} ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
echo "NODATA" >> ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
echo "REQUEST: ${DATASELECT_URL}" >> ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
echo "RET_CODE: ${RET_CODE}" >> ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
echo "HTTP_CODE: ${HTTP_CODE}" >> ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
fi
elif (( ${HTTP_CODE} == 403 )); then
echo " FORBIDDEN (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
if [ -f ${FILE_OUTPUT_MSEED} ]; then
mv ${FILE_OUTPUT_MSEED} ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
fi
elif (( ${HTTP_CODE} == 429 )); then
#echo " TOO MANY REQUEST (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
echo ${DATASELECT_URL} > ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).tooManyRequest
elif (( ${HTTP_CODE} == 503 )); then
#echo " SERVICE UNAVAILABLE (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
echo ${DATASELECT_URL} > ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).serviceUnavailable
else
echo " UNKNOWN (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
fi
else
echo " ERROR (for ${INPUT_STRING}) - retrieving \"${DATASELECT_URL}\". RET_CODE=${RET_CODE}, HTTP_CODE=${HTTP_CODE}"
if [ -f ${FILE_OUTPUT_MSEED} ]; then
mv ${FILE_OUTPUT_MSEED} ${DIR_MSEED_LOG}/$( basename ${FILE_OUTPUT_MSEED} ).log
fi
fi
#
if [ -f ${FILE_OUTPUT_MSEED_HTTPCODE_LOG} ]; then
rm -f ${FILE_OUTPUT_MSEED_HTTPCODE_LOG}
fi