!----- AGPL --------------------------------------------------------------------
!                                                                               
!  Copyright (C)  Stichting Deltares, 2015.                                     
!                                                                               
!  This file is part of Delft3D (D-Flow Flexible Mesh component).               
!                                                                               
!  Delft3D is free software: you can redistribute it and/or modify              
!  it under the terms of the GNU Affero General Public License as               
!  published by the Free Software Foundation version 3.                         
!                                                                               
!  Delft3D  is distributed in the hope that it will be useful,                  
!  but WITHOUT ANY WARRANTY; without even the implied warranty of               
!  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the                
!  GNU Affero General Public License for more details.                          
!                                                                               
!  You should have received a copy of the GNU Affero General Public License     
!  along with Delft3D.  If not, see <http://www.gnu.org/licenses/>.             
!                                                                               
!  contact: delft3d.support@deltares.nl                                         
!  Stichting Deltares                                                           
!  P.O. Box 177                                                                 
!  2600 MH Delft, The Netherlands                                               
!                                                                               
!  All indications and logos of, and references to, "Delft3D",                  
!  "D-Flow Flexible Mesh" and "Deltares" are registered trademarks of Stichting 
!  Deltares, and remain the property of Stichting Deltares. All rights reserved.
!                                                                               
!-------------------------------------------------------------------------------

! $Id: unstruc_display.F90 42642 2015-10-21 11:34:20Z dam_ar $
! $HeadURL: https://repos.deltares.nl/repos/ds/trunk/additional/unstruc/src/unstruc_display.F90 $

module m_WEARELT
   double precision :: XMIN,YMIN,XMAX,YMAX,X1,Y1,X2,Y2,RCIR,CR,DSIX
END module m_WEARELT

MODULE M_DEVICES
  INTEGER           :: NPX,NPY,NCOLR,NDEV,NOPSYS,IWS,IHS
END MODULE M_DEVICES

module m_textlines
  double precision  :: txsize
  double precision  :: txxpos
  double precision  :: txypos
  character (len=60):: txlin(3)
end module m_textlines

module unstruc_colors

   USE m_WEARELT
   USE M_DEVICES
   use m_textlines
!! Centralizes color definitions for unstruc.
!! Color specifications are based on Interactor.

implicit none

    integer :: klvec=4, klaxs=31, klscl=221, kltex=3, klfra=31, klobs=227, klsam=31, klzm=31, klank=31, klprof=221

    integer :: ncoldn    = 3   !< Design net
    integer :: ncolrn    = 211 !< Previous state net
    integer :: ncolnn    = 89  ! 203 !< Net node dots
    integer :: ncoldg    = 31  !< Design grid
    integer :: ncolrg    = 212 !< Previous state grid
    integer :: ncolln    = 120 !< Land boundary
    integer :: ncolsp    = 204 !< Splines
    integer :: ncoltx    = 210 !< Polygons
    integer :: ncolcrs   = 230 !< Cross sections
    integer :: ncolthd   = 231 !< Thin dams
    integer :: ncolfxw   = 232 !< Fixed weirs
    integer :: ncolmh    = 191 !< Fixed weirs
    integer :: ncolwarn1 = 191 ! warning1
    integer :: ncolwarn2 =  31 ! warning2
    integer :: ncolwarn3 =  22 ! warning3
    integer :: ncolhl    =  31 ! Highlight nodes/links
    integer :: ncolANA   =  63 ! 180! ANALYTIC SOLOUTIONS

    integer :: ncolblack = 3

    ! colors in text screens
    ! 0 : Black       4 : Cyan
    ! 1 : Red         5 : Blue
    ! 2 : Yellow      6 : Magenta
    ! 3 : Green       7 : White

    INTEGER :: STDFOR = 0, STDBCK = 5,    &   !   std
               MNUFOR = 0, MNUBCK = 4,    &   !   choice menu's
               INPFOR = 0, INPBCK = 4,    &   !   input menu's
               ERRFOR = 1, ERRBCK = 7,    &   !   error messages
               LBLFOR = 7, LBLBCK = 5,    &   !   menu names
               LINFOR = 0, LINBCK = 4,    &   !   lines
               TOPFOR = 1, TOPBCK = 7,    &   !   top line
               HLPFOR = 7, HLPBCK = 5,    &   !   help window
               BOTFOR = 7, BOTBCK = 5,    &   !   page indication
               KEYFOR = 1, KEYBCK = 4,    &   !   key indication
               WNDFOR = 0, WNDBCK = 4,    &   !   menu indication, POPUP WINDOW HELP
               SHAFOR = 7, SHABCK = 0         !   menu indication, shadow behind input forms

    integer :: nbluep
    integer :: nblues
    integer :: ngreenp
    integer :: ngreens
    integer :: nredp
    integer :: nreds

    character(len=255) :: coltabfile  = ' '
    character(len=255) :: coltabfile2 = ' '


end module unstruc_colors

module unstruc_display
!! Handles all display settings and screen plotting for Unstruc
!! (Not yet, a lot is still in REST.F90 [AvD])

! $Id: unstruc_display.F90 42642 2015-10-21 11:34:20Z dam_ar $

use unstruc_colors
implicit none
#ifndef HAVE_DISPLAY
#define HAVE_DISPLAY 1
#endif
#if HAVE_DISPLAY==1
    integer :: jaGUI      = 1          !< GUI (1) or not (0)
#else
    integer :: jaGUI      = 0          !< GUI (1) or not (0)
#endif

    integer :: ntek       = 0
    integer :: plottofile = 0
    integer :: jadatetime = 0



    ! Highlight certain net/flow node/link numbers
    integer :: nhlNetNode  = 0          !< Number of netnode  to be highlighted
    integer :: nhlNetLink  = 0          !< Number of netlink  to be highlighted
    integer :: nhlFlowNode = 0          !< Number of flownode to be highlighted
    integer :: nhlFlowLink = 0          !< Number of flowlink to be highlighted
    INTEGER :: NPOS(4)                  !< Size + position of HELP text screen
    integer :: jaHighlight =  0         !< Completely enable/disable highlighting.
 
    integer :: ndrawPol           = 2   !< Polygon, 1=No, 2=Regular, 3=plus numbers ZPL, 4=plus isocolour ZPL
    integer :: ndrawObs           = 2   !< Observationstation : 1='NO, 2=Cross, 3=Cross + name4=Polyfil,5='Polyfil + name,6=Cross+waterlevel,7=Cross+velocity magnitudes
    integer :: ndrawCrossSections = 5   !< how draw cross sections
    integer :: ndrawThinDams      = 2   !< show thin dams  0=no, 1=polylines, 2=net links
    integer :: ndrawFixedWeirs     = 1   !< show fixed weirs 0=no, 1=polylines, 2=flow links
    integer :: ndrawManholes      = 2   !< how draw manholes
  
    integer :: numzoomshift       = 250 !< nr of steps in zoomshift
    double precision :: wetplot   = 0.001 !< only show wet waterlevel points if (hs>wetplot)
    double precision :: yfac      = 0.0   !< cheap perspective
    integer :: jafullbottomline   = 0     !<larger bottomline with more complete description in screen 

contains

subroutine load_displaysettings(filename)
    use properties
    use unstruc_messages
    use unstruc_version_module
    use m_missing
    use M_RAAITEK
    use m_wearelt
    use M_isoscaleunit
    use m_transport, only: iconst_cur
    USE M_FLOW, only: kplot, nplot
!   use unstruc_opengl    ! circular dependency
 

    character(len=*), intent(in) :: filename
    
    type(tree_data), pointer :: dis_ptr
    character(len=2) :: nrstring
    integer :: istat, numdraw, i
    logical :: success, jawel
    integer :: JQN, jaeps, jaland
    
    integer          :: ndraw
    COMMON /DRAWTHIS/   NDRAW(40)
   
    integer          :: ncols, nv,nis,nie,jaAUTO
    double precision :: VMAX,VMIN,DV,VAL
    COMMON /DEPMAX/     VMAX,VMIN,DV,VAL(256),NCOLS(256),NV,NIS,NIE,JAAUTO
  
    integer          :: ncols2, nv2,nis2,nie2,jaAUTO2
    double precision :: VMAX2,VMIN2,DV2,VAL2
    COMMON /DEPMAX2/    VMAX2,VMIN2,DV2,VAL2(256),NCOLS2(256),NV2,NIS2,NIE2,JAAUTO2
 
    double precision :: DX, DY, X0, Y0, RMISS, DXSHOW, XD  
    COMMON /ARCINFO/    DX, DY, X0, Y0, RMISS, DXSHOW, XD

    integer          :: nvec 
    double precision :: vfac, vfacforce
    COMMON /VFAC/       VFAC,VFACFORCE,NVEC
   
    integer          :: jaxis 
    double precision :: xleft, ybot
    COMMON /SCREENAREA/ XLEFT, YBOT,JAXIS
   
    integer          :: ndec 
    double precision :: xsc, ysc, scalesize
    COMMON /SCALEPOS/   XSC ,YSC, SCALESIZE,NDEC
   
    integer          :: NHCDEV, NUMHCOPTS, IHCOPTS
    COMMON /HARDCOPY/   NHCDEV,NUMHCOPTS,IHCOPTS(2,20)
   
    double precision :: tsize
    COMMON /TEXTSIZE/   TSIZE
     
    integer :: jaopengl_loc


    ! Put .dis file into a property tree
    call tree_create(trim(filename), dis_ptr)
    call prop_file('ini',trim(filename), dis_ptr,istat)
    if (istat /= 0) then
        call mess(LEVEL_ERROR, 'Display settings file '''//trim(filename)//''' not found. Code: ', istat)
        return
    else
        call mess(LEVEL_DEBUG, 'Opened display settings file : ', trim(filename) )
    endif


    numdraw = 37
    do i=1,numdraw
        write (nrstring, '(I2)') i
        call prop_get_integer(dis_ptr, '*', 'ndraw('//trim(adjustl(nrstring))//')', ndraw(i), success)
    end do
    
!   load active constituent number
    call prop_get_integer(dis_ptr, '*', 'ICONST', iconst_cur, success)

    call prop_Get_integer(dis_ptr, '*', 'ndrawpol           '  , ndrawpol        , success)
    call prop_Get_integer(dis_ptr, '*', 'ndrawobs           '  , ndrawobs        , success)
    call prop_Get_integer(dis_ptr, '*', 'ndrawcrosssections '  , ndrawcrosssections , success)


    call prop_Get_integer(dis_ptr, '*', 'NHCDEV       '  , NHCDEV          , success)
    call prop_Get_integer(dis_ptr, '*', 'JAEPS        '  , JAEPS           , success)
    call prop_Get_integer(dis_ptr, '*', 'JALAND       '  , JALAND          , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'CR           '  , CR              , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'TSIZE        '  , TSIZE           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'dmiss        '  , dmiss           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'XLEFT        '  , XLEFT           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'YBOT         '  , YBOT            , success)
    call prop_Get_integer(dis_ptr, '*', 'JAXIS        '  , JAXIS           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'VFAC         '  , VFAC            , success)
    call prop_Get_integer(dis_ptr, '*', 'nvec         '  , nvec            , success)
    call prop_Get_integer(dis_ptr, '*', 'NTEK         '  , NTEK            , success)
    call prop_Get_integer(dis_ptr, '*', 'PLOTTOFILE   '  , PLOTTOFILE      , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'ZMINrai      '  , ZMINrai         , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'ZMAXrai      '  , ZMAXrai         , success)
    call prop_Get_integer(dis_ptr, '*', 'jtextflow    '  , jtextflow       , success)
    call prop_Get_integer(dis_ptr, '*', 'numzoomshift '  , numzoomshift    , success)
    call prop_Get_integer(dis_ptr, '*', 'jaHighlight  '  , jaHighlight     , success)
    call prop_Get_integer(dis_ptr, '*', 'nhlNetNode   '  , nhlNetNode      , success)
    call prop_Get_integer(dis_ptr, '*', 'nhlNetLink   '  , nhlNetLink      , success)
    call prop_Get_integer(dis_ptr, '*', 'nhlFlowNode  '  , nhlFlowNode     , success)
    call prop_Get_integer(dis_ptr, '*', 'nhlFlowLink  '  , nhlFlowLink     , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'wetplot      '  , wetplot         , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'yfac         '  , yfac            , success)

    call prop_Get_integer(dis_ptr, '*', 'JAAUTO       '  , JAAUTO          , success)
    call prop_Get_integer(dis_ptr, '*', 'NV           '  , NV              , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'VMIN,        '  , VMIN            , success) ! remove muuuchch later
    call prop_Get_DOUBLE (dis_ptr, '*', 'VMIN         '  , VMIN            , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'VMAX         '  , VMAX            , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'DV           '  , DV              , success)

    call prop_Get_integer(dis_ptr, '*', 'JAAUTO2      '  , JAAUTO2         , success)
    call prop_Get_integer(dis_ptr, '*', 'NV2          '  , NV2             , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'VMIN2        '  , VMIN2           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'VMAX2        '  , VMAX2           , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'DV2          '  , DV2             , success)

    
    call prop_Get_string (dis_ptr, '*', 'UNIT(1)      '  , UNIT(1)         , success)
    call prop_Get_string (dis_ptr, '*', 'PARAMTEX(1)  '  , PARAMTEX(1)     , success)
    call prop_Get_string (dis_ptr, '*', 'UNIT(2)      '  , UNIT(2)         , success)
    call prop_Get_string (dis_ptr, '*', 'PARAMTEX(2)  '  , PARAMTEX(2)     , success)

    call prop_Get_DOUBLE (dis_ptr, '*', 'X1           '  , X1              , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'Y1           '  , Y1              , success)
    call prop_Get_DOUBLE (dis_ptr, '*', 'X2           '  , X2              , success)


    ! Color scheme isolines
    call prop_get_string(dis_ptr, 'isocol', 'COLTABFILE', coltabfile)
    inquire(file = trim(coltabfile), exist = jawel)
    if (jawel) then
        call SETCOLTABFILE(coltabfile ,0)
    end if

    call prop_get_string(dis_ptr, 'isocol2', 'COLTABFILE2', coltabfile2)
    inquire(file = trim(coltabfile2), exist = jawel)
    if (jawel) then
        call SETCOLTABFILE(coltabfile2 ,1)
    end if

    jaopengl_loc = -1
    call prop_Get_integer(dis_ptr, '*', 'jaopengl'  , jaopengl_loc    , success)
    if (jaopengl_loc .ne. -1 ) then
       jaopengl_loc = 1 
       call iset_jaopengl(jaopengl_loc)
    end if
    
    call prop_get_integer(dis_ptr, '*', 'kplot', kplot, success)
    call prop_get_integer(dis_ptr, '*', 'nplot', nplot, success)

    RCIR  = CR*(X2 - X1)
    VFAC  = MAX( 0d0, VFAC)
    VFACFORCE  = MAX( 0d0, VFACFORCE)
    XLEFT = MAX( 0d0,(MIN(XLEFT,0.25d0) ) )
    YBOT  = MAX( 0d0,(MIN(YBOT ,0.25d0) ) )
    JAXIS = MIN( 1  ,(MAX(JAXIS,0) ) )
    IF (JAXIS .EQ. 1) THEN
       IF (XLEFT .EQ. 0) XLEFT = .15
       IF (YBOT  .EQ. 0) YBOT  = .10
    ENDIF
    !CALL NEWWORLD()
    call setwY(x1,y1,x2,y2)

    DO I = 1,NUMHCOPTS
       IF (IHCOPTS(1,I) .EQ. 22) IHCOPTS(2,I) = JAEPS
       IF (IHCOPTS(1,I) .EQ. 5)  IHCOPTS(2,I) = JALAND
    ENDDO

    CALL SETTEXTSIZE()
    if (plottofile == 1) then
        ndraw(10) = 1
    end if

end subroutine load_displaysettings


subroutine save_displaysettings(filename)
    use properties
    use unstruc_messages
    use unstruc_version_module
    use m_missing
    use M_RAAITEK
    use m_wearelt
    use M_isoscaleunit
    use m_transport, only: iconst_cur
    use m_flow
!   use unstruc_opengl    ! circular dependency


    character(len=*), intent(in) :: filename

    type(tree_data), pointer     :: dis_ptr
    character(len=2)             :: nrstring
    character(len=20)            :: rundat
    integer                      :: mfil, istat, numdraw, i
    
    
    integer          :: ndraw
    COMMON /DRAWTHIS/   NDRAW(40)
   
    integer          :: ncols, nv,nis,nie,jaAUTO
    double precision :: VMAX,VMIN,DV,VAL
    COMMON /DEPMAX/     VMAX,VMIN,DV,VAL(256),NCOLS(256),NV,NIS,NIE,JAAUTO
  
    integer          :: ncols2, nv2,nis2,nie2,jaAUTO2
    double precision :: VMAX2,VMIN2,DV2,VAL2
    COMMON /DEPMAX2/    VMAX2,VMIN2,DV2,VAL2(256),NCOLS2(256),NV2,NIS2,NIE2,JAAUTO2
 
    double precision :: DX, DY, X0, Y0, RMISS, DXSHOW, XD  
    COMMON /ARCINFO/    DX, DY, X0, Y0, RMISS, DXSHOW, XD

    integer          :: nvec 
    double precision :: vfac, vfacforce
    COMMON /VFAC/       VFAC,VFACFORCE,NVEC
   
    integer          :: jaxis 
    double precision :: xleft, ybot
    COMMON /SCREENAREA/ XLEFT, YBOT,JAXIS
   
    integer          :: ndec 
    double precision :: xsc, ysc, scalesize
    COMMON /SCALEPOS/   XSC ,YSC, SCALESIZE,NDEC
   
    integer          :: NHCDEV, NUMHCOPTS, IHCOPTS
    COMMON /HARDCOPY/   NHCDEV,NUMHCOPTS,IHCOPTS(2,20)
   
    double precision :: tsize
    COMMON /TEXTSIZE/   TSIZE
      
    integer :: JQN, jaeps, jaland
   
    integer, external :: iget_jaopengl

    call newfil(mfil, filename)

    ! Put .dis file into a property tree
    call tree_create('displaysettings', dis_ptr)

    call prop_set(dis_ptr, '*', 'ndraw(1) ' , ndraw(1) , ' ! clear screen yes/no                                                                                              ')
    call prop_set(dis_ptr, '*', 'ndraw(2) ' , ndraw(2) , ' ! display network                                                                                                  ')
    call prop_set(dis_ptr, '*', 'ndraw(3) ' , ndraw(3) , ' ! display landboundary                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(4) ' , ndraw(4) , ' ! void                                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(5) ' , ndraw(5) , ' ! void , this comment is copied from inidat, should be a copy                                                      ')
    call prop_set(dis_ptr, '*', 'ndraw(6) ' , ndraw(6) , ' ! void                                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(7) ' , ndraw(7) , ' ! WHICH netLINk VALue, 1 = NO  2=link nr, 3=node nrs based on links                                                ')
    call prop_set(dis_ptr, '*', 'ndraw(8) ' , ndraw(8) , ' ! WHICH netNODe VALue, 1 = NO, 2=node nr, 3=nr of links, 4=link nrs based on nodes                                 ')
    call prop_set(dis_ptr, '*', 'ndraw(9) ' , ndraw(9) , ' ! NORMAL OR PERSPECTIVE                                                                                            ')
    call prop_set(dis_ptr, '*', 'ndraw(10)' , ndraw(10), ' ! also plot to file yes/no                                                                                         ')
    call prop_set(dis_ptr, '*', 'ndraw(11)' , ndraw(11), ' ! HOW HISPLAY LINVAL, ,Linkhow: 1=no, 2=numbers, 3=isofil smooth, 4 = isofil, 5=dots                               ')
    call prop_set(dis_ptr, '*', 'ndraw(12)' , ndraw(12), ' ! 1=nodes, 2=links, 3=nodes+links, 4 = no isoscale                                                                 ')
    call prop_set(dis_ptr, '*', 'ndraw(13)' , ndraw(13), ' ! 1 !0  ! SHOW UV VECTORS AND FLOW FORCES                                                                          ')
    call prop_set(dis_ptr, '*', 'ndraw(14)' , ndraw(14), ' ! reserved for disval                                                                                              ')
    call prop_set(dis_ptr, '*', 'ndraw(15)' , ndraw(15), ' ! display splines                                                                                                  ')
    call prop_set(dis_ptr, '*', 'ndraw(16)' , ndraw(16), ' ! display PREVIOUS NETWORK                                                                                         ')
    call prop_set(dis_ptr, '*', 'ndraw(17)' , ndraw(17), ' ! void                                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(18)' , ndraw(18), ' ! Sideview, 1 = no, 2=small, 3=larger, 4=largest                                                                                                  ')
    call prop_set(dis_ptr, '*', 'ndraw(19)' , ndraw(19), ' ! HOW DISPLAY NODEVAL, Nodehow: 1=no, 2=numbers, 3=isofil smooth, 4 = isofil, 5=dots                               ')
    call prop_set(dis_ptr, '*', 'ndraw(20)' , ndraw(20), ' ! void                                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(21)' , ndraw(21), ' ! void                                                                                                             ')
    call prop_set(dis_ptr, '*', 'ndraw(22)' , ndraw(22), ' ! SHOW netcell types tri, quads penta, hexa                                                                        ')
    call prop_set(dis_ptr, '*', 'ndraw(26)' , ndraw(26), ' ! 1 = BITMAP                                                                                                       ')
    call prop_set(dis_ptr, '*', 'ndraw(27)' , ndraw(27), ' ! 1 = nothing, 2 = both surf and bot, 3 just bot, 4 just surf                                                      ')
    call prop_set(dis_ptr, '*', 'ndraw(28)' , ndraw(28), ' ! values at Flow Nodes 1=no, 2=s1, 3=bl, 4=ba, 5=v1    , nodewhat                                                  ')
    call prop_set(dis_ptr, '*', 'ndraw(29)' , ndraw(29), ' ! values at Flow Links 1=no, 2=u1, 3=q1, 4=au,         , linkwhat                                                  ')
    call prop_set(dis_ptr, '*', 'ndraw(30)' , ndraw(30), ' ! do not show all flow links                                                                                       ')
    call prop_set(dis_ptr, '*', 'ndraw(31)' , ndraw(31), ' ! do not show show values at cell corners, 2=ucnx, 3=ucny                                                          ')
    call prop_set(dis_ptr, '*', 'ndraw(32)' , ndraw(32), ' ! show samples coloured dot                                                                                        ')
    call prop_set(dis_ptr, '*', 'ndraw(33)' , ndraw(33), ' ! show values at net cells, 2=number, 3=aspect ratio, 4=orientation vectors, 5=aspect ratio and orientation vectors')
    call prop_set(dis_ptr, '*', 'ndraw(34)' , ndraw(34), ' ! Banf, 0=no, 1 =seqtr, 3=(seqtr-s), tekbanf                                                                       ')
    call prop_set(dis_ptr, '*', 'ndraw(35)' , ndraw(35), ' ! yes draw reference profiles (only for 3D)                                                                        ')
    call prop_set(dis_ptr, '*', 'ndraw(36)' , ndraw(36), ' ! values on flow nodes minus plotlin default 0                                                                     ')
    call prop_set(dis_ptr, '*', 'ndraw(37)' , ndraw(37), ' ! 0 = no, 1 = probe, no boundaryvalues, 2=probe+boundaryvalues                                                     ')

    call prop_set(dis_ptr, '*', 'ndrawpol'  , ndrawpol , ' ! Polygon, 1=No, 2=Regular, 3=plus numbers ZPL, 4=plus isocolour ZPL                                               ')
    call prop_set(dis_ptr, '*', 'ndrawobs'  , ndrawobs , ' ! 1=NO, 2=Cross, 3=Cross + name4=Polyfil,5=Polyfil + name,6=Cross+waterlevel,7=Cross+velocity magnitudes           ')
    call prop_set(dis_ptr, '*', 'ndrawcrosssections'  , ndrawcrosssections , ' ! 1=NO, etc                                                                                      ')
    
!   save active constituent number
    call prop_set_integer(dis_ptr, '*', 'ICONST', iconst_cur, ' ! active constituent number')

    DO I = 1,NUMHCOPTs
      IF (IHCOPTS(1,I) .EQ. 22) JAEPS  = IHCOPTS(2,I)
      IF (IHCOPTS(1,I) .EQ. 5)  JALAND = IHCOPTS(2,I)
    ENDDO

    call prop_set(dis_ptr, '*', 'NHCDEV       '  , NHCDEV          )
    call prop_set(dis_ptr, '*', 'JAEPS        '  , JAEPS           )
    call prop_set(dis_ptr, '*', 'JALAND       '  , JALAND          )
    call prop_set(dis_ptr, '*', 'CR           '  , CR              )
    call prop_set(dis_ptr, '*', 'TSIZE        '  , TSIZE           )
    call prop_set(dis_ptr, '*', 'dmiss        '  , dmiss           )
    call prop_set(dis_ptr, '*', 'XLEFT        '  , XLEFT           )
    call prop_set(dis_ptr, '*', 'YBOT         '  , YBOT            )
    call prop_set(dis_ptr, '*', 'JAXIS        '  , JAXIS           )
    call prop_set(dis_ptr, '*', 'VFAC         '  , VFAC            )
    call prop_set(dis_ptr, '*', 'nvec         '  , nvec            )
    call prop_set(dis_ptr, '*', 'NTEK         '  , NTEK            )
    call prop_set(dis_ptr, '*', 'PLOTTOFILE   '  , PLOTTOFILE      )
    call prop_set(dis_ptr, '*', 'ZMINrai      '  , ZMINrai         )
    call prop_set(dis_ptr, '*', 'ZMAXrai      '  , ZMAXrai         )
    call prop_set(dis_ptr, '*', 'jtextflow    '  , jtextflow       )
    call prop_set(dis_ptr, '*', 'numzoomshift '  , numzoomshift    )
    call prop_set(dis_ptr, '*', 'jaHighlight  '  , jaHighlight     )
    call prop_set(dis_ptr, '*', 'nhlNetNode   '  , nhlNetNode      )
    call prop_set(dis_ptr, '*', 'nhlNetLink   '  , nhlNetLink      )
    call prop_set(dis_ptr, '*', 'nhlFlowNode  '  , nhlFlowNode     )
    call prop_set(dis_ptr, '*', 'nhlFlowLink  '  , nhlFlowLink     )
    call prop_set(dis_ptr, '*', 'wetplot      '  , wetplot         )
    call prop_set(dis_ptr, '*', 'yfac         '  , yfac            )

    call prop_set(dis_ptr, '*', 'JAAUTO       '  , JAAUTO          )
    call prop_set(dis_ptr, '*', 'NV           '  , NV              )
    call prop_set(dis_ptr, '*', 'VMIN,        '  , VMIN            )
    call prop_set(dis_ptr, '*', 'VMAX         '  , VMAX            )
    call prop_set(dis_ptr, '*', 'DV           '  , DV              )

    call prop_set(dis_ptr, '*', 'JAAUTO2      '  , JAAUTO2         )
    call prop_set(dis_ptr, '*', 'NV2          '  , NV2             )
    call prop_set(dis_ptr, '*', 'VMIN2        '  , VMIN2           )
    call prop_set(dis_ptr, '*', 'VMAX2        '  , VMAX2           )
    call prop_set(dis_ptr, '*', 'DV2          '  , DV2             )

    call prop_set(dis_ptr, '*', 'UNIT(1)      '  , UNIT(1)         )
    call prop_set(dis_ptr, '*', 'PARAMTEX(1)  '  , PARAMTEX(1)     )
    call prop_set(dis_ptr, '*', 'UNIT(2)      '  , UNIT(2)         )
    call prop_set(dis_ptr, '*', 'PARAMTEX(2)  '  , PARAMTEX(2)     )

    call prop_set(dis_ptr, '*', 'X1           '  , X1              )
    call prop_set(dis_ptr, '*', 'Y1           '  , Y1              )
    call prop_set(dis_ptr, '*', 'X2           '  , X2              )

   !  call prop_set_string(dis_ptr, '*', 'isocol', 'COLTABFILE', coltabfile)
    call prop_set_string(dis_ptr, '*', 'COLTABFILE',  coltabfile)
    call prop_set_string(dis_ptr, '*', 'COLTABFILE2', coltabfile2)

    call prop_set(dis_ptr, '*', 'jaopengl', iget_jaopengl())

    call prop_set(dis_ptr, '*', 'kplot', kplot)
    call prop_set(dis_ptr, '*', 'nplot', nplot)
    
    call datum(rundat)
    write(mfil, '(a,a)') '# Generated on ', trim(rundat)
    write(mfil, '(a,a)') '# ', trim(unstruc_version_full)
    call prop_write_inifile(mfil, dis_ptr, istat)


   call doclose(mfil)
end subroutine save_displaysettings

!> Plots all observation points in the current viewport
subroutine plotObservations()
    use m_observations
    USE M_FLOWGEOM
    use m_flow
    integer      :: n, NN, K
    character*40 :: tex
    logical, external :: inview

    if (ndrawobs == 1 ) return

    call IGrCharJustify('L')
    call settextsizefac(1.0d0)

    do n  = 1, numobs+nummovobs
       if (.not. inview(xobs(n), yobs(n))) cycle

       call setcol(klobs)

       if (ndrawobs .ne. 4 .and. ndrawobs .ne. 5 ) then
            if (n > numobs) then
                ! It is a moving obs:
                call plotDiamond(xobs(n), yobs(n))
            else
                call plotCross(xobs(n), yobs(n))
            end if
       endif
       if (ndrawobs == 3) then
           call settextsizefac(1.5d0)
           call igrcharfont(7)
           call gtext(' '//trim(namobs(n)), xobs(n), yobs(n), klobs)
           call igrcharfont(1)
       endif

       K  = KOBS(N)
       IF (K > 0) THEN
          if (ndrawobs == 4 .or. ndrawobs == 5 ) then
             nn = size( nd(K)%x )
             call PFILLER(nd(k)%x, nd(k)%y, nn,klobs,klobs)
             if (ndrawobs == 5) then
                call gtext(' '//trim(namobs(n)), xobs(n), yobs(n), 221)
             endif
          else if (ndrawobs == 6) then
             write (tex,'(f10.4)') s1(k)
             call gtext(tex(1:14), xobs(n), yobs(n), ncolblack)
          else if (ndrawobs == 7) then
                       write (tex,'(f10.4)') sqrt(ucx(k)*ucx(k) + ucy(k)*ucy(k) )
             call gtext(tex(1:14), xobs(n), yobs(n), ncolblack)
          endif
       ENDIF
    end do

end subroutine plotObservations


!> Plots all manholes in the current viewport
subroutine plotManholes()
    use m_manholes
    USE M_FLOWGEOM
    use m_flow
    integer      :: n, NN, K
    character*40 :: tex
    logical, external :: inview
    if (ndrawmanholes == 1 ) return

    call setcol(klobs)
    call IGrCharJustify('L')


    !call add_manhole(.5d0*(x1+x2), .5d0*(x1+x2), "bla", MANHOLE_OPEN_MOMENTUM) ! AvD: TMP

    do n  = 1, nummh
       if (.not. inview(manholes(n)%x, manholes(n)%y)) cycle

       if (ndrawmanholes .ne. 4 .and. ndrawmanholes .ne. 5 ) then
            call movabs(manholes(n)%x, manholes(n)%y)
            call hlcir2(1.2*rcir, ncolmh, ncolblack)
       endif
!       if (ndrawobs == 3) then
!           call settextsizefac(1.5d0)
!           call igrcharfont(7)
!           call gtext(' '//trim(namobs(n)), xobs(n), yobs(n), klobs)
!           call igrcharfont(1)
!       endif

!       K  = KOBS(N)
!       IF (K > 0) THEN
!          if (ndrawobs == 4 .or. ndrawobs == 5 ) then
!             nn = size( nd(K)%x )
!             call PFILLER(nd(k)%x, nd(k)%y, nn,klobs,klobs)
!             if (ndrawobs == 5) then
!                call gtext(' '//trim(namobs(n)), xobs(n), yobs(n), 221)
!             endif
!          else if (ndrawobs == 6) then
!             write (tex,'(f10.4)') s1(k)
!             call gtext(tex(1:14), xobs(n), yobs(n), klobs)
!          else if (ndrawobs == 7) then
!                       write (tex,'(f10.4)') sqrt(ucx(k)*ucx(k) + ucy(k)*ucy(k) )
!             call gtext(tex(1:14), xobs(n), yobs(n), klobs)
!          endif
!       ENDIF
    end do

    !call init_manholes() ! AvD: TMP
end subroutine plotManholes

subroutine plotSplines(m1, m2, ncol)
    USE M_SPLINES
    use m_alloc

    implicit none

    integer, intent(in), optional :: m1
    integer, intent(in), optional :: m2
    integer, intent(in), optional :: ncol

    integer :: m1_, m2_, ncol_
    integer :: m, n2, numpi, numnew

    double precision, allocatable, dimension(:) :: xlist, ylist

!   allocate
    allocate(xlist(1), ylist(1))

    if (present(m1)) then
        m1_ = m1
    else
        m1_ = 1
    endif
    if (present(m2)) then
        m2_ = m2
    else
        m2_ = mcs
    endif
    if (present(ncol)) then
        ncol_ = ncol
    else
        ncol_ = ncolsp
    endif


    n2 = 0
    if (m1_ > 0) then
        do m = m1_,m2_
            CALL NUMP(m,NUMPI)

!           reallocate if necessary
            if ( numpi.gt.ubound(xlist,1) ) then
               numnew = int(1.2d0*dble(numpi))+1
               call realloc(xlist,numnew)
               call realloc(ylist,numnew)
            end if
            xlist(1:numpi) = xsp(m,1:numpi)
            ylist(1:numpi) = ysp(m,1:numpi)

            call plotSpline(xlist, ylist, numpi, ncol_)
        enddo
    endif

!   deallocate
    deallocate(xlist, ylist)

    return
end subroutine plotSplines


subroutine plotSpline(xh, yh, numpi, ncol)
    use m_wearelt
    double precision, dimension(numpi), intent(in) :: xh, yh
    integer, intent(in) :: numpi
    integer, intent(in) :: ncol

    !integer :: imax = 500 ! TODO: uit DIMENS [AvD]
!    double precision :: XH2(1000), YH2(1000)
    double precision, allocatable, dimension(:) :: xh2, yh2
    double precision :: xk, yk, tn


    integer :: ndraw
    COMMON /DRAWTHIS/ NDRAW(40)
    integer :: i, met, k, numk
    MET = NDRAW(15)

    if (met == 0) return

!   allocate
    allocate(xh2(numpi), yh2(numpi))

    NUMK  = 20
    CALL SETCOL(NCOL)

    IF (NUMPI .EQ. 1) THEN
        CALL MOVABS(XH(1),YH(1))
        IF (MET .LE. 2) CALL CIR(1.4*RCIR)
    ELSE IF (NUMPI .GT. 1) THEN
        CALL MOVABS(XH(1),YH(1))
        IF (MET .LE. 2) CALL CIR(1.4*RCIR)
        CALL SPLINE(XH,NUMPI,XH2)
        CALL SPLINE(YH,NUMPI,YH2)

         DO 10 I = 1,NUMPI-1
            DO 20 K = 1,NUMK
               TN = (I - 1) + dble(K) / dble(NUMK)
               CALL SPLINT(XH,XH2,NUMPI,TN,XK)
               CALL SPLINT(YH,YH2,NUMPI,TN,YK)
               CALL LNABS(XK,YK)
         20 CONTINUE
            IF (MET .LE. 2) CALL CIR(RCIR)
      10 CONTINUE
    ENDIF

!   allocate
    deallocate(xh2, yh2)

    RETURN
    ! TODO: M,N numbers (tekadministratie) [AvD]
end subroutine plotSpline


subroutine plotCrossSections()
    use m_crosssections

    integer :: i, met, jaArrow
    character :: tex*40

    met = ndrawCrosssections
    if (met == 1) return

    if (met >= 3) then
        jaArrow = 1
    else
        jaArrow = 0
    end if

    call thicklinetexcol(ncolcrs) 

    do i=1,ncrs   ! ToDo: writing labels only after first time step
        tex = ' '
        if (met == 4) then
            tex = trim(crs(i)%name)
        else if (met == 5) then
            tex = '1234567890 m3/s'
            call write_num_label(10, 3, crs(i)%sumvalcur(IPNT_Q1C)) ! discharge
        else if (met == 6) then
            tex = '1234567890 m2'
            call write_num_label(10, 3, crs(i)%sumvalcur(IPNT_AUC)) ! area
        else if (met == 7) then
            tex = '1234567890 m/s'
            call write_num_label(10, 3, crs(i)%sumvalcur(IPNT_U1A)) ! ave velocity
        else if (met == 8) then
            tex = '1234567890 m'
            call write_num_label(10, 4, crs(i)%sumvalcur(IPNT_S1A)) ! ave. waterlevel
        else if (met == 9) then
            tex = '1234567890 m'
            call write_num_label(10, 4, crs(i)%sumvalcur(IPNT_HUA)) ! ave. waterdepth 
        endif

        call plotCrossSectionPath(crs(i)%path, 2, ncolcrs, jaArrow, tex)
    end do
    
    call resetlinesizesetc()
    
contains
    !> Plot on TEX label, with certain nr of digits. UNLESS too big, then as integer.
    !! Make sure to keep TEX string length in line with numw variable (10?).
    subroutine write_num_label(numw, numd, val)
        implicit none
        integer, intent(in) :: numw !< Available chars (10?)
        integer, intent(in) :: numd !< Num digits preferred (3?)
        double precision, intent(in) :: val !< Value to be printed

        character(len=7) :: fmt
        fmt = '(f10.3)'
        if (val > -1d0 * 10**(numw-numd-2) .and. val < 10**(numw-numd-1)) then
            write(fmt(3:4), '(i2)') numw
            write(fmt(6:6), '(i1)') numd
            write(tex(1:numw), fmt) val
        else
            if (numw < 10) then
                fmt = '(i1)   '
                write(fmt(3:3), '(i1)') numw
            else
                fmt = '(i10)  '
                write(fmt(3:4), '(i2)') numw
            end if
            write(tex(1:numw),trim(fmt)) int(val,selected_int_kind(15)) ! discharge
        end if
    end subroutine write_num_label
end subroutine plotCrossSections


subroutine plotThinDams()
    use m_thindams

    integer :: i

    if (ndrawThinDams == 0 .or. nthd == 0) return

    call thicklinetexcol(ncolthd) 
    
    do i=1,nthd
        call plotCrossSectionPath(thd(i), ndrawThinDams, ncolthd, 0, ' ')
    end do
        
    call resetlinesizesetc()
end subroutine plotThinDams



subroutine plotFixedWeirs()
   use m_fixedweirs 
   use m_flowgeom, only : lnx, lncn, bob
   use m_flow,     only : hu, isimplefixedweirs  
   use m_netw,     only : xk, yk 

   integer             :: i, L, k3, k4, ncol
   double precision    :: xu, yu

   if (ndrawFixedWeirs == 0 .or. nfxw == 0 ) return

   call thicklinetexcol(ncolfxw) 


   if (isimplefixedweirs == 0) then 
       

      do i=1,nfxw
          call plotCrossSectionPath(fxw(i), ndrawFixedWeirs, ncolfxw, 0, ' ')
      end do

      
   else

  
       if (Lnx == 0) then 
          return
       endif

       do i = 1,nfxw
          L = lnfxw(i) 
          
          if (L > 0) then  
        
             call setcol(ncolfxw)

             k3=lncn(1,L) ; k4=lncn(2,L)  
         
             if (ndrawfixedweirs == 3 .or. ndrawfixedweirs == 4) then 
                call isocol(bob(1,L),ncol) 
             endif
             if (ndrawfixedweirs == 5) then 
                if (hu(L) >  0 ) then 
                   cycle
                endif 
             endif  
             call movabs( xk(k3), yk(k3) )
             call lnabs(  xk(k4), yk(k4) )
          endif      
     
       enddo


       if (ndrawFixedWeirs == 2 .or. ndrawFixedWeirs == 4) then 
          call setcol(ncolblack)
          do i  = 1,nfxw
             L  = lnfxw(i) ; k3=lncn(1,L) ; k4=lncn(2,L)  
             xu =  0.5d0*( xk(k3) +  xk(k4) ) 
             yu =  0.5d0*( yk(k3) +  yk(k4) ) 
             if ( ndrawFixedWeirs == 4) then 
                call isocol(bob(1,L),ncol) 
             endif          
             call htext(bob(1,L) , xu, yu) 
          enddo
       endif

   endif
    
   call resetlinesizesetc()    

end subroutine plotFixedWeirs


!> Plots a cross section path on the screen.
!! Prior to a 'geominit' the original polyline path is shown, for an
!! initialized model the crossed flow links are highlighted.
subroutine plotCrossSectionPath(path, met, ncol, jaArrow, label)
    use m_crspath
    use m_wearelt
    type(tcrspath),   intent(in) :: path    !< Path definition
    integer,          intent(in) :: met     !< Method: 1=plot polyline, 2=plot crossed net/flow links (as stored in path%xk)
    integer,          intent(in) :: ncol    !< Drawing color
    character(len=*), intent(in) :: label   !< Text label to be displayed.
    integer,          intent(in) :: jaArrow !< Whether or not (1/0) to draw an outgoing arrow.

    integer :: i, j, jj, jmin, jmax, ncolb = 3
    double precision :: xjmin, xt, yt, rn, rt, dis, xx1, yy1, xx2, yy2, xx, yy
    logical :: inview
    character tex*40


    call setcol(ncol)

    ! If crs is not yet placed on flow links, just plot the coarse polyline:
    ! (or flow links are already known, but plotmethod==1)
    if (path%np > 0 .and. (path%lnx <= 0 .or. met==1)) then
        call movabs(path%xp(1), path%yp(1))
        !call cir(.4d0*rcir)
        jmax = 1 ! jmax is the last visible point in coarse polyline.
                 ! Only #1 is not checked (so user should zoom out when even that one is not visible)
        do j=2,path%np
            call lnabs(path%xp(j), path%yp(j))

            if (inview(path%xp(j), path%yp(j))) then  ! find first and last j in viewing area
                jmax = j
            end if
            !call cir(.4d0*rcir)
        end do
        xx2 = path%xp(jmax)
        yy2 = path%yp(jmax)
    ! Else, default: plot all crossed flow links in crs.
    else if (path%lnx > 0 .and. met==2)  then
        jmin  = 0
        jmax  = 0

        do j=1,path%lnx
            call movabs(path%xk(1,j), path%yk(1,j))
            call cir(.2d0*rcir)
            xx = path%xk(2,j) ; yy = path%yk(2,j)
            call lnabs(xx,yy)
            call cir(.2d0*rcir)

            if (inview(xx,yy)) then  ! find first and last j in viewing area
                jmax = j
                if (jmin == 0) jmin = j
            end if
        end do
        if (jmax == 0 .and. jmin == 0) then
           call LINEWIDTH(1)
           return
        else
           if ( path%xk(1,jmin) > path%xk(1,jmax) ) then
              jj = jmin ; jmin = jmax ; jmax = jj
           endif
           xx1 = path%xk(1,jmin); yy1 = path%yk(1,jmin)
           xx2 = path%xk(1,jmax); yy2 = path%yk(1,jmax)
        endif

        ! For a monitoring cross section, plot the positive direction
        ! as an arrow in view area, and show discharge or other quant.
        if (jaArrow == 1) then
            xt = .5d0*(path%xk(1,jmin) + path%xk(2,jmin))
            yt = .5d0*(path%yk(1,jmin) + path%yk(2,jmin))
            call normalout(path%xk(1,jmin), path%yk(1,jmin), path%xk(2,jmin), path%yk(2,jmin), rn, rt)
            call arrowsxy(xt,yt,rn,rt,4d0*rcir)
        endif

    end if ! path%lnx > 0
    if (len_trim(label) > 0) then
        call igrcharfont(7)
        xt = xx2 ; yt = yy2
        if (xt > x2 - dsix) xt = x2-dsix
        call gtext(trim(label), xt, yt, kltex)
    endif



end subroutine plotCrossSectionPath

subroutine thicklinetexcol(ncol)
integer :: ncol
call settextsizefac(1.5d0)
call LINEWIDTH(2)
call setcol(ncol)
end subroutine thicklinetexcol


subroutine resetlinesizesetc()
call settextsize()
CALL LINEWIDTH(1)
call igrcharfont(1)
end subroutine resetlinesizesetc

SUBROUTINE MINMXNS()
      USE M_BITMAP
      use network_data
      USE M_SAMPLES
      USE M_grid
      USE M_SPLINES
      implicit none
      double precision :: aspect
      double precision :: dx
      double precision :: dy
      integer :: n
      integer :: ndraw
      double precision :: xcmax, xcmin, xkmax, xkmin, xlmax, xlmin, xpmax, xpmin, xsmax, xsmin, xspmax, xspmin
      double precision :: ycmax, ycmin, ykmax, ykmin, ylmax, ylmin, ypmax, ypmin, ysmax, ysmin, yspmax, yspmin
      double precision :: xm, ym
      double precision ::  XH(10), YH(10)
      COMMON /DRAWTHIS/ NDRAW(40)
      CALL DMINMAX(   XLAN,  MXLAN,  XLMIN,   XLMAX, MAXLAN)
      CALL DMINMAX(   YLAN,  MXLAN,  YLMIN,   YLMAX, MAXLAN)

      CALL DMINMAX(   XK  ,  NUMK ,  XKMIN,   XKMAX, KMAX  )
      CALL DMINMAX(   YK  ,  NUMK ,  YKMIN,   YKMAX, KMAX  )

      CALL DMINMAX(   Xc  ,  mc*nc,  XCMIN,   XCMAX, mc*nc )
      CALL DMINMAX(   Yc  ,  mc*nc,  YCMIN,   YCMAX, mc*nc )

      CALL DMINMAX(  XSP  ,   mcS*MAXSPLEN,  XSPMIN, XSPMAX,  mcS*MAXSPLEN) ! SPLINES
      CALL DMINMAX(  YSP  ,   mcS*MAXSPLEN,  YSPMIN, YSPMAX,  mcS*MAXSPLEN)

      CALL  DMINMAX(  XPL ,  NPL  ,  XPMIN,   XPMAX, MAXPOL)
      CALL  DMINMAX(  YPL ,  NPL  ,  YPMIN,   YPMAX, MAXPOL)

      if ( NS.gt.0 ) then
         CALL  DMINMAX(  XS  ,  NS   ,  XSMIN,   XSMAX, NS    )
         CALL  DMINMAX(  YS  ,  NS   ,  YSMIN,   YSMAX, NS    )
      else
         xsmin =  huge(0d0)
         xsmax =  huge(0d0)
         ysmin =  huge(0d0)
         ysmax =  huge(0d0)
      end if

      IF (NDRAW(26) .EQ. 1) THEN
         XMIN = MIN(XMIN,XP(1))
         XMAX = MAX(XMAX,XP(2))
         YMIN = MIN(YMIN,YP(1))
         YMAX = MAX(YMAX,YP(4))
      ENDIF

      N = 0
      IF (XKMAX .NE. XKMIN .OR. YKMAX .NE. YKMIN) THEN
         N     = N+1
         XH(N) = XKMAX
         YH(N) = YKMAX
         N     = N+1
         XH(N) = XKMIN
         YH(N) = YKMIN
      ENDIF

      IF (XLMAX .NE. XLMIN .OR. YLMAX .NE. YLMIN) THEN
         N     = N+1
         XH(N) = XLMAX
         YH(N) = YLMAX
         N     = N+1
         XH(N) = XLMIN
         YH(N) = YLMIN
      ENDIF

      IF (XPMAX .NE. XPMIN .OR. YPMAX .NE. YPMIN) THEN
         N     = N+1
         XH(N) = XPMAX
         YH(N) = YPMAX
         N     = N+1
         XH(N) = XPMIN
         YH(N) = YPMIN
      ENDIF

      IF (XSMAX .NE. XSMIN .OR. YSMAX .NE. YSMIN) THEN
         N     = N+1
         XH(N) = XSMAX
         YH(N) = YSMAX
         N     = N+1
         XH(N) = XSMIN
         YH(N) = YSMIN
      ENDIF

      IF (XCMAX .NE. XCMIN .OR. YCMAX .NE. YCMIN) THEN
         N     = N+1
         XH(N) = XCMAX
         YH(N) = YCMAX
         N     = N+1
         XH(N) = XCMIN
         YH(N) = YCMIN
      ENDIF

      IF (XSPMAX .NE. XSPMIN .OR. YSPMAX .NE. YSPMIN) THEN
         N     = N+1
         XH(N) = XSPMAX
         YH(N) = YSPMAX
         N     = N+1
         XH(N) = XSPMIN
         YH(N) = YSPMIN
      ENDIF

      CALL DMINMAX(     XH,      N,   XMIN,   XMAX, 10)
      CALL DMINMAX(     YH,      N,   YMIN,   YMAX, 10)

      CALL INQASP(ASPECT)

      IF (XMAX .EQ. XMIN .and. YMAX .EQ. YMIN) THEN
         XMIN = 0d0     ; YMIN = 0d0
         XMAX = 1000d0  ; Ymax = aspect*1000d0
      ENDIF

      DX   =  XMAX - XMIN
      DY   =  YMAX - YMIN
      XM   =  XMIN + DX/2
      YM   =  YMIN + DY/2

      IF (DY .LT. ASPECT*DX) THEN
         XMIN = XM - 0.55d0*DX
         XMAX = XM + 0.65d0*DX
         YMIN = YM - 0.5d0*DY - 0.05d0*DX
      ELSE
         YMIN = YM - 0.6d0*DY
         XMIN = XM - 0.6d0*DY/ASPECT
         XMAX = XM + 0.6d0*DY/ASPECT
      ENDIF

      CALL WEAREL()

      RETURN
      END subroutine minmxns

end module unstruc_display

 subroutine zoomshift(nshift) ! based on polygon
 use unstruc_display
 use m_flowtimes
 use m_polygon
 implicit none
 integer :: nshift, ndraw, i1
 double precision :: dr, x00, y00, dxw, dyw, rshift, dshift

 COMMON /DRAWTHIS/ NDRAW(40)

 nshift = nshift + 1
 rshift = dble(nshift)/dble(numzoomshift)
 i1     = int(rshift) + 1
 i1     = min(i1,npl-1)
 dr     = rshift - i1 + 1
 x00    = (1d0 - dr)*xpl(i1) + dr*xpl(i1+1)
 y00    = (1d0 - dr)*ypl(i1) + dr*ypl(i1+1)
 dxw    = 0.5d0*(x2-x1)
 dyw    = 0.5d0*(y2-y1)
 x1     = x00 - dxw
 x2     = x00 + dxw
 y1     = y00 - dyw
 y2     = y00 + dyw
 call setwor(x1,y1,x2,y2)
 ndraw(10) = 1 ! wel plotten
 end subroutine zoomshift


subroutine tekship()
 use m_ship
 implicit none
 double precision :: sx1,sx2,sy1,sy2,css,sns, rr, cr, sr, snum
 integer :: n
 if (iniship == 0) return

 call setcol(4)

 do n = 1,nshiptxy
    css = cos(shi(n))     ; sns = sin(shi(n))

    call smovabs(n,  1.0d0 ,  0.0d0)
    call slnabs (n,  0.9d0 , -1.0d0)
    call slnabs (n, -1.0d0 , -1.0d0)
    call slnabs (n, -1.0d0 ,  1.0d0)
    call slnabs (n,  0.9d0 ,  1.0d0)
    call slnabs (n,  1.0d0 ,  0.0d0)

    snum =  css*fx2(n) + sns*fy2(n)               ! pressure force in shipL dir
    call shtext(n, snum*1d-4, -1.3d0, 0d0)

    snum = -sns*fx2(n) + css*fy2(n)               ! pressure force in shipB dir
    call shtext(n, snum*1d-4, -1.3d0, 1d0)

    snum = fm2(n)/shL(n)                          ! pressure mom vertical ax
    call shtext(n, snum*1d-4, -1.3d0,-1d0)

    snum =  css*fricx(n) + sns*fricy(n)           ! fric force in shipL dir
    call shtext(n, snum*1d-4,  1.3d0, 0d0)

    snum = -sns*fricx(n) + css*fricy(n)           ! fric force in shipB dir
    call shtext(n, snum*1d-4,  1.3d0, 1d0)

    snum = fricm(n)/shL(n)                        ! fric mom vertical ax
    call shtext(n, snum*1d-4,  1.3d0,-1d0)

    snum =  css*stuwx(n) + sns*stuwy(n)           ! stuwforce in shipL dir
    call shtext(n, snum*1d-4, -0.8d0, 0d0)

    snum = -sns*stuwx(n) + css*stuwy(n)           ! stuwforce in shipB dir
    call shtext(n, snum*1d-4, -0.8d0, 1d0)

    snum = stuwm(n)/shL(n)                        ! stuwmom vertical ax normalised by half length
    call shtext(n, snum*1d-4, -0.8d0,-1d0)

    snum =  css*shu(n) + sns*shv(n)               ! snelheid in shipL dir
    call shtext(n, snum, -0.d0, 0d0)

    snum = -sns*shu(n) + css*shv(n)               ! snelheid in shipB dir
    call shtext(n, snum, -0.0d0, 1.1d0)

    snum = sho(n)                                 ! ronjes/minuut vertical ax
    call shtext(n, snum*60d0/6.28d0, -0.d0,-1.1d0)



    sx2 = shx(n) - shL(n)*css                  ! rudder
    sy2 = shy(n) - shL(n)*sns
    call movabs(sx2, sy2)
    rr = 0.4d0*shb(n) ; cr = cos(shi(n) + roer(n)) ; sr = sin(shi(n) + roer(n))
    call lnabs(sx2 - rr*cr, sy2 - rr*sr)



 enddo
 end subroutine tekship


 
subroutine tekwindvector()
 USE m_wind
 use m_wearelt
 use unstruc_display
 use m_flow, only : qinrain 
 implicit none
 double precision :: xp, yp, vfw, ws
 
 if (jawind > 0 ) then 

    xp = 0.90*x1 + 0.10*x2
    yp = 0.15*y1 + 0.85*y2

    vfw = 0.1d0*(x2-x1)/10d0   ! 10 m/s is 0.1*screen

    call thicklinetexcol(ncolln)  
    call arrowsxy( xp, yp, windxav, windyav, vfw)
    ws = sqrt(windxav*windxav + windyav*windyav)   
    xp = 0.97*x1 + 0.03*x2
    yp = 0.25*y1 + 0.75*y2
    call DHTEXT(ws,xp, yp, 0d0)
    call resetlinesizesetc()
 endif
 
 if (jarain > 0) then   
    xp = 0.97*x1 + 0.03*x2    
    yp = 0.25*y1 + 0.71*y2
    call thicklinetexcol(ncolln)  
    call DHTEXT(qinrain ,xp, yp, 0d0)
    call resetlinesizesetc()
 endif
 
end subroutine tekwindvector