Skip to content

README.md

Latest commit

 

History

History
78 lines (64 loc) · 2.44 KB

README.md

File metadata and controls

78 lines (64 loc) · 2.44 KB
Raw

PeakFQ for Linux

This program is an unofficial personal version of PeakFQ for Linux, which was modified from the official USGS PeakFQ for Windows. Their source code was missing some files and not compilable on Linux. Other than fixing compilation issues, the previous default computation of generalized skew is reenabled for my research. I used the following methods to compile it on Linux:

How to build

First, install the Intel Fortran Compiler on Linux.

Use make:

git clone https://github.com/HuidaeCho/peakfq-linux.git
cd peakfq-linux
make

# copy peakfq to your PATH

Or you can directly download the source code and other files, and compile it:

wget -O compile_peakfq.sh https://clawiki.isnew.info/_export/code/howtos/how_to_compile_peakfq_on_linux?codeblock=0
chmod a+x compile_peakfq.sh
./compile_peakfq.sh

# copy peakfq to your PATH

How to test

mkdir test
cd test

wget -O NM.INP "https://nwis.waterdata.usgs.gov/nwis/peak?state_cd=nm&group_key=NONE&sitefile_output_format=html_table&column_name=agency_cd&column_name=site_no&column_name=station_nm&set_logscale_y=1&date_format=YYYY-MM-DD&rdb_compression=file&format=hn2&hn2_compression=file&list_of_search_criteria=state_cd"

for i in `grep ^Z NM.INP  | sed 's/ .*//;s/^Z//'`; do
	cat << EOT > sta$i.psf
I ASCI STA$i.INP
O FILE STA$i.OUT
O DEBUG YES
O ADDITIONAL BOTH
STATION $i
  ANALYZE EMA
EOT
	grep "^.$i " NM.INP > STA$i.INP
	../peakfq sta$i.psf
done

# extract EMA estimates with generalized skew
(echo "station,probability,mean"
for i in *.OUT; do
	awk '
	/^ *Station -/{
		station = $3
	}
	/^PROBABILITY/{
		start = 1
		next
	}
	/^$/{
		if(start == 1)
			start = 2
		else if(start == 2)
			start = 0
		next
	}
	start == 2{
		prob = substr($0, 1, 9) + 0
		if(prob == 0.002 || prob == 0.01 || prob == 0.02 || prob == 0.04 || prob == 0.1 || prob == 0.2 || prob == 0.5)
			printf "%s,%f,%f\n", station, prob, substr($0, 10, 10)
	}' $i
done) > nm_ema.csv

Disclaimer

This program is my personal work derived from the USGS PeakFQ and I am not responsible or liable for any loss or damage resulting from using it.