# guess
respect --scf --inp=h2o+_dkh2

# mDKS-RI
ln -sf  h2o+_dkh2.50  h2o+_X.50
ln -sf  h2o+_dkh2.50  h2o+_Y.50
ln -sf  h2o+_dkh2.50  h2o+_Z.50
respect --scf --inp=h2o+_X --restart
respect --scf --inp=h2o+_Y --restart
respect --scf --inp=h2o+_Z --restart

# EPR parameters
respect --gt   --start-data="h2o+_X h2o+_Y h2o+_Z" --inp=GT
respect --hfcc --start-data="h2o+_X h2o+_Y h2o+_Z" --inp=HFS

# test paramagnetic part of pNMR only
util/pNMR-cs.py --quiet --t=293 --m=2 --select="1 2" --label=Hp --gt=h2o+-GT.out_gt --hfcc=h2o+-HFS.out_hfcc --skip-orb
util/pNMR-cs.py --quiet --t=293 --m=2 --select=3 --label=Op --gt=h2o+-GT.out_gt --hfcc=h2o+-HFS.out_hfcc --skip-orb

generic_test h2o+-GT.out_pnmr-cs_Hp
generic_test h2o+-GT.out_pnmr-cs_Op

# NMR
respect --cs --start-data=h2o+_X --inp=GIAO_RI-J_X
respect --cs --start-data=h2o+_Y --inp=GIAO_RI-J_Y
respect --cs --start-data=h2o+_Z --inp=GIAO_RI-J_Z

respect --cs --start-data=h2o+_X --inp=CGO-X
respect --cs --start-data=h2o+_Y --inp=CGO-Y
respect --cs --start-data=h2o+_Z --inp=CGO-Z

# test pNMR script
util/pNMR-cs.py --quiet --t=293 --m=2 --select="1 2" --label=H --gt=h2o+-GT.out_gt --hfcc=h2o+-HFS.out_hfcc --cs="h2o+_X-GIAO_RI-J_X.out_cs h2o+_Y-GIAO_RI-J_Y.out_cs h2o+_Z-GIAO_RI-J_Z.out_cs"
util/pNMR-cs.py --quiet --t=293 --m=2 --select=3 --label=O --gt=h2o+-GT.out_gt --hfcc=h2o+-HFS.out_hfcc --cs="h2o+_X-GIAO_RI-J_X.out_cs h2o+_Y-GIAO_RI-J_Y.out_cs h2o+_Z-GIAO_RI-J_Z.out_cs"

generic_test h2o+-GT.out_pnmr-cs_H
generic_test h2o+-GT.out_pnmr-cs_O

# NMR --> no RI
skip

respect --cs --start-data=h2o+_X --inp=GIAO_noRI_X
respect --cs --start-data=h2o+_Y --inp=GIAO_noRI_Y
respect --cs --start-data=h2o+_Z --inp=GIAO_noRI_Z

# clean
rm -rf *.50*
