#!/bin/sh

"$AMSBIN/ams" << eor

Task SinglePoint
System
    Atoms
        Tc 11.75071499547028 6.407900997529855 5.454960997897198 adf.R=1.967
        Tc 9.859737996199225 9.227135996443083 4.191235998384342 adf.R=1.967
        Tc 13.81217999467562 5.310298997952962 2.795978998922194 adf.R=1.967
        Tc 10.11900199609928 3.446142998671566 4.804797998147825 adf.R=1.967
        Tc 13.6988519947193 3.897665998497511 7.264892997199498 adf.R=1.967
        Tc 13.9152179946359 9.00725899652784 6.78883599738301 adf.R=1.967
        O 13.6420469947412 7.825387996983435 10.64965299589472 adf.R=1.72
        O 11.25048399566312 7.673856997041848 9.297477996415966 adf.R=1.72
        Tc 12.07871699534384 11.9885799953786 4.286675998347553 adf.R=1.967
        Tc 8.541336996707448 11.64838399550973 1.92812099925674 adf.R=1.967
        Tc 11.55529199554562 8.663581996660325 0.9649499996280269 adf.R=1.967
        O 12.3553739952372 4.903300998109853 11.58499299553416 adf.R=1.72
        Tc 7.644525997053155 7.224437997215093 2.086721999195602 adf.R=1.967
        O 8.95327799654865 6.558393997471843 10.31033799602552 adf.R=1.72
        O 9.025173996520936 6.705002997415326 13.08077299495756 adf.R=1.72
        O 10.82043799582888 7.818700996986013 4.789801998153607 adf.R=1.72
        O 8.942573996552778 9.369190996388323 5.74547499778521 adf.R=1.72
        O 12.97297499499912 6.456952997510947 3.863358998510736 adf.R=1.72
        O 14.81402699428942 6.157344997626439 1.757573999322482 adf.R=1.72
        O 14.75151999431352 4.266542998355314 3.68416999857981 adf.R=1.72
        O 12.71188799509976 4.42500199829423 1.931678999255368 adf.R=1.72
        O 10.58665999591901 5.133408998021152 4.449216998284896 adf.R=1.72
        O 10.11591499610048 3.182927998773032 6.460958997509402 adf.R=1.72
        O 11.20697699567988 2.389796999078771 4.070505998430884 adf.R=1.72
        O 8.58226299669167 3.181651998773523 4.186151998386304 adf.R=1.72
        O 12.7805499950733 4.792852998152431 6.031413997674984 adf.R=1.72
        O 15.22099599413254 3.505722998648598 6.68081499742465 adf.R=1.72
        O 12.87722299503604 2.473278999046589 7.650466997050864 adf.R=1.72
        O 13.85561799465888 4.815439998143723 8.673493996656502 adf.R=1.72
        O 13.0698039949618 7.572688997080847 6.406312997530467 adf.R=1.72
        O 15.03005299420616 8.714255996640793 7.930051996943089 adf.R=1.72
        O 14.54987899439126 9.553629996317227 5.329909997945403 adf.R=1.72
        O 12.93741699501282 10.26256799604394 7.283193997192443 adf.R=1.72
        O 11.29431599564622 10.48633299595768 4.796155998151157 adf.R=1.72
        O 13.4106649948304 11.66704899550254 3.3435839987111 adf.R=1.72
        O 12.5800259951506 12.86341199504136 5.6438079978244 adf.R=1.72
        O 10.99434399576186 12.93169999501504 3.393146998691994 adf.R=1.72
        O 8.892330996572147 10.82003299582904 3.45414799866848 adf.R=1.72
        O 9.919073996176351 11.68599499549524 0.9607559996296436 adf.R=1.72
        O 7.334056997172836 10.81534499583086 1.123423999566938 adf.R=1.72
        O 8.01835699690905 13.2218059949032 2.245576999134366 adf.R=1.72
        O 10.73615199586138 9.160251996468865 2.398077999075579 adf.R=1.72
        O 12.9937869949911 9.50324099633665 0.9061099996507087 adf.R=1.72
        O 10.79709999583789 4.898612998111661 9.254268996432625 adf.R=1.72
        O 11.8648159954263 7.034959997288134 0.9887149996188659 adf.R=1.72
        O 9.770355996233681 4.032650998445477 11.74765999547146 adf.R=1.72
        O 10.6642559958891 9.013642996525382 -0.3723569998564623 adf.R=1.72
        O 8.438946996746919 8.044788996898861 3.410938998685138 adf.R=1.72
        O 6.047383997668828 7.893257996957272 2.058762999206378 adf.R=1.72
        O 8.319316996793033 7.600873997069981 0.6849839997359495 adf.R=1.72
        O 7.462424997123352 5.635628997827554 2.368848999086846 adf.R=1.72
        Tc 7.93871499693975 9.310744996410856 7.253455997203907 adf.R=1.967
        Tc 9.829692996210808 6.491508997497625 8.517180996716762 adf.R=1.967
        Tc 5.877250997734412 10.40834599598774 9.912437996178911 adf.R=1.967
        Tc 9.570427996310748 12.27250199526914 7.903617996953278 adf.R=1.967
        Tc 5.990577997690727 11.8209799954432 5.443522997901607 adf.R=1.967
        Tc 5.774211997774132 6.711386997412865 5.919579997718095 adf.R=1.967
        O 11.370112995617 8.117771996870728 12.02343299536516 adf.R=1.72
        O 12.22700499528668 10.08301699611316 10.33956799601426 adf.R=1.72
        Tc 7.610713997066189 3.730064998562118 8.421740996753552 adf.R=1.967
        Tc 11.14809299570258 4.070261998430976 10.78029599584436 adf.R=1.967
        Tc 8.134138996864417 7.055063997280384 11.74346599547308 adf.R=1.967
        O 6.695643997418934 6.215404997604058 11.8023059954504 adf.R=1.72
        Tc 12.04490399535688 8.494206996725616 10.6216949959055 adf.R=1.967
        O 7.824613996983734 8.683685996652573 11.71970199548224 adf.R=1.72
        O 11.67107399550098 2.496839999037508 10.46283899596674 adf.R=1.72
        O 8.868991996581142 7.899944996954694 7.918613996947498 adf.R=1.72
        O 10.74685599585726 6.349453997552384 6.962940997315894 adf.R=1.72
        O 6.716454997410912 9.261691996429764 8.845057996590368 adf.R=1.72
        O 4.875402998120608 9.561300996314268 10.95084299577862 adf.R=1.72
        O 4.937909998096512 11.4521029955854 9.024246996521295 adf.R=1.72
        O 6.977541997310267 11.29364299564648 10.77673699584574 adf.R=1.72
        O 9.102769996491025 10.58523599591956 8.259199996816207 adf.R=1.72
        O 9.57351499630956 12.53571799516768 6.247457997591702 adf.R=1.72
        O 8.482453996730147 13.32884899486194 8.63791099667022 adf.R=1.72
        O 11.10716799571836 12.53699399516718 8.5222639967148 adf.R=1.72
        O 6.908880997336733 10.92579199578828 6.677001997426121 adf.R=1.72
        O 4.468434998277488 12.21292199529212 6.027601997676453 adf.R=1.72
        O 6.812207997374001 13.24536699489412 5.05794999805024 adf.R=1.72
        O 5.833811997751158 10.90320499579698 4.034921998444601 adf.R=1.72
        O 6.619625997448237 8.145955996859863 6.302103997570638 adf.R=1.72
        O 4.659377998203883 7.004389997299918 4.778364998158016 adf.R=1.72
        O 5.139551998018783 6.165015997623483 7.378506997155702 adf.R=1.72
        O 6.752013997397206 5.456077997896768 5.425222997908661 adf.R=1.72
        O 8.395114996763816 5.232312997983026 7.912259996949948 adf.R=1.72
        O 6.278764997579634 4.051596998438173 9.364831996390006 adf.R=1.72
        O 7.109404997259436 2.855232998899352 7.064608997276704 adf.R=1.72
        O 8.695086996648179 2.786944998925676 9.315268996409111 adf.R=1.72
    End
    Charge -4
    BondOrders
         1 16 1.0
         1 18 1.0
         1 22 1.0
         1 26 1.0
         1 30 1.0
         1 68 1.0
         2 16 1.0
         2 17 1.0
         2 34 1.0
         2 38 1.0
         2 42 1.0
         2 48 1.0
         3 18 1.0
         3 19 1.0
         3 20 1.0
         3 21 1.0
         4 22 1.0
         4 23 1.0
         4 24 1.0
         4 25 1.0
         5 26 1.0
         5 27 1.0
         5 28 1.0
         5 29 1.0
         6 30 1.0
         6 31 1.0
         6 32 1.0
         6 33 1.0
         7 64 1.0
         8 53 1.0
         8 64 1.0
         9 34 1.0
         9 35 1.0
         9 36 1.0
         9 37 1.0
         10 38 1.0
         10 39 1.0
         10 40 1.0
         10 41 1.0
         11 42 1.0
         11 43 1.0
         11 45 1.0
         11 47 1.0
         12 61 1.0
         13 48 1.0
         13 49 1.0
         13 50 1.0
         13 51 1.0
         14 53 1.0
         14 62 1.0
         15 62 1.0
         17 52 1.0
         44 53 1.0
         44 61 1.0
         46 61 1.0
         52 67 1.0
         52 69 1.0
         52 73 1.0
         52 77 1.0
         52 81 1.0
         53 67 1.0
         53 68 1.0
         53 85 1.0
         54 69 1.0
         54 70 1.0
         54 71 1.0
         54 72 1.0
         55 73 1.0
         55 74 1.0
         55 75 1.0
         55 76 1.0
         56 77 1.0
         56 78 1.0
         56 79 1.0
         56 80 1.0
         57 81 1.0
         57 82 1.0
         57 83 1.0
         57 84 1.0
         58 64 1.0
         59 64 1.0
         60 85 1.0
         60 86 1.0
         60 87 1.0
         60 88 1.0
         61 66 1.0
         62 63 1.0
         62 65 1.0
    End
End

Engine ADF
    Basis
        Type TZ2P
        Core None
    End
    Save TAPE10
    XC
        GGA OPBE
    End
    Relativity
        Level Spin-Orbit
    End
    Symmetry NOSYM
    Solvation
        Surf Delley
        Solv name=Water cav0=0.0 cav1=0.0067639
        Charged method=CONJ
        C-Mat POT
        SCF VAR ALL 
        CSMRSP
    End
    NumericalQuality VeryGood
EndEngine
eor

# =======
# NMR NMR
# =======

cp "${AMS_RESULTSDIR-${AMS_JOBNAME-ams}.results}/adf.rkf" TAPE21
cp "${AMS_RESULTSDIR-${AMS_JOBNAME-ams}.results}/TAPE10" TAPE10

"$AMSBIN/nmr" << eor
NMR
    out    iso    tens
    adfgui
    Atoms 1 2 11 13 3 4 5 6 52 53 
    Atoms 54 55 56 57 60 61 62 64 9 10 
    u1k best
    calc all
End

eor

mv TAPE21 "${AMS_RESULTSDIR-${AMS_JOBNAME-ams}.results}/adf.rkf"
rm -f TAPE10 TAPE15


