# 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 --hfcc --start-data="h2o+_X h2o+_Y h2o+_Z" --inp=HFS

# pSSCC
respect --sscc --start-data=h2o+_X --inp=JJ13_X
respect --sscc --start-data=h2o+_Y --inp=JJ13_Y
respect --sscc --start-data=h2o+_Z --inp=JJ13_Z
respect --sscc --start-data=h2o+_X --inp=JJ23_X
respect --sscc --start-data=h2o+_Y --inp=JJ23_Y
respect --sscc --start-data=h2o+_Z --inp=JJ23_Z

util/pNMR-sscc.py --quiet --t=293 --m=2 --coupling="1 3,2 3" --hfcc=h2o+-HFS.out_hfcc --sscc="h2o+_X-JJ13_X.out_sscc h2o+_Y-JJ13_Y.out_sscc h2o+_Z-JJ13_Z.out_sscc h2o+_X-JJ23_X.out_sscc h2o+_Y-JJ23_Y.out_sscc h2o+_Z-JJ23_Z.out_sscc" --label=JJ_H-O

generic_test h2o+-HFS.out_pnmr-sscc_JJ_H-O

# clean
rm -rf *.50*
