!---------------------------------------------------------------------- 
!---------------------------------------------------------------------- 
!  Muensteranian Torturials on Nonlinear Science
!  edited by Uwe Thiele, Oliver Kamps, Svetlana Gurevich
!  (Center for Nonlinear Science, Universitaet Muenster)
!----------------------------------------------------------------------
!  Section: Continuation with auto07p
!  Tutorial: SLIDROP: sliding drops on an inclined homogeneous substrate
!  by Uwe Thiele (www.uwethiele.de), supported by Christian Schelte,
!  Frank Ehebrecht, Thomas Seidel, Simon Hartmann
!  Version 3, Jan 2022
!  for complete set of files and up to date version see
!  http://dx.doi.org/10.5281/zenodo.4546381
!----------------------------------------------------------------------
!  Files: README_slidrop - Further details and screen output for runs 
!  with auto07p of slidrop.f90 producing sliding drops modelled by a thin
!  film equation 
!----------------------------------------------------------------------
!----------------------------------------------------------------------

General
-------
Here the shell command line commands and screen output are given.
See accompanying SLIDROP.pdf for alternative Python interface commands and 
detailed description of the tutorial

Commands and screen ouput
-------------------------

run 1
-----
# For convenience we first start with horizontal substrate.
# Parameters: hmean(H0)=5.0  (change in stpnt in slidrop.f90)
# Continuation in: parameters PAR(5) (period, per), PAR(6) (integr. constant, c0), PAR(42) (velocity)

@cl
@@R slidrop 1
@sv slidrop1
OUTPUT ON SCREEN:
  BR    PT  TY  LAB    PAR(5)      L2-NORM U(1)    MAX U(1)      MAX U(2)      MAX U(3)      PAR(6)        PAR(42)    
   1     1  EP    1   9.14243E+01   7.07107E-04   1.00000E-03   6.28319E-03   3.94784E-02   0.00000E+00   0.00000E+00
   1   225  UZ    2   5.00000E+01   1.83935E+00   2.19916E+00   3.61628E-01   1.25657E-01  -9.60618E-16   6.62713E-16
   1   389  UZ    3   5.00000E+01   2.73396E+00   3.42059E+00   5.77451E-01   2.10490E-01   8.23407E-15  -4.52913E-14
   1   500        4   8.73508E+01   3.92818E+00   5.90521E+00   6.10225E-01   2.20107E-01   5.14375E-13  -5.38978E-13
   1   538  UZ    5   1.00000E+02   4.20848E+00   6.61731E+00   6.17287E-01   2.22075E-01   1.01270E-12  -1.02066E-12
   1   736  UZ    6   2.00000E+02   5.66756E+00   1.11119E+01   6.48965E-01   2.30252E-01   1.22555E-12  -1.23180E-12
   1  1000  UZ    7   4.00000E+02   7.24113E+00   1.74554E+01   6.74175E-01   2.36042E-01   3.11096E-12  -3.10738E-12
-> NOTE, that the primary bifurcation is  subcritical here!


run 11
------
# Starting from label 7 of run 1 we incline the plane
#  Continuation in: PAR(41) (inclination, alp), PAR(6) (c0), PAR(42) (velocity ,vv)
# We also plot PAR(46), PAR(47) that correspond to advancing and receding
# cont. angle ,respectively.

@@R slidrop 11 slidrop1
@sv slidrop11
  BR    PT  TY  LAB    PAR(41)       L2-NORM       MAX U(1)      MAX U(2)      PAR(6)        PAR(42)       PAR(46)       PAR(47)    
   1   500        8   6.98697E-04   6.50227E+00   1.50682E+01   4.44993E-01  -1.21769E-01   1.16768E-01   7.71151E-01   4.44993E-01
   1  1000        9   8.50474E-04   5.71675E+00   1.28672E+01   4.39318E-01  -1.14577E-01   1.10107E-01   7.48403E-01   4.39318E-01
   1  1331  UZ   10   1.00000E-03   5.23240E+00   1.15514E+01   4.41132E-01  -1.17529E-01   1.13079E-01   7.39510E-01   4.41132E-01
   1  1500       11   1.17221E-03   4.87648E+00   1.05388E+01   4.31678E-01  -1.25059E-01   1.20133E-01   7.34966E-01   4.31678E-01
   1  2000       12   1.81196E-03   3.89983E+00   7.85644E+00   4.00102E-01  -1.42984E-01   1.36798E-01   7.14532E-01   4.00102E-01
   1  2376  UZ   13   2.50000E-03   3.21854E+00   6.17513E+00   3.77051E-01  -1.58247E-01   1.51002E-01   6.96917E-01   3.77051E-01
   1  2500       14   2.92228E-03   2.88886E+00   5.44012E+00   3.65347E-01  -1.66414E-01   1.58595E-01   6.87545E-01   3.65347E-01
   1  2945  UZ   15   5.00000E-03   1.63740E+00   3.26752E+00   3.21585E-01  -1.99901E-01   1.89451E-01   6.50894E-01   3.21585E-01
   1  3000       16   5.34155E-03   1.45328E+00   3.03804E+00   3.16525E-01  -2.05549E-01   1.94252E-01   6.45938E-01   3.16525E-01
   1  3153  UZ   17   7.50000E-03   8.01395E-01   2.26699E+00   3.27928E-01  -3.04406E-01   2.54171E-01   6.25255E-01   3.27928E-01
   1  3184  UZ   18   7.49994E-03   6.19498E-01   1.49420E+00   2.56327E-01  -5.48139E-01   3.02198E-01   4.21704E-01   2.56327E-01
   1  3196  UZ   19   5.00000E-03   5.52512E-01   1.03027E+00   1.30502E-01  -5.87432E-01   2.46015E-01   2.25182E-01   1.30502E-01
   1  3214  UZ   20   5.00000E-03   3.17930E-01   4.68631E-01   3.75615E-02  -9.63515E-01   3.19113E-01   7.19838E-02   3.75615E-02
   1  3231  UZ   21   7.49999E-03   1.83023E-01   2.56250E-01   1.76720E-02  -1.65362E+00   5.18952E-01   3.52736E-02   1.76720E-02
   1  3242  UZ   22   9.99996E-03   1.31356E-01   1.81331E-01   1.19144E-02  -2.29626E+00   7.09756E-01   2.41142E-02   1.19144E-02
-> you can nicely see the transition from sliding drops to surface waves as you increase the inclination


