PARSER = nm64file
include ../irena/libirena.make

%.ES: %.dat irenafile
	./irenafile < $< > $@

%.hist: %.ES
	./NM64.awk '{do_HISTS()};END{END_HIST()}' $(HISTFLAGS)  $< > $@

%.MIN: %.ES
	./NM64.awk '{Counts()};END{ENDCounts()}' $< > $@

.PRECIOUS: %.ES


%.HK: %.dat ./nm64file
	./nm64file < $< | grep ^H > $@

%.HK: %.dat.xz ./nm64file
	xzcat $< | ./nm64file | grep ^H > $@

%.NM: %.b64 ./nm64file
	./nm64file --base64 < $< > $@

%.NM: %.dat ./nm64file
	./nm64file < $< > $@

%.NM: %.dat.xz ./nm64file
	xzcat $< | ./nm64file > $@

%.NMnE: %.dat ./nm64file
	./nm64file < $< | grep -v ^E > $@

%.NMnE: %.dat.xz ./nm64file
	xzcat $< | ./nm64file | grep -v ^E > $@

%.P: %.dat ./nm64file NM64.awk ../pressure.awk
	./nm64file < $< | ./NM64.awk '{doPressure()}' > $@

%.P: %.dat.xz ./nm64file NM64.awk ../pressure.awk
	xzcat $< | ./nm64file | ./NM64.awk '{doPressure()}' > $@

%.nm_raw: %.dat # ./nm64file NM64.awk ../pressure.awk
	./nm64file < $< | ./NM64.awk '{nm_raw()}' > $@

%.nm_raw: %.dat.xz # ./nm64file NM64.awk ../pressure.awk
	xzcat $< | ./nm64file | ./NM64.awk '{nm_raw()}' > $@

%.pclk: %.dat
	./nm64file < $< | ./NM64.awk 'isBATE(){};/rtc.f=.*pclk=[56][0-9]{7}$$/{print $$0, TEMP}' | tr = \  >$@

nm64/nm64_kiel_2019-05-2x.thr_scan: nm64/nm64_kiel_2019-05-2[5-9]*.NM
	awk '/^H/{if(/nm64_thr/)thr=$$NF;T=$$2};/^C/&&thr>0{print T, thr+0, $$0}' \
            nm64/nm64_kiel_2019-05-2[5-9]*.NM \
      | sort -u -k1,1n \
      > $@

IGY_DIR=nm64/jfj/igy
$(IGY_DIR)/%.nm_avg: $(IGY_DIR)/%.nm_raw
	./nm-hour-rate.py --igy -s 600 -r 300 $< > $@

$(IGY_DIR)/%.1h_slope: $(IGY_DIR)/%.nm_avg
	awk '$$5>3000 && $$7>3000 && $$9>3000{print $$1, $$2, $$5+$$7+$$9}' $< \
      | $(HOME)/bin/derive.py -x1 -y2,3 -s3600 -p2 \
      > $@

YEAR:=$(shell date +%Y)
$(IGY_DIR)/JFJ_IGY_$(YEAR).nm_raw: $(IGY_DIR)/$(YEAR)/*Z.nm_raw
	find $(IGY_DIR)/$(YEAR)/ -maxdepth 1 -name '*Z.nm_raw' | xargs cat | sort -u > $@	

JFJ_IGY_2020_11_thrscan_FILES:=$(patsubst %.dat,%.NM,$(wildcard \
 nm64/jfj/igy/JFJ_IGY_2020-11-1[1-9]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-2[01]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-22T0[0-7]*Z.dat \
))
nm64/jfj/igy/JFJ_IGY_2020-11-thrscan.2dhist:%: $(JFJ_IGY_2020_11_thrscan_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(0)}' maxH=1000 T1=1605128400 $^ >$@

nm64/jfj/igy/JFJ_IGY_2020-11-thrscan.CH:%: $(JFJ_IGY_2020_11_thrscan_FILES)
	./NM64.awk '!OLF(){next};/^[CH]/;count_E()&&Time>T1{do_count_E()}' T1=1605128400 $^ > $@

JFJ_IGY_2020_11_22_thrscan_FILES:=$(patsubst %.dat,%.NM,$(wildcard \
 nm64/jfj/igy/JFJ_IGY_2020-11-22T0[7-9]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-22T1[0-7]*Z.dat \
))
nm64/jfj/igy/JFJ_IGY_2020-11-22-thrscan.2dhist:%: $(JFJ_IGY_2020_11_22_thrscan_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(0)}' maxH=2000 T2=1606065900 $^ >$@

nm64/jfj/igy/JFJ_IGY_2020-11-22-thrscan.CH:%: $(JFJ_IGY_2020_11_22_thrscan_FILES)
	./NM64.awk '!OLF(){next};/^[CH]/;count_E()&&Time>T1{do_count_E()}' $^ > $@

JFJ_IGY_2020_11_thrscan_2_FILES:=$(patsubst %.dat,%.NM,$(wildcard \
 nm64/jfj/igy/JFJ_IGY_2020-11-22T19*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-22T2*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-2[3-6]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-27T00*Z.dat \
))
nm64/jfj/igy/JFJ_IGY_2020-11-thrscan-2.2dhist:%: $(JFJ_IGY_2020_11_thrscan_2_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(0)}' maxH=2000 T1=1606071600 T2=1606438810 $^ >$@

nm64/jfj/igy/JFJ_IGY_2020-11-thrscan-2-r2.2dhist:%: $(JFJ_IGY_2020_11_thrscan_2_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(2)}' maxH=2000 T1=1606071600 T2=1606438810 $^ >$@

nm64/jfj/igy/JFJ_IGY_2020-11-thrscan-2.CH:%: $(JFJ_IGY_2020_11_thrscan_2_FILES)
	./NM64.awk '!OLF(){next};/^[CH]/;count_E()&&Time>=T1&&Time<=T2{do_count_E()}' \
             T1=1606071600 T2=1606438810 $^ > $@

JFJ_IGY_2020_11_thrscan_sec2_FILES:=$(patsubst %.dat,%.NM,$(wildcard \
 nm64/jfj/igy/JFJ_IGY_2020-11-29T[12]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-11-30T*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-0[1-3]T*Z.dat \
))
nm64/jfj/igy/JFJ_IGY_2020-11-thrscan-sec2.2dhist:%: $(JFJ_IGY_2020_11_thrscan_sec2_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(-2)}' maxH=2000 T2=1607035800 $^ >$@

JFJ_IGY_2020_12_thrscan_sec3_FILES:=$(patsubst %.dat,%.NM,$(wildcard \
 nm64/jfj/igy/JFJ_IGY_2020-12-08T0[7-9]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-08T[12]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-09T*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-1[01]*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-12T0*Z.dat \
 nm64/jfj/igy/JFJ_IGY_2020-12-12T1[0-6]*Z.dat \
))
nm64/jfj/igy/JFJ_IGY_2020-12-thrscan-sec3.2dhist:%: $(JFJ_IGY_2020_12_thrscan_sec3_FILES)
	./NM64.awk 'is_THR2D(){do_THR2D(-3)}' maxH=2400 $^ >$@

%.scanspec: %.2dhist
	awk 'NF{R+=$$5;D=$$3;T=-$$2};/^$$/&&T>700{print T,R/D,sqrt(R)/D;R=0}' $< \
      | derive.py -x 1 -y 2,3 -s 800 -p 3 -N 2 >$@

SECTION=1
MONTH:=$(shell date -d -2days +"%Y-%m")
DYNSPEC: nm64/jfj/igy/JFJ_IGY_$(MONTH)_s$(SECTION).dynspec

J=-j8
NM:
	find nm64/jfj/igy -name 'JFJ_IGY_$(MONTH)*.dat' | sed s,.dat,.NM, | sort | xargs make $J

nm64/jfj/igy/JFJ_IGY_$(MONTH)_s$(SECTION).dynspec: 
	find nm64/jfj/igy/ -path nm64/jfj/igy/JFJ_IGY_$(MONTH)\*.NM | sort | xargs cat | ./NM64.awk '{dynspec()}' ds_ch=$(SECTION) > $@