$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$   ----------    Mooring at a Quay   --------
$
$    This examples shows how to moor a barge bow to a quay. 
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      set basics
$
&device -g_default device -oecho yes
&dimen -dimen meters k-nts
$$ &device -g_default screen
$
$ ***** variables
$
&set use_obs = .false. $ set to true to account for interaction between barge and quay
$
$*********************************************     read model
$
inmodel
$
$*********************************************     set condition
$
$ set barge draft to 10m so fender point align
&instate -loc barge 105   0 -10 0 0 0
&weight -compute barge 5 20 20 20
&status
&status force
$
$ **** verify / check location of fender pts
$
chkpts
$
$*********************************************     interest points
$
&select :interest -select *b@ *q@
&describe interest -associate :interest
$
$*********************************************     define connectors
$
medit
    ~fend gspr  compression x 1000 2000 y 0.001 2000 z 0.001 2000
    CONNECTOR f1 ~fend *b001 *q001 -euler -x
    CONNECTOR f2 ~fend *b002 *q002 -euler -x
    CONNECTOR f3 ~fend *b003 *q003 -euler -x
    CONNECTOR f4 ~fend *b004 *q004 -euler -x
    CONNECTOR f5 ~fend *b005 *q005 -euler -x
$
$ **** barge to quay
    ~WIRE1 h_cat   25  -len 2000 -buoy 0 -wtpl .0125
    connector mq1 ~wire1 *ma *qm1
    connector mq2 ~wire1 *mc *qm2
$
$ **** barge to anchor
    ~WIRE  b_cat  100  -dep 40 -len 2000 -buoy 0 -wtpl .0125
    connector ma3 ~wire *mb -anc  -90 10
    connector ma4 ~wire *mb -anc  -45 10
    connector ma5 ~wire *md -anc   45 10
    connector ma6 ~wire *md -anc   90 10
end
$
$ **** verify / check location of fender pts
$
chkpts
$
$*********************************************     set tensions 
$
&connector mq/ -l_tension 80
&connector ma@ -a_tension 80
&status f_connect
&status g_connect
&picture iso -render gl -conn no
&picture top -render gl -conn no
&dcptime Time to Connect
$
$*********************************************     check tensions 
$ demonstration of using macros for repetitive reporting
$ macro 'repos' defined in the .dat file
$
repos  1 $ move to 1m away from quay
repos  1 $ move to 2m away from quay
repos -2 $ move back to initial position
repos -1 $ move to 1m away from quay
repos -2 $ move to 3m away from quay
repos  3 $ move back to initial position
$
$*********************************************      Frequency Domain
$
&parameter -m_distance 10
hydrodynamics
    g_pressure barge -heading 90
end
$
$*********************************************     define environment
$
&ENV TEST -sea issc -90 5  7.5 -TIME 100 1 
&STATUS F_CONNECT
&STATUS Force
&equi -omega .4
&dcptime Time For Equilibrium
$
$*********************************************     Define Report Points
$
&describe interest -associate   *M@
$
$*********************************************     Time Domain
$
TDOM   -newmark
&DCPTIME Time For Time Domain
$
$*********************************************     Time post-processing
$
prcpost
&subt fender forces
    conforce
       report
       vlist
       plot 1 8 16 24 32 40 -no
    end
&subt relative motion *b001 *q001
    rel_motion *b001 *q001
       report
       vlist
       extreme 1 2 3 4 5 
       store 1 2 3 4 5
    end
&subt relative motion *b002 *q002
    rel_motion *b002 *q002
       report
       vlist
       extreme 1 2 3 4 5 
       store 1 2 3 4 5
    end
&subt relative motion *b003 *q003
    rel_motion *b003 *q003
       report
       vlist
       extreme 1 2 3 4 5 
       store 1 2 3 4 5
    end
&subt relative motion *b004 *q004
    rel_motion *b004 *q004
       report
       vlist
       extreme 1 2 3 4 5 
       store 1 2 3 4 5
    end
&subt relative motion *b005 *q005
    rel_motion *b005 *q005
       report
       vlist
       extreme 1 2 3 4 5 
       store 1 2 3 4 5
    end
end
$
$*********************************************     all done
$
&finish