#!/bin/csh

# -------------------------------------
# Calculate PGVPGA Values for Grid Data
# -------------------------------------

# Set Parameters
# --------------
  set NRUN       = 11C
  set DIR1       = /data4/s1906/runs/run$NRUN/raw
  set DIR2       = /data4/s1906/runs/run$NRUN/data
  set GETPGV     = /u0/shawn/src/s1906/src/getpgv
  set GETPGA     = /u0/shawn/src/s1906/src/getpga
  set FILTER     = /u0/shawn/src/misc/filter/filterfile
  set NAME2      = $DIR2/PGAT
  set NAME3      = $DIR2/PGVT
  set PGVPAR     = "202551  750     "               # for run11E 250 m
  set PGAPAR     = "202551  750 0.30"               # for run11E 250 m
  set PGVPAR     = "202551  500     "               # for run11C 250 m
  set PGAPAR     = "202551  500 0.30"               # for run11C 250 m
  set PGVPAR     = "202551  750     "               # for run14C 125 m
  set PGAPAR     = "202551  750 0.30"               # for run14C 125 m
  set FLTPAR     = "202551  500 0.30 0.00 0.40 LP"
  set FRQ        = "0.40"

# Preprocess Files
# ----------------
  if (0) then
  getdtraces < $DIR1/T.0.TVx > $DIR2/T.X.demul
  getdtraces < $DIR1/T.0.TVy > $DIR2/T.Y.demul
  endif

# Filter Data
# -----------
  if (0) then
  $FILTER $FLTPAR < $DIR2/T.X.demul > $DIR2/T.X.demul_$FRQ
  $FILTER $FLTPAR < $DIR2/T.Y.demul > $DIR2/T.Y.demul_$FRQ
  endif

# Compute PGV and PGA for Grids 
# -----------------------------
  if (1) then
  $GETPGV $PGVPAR $DIR2/T.X.demul      $DIR2/T.Y.demul      > $DIR2/PGVT
  $GETPGA $PGAPAR $DIR2/T.X.demul      $DIR2/T.Y.demul      > $DIR2/PGAT
 #$GETPGV $PGVPAR $DIR2/T.X.demul_$FRQ $DIR2/T.Y.demul_$FRQ > $DIR2/PGVT_$FRQ
 #$GETPGA $PGAPAR $DIR2/T.X.demul_$FRQ $DIR2/T.Y.demul_$FRQ > $DIR2/PGAT_$FRQ
  endif

