diff --git a/LINELABELS.xlsx b/LINELABELS.xlsx
new file mode 100644
index 0000000..a2d9fca
Binary files /dev/null and b/LINELABELS.xlsx differ
diff --git a/RMA10PRP/DUMMY.R10 b/RMA10PRP/DUMMY.R10
index e58b7e5..f8da133 100644
--- a/RMA10PRP/DUMMY.R10
+++ b/RMA10PRP/DUMMY.R10
@@ -1,100 +1,110 @@
INBNGEO MESH.GEO & Input binary file with geometric data
OUTFIL RESULTS.OUT & Output ASCII results
OUTBNRMA BRESULTS.RMA & Output binary file with new format results data
+OUTBNRST BRESULTS.RST & Output binary file with restart data
+OUTCON CONTIN.DAT & Output ASCII file with continuity line flow data
+OUTASXTR EXTRACT.DAT & Output ASCII file listing extracted results for specified nodes
ENDFIL
ENDLIMIT
-TI
-STARTIM & Starting date and time
-ENDTIM & Ending date and time
+TI TEST ESTUARY
+3DTYPLD2 & 3-d model with type LD2 input
+1DTYP 0 & 1-D cross section type (see users manual for more details)
+ELEV 0.000 & Reference level for vertical transformation / Starting water level
+STARTIM 01/01/2019 00.00 & Starting date and time
+OVERLAYT & Overlay input time/date onto restart file
+SSITN 5 & Number of steady state iteration cycles to be applied
+TRITN 5 & Number of time transient iteration cycles to be applied per time step
+TRSTEPS 4 & Number of transient steps to be simulated
+TIMPROJ 1 & Time projection switch 1= no projection, 0 = use time derivative (least stable), 2 = project u
+SOLVER 5 & Equation solver type 0=FRONT > 0 PARDISO SOLVER
+VERTEQSW 20000 & Number of vertical velocity equations to switch to PARDISO (Default=20000)
+SLVERDSK 0 & Switch to initiate out of local memory (disk) version of PARDISO for very large prblems
VERTTURB 0 & Vertical turbulence option
-MARSHOPT 0 & Marsh option 0=inactive -2 = input marsh parameters
EDDYTYPE 1 & Eddy type 0 = fixed, 1= scaled, 2 = Smagorinsky
TBFACT 0.200 & Turbulence factor when Smagorinsky is active
TBMIN 1.000 & Turbulence minimum when Smagorinsky is active
-TIMPROJ 1 & Time projection switch 1= no projection, 0 = use time derivative (least stable), 2 = project u
-2DMODEL & 2-d model approximation
-GRAVITY 1 & Units switch 1 = metric 0 = english
-ZEROBED 0 & Force zero bed velocity 1 = at all bed levels 2 = limited by bed level input
-ZERBELLV 0.000 & Force zero velocity on water surface
-PASSSAL 0 & Bed level below which zero bed level in applied
-PASSTMP 0 & Ignore influence of salinity on water density
-PASSSED 0 & Ignore influence of temperature on water density
-ZEROSURF 0 & Ignore influence of sediment on water density
-SAVITR 0 & Save binary file for all iterations 0 = ignore, 1 =save
REVRTDEN 0.000 & Option for vertical density model 0 = standard 1 = revised method
+SURFMIX 0.400 & Depth of surface mixing (associated with free water surface local turbulence
+TRANSIT 0.50000 & Transition depth for collapse from 3-D to 2-D approximation
+GRAVITY 1 & Units switch 1 = metric 0 = english
OMEGA 0.000 & Latitude (degrees) positive in northern hemisphere
-ELEV 0.000 & Reference level for vertical transformation / Starting water level
-X-SCALE 1.000 & x scale factor (0.0 equivalent to 1.0)
-Y-SCALE 1.000 & y scale factor (0.0 equivalent to 1.0)
-Z-SCALE 1.000 & z sacle factor (0.0 equivalent to 1.0)
-BOUNDCMN 1.000 & Defines constant associated with vertical constituent distribution at boundary
-BOUNDCPW 1.000 & Defines power associated with vertical constituent distribution at boundary
-UNOM 0.000 & Nominal water velocity at startup
-UDIR 0.000 & Nominal current direction in radians c-clockwise from horzontal
-HMIN 0.000 & Minimum depth as startup (ignored if set = 0.0) Use with care can cause initial instabilty
-DSET 0.200 & Depth at which drying (element removal) is initiated when drying
-DSETD 0.300 & Depth at which rewetting (element return) occurs when flooding
-FIXBC 0 & Switch that forces boundary conditions to be applied regardless of direction
-VERTEQSW 20000 & Number of vertical velocity equations to switch to PARDISO (Default=20000)
+PASSSAL 0 & Ignore influence of salinity on water density
+PASSTMP 0 & Ignore influence of temperature on water density
+PASSSED 0 & Ignore influence of sediment on water density
+ZEROBED 0 & Force zero bed velocity 1 = at all bed levels 2 = limited by bed level input
+ZERBELLV 0.000 & Bed level below which zero bed level in applied
+ZEROSURF 0 & Force zero velocity on water surface
+INITWSRF 0.000 & Initial water surface elevation, Required if different from tra nsformation reference level
INITSAL 0.000 & Initial salinity when no restart file specified
INITTEMP 20.000 & Initial temperature when no restart file specified
INITSED 0.000 & Initial sediment concentration when no restart file specified
-INITUBED 0.100 & Bed x-velocity used to compute resistance to flow for initial conditions regardless of initial
-INITVBED 0.100 & Bed y-velocity used to compute resistance to flow for initial conditions regardless of initial
+INITUBED 0.000 & Bed x-velocity used to compute resistance to flow for initial conditions regardless of initial
+INITVBED 0.000 & Bed y-velocity used to compute resistance to flow for initial conditions regardless of initial
RETRNPCT 0.000 & Percent of outflowing average concentration returned on next tidal inflow
-SURFMIX 0.500 & Depth of surface mixing (associated with free water surface local turbulence
-SSITN 20 & Number of steady state iteration cycles to be applied
-TRITN 20 & Number of time transient iteration cycles to be applied per time step
-TRSTEPS 12 & Number of transient steps to be simulated
+CONV-VEL 0.01000 & Convergence limit for velocities
+CONV-DEP 0.00100 & Convergence limit for depth
+CONV-SAL 0.01000 & Convergence limit for salinity
+CONV-TMP 0.01000 & Convergence limit for temperature
+CONV-SED 0.01000 & Convergence limit for sediment concentration
ECHSVALL & Show element and node data in echo print
ITERINTV 0 & Iteration frequency for ASCII output. Skip if = 0
ASCSVFRQ 1 & Time step frequency for ASCII output
STBIN 1 & Save startup to binary results file
-DRYSW 0 & Switch to initiate element elimination
FRQBIN 1 & Time step frequency for binary output
FRQRST 10 & Time step frequency for individual binary restart file
RWDASFRQ 999999 & Time step frequency for restarting ASCII results file (saves file size for very large number o
-1DTYP 1 & 1-D cross section type (see users manual for more details)
-SOLVER 4 & Equation solver type 0=FRONT > 0 PARDISO SOLVER
-SLVERDSK 0 & Switch to initiate out of local memory (disk) version of PARDISO for very large prblems
MESSRWND 999999 & Time step frequency for restarting ASCII message file (saves file size for very large number o
-INITWSRF 0.000 & Initial water surface elevation, Required if different from tra nsformation reference level
-CONV-VEL 0.01000 & Convergence limit for velocities
-CONV-DEP 0.00100 & Convergence limit for depth
-CONV-SAL 0.01000 & Convergence limit for salinity
-CONV-TMP 0.01000 & Convergence limit for temperature
-CONV-SED 0.01000 & Convergence limit for sediment concentration
- Continuity lines for selective output
-CCLINOUT 1 4 5
+SAVITR 0 & Save binary file for all iterations 0 = ignore, 1 =save
+ & Nodes for extraction of all constiuents to a time series file
+SPR-NOD 23 45 56
+ & Continuity lines for selective output
+CCLINOUT 1
CONT-FLW & Save flows in continuity line output
+CONT-DEP & Save depths in continuity line output
CONT-ELV & Save average water surface elevations in continuity line output
+CONT-SAL & Save average salinity in continuity line output
CONT-TMP & Save average temperature in continuity line output
CONT-SED & Save average sediment concentration in continuity line output
-EDDY-MAT 1 -0.20000 & Horizontal turbulent eddy coefficient
-EDDY-VRT 1 0.000E+00 & Vertical turbulent eddy coefficient
-DIFF-MAT 1 0.00000 & Horizontal turbulent diffusion coefficient
-DIFF-VRT 1 0.000E+00 & Vertical turbulent diffusion coefficient
+X-SCALE 1.000 & x scale factor (0.0 equivalent to 1.0)
+Y-SCALE 1.000 & y scale factor (0.0 equivalent to 1.0)
+Z-SCALE 1.000 & z sacle factor (0.0 equivalent to 1.0)
+DRYSW 0 & Switch to initiate element elimination
+DSET 0.200 & Depth at which drying (element removal) is initiated when drying
+DSETD 0.300 & Depth at which rewetting (element return) occurs when flooding
+FIXBC 0 & Switch that forces boundary conditions to be applied regardless of direction
+UNOM 0.250 & Nominal water velocity at startup
+UDIR 0.000 & Nominal current direction in radians c-clockwise from horzontal
+HMIN 0.000 & Minimum depth as startup (ignored if set = 0.0) Use with care can cause initial instabilty
+EDDY-MAT 1 0.20000 & Horizontal turbulent eddy coefficient
+EDDY-VRT 1 1.000E-01 & Vertical turbulent eddy coefficient
+DIFF-MAT 1 0.10000 & Horizontal turbulent diffusion coefficient
+DIFF-VRT 1 1.000E-03 & Vertical turbulent diffusion coefficient
MANN-MAT 1 0.02000 & Manning/Chezy bed friction coefficient >1.0 is Chezy coefficient
MANN-BNK 1 0.00000 & Bank Manning coefficient if appropriate
MANN-SRF 1 0.00000 & Water surface Manning coefficient if appropriate
-MARSH-FR 1 0.00000 & Marsh factor when drying applied to friction coefficient
+MARSH-FR 1 10.00000 & Marsh factor when drying applied to friction coefficient
ED-VTVAR 1 1.000 0.000 0.000 & Eqn for vertical distribution factor for horizontal eddy coefficient F=a+z*(b+c*z)
+MARSHOPT -2 & Marsh option 0=inactive -2 = input marsh parameters
+MARSH-SH 2.000 & Marsh coefficient depth shift
+MARSH-RG 0.500 & Marsh coefficient range
+MARSH-PR 0.010 & Marsh coefficient porosity
+MARSH-LM 0.000 & Marsh coefficient limit
+ & Nodal layer data type LD2, 0 = apply to all nodes, number of layers Proportional factors
+LD2 0 4 1.000 2.000 3.000 4.000
+BOUNDCMN 1.000 & Defines constant associated with vertical constituent distribution at boundary
+BOUNDCPW 1.000 & Defines power associated with vertical constituent distribution at boundary
BN-V-MIN 1.000 & Boundary velocity vertical distribution constant
BN-V-PWR 1.000 & Boundary velocity vertical distribution power
TR-V-MIN 1.000 & 2-D to 3-D velocity vertical distribution constant
TR-V-PWR 1.000 & 2-D to 3-D velocity vertical distribution power
-MARSH-SH 1.500
-MARSH-RG 0.670
-MARSH-PR 0.040
-MARSH-LM 0.000
-TAB-END
-ENDGEO
+TAB-END & Indicator of end of tabular data set
+ENDGEO & Indicator of end of end control and geometry data
DELTA-HR 0.0000 & Time step in minutes
CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
ENDSTEP
DELTA-MN 15.0000 & Time step in hours
-DATE-END 31/01/2018 12.00 & Ending time for the time step block
+DATE-END 31/01/2019 12.00 & Ending time for the time step block
CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
ENDSTEP
diff --git a/RMA10PRP/New Text Document.txt b/RMA10PRP/New Text Document.txt
new file mode 100644
index 0000000..387949a
--- /dev/null
+++ b/RMA10PRP/New Text Document.txt
@@ -0,0 +1,4 @@
+SED SETTLING
+IOV
+TRANSIT
+
diff --git a/RMA10PRP/RMA10PRP.u2d b/RMA10PRP/RMA10PRP.u2d
index e69de29..d45a425 100644
Binary files a/RMA10PRP/RMA10PRP.u2d and b/RMA10PRP/RMA10PRP.u2d differ
diff --git a/RMA10PRP/RMA10PRP.vfproj b/RMA10PRP/RMA10PRP.vfproj
index 62f017f..85c8e30 100644
--- a/RMA10PRP/RMA10PRP.vfproj
+++ b/RMA10PRP/RMA10PRP.vfproj
@@ -37,7 +37,7 @@
-
+
diff --git a/RMA10PRP/SRC/BLKRM10.F90 b/RMA10PRP/SRC/BLKRM10.F90
index a313e23..da8c2f8 100644
--- a/RMA10PRP/SRC/BLKRM10.F90
+++ b/RMA10PRP/SRC/BLKRM10.F90
@@ -13,13 +13,13 @@
INTEGER IWIND(25),IWUNITS(25), NRITR(25)
INTEGER ITYP(1000),INTG(100),IDEFLT(100),IQLIN(20,25),IQLIN1(20,25),IHLIN(20,25),IELIN(500,25),LAYNUMH(20,25),LAYNUMQ(20,25),LAYNUME(20,25),IITYPE(25)
INTEGER qcklin(20,25),hcklin(20,25),ecklin(500,25),ecklin1(500,25),ecklin2(500,25),IETP(500,20),hcklin1(20,25),IELIN1(500,25)
- INTEGER ICFLOWS(20),ICELEV(20),ICDEP(20),ISPRTN(20),NCFLOWS,NCELEV,NCDEP,NSPRTN
+ INTEGER ICFLOWS(20),ICELEV(20),ICDEP(20),NCFLOWS,NCELEV,NCDEP,NSPRTN
INTEGER ISFLIN(20,25),NSFLIN(25),IFCLIN(20,25),IFCLIN1(20,25),NFCLIN(25)
INTEGER LIMITVAL(13),LIMITDFT(13),IWTYP(1000),NWR,NHOLDC,ICPLIN(27)
INTEGER IELTOPIP(10),IPIPTOEL(10),JET,JGATE,L3DCNT,L3DTYP
INTEGER IBF1(1000),IEXTYP(1000),IGATE(1000,2)
INTEGER JND(1000),KKLAY(1000)
- INTEGER NCSPLPT(30)
+ INTEGER NCSPLPT(30),NSPLPT(10)
REAL L3DDAT(1000,40),BCSND(500,25,6)
INTEGER IITYP(1000)
INTEGER KHCN(25),IHCN(20,25)
diff --git a/RMA10PRP/SRC/CONVERT.F90 b/RMA10PRP/SRC/CONVERT.F90
index 3ad45f7..8b459d8 100644
--- a/RMA10PRP/SRC/CONVERT.F90
+++ b/RMA10PRP/SRC/CONVERT.F90
@@ -4,7 +4,7 @@
USE BLKMAT
USE BLKRM10
integer inrm1,k
- CHARACTER*8 LABELF(43),labelin,ID
+ CHARACTER*8 LABELF(44),labelin,ID
character*72 namin,DLIN
CHARACTER*8 DLINEXTRA
COMMON /DLINF/ DLINEXTRA
@@ -16,7 +16,7 @@
INTEGER IMTS(12)
DATA IMTS/0,31,59,90,120,151,181,212,243,273,304,334/
DATA LABELF/&
- 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELFL ','INBNELF ','INHARM ','METFIL '&
+ 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELTFL ','INBNELF ','INHARM ','INLAYDAT','METFIL '&
,'OUTFIL ','OUTBNRMA','OUTBNRST','OUTBN3GE','OUTCON ','OUTMET ','TIMFIL ','BWINDIN ','AWINDIN ','OUTBNXTR'&
,'OUTASXTR','OUTBNELF','INSRCORD','INBNWGT ','INBNSTRS','OUTBNWGT','OUTWGT ','GROUPNUM','STFLFIL ','INASTRAT'&
,'BCFIL ','VOLFIL ','OUTBNRES','VELBNFIL','IN3DBNGE','OUTBN2GE','INSMSGN ','OUTSMS '&
@@ -34,7 +34,7 @@
endif
100 CONTINUE
- DO J=1,43
+ DO J=1,45
DO K=1,72
FNAME(J)(K:k)=' '
ENDDO
@@ -94,7 +94,7 @@
OPEN(INRM2,FILE=NAMEFL,FORM='FORMATTED')
! Now read files
- DO J=1,44
+ DO J=1,45
IF(LABELIN(1:6) .EQ. 'ENDFIL') THEN
NFIL=J-1
EXIT
@@ -195,6 +195,9 @@
ELSEIF(ID(1:2) .EQ. 'C2') THEN
!IPK DEC07 REPLACE OMEGA BY OMEGA1
READ(DLIN,5020) OMEGA1,ELEV,XSCALE,YSCALE,ZSCALE
+ IF(XSCALE .EQ. 0.) XSCALE=1.0
+ IF(YSCALE .EQ. 0.) YSCALE=1.0
+ IF(ZSCALE .EQ. 0.) ZSCALE=1.0
ELEV1=ELEV
5020 FORMAT(7F8.0,2I8)
write(75,*) 'read c2'
@@ -204,7 +207,10 @@
READ(DLIN,5020) CMIN,CPR,UNOM,UDIR,HMNN,DSET,DSETD,IFXBC,IEQSWT
!IPK NOV97 READ(LIN,7000) ID,DLIN
+ IF(CMIN .EQ. 0.) CMIN=1.0
+ IF(CPR .EQ. 0.) CPR=1.0
IF(UNOM .EQ. 0.) UNOM=0.25
+ IF(IEQSWT .EQ. 0) IEQSWT=20000
GO TO 250
ELSEIF(ID(1:2) .EQ. 'C4') THEN
@@ -217,7 +223,8 @@
ELSEIF(ID(1:2) .EQ. 'C5') THEN
READ(DLIN,5011) NITI,NITN,TSTART,NCYC,IPRT,NPRTI,NPRTF,IRSAV,IDSWT
- 5011 FORMAT( 2I8,F8.0,6I8)
+5011 FORMAT( 2I8,F8.0,6I8)
+ IF(IRSAV .EQ. 0) IRSAV=1
write(75,*) 'read c5'
@@ -227,6 +234,10 @@
!IPK APR11
READ(DLIN,'(8I8)') IOUTFREQ,IOUTRST,IREWOUT,ID1DND,ICPU,IOOC,IREWMES,IDEBUG1
+ IF(IOUTRST .EQ. 0) IOUTRST=10
+ IF(IOUTFREQ .EQ. 0) IOUTFREQ=1
+ IF(IREWOUT .EQ. 0) IREWOUT=999999
+ IF(IREWMES .EQ. 0) IREWMES=999999
!IPK JAN09
GO TO 250
@@ -370,17 +381,17 @@
ENDDO
GO TO 250
- ! ELSEIF(ID(1:3) .EQ. 'SP ') THEN
- ! READ(DLIN,5010) (NSPLPT(J),J= 1,9)
- ! READ(DLINEXTRA,'(I8)') NSPLPT(10)
- ! DO J=10,1,-1
- ! IF(NSPLPT(J) .GT. 0) THEN
- ! JSPLPT=J
- ! GO TO 22
- ! ENDIF
- ! ENDDO
- !22 CONTINUE
- ! GO TO 250
+ ELSEIF(ID(1:3) .EQ. 'SP ') THEN
+ READ(DLIN,5010) (NSPLPT(J),J= 1,9)
+ READ(DLINEXTRA,'(I8)') NSPLPT(10)
+ DO J=10,1,-1
+ IF(NSPLPT(J) .GT. 0) THEN
+ JSPLPT=J
+ EXIT
+ ENDIF
+ ENDDO
+ CONTINUE
+ GO TO 250
ELSEIF(ID(1:4) .EQ. 'CSPC') THEN
READ(DLIN,'(9A8)') (ACSPLPT(J),J= 1,8)
@@ -426,6 +437,7 @@
5031 FORMAT(8X,8F8.0)
5032 FORMAT(8X,I8,2F8.0)
READ(DLIN,5030) J,(ORT(J,K),K=1,7)
+ ITYP(J)=J
write(75,*) 'read ed1'
IF(NMAT .LT. J) NMAT=J
diff --git a/RMA10PRP/SRC/DEFAULTS.F90 b/RMA10PRP/SRC/DEFAULTS.F90
index b632aa1..eccb4e5 100644
--- a/RMA10PRP/SRC/DEFAULTS.F90
+++ b/RMA10PRP/SRC/DEFAULTS.F90
@@ -67,8 +67,8 @@
CPR=1.0
CMINDF=1.0
CPRDF=1.0
- UNOM=0.
- UNOMDF=0.
+ UNOM=0.25
+ UNOMDF=0.25
UDIR=0.
UDIRDF=0.
HMIN=0.0
@@ -87,10 +87,10 @@
TEMPIDF=20.
SEDI=0.
SEDIDF=0.
- UINP=0.1
- UINPDF=0.1
- VINP=0.1
- VINPDF=0.1
+ UINP=0.0
+ UINPDF=0.0
+ VINP=0.0
+ VINPDF=0.0
prcnt=0.
prcntDF=0.
DMIX=0.5
@@ -119,10 +119,10 @@
IOUTRSTDF=10
IREWOUT=999999
IREWOUTDF=999999
- ID1DND=1
- ID1DNDDF=1
- ICPU=4
- ICPUDF=4
+ ID1DND=0
+ ID1DNDDF=0
+ ICPU=5
+ ICPUDF=5
IOOC=0
IOOCDF=0
IREWMES=999999
@@ -136,6 +136,12 @@
CONV(4)=0.01
CONV(5)=0.01
CONV(6)=0.01
+ CONVDF(1)=0.01
+ CONVDF(2)=0.01
+ CONVDF(3)=0.001
+ CONVDF(4)=0.01
+ CONVDF(5)=0.01
+ CONVDF(6)=0.01
IDRPT=0
DRFACT=0.1
IOV=0
diff --git a/RMA10PRP/SRC/FORMFILE - Copy.F90 b/RMA10PRP/SRC/FORMFILE - Copy.F90
new file mode 100644
index 0000000..f1f67ae
--- /dev/null
+++ b/RMA10PRP/SRC/FORMFILE - Copy.F90
@@ -0,0 +1,195 @@
+
+
+
+
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6091) helpint(95)
+6091 FORMAT('TAB-END',t70,a96)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6092) helpint(96)
+6092 FORMAT('ENDGEO',t70,a96)
+
+ DO J=1,JT
+ LINE=LINE+1
+ IF(DELT(J) .EQ. 0 .AND. DELTM(J) .EQ. 0. .OR. DELT(J) .GT. 0.) THEN
+ WRITE(DATALIN(LINE),6100) DELT(J),HELPTIME(1)
+ ELSE
+ WRITE(DATALIN(LINE),6101) DELTM(J),HELPTIME(2)
+ ENDIF
+6100 FORMAT('DELTA-HR',F12.4,T70,A162)
+6101 FORMAT('DELTA-MN',F12.4,T70,A162)
+ IF(DATENDST(J) .NE. ' ' .AND. J .GT. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6102) DATENDST(J),HELPTIME(3)
+6102 FORMAT('DATE-END',4X,A16,T70,A162)
+ ENDIF
+ IF(IITYPE(J) .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6103),HELPTIME(4)
+6103 FORMAT('Q-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6104),HELPTIME(5)
+6104 FORMAT('QSL-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 3) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6105),HELPTIME(6)
+6105 FORMAT('QTM-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 4) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6106),HELPTIME(7)
+6106 FORMAT('QSD-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 5) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6107),HELPTIME(8)
+6107 FORMAT('QST-ITR',T70,A162)
+ ENDIF
+ KH=NHLIN(J)
+ KQ=NQLIN(J)
+ KE=NELIN(J)
+ DO K=1,KH
+ LINE=LINE+1
+ IF(HLIN1(K,J) .GT. -9000.) THEN
+ IF(HCKLIN(K,J) .EQ. 1 ) THEN
+ WRITE(DATALIN(LINE),6108) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),HLIN1(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(9)
+6108 FORMAT('FL-ELEVI',2I8,2F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6109) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),HLIN1(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(10)
+6109 FORMAT('CN-ELEVI',2I8,2F8.4,3F8.3,T70,A162)
+ ENDIF
+ ELSE
+ IF(HCKLIN(K,J) .EQ. 1 ) THEN
+ WRITE(DATALIN(LINE),6110) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(11)
+6110 FORMAT('FL-ELEVC',2I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6111) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(12)
+6111 FORMAT('CN-ELEVC',2I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDIF
+ ENDDO
+
+ IF(KHCN(J) .GT. 0) THEN
+ DO K=1,KHCN(J)
+ L1=MIN(IHCN(K,J),9)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=1,L1)
+61111 FORMAT('HCN',5X,9F8.3)
+ IF(IHCN(K,J) .LE. 9) CYCLE
+ L1=MIN(IHCN(K,J),18)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=10,L1)
+ IF(IHCN(K,J) .LE. 18) CYCLE
+ L1=MIN(IHCN(K,J),27)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=19,L1)
+ IF(IHCN(K,J) .LE. 27) CYCLE
+ L1=MIN(IHCN(K,J),36)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=28,L1)
+ ENDDO
+ ENDIF
+
+ IF(WVEL(J) .NE. 0.) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61112) WVEL(J),WDIR(J)
+61112 FORMAT('GBWIND ',1X,2F8.3)
+ ENDIF
+
+ DO K=1,KQ
+ LINE=LINE+1
+ IF(IQLIN1(K,J) .EQ. 1) ALAB='C'
+ IF(IQLIN1(K,J) .EQ. 2) ALAB='F'
+ IF(IQLIN1(K,J) .EQ. 3) ALAB='H'
+ IF(QCKLIN(K,J) .EQ. 0) THEN
+ BLAB='CN'
+ ELSE
+ BLAB='FL'
+ ENDIF
+ IF(ALAB .EQ. 'C') THEN
+ IHP=13
+ ELSEIF(ALAB .EQ. 'F') THEN
+ IHP=14
+ ELSE
+ IHP=15
+ ENDIF
+ WRITE(DATALIN(LINE),6112) BLAB,ALAB,IQLIN(K,J),LAYNUMQ(K,J),QLIN(K,J),QDLIN(K,J),(QUALKQ(K,J,KK),KK=1,3),HELPTIME(IHP)
+6112 FORMAT(A2,'-Q',A1,3X,2I8,2F10.2,3F8.3,T70,A162)
+ ENDDO
+
+ DO K=1,KE
+ IF(ECKLIN(K,J) .EQ. 1) THEN
+ BLAB='FL'
+ ELSE
+ BLAB='CN'
+ ENDIF
+ LINE=LINE+1
+ IF(IETP(K,J) .EQ. 4) THEN
+ WRITE(DATALIN(LINE),6113) BLAB,ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(16)
+6113 FORMAT(A2,'-RAIN ',F8.4,3F8.3,T70,A162)
+ ELSEIF(IETP(K,J) .EQ. 1) THEN
+ IF(IELIN(K,J) .GT. 0) THEN
+ IF(ECKLIN1(K,J) .EQ. 1) THEN
+ WRITE(DATALIN(LINE),6116) BLAB,IELIN(K,J),LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(19)
+6116 FORMAT(A2,'-ELMET',2I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6117) BLAB,IELIN(K,J),LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(20)
+6117 FORMAT(A2,'-ELMED',2I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ELSE
+ IF(ECKLIN1(K,J) .EQ. 1) THEN
+ WRITE(DATALIN(LINE),6115) BLAB,LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(18)
+6115 FORMAT(A2,'-ELMGT',I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6114) BLAB,LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(17)
+6114 FORMAT(A2,'-ELMGD',I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDIF
+ ELSEIF(IETP(K,J) .EQ. 6) THEN
+ WRITE(DATALIN(LINE),6118) BLAB,IELIN(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(21)
+6118 FORMAT(A2,'-GRAIN',I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDDO
+
+ IF(NSND(J) .GT. 0) THEN
+ DO K=1,NSND(J)
+ IF(NODSND(K,J) .GT. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6119) NODSND(K,J),NCODE(K,J),(BCSND(K,J,L),L=1,6)
+6119 FORMAT('SN',8X,2I7,2f10.5,4F8.3)
+ ELSE
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61191) NCODE(K,J),(BCSND(K,J,L),L=1,6)
+61191 FORMAT('SD'15X,I7,2F10.5,4F8.2)
+ ENDIF
+ ENDDO
+ ENDIF
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6200)
+6200 FORMAT('ENDSTEP')
+
+ ENDDO
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6999)
+6999 FORMAT('ENDDATA')
+
+ filter="r10 file *.r10|*.R10|"
+ IOUT=66
+
+ CALL WSelectFile(Filter,SaveDialog+PromptOn+AppendExt+DirChange,NAMEFL,'Save R10 File for Execution')
+
+ IF (WInfoDialog(ExitButtonCommon).EQ.CommonOpen) THEN
+ go to 500
+ else
+ RETURN
+ endif
+
+500 CONTINUE
+ OPEN(IOUT,FILE=NAMEFL,FORM='FORMATTED')
+ DO K=1,LINE
+ WRITE(IOUT,'(A230)') DATALIN(K)
+ ENDDO
+ CLOSE(IOUT)
+ RETURN
+ END
\ No newline at end of file
diff --git a/RMA10PRP/SRC/FORMFILE.F90 b/RMA10PRP/SRC/FORMFILE.F90
index 6efaf5e..e34bf4d 100644
--- a/RMA10PRP/SRC/FORMFILE.F90
+++ b/RMA10PRP/SRC/FORMFILE.F90
@@ -5,11 +5,11 @@
USE RESID
INTEGER ICROSREF(45)
- CHARACTER*8 LABELF(43),LABELC(26),LABELFL(21),LABELCE,LABELFE,LIMITL(8),LABELNDL
+ CHARACTER*8 LABELF(44),LABELC(26),LABELFL(21),LABELCE,LABELFE,LIMITL(8),LABELNDL
CHARACTER(LEN=255) :: filter
- character*137 helpint(110)
+ character*137 helpint(97)
character*162 helptime(43)
- character*89 helpfile(44)
+ character*89 helpfile(45)
CHARACTER*1 ILAB(4),ALAB
CHARACTER*2 BLAB
@@ -19,7 +19,7 @@
,59,36,60,61,62/
! DATA LABELF/'INBNGEO ','INBNRST ','INRM1 ','INELEV ','INHYD ','INELTFL ','OUTFIL ','OUTBNRMA','OUTBNRST','TIMFIL ','OUTSMS ','OUTSPL ','OUTCON ','INCROS ','INSRCORD','INSTRESS','INBNSTRS','INWGT ','INBNWGT ','OUTWGT ','OUTBNWGT','OUTBNSTR','WEIRDATA','AWINDIN ','BWINDIN ','OUTFLOW ','INHARM ','OUTBN3GE'/
DATA LABELF/&
- 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELFL ','INBNELF ','INHARM ','METFIL '&
+ 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELTFL ','INBNELF ','INHARM ','INLAYDAT','METFIL '&
,'OUTFIL ','OUTBNRMA','OUTBNRST','OUTBN3GE','OUTCON ','OUTMET ','TIMFIL ','BWINDIN ','AWINDIN ','OUTBNXTR'&
,'OUTASXTR','OUTBNELF','INSRCORD','INBNWGT ','INBNSTRS','OUTBNWGT','OUTWGT ','GROUPNUM','STFLFIL ','INASTRAT'&
,'BCFIL ','VOLFIL ','OUTBNRES','VELBNFIL','IN3DBNGE','OUTBN2GE','INSMSGN ','OUTSMS '&
@@ -103,7 +103,7 @@ DATA HELPINT/&
'& Convergence limit for salinity ',&
'& Convergence limit for temperature ',&
'& Convergence limit for sediment concentration ',&
- '& Transition depth forl collapse from 3-D to 2-D approximation ',&
+ '& Transition depth for collapse from 3-D to 2-D approximation ',&
'& Save flows in continuity line output ',&
'& Save depths in continuity line output ',&
'& Save average water surface elevations in continuity line output ',&
@@ -123,25 +123,12 @@ DATA HELPINT/&
'& Boundary velocity vertical distribution power ',&
'& 2-D to 3-D velocity vertical distribution constant ',&
'& 2-D to 3-D velocity vertical distribution power ',&
- 'Pipe diameter for this element type ',&
- 'Pipe Manning coefficient at nominal water levels for this element type ',&
- 'Pipe Manning coefficient when surcharging starts for this element type ',&
- 'Pipe entrance loss factor for this element type ',&
- 'Pipe exit losses factor for this element type ',&
- 'Number of pipes for this element type ',&
- 'Culvert height for this element type ',&
- 'Culvert Manning coefficient at nominal water level for this element type ',&
- 'Culvert Manning coefficient when surcharging starts for this element type ',&
- 'Culvert entrance loss factor for this element type ',&
- 'Culvert exit losses factor for this element type ',&
- 'Number of Culverts for this element type ',&
- 'Variable Manning coefficient transition lower level for this element type ',&
- 'Variable Manning coefficient transition lower level Manning coefficient for this element type ',&
- 'Variable Manning coefficient transition upper level for this element type ',&
- 'Variable Manning coefficient transition upper level Manning coefficient for this element type ',&
- 'Variable Manning coefficient slope of Manning coefficient above upper level ',&
- 'Gate reference node 1 for gate opening decisions for this element type ',&
- 'Gate reference node 2 for gate opening decisions for this element type Gate closed when elev at node 1 ,elev at node2 ',&
+ '& Marsh coefficient depth shift ',&
+ '& Marsh coefficient range ',&
+ '& Marsh coefficient porosity ',&
+ '& Marsh coefficient limit ',&
+ '& Indicator of end of tabular data set ',&
+ '& Indicator of end of end control and geometry data ',&
' '/
DATA HELPTIME/&
'& Time step in minutes ',&
@@ -198,6 +185,7 @@ DATA HELPINT/&
'& Input ASCII file with element inflow data',&
'& Input binary element inflow data',&
'& Input ASCII file with elevation harmonic boundary data',&
+ '& Input ASCII file with layer data',&
'& Input ASCII file with meteorological data',&
'& Output ASCII results',&
'& Output binary file with new format results data',&
@@ -243,7 +231,7 @@ DATA HELPINT/&
ISHOWSW=0
ENDIF
- DO K=1,43
+ DO K=1,44
IF(FNAME(K)(1:16) .NE. ' ') THEN
LINE=LINE+1
DO KK=1,43
@@ -285,8 +273,10 @@ DATA HELPINT/&
WRITE(DATALIN(LINE),60002) LABELC(1),DATEST,helpint(1)
60002 FORMAT(A8,3X,A16,t70,a96)
IF(JT .GT. 1 .OR. DELT(1) .GT. 0. .OR. DELTM(1) .GT. 0. .or. DELT(2) .GT. 0. .OR. DELTM(2) .GT. 0.) THEN
- LINE=LINE+1
- WRITE(DATALIN(LINE),60002) LABELC(2),DATEND,helpint(2)
+ IF(DATEND .NE. ' ') THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60002) LABELC(2),DATEND,helpint(2)
+ ENDIF
ENDIF
! IF(IYRR .NE. IYRRDF .OR. ISHOWSW .EQ. 1) THEN
@@ -366,27 +356,27 @@ DATA HELPINT/&
ENDIF
IF(ZERBED .NE. ZERBEDDF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6022) ZERBED,helpint(24)
+ WRITE(DATALIN(LINE),6022) ZERBED,helpint(20)
6022 FORMAT('ZERBELLV',F8.3,t70,a96)
ENDIF
IF(IPASS1 .NE. IPASS1DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6018) IPASS1,helpint(20)
+ WRITE(DATALIN(LINE),6018) IPASS1,helpint(21)
6018 FORMAT('PASSSAL ',I8,t70,a96)
ENDIF
IF(IPASS2 .NE. IPASS2DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6019) IPASS2,helpint(21)
+ WRITE(DATALIN(LINE),6019) IPASS2,helpint(22)
6019 FORMAT('PASSTMP ',I8,t70,a96)
ENDIF
IF(IPASS3 .NE. IPASS3DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6020) IPASS3,helpint(22)
+ WRITE(DATALIN(LINE),6020) IPASS3,helpint(23)
6020 FORMAT('PASSSED ',I8,t70,a96)
ENDIF
IF(IZERS .NE. IZERSDF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6021) IZERS,helpint(23)
+ WRITE(DATALIN(LINE),6021) IZERS,helpint(24)
6021 FORMAT('ZEROSURF',I8,t70,a96)
ENDIF
IF(ISAVITR .NE. ISAVITRDF .OR. ISHOWSW .EQ. 1) THEN
@@ -620,7 +610,7 @@ DATA HELPINT/&
IF(ITRANSIT .EQ. 1) THEN
LINE=LINE+1
WRITE(DATALIN(LINE),60651) TRANSIT,helpint(71)
-60651 FORMAT('TRANSIT',3X,F8.5,t70,a96)
+60651 FORMAT('TRANSIT',1X,F8.5,t70,a96)
ENDIF
IF(ICTOTC .GT. 0) THEN
LINE=LINE+1
@@ -721,45 +711,58 @@ DATA HELPINT/&
6084 FORMAT('TR-V-PWR',F8.3,t70,a96)
ENDIF
IF(L3DTYP .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60841)
+60841 FORMAT(30X,'& Nodal layer data type LD2, 0 = apply to all nodes, number of layers Proportional factors')
DO K=1,L3DCNT
LINE=LINE+1
WRITE(DATALIN(LINE),6085) JND(K),KKLAY(K),(L3DDAT(K,J),J=1,KKLAY(K))
6085 FORMAT('LD2 ',2I8,20F8.3)
ENDDO
ELSEIF(L3DTYP .EQ. 3) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60851)
+60851 FORMAT(30X,'& Nodal layer data type LD3, 0 = apply to all nodes, number of layers Layer elevations')
LINE=LINE+1
DO K=1,L3DCNT
WRITE(DATALIN(LINE),6086) JND(K),KKLAY(K),(L3DDAT(K,J),J=1,KKLAY(K))
6086 FORMAT('LD3 ',2I8,20F8.3)
ENDDO
+ ELSEIF(JSPLPT .GT. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60861)
+60861 FORMAT(30X,'& Nodes for extraction of all constiuents to a time series file')
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60862) (NSPLPT(J),J=1,JSPLPT)
+60862 FORMAT('SPR-NOD ',1X,10I7)
ENDIF
IF(ACT1 .NE. ACT1DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6087) ACT1
-6087 FORMAT('MARSH-SH',F8.3)
+ WRITE(DATALIN(LINE),6087) ACT1,helpint(91)
+6087 FORMAT('MARSH-SH',F8.3,t70,a96)
ENDIF
IF(ACT2 .NE. ACT2DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6088) ACT2
-6088 FORMAT('MARSH-RG',F8.3)
+ WRITE(DATALIN(LINE),6088) ACT2,helpint(92)
+6088 FORMAT('MARSH-RG',F8.3,t70,a96)
ENDIF
IF(ACT3 .NE. ACT3DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6089) ACT3
-6089 FORMAT('MARSH-PR',F8.3)
+ WRITE(DATALIN(LINE),6089) ACT3,helpint(93)
+6089 FORMAT('MARSH-PR',F8.3,t70,a96)
ENDIF
IF(ACT4 .NE. ACT4DF .OR. ISHOWSW .EQ. 1) THEN
LINE=LINE+1
- WRITE(DATALIN(LINE),6090) ACT4
-6090 FORMAT('MARSH-LM',F8.3)
+ WRITE(DATALIN(LINE),6090) ACT4,helpint(94)
+6090 FORMAT('MARSH-LM',F8.3,t70,a96)
ENDIF
LINE=LINE+1
- WRITE(DATALIN(LINE),6091)
-6091 FORMAT('TAB-END')
+ WRITE(DATALIN(LINE),6091) helpint(95)
+6091 FORMAT('TAB-END',t70,a96)
LINE=LINE+1
- WRITE(DATALIN(LINE),6092)
-6092 FORMAT('ENDGEO')
+ WRITE(DATALIN(LINE),6092) helpint(96)
+6092 FORMAT('ENDGEO',t70,a96)
DO J=1,JT
LINE=LINE+1
diff --git a/RMA10PRP/SRC/FORMFILEV2.F90 b/RMA10PRP/SRC/FORMFILEV2.F90
new file mode 100644
index 0000000..690d0fc
--- /dev/null
+++ b/RMA10PRP/SRC/FORMFILEV2.F90
@@ -0,0 +1,937 @@
+ SUBROUTINE FORMFILE
+
+ use winteracter
+ USE BLKRM10
+
+ USE RESID
+ INTEGER ICROSREF(45)
+ CHARACTER*8 LABELF(44),LABELC(26),LABELFL(21),LABELCE,LABELFE,LIMITL(8),LABELNDL
+ CHARACTER(LEN=255) :: filter
+ character*137 helpint(98)
+ character*162 helptime(43)
+ character*89 helpfile(45)
+ CHARACTER*1 ILAB(4),ALAB
+ CHARACTER*2 BLAB
+
+ DATA ILAB/'C','F','G','H'/
+ DATA ICROSREF/21,17,18,52,13,22,16,45, 1, 3, 4,14, 2,42,10,25,57,58,44,37&
+ ,38,39,40,15,31,30,19,20,53,54,55,56,11,12, 9,32,33,34,26,27&
+ ,59,36,60,61,62/
+! DATA LABELF/'INBNGEO ','INBNRST ','INRM1 ','INELEV ','INHYD ','INELTFL ','OUTFIL ','OUTBNRMA','OUTBNRST','TIMFIL ','OUTSMS ','OUTSPL ','OUTCON ','INCROS ','INSRCORD','INSTRESS','INBNSTRS','INWGT ','INBNWGT ','OUTWGT ','OUTBNWGT','OUTBNSTR','WEIRDATA','AWINDIN ','BWINDIN ','OUTFLOW ','INHARM ','OUTBN3GE'/
+ DATA LABELF/&
+ 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELTFL ','INBNELF ','INHARM ','INLAYDAT','METFIL '&
+ ,'OUTFIL ','OUTBNRMA','OUTBNRST','OUTBN3GE','OUTCON ','OUTMET ','TIMFIL ','BWINDIN ','AWINDIN ','OUTBNXTR'&
+ ,'OUTASXTR','OUTBNELF','INSRCORD','INBNWGT ','INBNSTRS','OUTBNWGT','OUTWGT ','GROUPNUM','STFLFIL ','INASTRAT'&
+ ,'BCFIL ','VOLFIL ','OUTBNRES','VELBNFIL','IN3DBNGE','OUTBN2GE','INSMSGN ','OUTSMS '&
+ ,'OUTSMS1 ','OUTFLOW ','INBC ','INBNBC ','OUTBNICE'/
+
+
+! DATA LIMITL/ 'MAXFRONT','BUFFSIZL','MAXSTEPS','MAXPBUFR','MAXCROSS','MAXCPTS ','MAXQINPT','MAXHINPT','MAXEINPT','MAXQPTS ','MAXHPTS ','MAXELPTS','BUFFSIZ '/
+ DATA LIMITL/ 'MAXFRONT','BUFFSIZL','MAXLAY ','MAXCROSS','MAXCPTS ','MAXQINPT','MAXHINPT','MAXEINPT'/
+
+ DATA LABELC/ 'STARTIM ','ENDTIM ','PRTOPT ','SSITN ','TRITN ','MARSHOPT','TIMPROJ ','ITERINTV','1DTYP ','MESSRWND','TRSTEPS ','HALF-TIM','HALF-HOL','GRAVITY ','OVERLAYT','WRTINIT ','EDDYTYPE','DRYSW ','HIGH-ORD','9-NODE ','RSTOPT ','STBIN ','FRQBIN ','FRQRST ','NUMRST ','SOLVER '/
+ DATA LABELFL/'ELEV ','OMEGA ','CONV-VEL','CONV-DEP','MARSH-SH','MARSH-RG','MARSH-PR','MARSH-LM','TBFACT ','TBMIN ','HALF-TOL','X-SCALE ','Y-SCALE ','V-SCALE ','DSET ','DSETD ','UNOM ','HMIN ',' ','DPCUTOF1','DPCUTOF2'/
+ DATA LABELCE/'ENDGEO '/,LABELFE/'ENDFIL '/,LABELNDL/'ENDLIMIT'/
+DATA HELPINT/&
+ '& Starting date and time ',&
+ '& Ending date and time ',&
+ '& Ending year ',&
+ '& Ending month ',&
+ '& Ending day ',&
+ '& Ending hour ',&
+ '& Vertical turbulence option ',&
+ '& Marsh option 0=inactive -2 = input marsh parameters ',&
+ '& Eddy type 0 = fixed, 1= scaled, 2 = Smagorinsky ',&
+ '& Turbulence factor when Smagorinsky is active ',&
+ '& Turbulence minimum when Smagorinsky is active ',&
+ '& Time projection switch 1= no projection, 0 = use time derivative (least stable), 2 = project using straight line ',&
+ '& 2-d model approximation ',&
+ '& 3-d model with a fixed number of layers, number of layers ',&
+ '& 3-d model with type LD1 input that is define number of layers for all nodes ',&
+ '& 3-d model with type LD2 input ',&
+ '& 3-d model with type LD3 input ',&
+ '& Units switch 1 = metric 0 = english ',&
+ '& Force zero bed velocity 1 = at all bed levels 2 = limited by bed level input ',&
+ '& Bed level below which zero bed level in applied ',&
+ '& Ignore influence of salinity on water density ',&
+ '& Ignore influence of temperature on water density ',&
+ '& Ignore influence of sediment on water density ',&
+ '& Force zero velocity on water surface ',&
+ '& Save binary file for all iterations 0 = ignore, 1 =save ',&
+ '& Option for vertical density model 0 = standard 1 = revised method ',&
+ '& Latitude (degrees) positive in northern hemisphere ',&
+ '& Reference level for vertical transformation / Starting water level ',&
+ '& x scale factor (0.0 equivalent to 1.0) ',&
+ '& y scale factor (0.0 equivalent to 1.0) ',&
+ '& z sacle factor (0.0 equivalent to 1.0) ',&
+ '& Defines constant associated with vertical constituent distribution at boundary ',&
+ '& Defines power associated with vertical constituent distribution at boundary ',&
+ '& Nominal water velocity at startup ',&
+ '& Nominal current direction in radians c-clockwise from horzontal ',&
+ '& Minimum depth as startup (ignored if set = 0.0) Use with care can cause initial instabilty ',&
+ '& Depth at which drying (element removal) is initiated when drying ',&
+ '& Depth at which rewetting (element return) occurs when flooding ',&
+ '& Switch that forces boundary conditions to be applied regardless of direction ',&
+ '& Number of vertical velocity equations to switch to PARDISO (Default=20000) ',&
+ '& Initial salinity when no restart file specified ',&
+ '& Initial temperature when no restart file specified ',&
+ '& Initial sediment concentration when no restart file specified ',&
+ '& Bed x-velocity used to compute resistance to flow for initial conditions regardless of initial velocity ',&
+ '& Bed y-velocity used to compute resistance to flow for initial conditions regardless of initial velocity ',&
+ '& Percent of outflowing average concentration returned on next tidal inflow ',&
+ '& Depth of surface mixing (associated with free water surface local turbulence ',&
+ '& Number of steady state iteration cycles to be applied ',&
+ '& Number of time transient iteration cycles to be applied per time step ',&
+ '& Number of transient steps to be simulated ',&
+ '& Skip element and node data and initial conditions in echo print ',&
+ '& Skip element and node data show initial conditions in echoprint ',&
+ '& Show element and node data in echo print ',&
+ '& Iteration frequency for ASCII output. Skip if = 0 ',&
+ '& Time step frequency for ASCII output ',&
+ '& Save startup to binary results file ',&
+ '& Switch to initiate element elimination ',&
+ '& Time step frequency for binary output ',&
+ '& Time step frequency for individual binary restart file ',&
+ '& Time step frequency for restarting ASCII results file (saves file size for very large number of steps) ',&
+ '& 1-D cross section type (see users manual for more details) ',&
+ '& Equation solver type 0=FRONT > 0 PARDISO SOLVER ',&
+ '& Switch to initiate out of local memory (disk) version of PARDISO for very large prblems ',&
+ '& Time step frequency for restarting ASCII message file (saves file size for very large number of steps) ',&
+ '& Initial water surface elevation, Required if different from tra nsformation reference level ',&
+ '& Convergence limit for velocities ',&
+ '& Convergence limit for depth ',&
+ '& Convergence limit for salinity ',&
+ '& Convergence limit for temperature ',&
+ '& Convergence limit for sediment concentration ',&
+ '& Transition depth for collapse from 3-D to 2-D approximation ',&
+ '& Save flows in continuity line output ',&
+ '& Save depths in continuity line output ',&
+ '& Save average water surface elevations in continuity line output ',&
+ '& Save average salinity in continuity line output ',&
+ '& Save average temperature in continuity line output ',&
+ '& Save average sediment concentration in continuity line output ',&
+ '& Horizontal turbulent eddy coefficient ',&
+ '& Vertical turbulent eddy coefficient ',&
+ '& Horizontal turbulent diffusion coefficient ',&
+ '& Vertical turbulent diffusion coefficient ',&
+ '& Manning/Chezy bed friction coefficient >1.0 is Chezy coefficient ',&
+ '& Bank Manning coefficient if appropriate ',&
+ '& Water surface Manning coefficient if appropriate ',&
+ '& Marsh factor when drying applied to friction coefficient ',&
+ '& Eqn for vertical distribution factor for horizontal eddy coefficient F=a+z*(b+c*z) ',&
+ '& Boundary velocity vertical distribution constant ',&
+ '& Boundary velocity vertical distribution power ',&
+ '& 2-D to 3-D velocity vertical distribution constant ',&
+ '& 2-D to 3-D velocity vertical distribution power ',&
+ '& Marsh coefficient depth shift ',&
+ '& Marsh coefficient range ',&
+ '& Marsh coefficient porosity ',&
+ '& Marsh coefficient limit ',&
+ '& Indicator of end of tabular data set ',&
+ '& Indicator of end of end control and geometry data ',&
+ '& Overlay input time/date onto restart file ',&
+ ' '/
+ DATA HELPTIME/&
+ '& Time step in minutes ',&
+ '& Time step in hours ',&
+ '& Ending time for the time step block ',&
+ '& Iterate on flow/depth alone ',&
+ '& iterate on flow/depth and salinity -2- flow/depth and -1- salinty iterations per cycle ',&
+ '& iterate on flow/depth and temperature -2- flow/depth and -1- salinty iterations per cycle ',&
+ '& iterate on flow/depth and sediment -2- flow/depth and- 1- salinty iterations per cycle ',&
+ '& iterate on flow/depth salinity and temperature -2- flow/depth and- 1- salinty -1- temperature iterations per cycle ',&
+ '& Continuity line elevation FL-ELEVI is elevation specified as tapered between two end values (subs values from file)- CCLINE, Elev1, Elev2,bc-sal,bc-temp,bc-sed',&
+ '& Continuity line elevation CN-ELEVI is elevation specified as tapered between two end values- CCLINE, Elev1, Elev2 ,bc-sal,bc-temp,bc-sed ',&
+ '& Continuity line elevation FL-ELEVC is elevation specified as constant along a line (subs values from file)- CCLINE, Elevation ,bc-sal,bc-temp,bc-sed ',&
+ '& Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed ',&
+ '& Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file ',&
+ '& Continuity line inflow CN-QF is inflow tapered to zero at end element, flow parallel allowed - CCLINE, Total Flow, dir,bc-sal ,bc-temp,bc-sed -FL means file ',&
+ '& Continuity line inflow CN-QH is specified inflow distribution - CCLINE, Total Flow, dir,bc-sal ,bc-temp,bc-sed -FL means subsequent values from file ',&
+ '& Element Inflow treated as a rainfall rate (in mm/day) over all surface elementa. flow rate,salinity,temp,sediment -FL means subsequent values from file ',&
+ '& Element Inflow treated as a total rate over all given layer of elements. Lay num,flow rate,salinity,temp,sediment -FL means subsequent values from file ',&
+ '& Element Inflow treated as a spcific rate over all given layer of elements. Lay num,flow rate,salinity,temp,sediment -FL means subsequent values from file ',&
+ '& Element Inflow listing element number and total inflow into element. Elt num, Lay num,flow rate,salinity,temp,sediment -FL means subsequent values from file ',&
+ '& Element Inflow listing element number and spec flow into element. Elt num, Lay num,flow rate,salinity,temp,sediment -FL means subsequent values from file ',&
+ '& Element Inflow listing element inflow for a given group and rainfall into element. Grp num,flow rate,salinity,temp,sed -FL means subsequent values from file ',&
+ '& Element Inflow listing element inflow for a given type and specific flow entering element. ',&
+ '& Element Inflow listing element inflow for a given group and specific flow entering element. ',&
+ '& Element Inflow listing element number and specific flow entering element. SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow listing element inflow for a given type and specific flow entering element. SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow listing element inflow for a given type treated as a rainfall rate (in mm/day) and rainfall entering element. (SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow listing element inflow for a given group and specific flow entering element. (SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow listing rainfall entering all elements (in mm/day) . Element number ',&
+ '& Continuity line elevation FL-ELVHM is elevation specified as constant along a line (Computes values from harmonic file)- CCLINE, Elevation ',&
+ '& Element Inflow listing element inflow for a given type treated as a rainfall rate (in mm/day) and rainfall entering element. ) ',&
+ '& Element Inflow listing rainfall entering a single elements (in mm/day) . Element number = 0 apply globally ',&
+ '& Element Inflow listing element inflow treated as a rainfall rate (in mm/day) and rainfall entering all elements. (SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow listing element inflow for a given group and rainfall entering element. (SUBSQT VALUES FROM A FILE) ',&
+ '& Stage-Flow boundary condition. See line type SQC of user document ',&
+ '& Flag to identify units of wind vel = 0, then values are given in miles/hour = 1, then values are given in meters/second ',&
+ '& Wind velocity applied at all nodes. ',&
+ '& The angle between the wind velocity and the x-axis for all nodes (degrees measured anti-clockwise). ',&
+ '& Flow controller identifier, Flow controller type AJ1, BJ1, CJ1, GAM1, QD1, DJ1, TRCCE in flow controller equation for line type FC in the user document ',&
+ '& Element Inflow and total flow entering element. SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow and specific flow entering element. SUBSQT VALUES FROM A FILE) ',&
+ '& Element Inflow and total flow entering element. ',&
+ '& Element Inflow and specific flow entering element. ',&
+ '& '/
+
+ DATA HELPFILE/&
+ '& Input binary file with geometric data',&
+ '& Input ASCII file with geometric data',&
+ '& Input binary file with restart data',&
+ '& Input ASCII file with cross-section data',&
+ '& Input ASCII file with boundary flow data',&
+ '& Input ASCII file with elevation boundary data',&
+ '& Input ASCII file with element inflow data',&
+ '& Input binary element inflow data',&
+ '& Input ASCII file with elevation harmonic boundary data',&
+ '& Input ASCII file with layer data',&
+ '& Input ASCII file with meteorological data',&
+ '& Output ASCII results',&
+ '& Output binary file with new format results data',&
+ '& Output binary file with restart data',&
+ '& Output binary file with 3-D geometry data',&
+ '& Output ASCII file with continuity line flow data',&
+ '& Output ASCII file with summary meteorological input',&
+ '& Output ASCII file with cpu timing data',&
+ '& Input binary file with wind data',&
+ '& Input ASCII file with wind data',&
+ '& Output binary file listing extracted results for specified nodes',&
+ '& Output ASCII file listing extracted results for specified nodes',&
+ '& Output binary element inflow data',&
+ '& Input coordinates for surface stresses',&
+ '& Input binary weighting file for surface stresses',&
+ '& Input binary file for surface stresses',&
+ '& Output binary weighting file for surface stresses',&
+ '& Output ASCII weighting file for surface stresses',&
+ '& Input ASCII file with group number data',&
+ '& Input ASCII file with stage flow data for continuity lines',&
+ '& Input ASCII temperature stratification data file',&
+ '& Input ASCII file with continuation transient step data',&
+ '& Output ASCII file with average salinity/total flow',&
+ '& Output binary file with old format results data',&
+ '& Input binary file with velocity/depth data',&
+ '& Input binary file with 3-D geometry data',&
+ '& Output binary file with 2-D geometry data',&
+ '& Input SMS format geometry file',&
+ '& Output SMS RMA-2 format file with horizontal velocity and depth data',&
+ '& Output SMS RMA4 format file with salinity, temperate and sediment data',&
+ '& Output binary RMA-11 format file listing 1-d flows as a pseudo WQ constituent',&
+ '& Input ASCII file containing nodal boundary condition data',&
+ '& Input binary file containing nodal boundary condition data',&
+ '& Output ice results',&
+ '& '/
+
+ LINE=0
+ CALL WMessageBox(YesNo,QuestionIcon,CommonYes, &
+ 'Do you wish to save all data lines(YES) or skip default values (NO) ?', 'OUTPUT DEFAULT VALUES?')
+ IF(WInfoDialog(4) .eq. 1) then
+ ISHOWSW=1
+ ELSE
+ ISHOWSW=0
+ ENDIF
+
+ DO K=1,44
+ IF(FNAME(K)(1:16) .NE. ' ') THEN
+ LINE=LINE+1
+ DO KK=1,43
+ IF(FNAME(K)(KK:KK) .NE. ' ') THEN
+ DO KKK=1,72-KK
+ FNAME(K)(KKK:KKK)=FNAME(K)(KK-1+KKK:KK-1+KKK)
+ ENDDO
+ EXIT
+ ENDIF
+ ENDDO
+ WRITE(DATALIN(LINE),6000) LABELF(K),FNAME(K)(1:69),helpfile(k)
+6000 FORMAT(A8,3X,A69,a89)
+ ENDIF
+ ENDDO
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6000) LABELFE
+
+ DO K=1,8
+ IF(LIMITVAL(K) .NE. LIMITDFT(K)) THEN
+ LINE=LINE+1
+ IF(K .EQ. 2) THEN
+ WRITE(DATALIN(LINE),'(A8,I16)') LIMITL(K),LIMITVAL(K)
+ ELSE
+ WRITE(DATALIN(LINE),'(A8,I8)') LIMITL(K),LIMITVAL(K)
+ ENDIF
+ ENDIF
+ ENDDO
+
+
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6000) LABELNDL
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60001) TITLEHEAD
+60001 FORMAT('TI',6X,A72)
+
+ IF(NLAYTYP .NE. NLAYTYPDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ IF(NLAYTYP .EQ. 0) WRITE(DATALIN(LINE),6011) helpint(13)
+6011 FORMAT('2DMODEL ',t70,a96)
+ IF(NLAYTYP .EQ. 1) WRITE(DATALIN(LINE),6012) helpint(14)
+6012 FORMAT('3DFIXLAY',t70,a96)
+ IF(NLAYTYP .EQ. 2) WRITE(DATALIN(LINE),6013) helpint(15)
+6013 FORMAT('3DTYPLD1',t70,a96)
+ IF(NLAYTYP .EQ. 3) WRITE(DATALIN(LINE),6014) helpint(16)
+6014 FORMAT('3DTYPLD2',t70,a96)
+ IF(NLAYTYP .EQ. 4) WRITE(DATALIN(LINE),6015) helpint(17)
+6015 FORMAT('3DTYPLD3',t70,a96)
+ ENDIF
+ IF(ID1DND .NE. ID1DNDDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6056) ID1DND,helpint(61)
+6056 FORMAT('1DTYP ',I8,t70,a96)
+ ENDIF
+ IF(ELEV .NE. ELEVDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6026) ELEV,helpint(28)
+6026 FORMAT('ELEV ',F8.3,t70,a96)
+ ENDIF
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60002) LABELC(1),DATEST,helpint(1)
+60002 FORMAT(A8,3X,A16,t70,a96)
+ IF(JT .GT. 1 .OR. DELT(1) .GT. 0. .OR. DELTM(1) .GT. 0. .or. DELT(2) .GT. 0. .OR. DELTM(2) .GT. 0.) THEN
+ IF(DATEND .NE. ' ') THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60002) LABELC(2),DATEND,helpint(2)
+ ENDIF
+ ENDIF
+ IF(IOV .NE. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60003) helpint(97)
+60003 FORMAT('OVERLAYT',t70,a96)
+ ENDIF
+
+ IF(NITI .NE. NITIDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6046) NITI,helpint(48)
+6046 FORMAT('SSITN ',I8,t70,a96)
+ ENDIF
+ IF(NITN .NE. NITNDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6047) NITN,helpint(49)
+6047 FORMAT('TRITN ',I8,t70,a96)
+ ENDIF
+ IF(NCYC .NE. NCYCDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6048) NCYC,helpint(50)
+6048 FORMAT('TRSTEPS ',I8,t70,a96)
+ ENDIF
+ IF(IPROJ .NE. IPROJDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6010) IPROJ,helpint(12)
+6010 FORMAT('TIMPROJ ',I8,t70,a96)
+ ENDIF
+ IF(ICPU .NE. ICPUDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6057) ICPU,helpint(62)
+6057 FORMAT('SOLVER ',I8,t70,a96)
+ ENDIF
+ IF(IEQSWT .NE. IEQSWTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6038) IEQSWT,helpint(40)
+6038 FORMAT('VERTEQSW',I8,t70,a96)
+ ENDIF
+ IF(IOOC .NE. IOOCDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6058) IOOC,helpint(63)
+6058 FORMAT('SLVERDSK',I8,t70,a96)
+ ENDIF
+ IF(IOPTZD .NE. IOPTZDDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6005) IOPTZD,helpint(7)
+6005 FORMAT('VERTTURB',I8,t70,a96)
+ ENDIF
+ IF(IEDSW .NE. IEDSWDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6007) IEDSW,helpint(9)
+6007 FORMAT('EDDYTYPE',I8,t70,a96)
+ ENDIF
+ IF(TBFACT .NE. TBFACTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6008) TBFACT,helpint(10)
+6008 FORMAT('TBFACT ',F8.3,t70,a96)
+ ENDIF
+ IF(TBMIN .NE. TBMINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6009) TBMIN,helpint(11)
+6009 FORMAT('TBMIN ',F8.3,t70,a96)
+ ENDIF
+ IF(IVERTDEN .NE. IVERTDENDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6024) IVERTDEN,helpint(26)
+6024 FORMAT('REVRTDEN',F8.3,t70,a96)
+ ENDIF
+ IF(DMIX .NE. DMIXDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6045) DMIX,helpint(47)
+6045 FORMAT('SURFMIX ',F8.3,t70,a96)
+ ENDIF
+ IF(ITRANSIT .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60651) TRANSIT,helpint(71)
+60651 FORMAT('TRANSIT',1X,F8.5,t70,a96)
+ ENDIF
+ IF(IGRV .NE. IGRVDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6016) IGRV,helpint(18)
+6016 FORMAT('GRAVITY ',I8,t70,a96)
+ ENDIF
+ IF(OMEGA1 .NE. OMEGA1DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6025) OMEGA1,helpint(27)
+6025 FORMAT('OMEGA ',F8.3,t70,a96)
+ ENDIF
+ IF(IPASS1 .NE. IPASS1DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6018) IPASS1,helpint(21)
+6018 FORMAT('PASSSAL ',I8,t70,a96)
+ ENDIF
+ IF(IPASS2 .NE. IPASS2DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6019) IPASS2,helpint(22)
+6019 FORMAT('PASSTMP ',I8,t70,a96)
+ ENDIF
+ IF(IPASS3 .NE. IPASS3DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6020) IPASS3,helpint(23)
+6020 FORMAT('PASSSED ',I8,t70,a96)
+ ENDIF
+ IF(IZB .NE. IZBDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6017) IZB,helpint(19)
+6017 FORMAT('ZEROBED ',I8,t70,a96)
+ ENDIF
+ IF(ZERBED .NE. ZERBEDDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6022) ZERBED,helpint(20)
+6022 FORMAT('ZERBELLV',F8.3,t70,a96)
+ ENDIF
+ IF(IZERS .NE. IZERSDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6021) IZERS,helpint(24)
+6021 FORMAT('ZEROSURF',I8,t70,a96)
+ ENDIF
+ IF(ELEV1 .NE. ELEV .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6060) ELEV1,helpint(65)
+6060 FORMAT('INITWSRF',F8.3,t70,a96)
+ ENDIF
+ IF(SALI .NE. SALIDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6039) SALI,helpint(41)
+6039 FORMAT('INITSAL ',F8.3,t70,a96)
+ ENDIF
+ IF(TEMPI .NE. TEMPIDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6040) TEMPI,helpint(42)
+6040 FORMAT('INITTEMP',F8.3,t70,a96)
+ ENDIF
+ IF(SEDI .NE. SEDIDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6041) SEDI,helpint(43)
+6041 FORMAT('INITSED ',F8.3,t70,a96)
+ ENDIF
+ IF(UINP .NE. UINPDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6042) UINP,helpint(44)
+6042 FORMAT('INITUBED',F8.3,t70,a96)
+ ENDIF
+ IF(VINP .NE. VINPDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6043) VINP,helpint(45)
+6043 FORMAT('INITVBED',F8.3,t70,a96)
+ ENDIF
+ IF(PRCNT .NE. PRCNTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6044) PRCNT,helpint(46)
+6044 FORMAT('RETRNPCT',F8.3,t70,a96)
+ ENDIF
+ IF(CONV(1) .NE. CONVDF(1) .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6061) CONV(1),helpint(66)
+6061 FORMAT('CONV-VEL',F8.5,t70,a96)
+ ENDIF
+ IF(CONV(3) .NE. CONVDF(3) .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6062) CONV(3),helpint(67)
+6062 FORMAT('CONV-DEP',F8.5,t70,a96)
+ ENDIF
+ IF(CONV(4) .NE. CONVDF(4) .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6063) CONV(4),helpint(68)
+6063 FORMAT('CONV-SAL',F8.5,t70,a96)
+ ENDIF
+ IF(CONV(5) .NE. CONVDF(5) .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6064) CONV(5),helpint(69)
+6064 FORMAT('CONV-TMP',F8.5,t70,a96)
+ ENDIF
+ IF(CONV(6) .NE. CONVDF(6) .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6065) CONV(6),helpint(70)
+6065 FORMAT('CONV-SED',F8.5,t70,a96)
+ ENDIF
+ IF(IPRT .EQ. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60481) helpint(51)
+60481 FORMAT('ECHSVNON',t70,a96)
+ ELSEIF(IPRT .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60482) helpint(52)
+60482 FORMAT('ECHSVINI',t70,a96)
+ ELSEIF(ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60483) helpint(53)
+60483 FORMAT('ECHSVALL',t70,a96)
+ ENDIF
+ IF(NPRTI .NE. NPRTIDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6049) NPRTI,helpint(54)
+6049 FORMAT('ITERINTV',I8,t70,a96)
+ ENDIF
+ IF(NPRTF .NE. NPRTFDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6050) NPRTF,helpint(55)
+6050 FORMAT('ASCSVFRQ',I8,t70,a96)
+ ENDIF
+ IF(IRSAV .NE. IRSAVDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6051) IRSAV,helpint(56)
+6051 FORMAT('STBIN ',I8,t70,a96)
+ ENDIF
+ IF(IOUTFREQ .NE. IOUTFREQDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6053) IOUTFREQ,helpint(58)
+6053 FORMAT('FRQBIN ',I8,t70,a96)
+ ENDIF
+ IF(IOUTRST .NE. IOUTRSTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6054) IOUTRST,helpint(59)
+6054 FORMAT('FRQRST ',I8,t70,a96)
+ ENDIF
+ IF(IREWOUT .NE. IREWOUTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6055) IREWOUT,helpint(60)
+6055 FORMAT('RWDASFRQ',I8,t70,a96)
+ ENDIF
+ IF(IREWMES .NE. IREWMESDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6059) IREWMES,helpint(64)
+6059 FORMAT('MESSRWND',I8,t70,a96)
+ ENDIF
+ IF(ISAVITR .NE. ISAVITRDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6023) ISAVITR,helpint(25)
+6023 FORMAT('SAVITR',I8,t70,a96)
+ ENDIF
+ IF(JSPLPT .GT. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60861)
+60861 FORMAT(30X,'& Nodes for extraction of all constiuents to a time series file')
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60862) (NSPLPT(J),J=1,JSPLPT)
+60862 FORMAT('SPR-NOD ',1X,10I7)
+ ENDIF
+ IF(ICTOTC .GT. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6066)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60661) (NCSPLPT(J),J=1,ICTOTC)
+6066 FORMAT(40X,'& Continuity lines for selective output')
+60661 FORMAT('CCLINOUT',20I7)
+ ENDIF
+
+ DO J=1,6
+ IF(ICCOMP(J) .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6067) helpint(72)
+6067 FORMAT('CONT-FLW',t70,a96)
+ ELSEIF(ICCOMP(J) .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6068) helpint(73)
+6068 FORMAT('CONT-DEP',t70,a96)
+ ELSEIF(ICCOMP(J) .EQ. 3) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6069) helpint(74)
+6069 FORMAT('CONT-ELV',t70,a96)
+ ELSEIF(ICCOMP(J) .EQ. 4) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6070) helpint(75)
+6070 FORMAT('CONT-SAL',t70,a96)
+ ELSEIF(ICCOMP(J) .EQ. 5) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6071) helpint(76)
+6071 FORMAT('CONT-TMP',t70,a96)
+ ELSEIF(ICCOMP(J) .EQ. 6) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6072) helpint(77)
+6072 FORMAT('CONT-SED',t70,a96)
+ ENDIF
+ ENDDO
+ IF(XSCALE .NE. XSCALEDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6027) XSCALE,helpint(29)
+6027 FORMAT('X-SCALE ',F8.3,t70,a96)
+ ENDIF
+ IF(YSCALE .NE. YSCALEDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6028) YSCALE,helpint(30)
+6028 FORMAT('Y-SCALE ',F8.3,t70,a96)
+ ENDIF
+ IF(ZSCALE .NE. ZSCALEDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6029) ZSCALE,helpint(31)
+6029 FORMAT('Z-SCALE ',F8.3,t70,a96)
+ ENDIF
+ IF(IDSWT .NE. IDSWTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6052) IDSWT,helpint(57)
+6052 FORMAT('DRYSW ',I8,t70,a96)
+ ENDIF
+ IF(DSET .NE. DSETDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6035) DSET,helpint(37)
+6035 FORMAT('DSET ',F8.3,t70,a96)
+ ENDIF
+ IF(DSETD .NE. DSETDDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6036) DSETD,helpint(38)
+6036 FORMAT('DSETD ',F8.3,t70,a96)
+ ENDIF
+ IF(IFXBC .NE. IFXBCDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6037) IFXBC,helpint(39)
+6037 FORMAT('FIXBC ',I8,t70,a96)
+ ENDIF
+ IF(UNOM .NE. UNOMDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6032) UNOM,helpint(34)
+6032 FORMAT('UNOM ',F8.3,t70,a96)
+ ENDIF
+ IF(UDIR .NE. UDIRDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6033) UDIR,helpint(35)
+6033 FORMAT('UDIR ',F8.3,t70,a96)
+ ENDIF
+ IF(HMIN .NE. HMINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6034) HMIN,helpint(36)
+6034 FORMAT('HMIN ',F8.3,t70,a96)
+ ENDIF
+ DO K=1,NMAT
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6073) K,ORT(K,1),helpint(78)
+6073 FORMAT('EDDY-MAT',I8,F10.5,t70,a96)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6074) K,ORT(K,6),helpint(79)
+6074 FORMAT('EDDY-VRT',I8,1PE12.3,t70,a96)
+ IF(ORT(K,8) .NE. 0. .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6075) K,ORT(K,8),helpint(80)
+6075 FORMAT('DIFF-MAT',I8,F10.5,t70,a96)
+ ENDIF
+ IF(ORT(K,10) .NE. 0. .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6076) K,ORT(K,10),helpint(81)
+6076 FORMAT('DIFF-VRT',I8,1PE12.3,t70,a96)
+ ENDIF
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6077) K,ORT(K,5),helpint(82)
+6077 FORMAT('MANN-MAT',I8,F10.5,t70,a96)
+ IF(ORT(K,11) .NE. 0. .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6078) K,ORT(K,11),helpint(83)
+6078 FORMAT('MANN-BNK',I8,F10.5,t70,a96)
+ ENDIF
+ IF(ORT(K,13) .NE. 0. .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6079) K,ORT(K,13),helpint(84)
+6079 FORMAT('MANN-SRF',I8,F10.5,t70,a96)
+ ENDIF
+ IF(ORT(K,12) .NE. 0. .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6080) K,ORT(K,12),helpint(85)
+6080 FORMAT('MARSH-FR',I8,F10.5,t70,a96)
+ ENDIF
+ IF(EDD1(K) .NE. 0. .OR. EDD2(K) .NE. 0. .OR. EDD3(K) .NE. 0 .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60811) K,EDD1(K),EDD2(K),EDD3(K),helpint(86)
+60811 FORMAT('ED-VTVAR',I8,3F8.3,t70,a96)
+ ENDIF
+ ENDDO
+ IF(IDNOPT .NE. IDNOPTDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6006) IDNOPT,helpint(8)
+6006 FORMAT('MARSHOPT',I8,t70,a96)
+ ENDIF
+ IF(ACT1 .NE. ACT1DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6087) ACT1,helpint(91)
+6087 FORMAT('MARSH-SH',F8.3,t70,a96)
+ ENDIF
+ IF(ACT2 .NE. ACT2DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6088) ACT2,helpint(92)
+6088 FORMAT('MARSH-RG',F8.3,t70,a96)
+ ENDIF
+ IF(ACT3 .NE. ACT3DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6089) ACT3,helpint(93)
+6089 FORMAT('MARSH-PR',F8.3,t70,a96)
+ ENDIF
+ IF(ACT4 .NE. ACT4DF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6090) ACT4,helpint(94)
+6090 FORMAT('MARSH-LM',F8.3,t70,a96)
+ ENDIF
+ IF(L3DTYP .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60841)
+60841 FORMAT(30X,'& Nodal layer data type LD2, 0 = apply to all nodes, number of layers Proportional factors')
+ DO K=1,L3DCNT
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6085) JND(K),KKLAY(K),(L3DDAT(K,J),J=1,KKLAY(K))
+6085 FORMAT('LD2 ',2I8,20F8.3)
+ ENDDO
+ ELSEIF(L3DTYP .EQ. 3) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),60851)
+60851 FORMAT(30X,'& Nodal layer data type LD3, 0 = apply to all nodes, number of layers Layer elevations')
+ LINE=LINE+1
+ DO K=1,L3DCNT
+ WRITE(DATALIN(LINE),6086) JND(K),KKLAY(K),(L3DDAT(K,J),J=1,KKLAY(K))
+6086 FORMAT('LD3 ',2I8,20F8.3)
+ ENDDO
+ endif
+ IF(CMIN .NE. CMINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6030) CMIN,helpint(32)
+6030 FORMAT('BOUNDCMN',F8.3,t70,a96)
+ ENDIF
+ IF(CPR .NE. CPRDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6031) CPR,helpint(33)
+6031 FORMAT('BOUNDCPW',F8.3,t70,a96)
+ ENDIF
+ IF(VMIN .NE. UMINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6081) VMIN,helpint(87)
+6081 FORMAT('BN-V-MIN',F8.3,t70,a96)
+ ENDIF
+ IF(POWER .NE. POWERDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6082) POWER,helpint(88)
+6082 FORMAT('BN-V-PWR',F8.3,t70,a96)
+ ENDIF
+ IF(UMIN .NE. UMINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6083) UMIN,helpint(89)
+6083 FORMAT('TR-V-MIN',F8.3,t70,a96)
+ ENDIF
+ IF(PWERIN .NE. PWERINDF .OR. ISHOWSW .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6084) PWERIN,helpint(90)
+6084 FORMAT('TR-V-PWR',F8.3,t70,a96)
+ ENDIF
+
+
+
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6091) helpint(95)
+6091 FORMAT('TAB-END',t70,a96)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6092) helpint(96)
+6092 FORMAT('ENDGEO',t70,a96)
+
+ DO J=1,JT
+ LINE=LINE+1
+ IF(DELT(J) .EQ. 0 .AND. DELTM(J) .EQ. 0. .OR. DELT(J) .GT. 0.) THEN
+ WRITE(DATALIN(LINE),6100) DELT(J),HELPTIME(1)
+ ELSE
+ WRITE(DATALIN(LINE),6101) DELTM(J),HELPTIME(2)
+ ENDIF
+6100 FORMAT('DELTA-HR',F12.4,T70,A162)
+6101 FORMAT('DELTA-MN',F12.4,T70,A162)
+ IF(DATENDST(J) .NE. ' ' .AND. J .GT. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6102) DATENDST(J),HELPTIME(3)
+6102 FORMAT('DATE-END',4X,A16,T70,A162)
+ ENDIF
+ IF(IITYPE(J) .EQ. 1) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6103),HELPTIME(4)
+6103 FORMAT('Q-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 2) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6104),HELPTIME(5)
+6104 FORMAT('QSL-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 3) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6105),HELPTIME(6)
+6105 FORMAT('QTM-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 4) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6106),HELPTIME(7)
+6106 FORMAT('QSD-ITR',T70,A162)
+ ELSEIF(IITYPE(J) .EQ. 5) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6107),HELPTIME(8)
+6107 FORMAT('QST-ITR',T70,A162)
+ ENDIF
+ KH=NHLIN(J)
+ KQ=NQLIN(J)
+ KE=NELIN(J)
+ DO K=1,KH
+ LINE=LINE+1
+ IF(HLIN1(K,J) .GT. -9000.) THEN
+ IF(HCKLIN(K,J) .EQ. 1 ) THEN
+ WRITE(DATALIN(LINE),6108) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),HLIN1(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(9)
+6108 FORMAT('FL-ELEVI',2I8,2F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6109) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),HLIN1(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(10)
+6109 FORMAT('CN-ELEVI',2I8,2F8.4,3F8.3,T70,A162)
+ ENDIF
+ ELSE
+ IF(HCKLIN(K,J) .EQ. 1 ) THEN
+ WRITE(DATALIN(LINE),6110) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(11)
+6110 FORMAT('FL-ELEVC',2I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6111) IHLIN(K,J),LAYNUMH(K,J),HLIN(K,J),(QUALKH(K,J,KK),KK=1,3),HELPTIME(12)
+6111 FORMAT('CN-ELEVC',2I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDIF
+ ENDDO
+
+ IF(KHCN(J) .GT. 0) THEN
+ DO K=1,KHCN(J)
+ L1=MIN(IHCN(K,J),9)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=1,L1)
+61111 FORMAT('HCN',5X,9F8.3)
+ IF(IHCN(K,J) .LE. 9) CYCLE
+ L1=MIN(IHCN(K,J),18)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=10,L1)
+ IF(IHCN(K,J) .LE. 18) CYCLE
+ L1=MIN(IHCN(K,J),27)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=19,L1)
+ IF(IHCN(K,J) .LE. 27) CYCLE
+ L1=MIN(IHCN(K,J),36)
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61111) (FHCN(K,J,L),L=28,L1)
+ ENDDO
+ ENDIF
+
+ IF(WVEL(J) .NE. 0.) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61112) WVEL(J),WDIR(J)
+61112 FORMAT('GBWIND ',1X,2F8.3)
+ ENDIF
+
+ DO K=1,KQ
+ LINE=LINE+1
+ IF(IQLIN1(K,J) .EQ. 1) ALAB='C'
+ IF(IQLIN1(K,J) .EQ. 2) ALAB='F'
+ IF(IQLIN1(K,J) .EQ. 3) ALAB='H'
+ IF(QCKLIN(K,J) .EQ. 0) THEN
+ BLAB='CN'
+ ELSE
+ BLAB='FL'
+ ENDIF
+ IF(ALAB .EQ. 'C') THEN
+ IHP=13
+ ELSEIF(ALAB .EQ. 'F') THEN
+ IHP=14
+ ELSE
+ IHP=15
+ ENDIF
+ WRITE(DATALIN(LINE),6112) BLAB,ALAB,IQLIN(K,J),LAYNUMQ(K,J),QLIN(K,J),QDLIN(K,J),(QUALKQ(K,J,KK),KK=1,3),HELPTIME(IHP)
+6112 FORMAT(A2,'-Q',A1,3X,2I8,2F10.2,3F8.3,T70,A162)
+ ENDDO
+
+ DO K=1,KE
+ IF(ECKLIN(K,J) .EQ. 1) THEN
+ BLAB='FL'
+ ELSE
+ BLAB='CN'
+ ENDIF
+ LINE=LINE+1
+ IF(IETP(K,J) .EQ. 4) THEN
+ WRITE(DATALIN(LINE),6113) BLAB,ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(16)
+6113 FORMAT(A2,'-RAIN ',F8.4,3F8.3,T70,A162)
+ ELSEIF(IETP(K,J) .EQ. 1) THEN
+ IF(IELIN(K,J) .GT. 0) THEN
+ IF(ECKLIN1(K,J) .EQ. 1) THEN
+ WRITE(DATALIN(LINE),6116) BLAB,IELIN(K,J),LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(19)
+6116 FORMAT(A2,'-ELMET',2I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6117) BLAB,IELIN(K,J),LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(20)
+6117 FORMAT(A2,'-ELMED',2I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ELSE
+ IF(ECKLIN1(K,J) .EQ. 1) THEN
+ WRITE(DATALIN(LINE),6115) BLAB,LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(18)
+6115 FORMAT(A2,'-ELMGT',I8,F8.4,3F8.3,T70,A162)
+ ELSE
+ WRITE(DATALIN(LINE),6114) BLAB,LAYNUME(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(17)
+6114 FORMAT(A2,'-ELMGD',I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDIF
+ ELSEIF(IETP(K,J) .EQ. 6) THEN
+ WRITE(DATALIN(LINE),6118) BLAB,IELIN(K,J),ELIN(K,J),(QUALKE(K,J,KK),KK=1,3),HELPTIME(21)
+6118 FORMAT(A2,'-GRAIN',I8,F8.4,3F8.3,T70,A162)
+ ENDIF
+ ENDDO
+
+ IF(NSND(J) .GT. 0) THEN
+ DO K=1,NSND(J)
+ IF(NODSND(K,J) .GT. 0) THEN
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6119) NODSND(K,J),NCODE(K,J),(BCSND(K,J,L),L=1,6)
+6119 FORMAT('SN',8X,2I7,2f10.5,4F8.3)
+ ELSE
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),61191) NCODE(K,J),(BCSND(K,J,L),L=1,6)
+61191 FORMAT('SD'15X,I7,2F10.5,4F8.2)
+ ENDIF
+ ENDDO
+ ENDIF
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6200)
+6200 FORMAT('ENDSTEP')
+
+ ENDDO
+
+ LINE=LINE+1
+ WRITE(DATALIN(LINE),6999)
+6999 FORMAT('ENDDATA')
+
+ filter="r10 file *.r10|*.R10|"
+ IOUT=66
+
+ CALL WSelectFile(Filter,SaveDialog+PromptOn+AppendExt+DirChange,NAMEFL,'Save R10 File for Execution')
+
+ IF (WInfoDialog(ExitButtonCommon).EQ.CommonOpen) THEN
+ go to 500
+ else
+ RETURN
+ endif
+
+500 CONTINUE
+ OPEN(IOUT,FILE=NAMEFL,FORM='FORMATTED')
+ DO K=1,LINE
+ WRITE(IOUT,'(A230)') DATALIN(K)
+ ENDDO
+ CLOSE(IOUT)
+ RETURN
+ END
\ No newline at end of file
diff --git a/RMA10PRP/SRC/READRM10.F90 b/RMA10PRP/SRC/READRM10.F90
index f9788f7..ed5dbc7 100644
--- a/RMA10PRP/SRC/READRM10.F90
+++ b/RMA10PRP/SRC/READRM10.F90
@@ -13,14 +13,14 @@
REAL FCOEF,EDD11
REAL AA
CHARACTER(LEN=255) :: filter
- CHARACTER*8 LABELF(43),LABELC(31),LABELFL(21),LABELCE,LABELFE,ID,LABELCF,LABELCEL,LABELCD,LABELSP,LIMITL(9)
+ CHARACTER*8 LABELF(44),LABELC(31),LABELFL(21),LABELCE,LABELFE,ID,LABELCF,LABELCEL,LABELCD,LABELSP,LIMITL(9)
CHARACTER*72 DLIN
CHARACTER*72 DLININ
DATA LIMITL/ 'MAXFRONT','BUFFSIZL','MAXLAY ','MAXCROSS','MAXCPTS ','MAXQINPT','MAXHINPT','MAXEINPT','BUFFSIZ '/
! DATA LABELF/'INBNGEO ','INBNRST ','INRM1 ','INELEV ','INHYD ','INELTFL ','OUTFIL ','OUTBNRMA','OUTBNRST','TIMFIL ','OUTSMS ','OUTSPL ','OUTCON ','INCROS ','INSRCORD','INSTRESS','INBNSTRS','INWGT ','INBNWGT ','OUTWGT ','OUTBNWGT','OUTBNSTR','WEIRDATA','AWINDIN ','BWINDIN ','OUTFLOW ','INHARM ','OUTBN3GE'/
DATA LABELF/&
- 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELFL ','INBNELF ','INHARM ','METFIL '&
+ 'INBNGEO ','INRM1 ','INBNRST ','INCROS ','INHYD ','INELEV ','INELTFL ','INBNELF ','INHARM ','INLAYDAT','METFIL '&
,'OUTFIL ','OUTBNRMA','OUTBNRST','OUTBN3GE','OUTCON ','OUTMET ','TIMFIL ','BWINDIN ','AWINDIN ','OUTBNXTR'&
,'OUTASXTR','OUTBNELF','INSRCORD','INBNWGT ','INBNSTRS','OUTBNWGT','OUTWGT ','GROUPNUM','STFLFIL ','INASTRAT'&
,'BCFIL ','VOLFIL ','OUTBNRES','VELBNFIL','IN3DBNGE','OUTBN2GE','INSMSGN ','OUTSMS '&
@@ -101,7 +101,7 @@
NQLIN=0
NHLIN=0
JT=0
- DO K=1,43
+ DO K=1,44
DO L=1,72
FNAME(K)(L:L)=' '
ENDDO
@@ -152,7 +152,7 @@
close(iin)
RETURN
ENDIF
- DO L=1,28
+ DO L=1,43
IF(ID .EQ. LABELF(L)) THEN
FNAME(L)=DLININ
GO TO 200
@@ -191,6 +191,8 @@
ENDDO
130 READ(DLININ(LIMIT:LIMIT+15),5010) DATEND
+ ELSEIF(ID(1:8) .EQ. 'OVERLAYT') THEN
+ IOV=1
ELSEIF(ID(1:6) .EQ. 'YEAR ' ) THEN
READ(DLIN,*) IYRR
IYINFO=2
@@ -268,8 +270,10 @@
IF(ZSCALE .EQ. 0.) ZSCALE=1.0
ELSEIF(ID(1:8) .EQ. 'BOUNDCMN') THEN
READ(DLIN,*) CMIN
+ IF(CMIN .EQ. 0.) CMIN=1.0
ELSEIF(ID(1:8) .EQ. 'BOUNDCPW' ) THEN
READ(DLIN,*) CPR
+ IF(CPR .EQ. 0.) CPR=1.0
ELSEIF(ID(1:4) .EQ. 'UNOM' ) THEN
READ(DLIN,*) UNOM
ELSEIF(ID(1:4) .EQ. 'UDIR' ) THEN
@@ -284,6 +288,8 @@
READ(DLIN,*) IFXBC
ELSEIF(ID(1:8) .EQ. 'VERTEQSW' ) THEN
READ(DLIN,*) IEQSWT
+ IF(IEQSWT .EQ. 0) IEQSWT=20000
+
ELSEIF(ID(1:8) .EQ. 'INITSAL ' ) THEN
READ(DLIN,*) SALI
ELSEIF(ID(1:8) .EQ. 'INITTEMP' ) THEN
@@ -319,14 +325,18 @@
READ(DLIN,*) NPRTF
ELSEIF(ID(1:5) .EQ. 'STBIN' ) THEN
READ(DLIN,*) IRSAV
+ IF(IRSAV .EQ. 0) IRSAV=1
ELSEIF(ID(1:5) .EQ. 'DRYSW' ) THEN
READ(DLIN,*) IDSWT
ELSEIF(ID(1:6) .EQ. 'FRQBIN' ) THEN
READ(DLIN,*) IOUTFREQ
+ IF(IOUTFREQ .EQ. 0) IOUTFREQ=1
ELSEIF(ID(1:6) .EQ. 'FRQRST' ) THEN
READ(DLIN,*) IOUTRST
+ IF(IOUTRST .EQ. 0) IOUTRST=10
ELSEIF(ID(1:8) .EQ. 'RWDASFRQ' ) THEN
READ(DLIN,*) IREWOUT
+ IF(IREWOUT .EQ. 0) IREWOUT=999999
ELSEIF(ID(1:5) .EQ. '1DTYP' ) THEN
READ(DLIN,*) ID1DND
ELSEIF(ID(1:6) .EQ. 'SOLVER') THEN
@@ -335,6 +345,7 @@
READ(DLIN,*) IOOC
ELSEIF(ID(1:8) .EQ. 'MESSRWND' ) THEN
READ(DLIN,*) IREWMES
+ IF(IREWMES .EQ. 0) IREWMES=999999
ELSEIF(ID(1:8) .EQ. 'INITWSRF' ) THEN
IINITWS=1
READ(DLIN,*) ELEV1
@@ -437,6 +448,10 @@
L3DCNT=L3DCNT+1
L3DTYP=2
READ(DLIN,*) JND(L3DCNT),KKLAY(L3DCNT),(L3DDAT(L3DCNT,K),K=1,KKLAY(L3DCNT))
+ ELSEIF(ID(1:8) .EQ. 'SPR-NOD ' ) THEN
+ READ(DLIN,*) (NSPLPT(J),J=JSPLPT+1,NCHARBL+JSPLPT)
+ JSPLPT=JSPLPT+NCHARBL
+
! ENDIF
diff --git a/RMA10PRP/SRC/RESOURCE10.res b/RMA10PRP/SRC/RESOURCE10.res
index 931bcb5..a46b840 100644
Binary files a/RMA10PRP/SRC/RESOURCE10.res and b/RMA10PRP/SRC/RESOURCE10.res differ
diff --git a/RMA10PRP/SRC/SHOWCNTRL.F90 b/RMA10PRP/SRC/SHOWCNTRL.F90
index 055b231..10a84c7 100644
--- a/RMA10PRP/SRC/SHOWCNTRL.F90
+++ b/RMA10PRP/SRC/SHOWCNTRL.F90
@@ -7,7 +7,7 @@
USE BLKRM10
- integer ilb(6)
+ integer ilb(6)
IF(ITEM .NE. 1) GO TO 400
300 CONTINUE
305 CONTINUE
@@ -100,7 +100,7 @@
CALL HELP10(IHSW)
GO TO 305
ELSEIF (WInfoDialog(ExitButton) .EQ. IDOK .OR. WInfoDialog(ExitButton) .EQ. IDF_PART2 .OR. WInfoDialog(ExitButton) .EQ. IDF_PART3 .OR. WInfoDialog(ExitButton) .EQ. IDF_PART4) THEN
- ITYP=0
+! ITYP=0
! ORT=0
CALL WGridGetCellInteger(IDF_GRID5,1,1,IOPTZD)
@@ -168,7 +168,9 @@
IF(NLAYTYP .GT. 0) THEN
CALL WGridGetCellReal(IDF_GRID8,1,1,DMIX)
ENDIF
-
+ IF(NLAYTYP .EQ. 4) L3DTYP=3
+ IF(NLAYTYP .EQ. 3) L3DTYP=2
+ IF(NLAYTYP .EQ. 2) L3DTYP=1
! CALL WDialogGetInteger(IDF_INTEGER14,INTG(26))
@@ -200,9 +202,10 @@
NN=0
do k=1,1000
- IF(ORT(K,1) .NE. 0.) THEN
+ IF(ITYP(K) .GT. 0) THEN
+! IF(ORT(K,1) .NE. 0. .AND. ORT(K,1) .GT. -900.) THEN
NN=NN+1
- CALL WGridPutCellInteger(IDF_GRID1,1,NN,K)
+ CALL WGridPutCellInteger(IDF_GRID1,1,NN,ITYP(K))
CALL WGridPutCellReal(IDF_GRID1,2,NN,ORT(K,1))
CALL WGridPutCellReal(IDF_GRID1,3,NN,ORT(K,6))
CALL WGridPutCellReal(IDF_GRID1,4,NN,ORT(K,8))
@@ -234,26 +237,26 @@
NMAT=0
NN=0
DO K=1,1000
- CALL WGridGetCellInteger(IDF_GRID1,1,K,KK)
-
- IF(KK .LE. 0) THEN
- EXIT
+ CALL WGridGetCellInteger(IDF_GRID1,1,K,ITYP(K))
+
+ IF(ITYP(K) .LE. 0) THEN
+ CYCLE
ELSE
- IF(KK .GT. NMAT) NMAT=KK
+ IF(ITYP(K) .GT. NMAT) NMAT=ITYP(K)
ENDIF
- CALL WGridGetCellReal(IDF_GRID1,2,K,ORT(KK,1))
- CALL WGridGetCellReal(IDF_GRID1,3,K,ORT(KK,6))
- CALL WGridGetCellReal(IDF_GRID1,4,K,ORT(KK,8))
- CALL WGridGetCellReal(IDF_GRID1,5,K,ORT(KK,10))
- CALL WGridGetCellReal(IDF_GRID1,6,K,ORT(KK,5))
- CALL WGridGetCellReal(IDF_GRID1,7,K,ORT(KK,12))
- CALL WGridGetCellReal(IDF_GRID1,8,K,ORT(KK,11))
- CALL WGridGetCellReal(IDF_GRID1,9,K,ORT(KK,13))
+ CALL WGridGetCellReal(IDF_GRID1,2,K,ORT(K,1))
+ CALL WGridGetCellReal(IDF_GRID1,3,K,ORT(K,6))
+ CALL WGridGetCellReal(IDF_GRID1,4,K,ORT(K,8))
+ CALL WGridGetCellReal(IDF_GRID1,5,K,ORT(K,10))
+ CALL WGridGetCellReal(IDF_GRID1,6,K,ORT(K,5))
+ CALL WGridGetCellReal(IDF_GRID1,7,K,ORT(K,12))
+ CALL WGridGetCellReal(IDF_GRID1,8,K,ORT(K,11))
+ CALL WGridGetCellReal(IDF_GRID1,9,K,ORT(K,13))
ENDDO
DO K=1,1000
CALL WGridGetCellInteger(IDF_GRID6,1,K,JND(K))
- IF(JND(K) .LE. 0) THEN
+ IF(JND(K) .LT. 0) THEN
L3DCNT=K-1
EXIT
ENDIF
@@ -314,10 +317,19 @@
CALL WGridPutCellInteger(IDF_GRID6,1,11,ISAVITR)
CALL WGridPutCellInteger(IDF_GRID6,1,12,IVERTDEN)
CALL WGridPutCellInteger(IDF_GRID6,1,13,IFXBC)
+
IF(IEQSWT .EQ. 0) IEQSWT=20000
CALL WGridPutCellInteger(IDF_GRID6,1,14,IEQSWT)
-
-
+
+ IF(JSPLPT .GT. 0) THEN
+ DO KK=1,5
+ IF(JSPLPT .LT. KK) EXIT
+ CALL WGridPutCellInteger(IDF_GRID4,1,KK,NSPLPT(KK))
+ IF(JSPLPT .GT. 5 .AND. KK+5 .LE. JSPLPT) THEN
+ CALL WGridPutCellInteger(IDF_GRID4,2,KK,NSPLPT(KK+5))
+ ENDIF
+ ENDDO
+ ENDIF
call wdialogPutcheckbox(IDF_CHECK1,IGRV)
call wdialogPutcheckbox(IDF_CHECK2,IOV)
call wdialogPutcheckbox(IDF_CHECK3,INTG(16))
@@ -352,12 +364,7 @@
DO KK=1,ICTOTC
CALL WGridPutCellInteger(IDF_CCLINE,1,KK,NCSPLPT(KK))
- ENDDO
-
- DO KK=1,20
- CALL WGridPutCellInteger(IDF_GRID4,1,KK,ISPRTN(KK))
- ENDDO
-
+ ENDDO
CALL WDialogShow(-1,-1,0,Modal)
@@ -404,6 +411,21 @@
CALL WGridGetCellInteger(IDF_GRID6,1,14,IEQSWT)
IF(IEQSWT .EQ. 0) IEQSWT=20000
+ DO KK=1,10
+ LL=1
+ KL=KK
+ IF(KK .GT. 5) THEN
+ KL=KK-5
+ LL=2
+ ENDIF
+ CALL WGridGetCellInteger(IDF_GRID4,LL,KL,NSPLPT(KK))
+ IF(NSPLPT(KK) .LE. 0) THEN
+ JSPLPT=KK-1
+ EXIT
+ ENDIF
+ ENDDO
+
+
call wdialogGetcheckbox(IDF_CHECK1,IGRV)
call wdialogGetcheckbox(IDF_CHECK2,IOV)
call wdialogGetcheckbox(IDF_CHECK3,INTG(16))
@@ -461,13 +483,6 @@
ENDIF
ENDDO
- DO KK=1,20
- CALL WGridGetCellInteger(IDF_GRID4,1,KK,ISPRTN(KK))
- IF(ISPRTN(KK) .EQ. 0) THEN
- NSPRTN=KK-1
- EXIT
- ENDIF
- ENDDO
IF (WInfoDialog(ExitButton) .EQ. IDF_PART1) GO TO 305
IF (WInfoDialog(ExitButton) .EQ. IDF_PART2) GO TO 405
IF (WInfoDialog(ExitButton) .EQ. IDF_PART4) GO TO 705
@@ -567,11 +582,11 @@
ENDDO
FNAME(1)='MESH.GEO'
- FNAME(11)='RESULTS.OUT'
- FNAME(12)='BRESULTS.RMA'
- DATEST=' '
+ FNAME(12)='RESULTS.OUT'
+ FNAME(13)='BRESULTS.RMA'
+ DATEST='01/01/2019 00.00'
DATEND=' '
- DATENDST='31/01/2018 12.00'
+ DATENDST='31/01/2019 12.00'
! MFWW=1000
LIMITVAL(1)=1000
@@ -763,17 +778,18 @@
-
- DO K=1,1000
- ITYP(K)=0
- ORT(K,1)=-0.2
- ORT(K,5)=0.02
- ORT(K,7)=10.
- ORTDFLT(K,1)=-0.2
- ORTDFLT(K,5)=0.02
- ORTDFLT(K,7)=10.
- ENDDO
- ITYP(1)=1
+ ITYP=0
+ ORT=0.
+ !DO K=1,1000
+ ! ITYP(K)=0
+ ! ORT(K,1)=-0.2
+ ! ORT(K,5)=0.02
+ ! ORT(K,7)=10.
+ ! ORTDFLT(K,1)=-0.2
+ ! ORTDFLT(K,5)=0.02
+ ! ORTDFLT(K,7)=10.
+ !ENDDO
+ !ITYP(1)=1
DELTM=0.
DELTH=0.
@@ -847,7 +863,6 @@
ICFLOWS=0
ICELEV=0
ICDEP=0
- ISPRTN=0
WDT=0.
IWTYP=0
diff --git a/RMA10PRP/SRC/SHOWFILES.f90 b/RMA10PRP/SRC/SHOWFILES.f90
index af2b9e7..a9c1fd8 100644
--- a/RMA10PRP/SRC/SHOWFILES.f90
+++ b/RMA10PRP/SRC/SHOWFILES.f90
@@ -24,41 +24,42 @@
CALL WGridPutCellString(IDF_GRID1,1, 7,'Input ASCII file with element inflow data')
CALL WGridPutCellString(IDF_GRID1,1, 8,'Input binary element inflow data')
CALL WGridPutCellString(IDF_GRID1,1, 9,'Input ASCII file with elevation harmonic boundary data')
- CALL WGridPutCellString(IDF_GRID1,1,10,'Input ASCII file with meteorological data')
- CALL WGridPutCellString(IDF_GRID1,1,11,'Output ASCII results')
- CALL WGridPutCellString(IDF_GRID1,1,12,'Output binary file with new format results data')
- CALL WGridPutCellString(IDF_GRID1,1,13,'Output binary file with restart data')
- CALL WGridPutCellString(IDF_GRID1,1,14,'Output binary file with 3-D geometry data')
- CALL WGridPutCellString(IDF_GRID1,1,15,'Output ASCII file with continuity line flow data')
- CALL WGridPutCellString(IDF_GRID1,1,16,'Output ASCII file with summary meteorological input')
- CALL WGridPutCellString(IDF_GRID1,1,17,'Output ASCII file with cpu timing data')
- CALL WGridPutCellString(IDF_GRID1,1,18,'Input binary file with wind data')
- CALL WGridPutCellString(IDF_GRID1,1,19,'Input ASCII file with wind data')
- CALL WGridPutCellString(IDF_GRID1,1,20,'Output binary file listing extracted results for specified nodes')
- CALL WGridPutCellString(IDF_GRID1,1,21,'Output ASCII file listing extracted results for specified nodes')
- CALL WGridPutCellString(IDF_GRID1,1,22,'Output binary element inflow data')
- CALL WGridPutCellString(IDF_GRID1,1,23,'Input coordinates for surface stresses')
- CALL WGridPutCellString(IDF_GRID1,1,24,'Input binary weighting file for surface stresses')
- CALL WGridPutCellString(IDF_GRID1,1,25,'Input binary file for surface stresses')
- CALL WGridPutCellString(IDF_GRID1,1,26,'Output binary weighting file for surface stresses')
- CALL WGridPutCellString(IDF_GRID1,1,27,'Output ASCII weighting file for surface stresses')
- CALL WGridPutCellString(IDF_GRID1,1,28,'Input ASCII file with group number data')
- CALL WGridPutCellString(IDF_GRID1,1,29,'Input ASCII file with stage flow data for continuity lines')
- CALL WGridPutCellString(IDF_GRID1,1,30,'Input ASCII temperature stratification data file')
- CALL WGridPutCellString(IDF_GRID1,1,31,'Input ASCII file with continuation transient step data')
- CALL WGridPutCellString(IDF_GRID1,1,32,'Output ASCII file with average salinity/total flow')
- CALL WGridPutCellString(IDF_GRID1,1,33,'Output binary file with old format results data')
- CALL WGridPutCellString(IDF_GRID1,1,34,'Input binary file with velocity/depth data')
- CALL WGridPutCellString(IDF_GRID1,1,35,'Input binary file with 3-D geometry data')
- CALL WGridPutCellString(IDF_GRID1,1,36,'Output binary file with 2-D geometry data')
- CALL WGridPutCellString(IDF_GRID1,1,37,'Input SMS format geometry file')
- CALL WGridPutCellString(IDF_GRID1,1,38,'Output SMS RMA-2 format file with horizontal velocity and depth data')
- CALL WGridPutCellString(IDF_GRID1,1,39,'Output SMS RMA4 format file with salinity, temperate and sediment data')
- CALL WGridPutCellString(IDF_GRID1,1,40,'Output binary RMA-11 format file listing 1-d flows as a pseudo WQ constituent')
- CALL WGridPutCellString(IDF_GRID1,1,41,'Input ASCII file containing nodal boundary condition data')
- CALL WGridPutCellString(IDF_GRID1,1,42,'Input binary file containing nodal boundary condition data')
- CALL WGridPutCellString(IDF_GRID1,1,43,'Output ice results')
- DO K=1,43
+ CALL WGridPutCellString(IDF_GRID1,1,10,'Input ASCII file with layer data')
+ CALL WGridPutCellString(IDF_GRID1,1,11,'Input ASCII file with meteorological data')
+ CALL WGridPutCellString(IDF_GRID1,1,12,'Output ASCII results')
+ CALL WGridPutCellString(IDF_GRID1,1,13,'Output binary file with new format results data')
+ CALL WGridPutCellString(IDF_GRID1,1,14,'Output binary file with restart data')
+ CALL WGridPutCellString(IDF_GRID1,1,15,'Output binary file with 3-D geometry data')
+ CALL WGridPutCellString(IDF_GRID1,1,16,'Output ASCII file with continuity line flow data')
+ CALL WGridPutCellString(IDF_GRID1,1,17,'Output ASCII file with summary meteorological input')
+ CALL WGridPutCellString(IDF_GRID1,1,18,'Output ASCII file with cpu timing data')
+ CALL WGridPutCellString(IDF_GRID1,1,19,'Input binary file with wind data')
+ CALL WGridPutCellString(IDF_GRID1,1,20,'Input ASCII file with wind data')
+ CALL WGridPutCellString(IDF_GRID1,1,21,'Output binary file listing extracted results for specified nodes')
+ CALL WGridPutCellString(IDF_GRID1,1,22,'Output ASCII file listing extracted results for specified nodes')
+ CALL WGridPutCellString(IDF_GRID1,1,23,'Output binary element inflow data')
+ CALL WGridPutCellString(IDF_GRID1,1,24,'Input coordinates for surface stresses')
+ CALL WGridPutCellString(IDF_GRID1,1,25,'Input binary weighting file for surface stresses')
+ CALL WGridPutCellString(IDF_GRID1,1,26,'Input binary file for surface stresses')
+ CALL WGridPutCellString(IDF_GRID1,1,27,'Output binary weighting file for surface stresses')
+ CALL WGridPutCellString(IDF_GRID1,1,28,'Output ASCII weighting file for surface stresses')
+ CALL WGridPutCellString(IDF_GRID1,1,29,'Input ASCII file with group number data')
+ CALL WGridPutCellString(IDF_GRID1,1,30,'Input ASCII file with stage flow data for continuity lines')
+ CALL WGridPutCellString(IDF_GRID1,1,31,'Input ASCII temperature stratification data file')
+ CALL WGridPutCellString(IDF_GRID1,1,32,'Input ASCII file with continuation transient step data')
+ CALL WGridPutCellString(IDF_GRID1,1,33,'Output ASCII file with average salinity/total flow')
+ CALL WGridPutCellString(IDF_GRID1,1,34,'Output binary file with old format results data')
+ CALL WGridPutCellString(IDF_GRID1,1,35,'Input binary file with velocity/depth data')
+ CALL WGridPutCellString(IDF_GRID1,1,36,'Input binary file with 3-D geometry data')
+ CALL WGridPutCellString(IDF_GRID1,1,37,'Output binary file with 2-D geometry data')
+ CALL WGridPutCellString(IDF_GRID1,1,38,'Input SMS format geometry file')
+ CALL WGridPutCellString(IDF_GRID1,1,39,'Output SMS RMA-2 format file with horizontal velocity and depth data')
+ CALL WGridPutCellString(IDF_GRID1,1,40,'Output SMS RMA4 format file with salinity, temperate and sediment data')
+ CALL WGridPutCellString(IDF_GRID1,1,41,'Output binary RMA-11 format file listing 1-d flows as a pseudo WQ constituent')
+ CALL WGridPutCellString(IDF_GRID1,1,42,'Input ASCII file containing nodal boundary condition data')
+ CALL WGridPutCellString(IDF_GRID1,1,43,'Input binary file containing nodal boundary condition data')
+ CALL WGridPutCellString(IDF_GRID1,1,44,'Output ice results')
+ DO K=1,44
CALL WGridPutCellString(IDF_GRID1,2,K,fname(K))
ENDDO
@@ -88,7 +89,7 @@
CALL WDialogGetString(idf_string1,TITLEHEAD)
- DO K=1,43
+ DO K=1,44
CALL WGridGetCellString(IDF_GRID1,2,K,fname(K))
ENDDO
diff --git a/RMA10PRP/SRC/resid.f90 b/RMA10PRP/SRC/resid.f90
index 09dd57d..59fadfd 100644
--- a/RMA10PRP/SRC/resid.f90
+++ b/RMA10PRP/SRC/resid.f90
@@ -1,4 +1,4 @@
-! Winteracter resource identifiers. Created : 10/Aug/2019 16:08:19
+! Winteracter resource identifiers. Created : 17/Aug/2019 17:31:52
!
! This file is generated by the Winteracter resource editor.
! It should not be edited manually. It is also not advisable to load this
@@ -213,4 +213,10 @@
INTEGER, PARAMETER :: ISS34 = 1102
INTEGER, PARAMETER :: ISS35 = 1103
INTEGER, PARAMETER :: IDD_TEMPLATE028 = 107
+ INTEGER, PARAMETER :: IDF_CHECK16 = 1064
+ INTEGER, PARAMETER :: IDF_CHECK17 = 1065
+ INTEGER, PARAMETER :: IDF_CHECK18 = 1066
+ INTEGER, PARAMETER :: IDF_CHECK19 = 1067
+ INTEGER, PARAMETER :: IDF_CHECK20 = 1069
+ INTEGER, PARAMETER :: IDF_CHECK21 = 1070
END MODULE RESID
diff --git a/RMA10PRP/SRC/resource10.rc b/RMA10PRP/SRC/resource10.rc
index 95b8ee3..a552b4f 100644
--- a/RMA10PRP/SRC/resource10.rc
+++ b/RMA10PRP/SRC/resource10.rc
@@ -8,7 +8,7 @@
//
// Winteracter resource script.
//
-// Modified : 10/Aug/2019 16:08:19
+// Modified : 17/Aug/2019 17:31:52
//
///////////////////////////////////////////////////
//
@@ -221,6 +221,12 @@
#define ISS34 1102
#define ISS35 1103
#define IDD_TEMPLATE028 107
+#define IDF_CHECK16 1064
+#define IDF_CHECK17 1065
+#define IDF_CHECK18 1066
+#define IDF_CHECK19 1067
+#define IDF_CHECK20 1069
+#define IDF_CHECK21 1070
///////////////////////////////////////////////////
//
@@ -231,10 +237,10 @@ STYLE WS_POPUP | WS_BORDER | WS_DLGFRAME | WS_SYSMENU | DS_3DLOOK | DS_MODALFRAM
FONT 8, "MS Sans Serif"
CAPTION "DEFINE FILES and LIMITS"
BEGIN
- CONTROL "",IDF_GRID1,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_VSCROLL | WS_HSCROLL | WS_GROUP | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS | GS_WANTRETURN, 8, 5, 528, 208
+ CONTROL "",IDF_GRID1,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_VSCROLL | WS_HSCROLL | WS_GROUP | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS | GS_WANTRETURN, 8, 5, 536, 208
CONTROL "OK",IDOK,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_DEFPUSHBUTTON | BS_TEXT, 345, 269, 40, 14
CONTROL "CANCEL",IDCANCEL,"BUTTON",WS_CHILD | WS_VISIBLE | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 498, 268, 40, 14
- CONTROL "",IDF_GRID2,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS, 550, 42, 172, 122
+ CONTROL "",IDF_GRID2,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS, 553, 42, 172, 122
CONTROL "HELP",IDF_HELP,"BUTTON",WS_CHILD | WS_VISIBLE | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 186, 268, 40, 14
CONTROL "",IDF_STRING1,"EDIT",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | ES_LEFT, 136, 228, 470, 14
CONTROL "TITLE",IDF_STRING2,"EDIT",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | ES_LEFT, 65, 227, 40, 14
@@ -450,29 +456,35 @@ BEGIN
CONTROL "OK",IDOK,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_DEFPUSHBUTTON | BS_TEXT, 249, 280, 40, 14
CONTROL "Cancel",IDCANCEL,"BUTTON",WS_CHILD | WS_VISIBLE | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 140, 280, 40, 14
CONTROL "",IDF_CCLINE,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_VSCROLL | WS_GROUP | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS | GS_WANTRETURN | GS_WRAP, 386, 100, 108, 69
- CONTROL "",IDF_GRID4,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_VSCROLL | WS_HSCROLL | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS | GS_WANTRETURN | GS_WRAP, 28, 196, 92, 29
+ CONTROL "",IDF_GRID4,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | GS_DEFROWLABELS | GS_COLUMNLABELS | GS_WANTRETURN | GS_WRAP, 358, 20, 82, 68
CONTROL "Metric Units",IDF_CHECK1,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 198, 169, 112, 14
CONTROL "Ovelay Time for Restart",IDF_CHECK2,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 198, 187, 112, 14
CONTROL "Write Initial Conditions",IDF_CHECK3,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 198, 206, 112, 14
CONTROL "Page 4",IDF_PART4,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 360, 260, 40, 14
CONTROL "Page 1",IDF_PART1,"BUTTON",WS_CHILD | WS_VISIBLE | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 140, 260, 40, 14
- CONTROL "Output Continuity Line Flows",IDF_CHECK4,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 170, 112, 14
- CONTROL "Output Ave Continuity Line Elevations",IDF_CHECK5,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 182, 112, 14
- CONTROL "Output Ave Continuity Line Depths",IDF_CHECK6,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 194, 112, 14
+ CONTROL "Output Continuity Line Flows",IDF_CHECK4,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 178, 112, 14
+ CONTROL "Output Ave Continuity Line Elevations",IDF_CHECK5,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 190, 112, 14
+ CONTROL "Output Ave Continuity Line Depths",IDF_CHECK6,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 202, 112, 14
CONTROL "HELP",IDF_HELP10,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 360, 280, 40, 14
CONTROL "",IDF_GRID3,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | GS_ROWLABELS | GS_WANTRETURN, 21, 12, 134, 69
CONTROL "",IDF_GRID6,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | GS_ROWLABELS | GS_WANTRETURN, 194, 8, 124, 160
- CONTROL "Zero Velocity on Bed",IDF_CHECK7,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 383, 13, 112, 14
- CONTROL "Passive Saliniity Density",IDF_CHECK8,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 43, 112, 14
- CONTROL "Passive Temperature Density",IDF_CHECK9,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 57, 112, 14
- CONTROL "Passive Sediment Density",IDF_CHECK10,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 71, 112, 14
- CONTROL "Zero Velocity at Water Surface",IDF_CHECK11,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 383, 25, 112, 14
- CONTROL "",IDF_GRID7,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | GS_ROWLABELS, 22, 101, 136, 78
+ CONTROL "Zero Velocity on Bed",IDF_CHECK7,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 24, 184, 112, 14
+ CONTROL "Passive Saliniity Density",IDF_CHECK8,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 25, 214, 112, 14
+ CONTROL "Passive Temperature Density",IDF_CHECK9,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 25, 228, 112, 14
+ CONTROL "Passive Sediment Density",IDF_CHECK10,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 25, 242, 112, 14
+ CONTROL "Zero Velocity at Water Surface",IDF_CHECK11,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 24, 196, 112, 14
+ CONTROL "",IDF_GRID7,"ISSGRID",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_GROUP | WS_TABSTOP | GS_ROWLABELS, 19, 93, 136, 78
CONTROL "Use Out of Core Pardiso (very large problems)",IDF_CHECK12,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 188, 224, 140, 14
- CONTROL "Output Ave Continuity Line Salinity",IDF_CHECK13,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 206, 112, 14
- CONTROL "Output Ave Continuity Line Temps",IDF_CHECK14,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 218, 112, 14
- CONTROL "Output Ave Continuity Line Sedment",IDF_CHECK15,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 231, 112, 14
+ CONTROL "Output Ave Continuity Line Salinity",IDF_CHECK13,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 214, 112, 14
+ CONTROL "Output Ave Continuity Line Temps",IDF_CHECK14,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 226, 112, 14
+ CONTROL "Output Ave Continuity Line Sedment",IDF_CHECK15,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 384, 239, 112, 14
CONTROL "Page 2",IDF_PART2,"BUTTON",WS_CHILD | WS_VISIBLE | WS_TABSTOP | BS_PUSHBUTTON | BS_TEXT, 249, 260, 40, 14
+ CONTROL "Extract Nodal Velocities",IDF_CHECK16,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 16, 80, 14
+ CONTROL "Extract Nodal Elevations",IDF_CHECK17,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 28, 80, 14
+ CONTROL "Extract Nodal Depths",IDF_CHECK18,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 40, 80, 14
+ CONTROL "Extract Nodal Salinity",IDF_CHECK19,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 52, 80, 14
+ CONTROL "Extract Nodal Temps",IDF_CHECK20,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 64, 80, 14
+ CONTROL "Extract Nodal Sedment",IDF_CHECK21,"BUTTON",WS_CHILD | WS_VISIBLE | WS_GROUP | WS_TABSTOP | BS_AUTOCHECKBOX | BS_LEFTTEXT | BS_TEXT, 450, 77, 80, 14
END
IDD_CONTROLS2 RCDATA
@@ -493,6 +505,12 @@ BEGIN
" 1056 0 \n"
" 1057 0 \n"
" 1063 0 \n"
+" 1064 0 \n"
+" 1065 0 \n"
+" 1066 0 \n"
+" 1067 0 \n"
+" 1069 0 \n"
+" 1070 0 \n"
"[Fonts] \n"
" 1 9 Arial, 700 0 \n"
" 2 9 Arial, 700 0 \n"
@@ -517,9 +535,15 @@ BEGIN
" 1057 9 Arial, 700 0 \n"
" 1063 9 Arial, 700 0 \n"
" 1041 10 MS Sans Serif, 700 0 \n"
+" 1064 9 Arial, 700 0 \n"
+" 1065 9 Arial, 700 0 \n"
+" 1066 9 Arial, 700 0 \n"
+" 1067 9 Arial, 700 0 \n"
+" 1069 9 Arial, 700 0 \n"
+" 1070 9 Arial, 700 0 \n"
"[Grids] \n"
" 1045 1 30 111 \n"
-" 1035 1 15 115 \n"
+" 1035 2 5 115 \n"
" 1034 1 6 103 \n"
" 100 \n"
"Turbulence Factor \n"
@@ -573,7 +597,8 @@ IDD_TEMPLATE010 DIALOG 0, 0, 1000, 16
STYLE DS_3DLOOK
FONT 9, "Arial"
BEGIN
- CONTROL "SPECIAL PRT NODES",ISS1,"INTEGEREDIT",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | ES_LEFT | ES_CENTER | ES_MULTILINE, 0, 0, 90, 14
+ CONTROL "EXTRACTION",ISS1,"INTEGEREDIT",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | ES_LEFT | ES_CENTER | ES_MULTILINE, 0, 0, 50, 14
+ CONTROL "NODES",ISS2,"INTEGEREDIT",WS_CHILD | WS_VISIBLE | WS_BORDER | WS_TABSTOP | ES_LEFT | ES_CENTER | ES_MULTILINE, 0, 0, 40, 14
END
IDD_TEMPLATE010 RCDATA
@@ -1548,4 +1573,4 @@ END
//*WI* FILENAME resid.f90
//*WI* FMODNAME RESID
//*WI* LASTTYPE 2
-//*WI* LASTRES 1
+//*WI* LASTRES 7
diff --git a/RMA10PRP/SUMMYSH.R10 b/RMA10PRP/SUMMYSH.R10
new file mode 100644
index 0000000..f895158
--- /dev/null
+++ b/RMA10PRP/SUMMYSH.R10
@@ -0,0 +1,52 @@
+INBNGEO MESH.GEO & Input binary file with geometric data
+OUTFIL RESULTS.OUT & Output ASCII results
+OUTBNRMA BRESULTS.RMA & Output binary file with new format results data
+OUTBNRST BRESULTS.RST & Output binary file with restart data
+OUTCON CONTIN.DAT & Output ASCII file with continuity line flow data
+OUTASXTR EXTRACT.DAT & Output ASCII file listing extracted results for specified nodes
+ENDFIL
+ENDLIMIT
+TI TEST ESTUARY
+3DTYPLD2 & 3-d model with type LD2 input
+STARTIM 01/01/2019 00.00 & Starting date and time
+OVERLAYT & Overlay input time/date onto restart file
+SSITN 5 & Number of steady state iteration cycles to be applied
+TRITN 5 & Number of time transient iteration cycles to be applied per time step
+TRSTEPS 4 & Number of transient steps to be simulated
+SURFMIX 0.400 & Depth of surface mixing (associated with free water surface local turbulence
+TRANSIT 0.50000 & Transition depth for collapse from 3-D to 2-D approximation
+ & Nodes for extraction of all constiuents to a time series file
+SPR-NOD 23 45 56
+ & Continuity lines for selective output
+CCLINOUT 1
+CONT-FLW & Save flows in continuity line output
+CONT-DEP & Save depths in continuity line output
+CONT-ELV & Save average water surface elevations in continuity line output
+CONT-SAL & Save average salinity in continuity line output
+CONT-TMP & Save average temperature in continuity line output
+CONT-SED & Save average sediment concentration in continuity line output
+EDDY-MAT 1 0.20000 & Horizontal turbulent eddy coefficient
+EDDY-VRT 1 1.000E-01 & Vertical turbulent eddy coefficient
+DIFF-MAT 1 0.10000 & Horizontal turbulent diffusion coefficient
+DIFF-VRT 1 1.000E-03 & Vertical turbulent diffusion coefficient
+MANN-MAT 1 0.02000 & Manning/Chezy bed friction coefficient >1.0 is Chezy coefficient
+MARSH-FR 1 10.00000 & Marsh factor when drying applied to friction coefficient
+ED-VTVAR 1 1.000 0.000 0.000 & Eqn for vertical distribution factor for horizontal eddy coefficient F=a+z*(b+c*z)
+MARSHOPT -2 & Marsh option 0=inactive -2 = input marsh parameters
+MARSH-SH 2.000 & Marsh coefficient depth shift
+MARSH-RG 0.500 & Marsh coefficient range
+MARSH-PR 0.010 & Marsh coefficient porosity
+ & Nodal layer data type LD2, 0 = apply to all nodes, number of layers Proportional factors
+LD2 0 4 1.000 2.000 3.000 4.000
+TAB-END & Indicator of end of tabular data set
+ENDGEO & Indicator of end of end control and geometry data
+DELTA-HR 0.0000 & Time step in minutes
+CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
+CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
+ENDSTEP
+DELTA-MN 15.0000 & Time step in hours
+DATE-END 31/01/2019 12.00 & Ending time for the time step block
+CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
+CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
+ENDSTEP
+ENDDATA
diff --git a/defaults.R10 b/defaults.R10
new file mode 100644
index 0000000..b4f988e
--- /dev/null
+++ b/defaults.R10
@@ -0,0 +1,94 @@
+INBNGEO MESH.GEO & Input binary file with geometric data
+METFIL RESULTS.OUT & Input ASCII file with meteorological data
+OUTFIL BRESULTS.RMA & Output ASCII results
+ENDFIL
+ENDLIMIT
+TI
+STARTIM & Starting date and time
+VERTTURB 0 & Vertical turbulence option
+MARSHOPT 0 & Marsh option 0=inactive -2 = input marsh parameters
+EDDYTYPE 1 & Eddy type 0 = fixed, 1= scaled, 2 = Smagorinsky
+TBFACT 0.200 & Turbulence factor when Smagorinsky is active
+TBMIN 1.000 & Turbulence minimum when Smagorinsky is active
+TIMPROJ 1 & Time projection switch 1= no projection, 0 = use time derivative (least stable), 2 = project u
+2DMODEL & 2-d model approximation
+GRAVITY 1 & Units switch 1 = metric 0 = english
+ZEROBED 0 & Force zero bed velocity 1 = at all bed levels 2 = limited by bed level input
+ZERBELLV 0.000 & Force zero velocity on water surface
+PASSSAL 0 & Bed level below which zero bed level in applied
+PASSTMP 0 & Ignore influence of salinity on water density
+PASSSED 0 & Ignore influence of temperature on water density
+ZEROSURF 0 & Ignore influence of sediment on water density
+SAVITR 0 & Save binary file for all iterations 0 = ignore, 1 =save
+REVRTDEN 0.000 & Option for vertical density model 0 = standard 1 = revised method
+OMEGA 0.000 & Latitude (degrees) positive in northern hemisphere
+ELEV 0.000 & Reference level for vertical transformation / Starting water level
+X-SCALE 1.000 & x scale factor (0.0 equivalent to 1.0)
+Y-SCALE 1.000 & y scale factor (0.0 equivalent to 1.0)
+Z-SCALE 1.000 & z sacle factor (0.0 equivalent to 1.0)
+BOUNDCMN 1.000 & Defines constant associated with vertical constituent distribution at boundary
+BOUNDCPW 1.000 & Defines power associated with vertical constituent distribution at boundary
+UNOM 0.250 & Nominal water velocity at startup
+UDIR 0.000 & Nominal current direction in radians c-clockwise from horzontal
+HMIN 0.000 & Minimum depth as startup (ignored if set = 0.0) Use with care can cause initial instabilty
+DSET 0.200 & Depth at which drying (element removal) is initiated when drying
+DSETD 0.300 & Depth at which rewetting (element return) occurs when flooding
+FIXBC 0 & Switch that forces boundary conditions to be applied regardless of direction
+VERTEQSW 20000 & Number of vertical velocity equations to switch to PARDISO (Default=20000)
+INITSAL 0.000 & Initial salinity when no restart file specified
+INITTEMP 20.000 & Initial temperature when no restart file specified
+INITSED 0.000 & Initial sediment concentration when no restart file specified
+INITUBED 0.000 & Bed x-velocity used to compute resistance to flow for initial conditions regardless of initial
+INITVBED 0.000 & Bed y-velocity used to compute resistance to flow for initial conditions regardless of initial
+RETRNPCT 0.000 & Percent of outflowing average concentration returned on next tidal inflow
+SURFMIX 0.500 & Depth of surface mixing (associated with free water surface local turbulence
+SSITN 20 & Number of steady state iteration cycles to be applied
+TRITN 20 & Number of time transient iteration cycles to be applied per time step
+TRSTEPS 12 & Number of transient steps to be simulated
+ECHSVALL & Show element and node data in echo print
+ITERINTV 0 & Iteration frequency for ASCII output. Skip if = 0
+ASCSVFRQ 1 & Time step frequency for ASCII output
+STBIN 1 & Save startup to binary results file
+DRYSW 0 & Switch to initiate element elimination
+FRQBIN 1 & Time step frequency for binary output
+FRQRST 10 & Time step frequency for individual binary restart file
+RWDASFRQ 999999 & Time step frequency for restarting ASCII results file (saves file size for very large number o
+1DTYP 0 & 1-D cross section type (see users manual for more details)
+SOLVER 4 & Equation solver type 0=FRONT > 0 PARDISO SOLVER
+SLVERDSK 0 & Switch to initiate out of local memory (disk) version of PARDISO for very large prblems
+MESSRWND 999999 & Time step frequency for restarting ASCII message file (saves file size for very large number o
+INITWSRF 0.000 & Initial water surface elevation, Required if different from tra nsformation reference level
+CONV-VEL 0.01000 & Convergence limit for velocities
+CONV-DEP 0.00100 & Convergence limit for depth
+CONV-SAL 0.01000 & Convergence limit for salinity
+CONV-TMP 0.01000 & Convergence limit for temperature
+CONV-SED 0.01000 & Convergence limit for sediment concentration
+EDDY-MAT 1 -0.20000 & Horizontal turbulent eddy coefficient
+EDDY-VRT 1 0.000E+00 & Vertical turbulent eddy coefficient
+DIFF-MAT 1 0.00000 & Horizontal turbulent diffusion coefficient
+DIFF-VRT 1 0.000E+00 & Vertical turbulent diffusion coefficient
+MANN-MAT 1 0.02000 & Manning/Chezy bed friction coefficient >1.0 is Chezy coefficient
+MANN-BNK 1 0.00000 & Bank Manning coefficient if appropriate
+MANN-SRF 1 0.00000 & Water surface Manning coefficient if appropriate
+MARSH-FR 1 0.00000 & Marsh factor when drying applied to friction coefficient
+ED-VTVAR 1 1.000 0.000 0.000 & Eqn for vertical distribution factor for horizontal eddy coefficient F=a+z*(b+c*z)
+BN-V-MIN 1.000 & Boundary velocity vertical distribution constant
+BN-V-PWR 1.000 & Boundary velocity vertical distribution power
+TR-V-MIN 1.000 & 2-D to 3-D velocity vertical distribution constant
+TR-V-PWR 1.000 & 2-D to 3-D velocity vertical distribution power
+MARSH-SH 1.500 & Marsh coefficient depth shift
+MARSH-RG 0.670 & Marsh coefficient range
+MARSH-PR 0.040 & Marsh coefficient porosity
+MARSH-LM 0.000 & Marsh coefficient limit
+TAB-END & Indicator of end of tabular data set
+ENDGEO & Indicator of end of end control and geometry data
+DELTA-HR 0.0000 & Time step in minutes
+CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
+CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
+ENDSTEP
+DELTA-MN 15.0000 & Time step in hours
+DATE-END 31/01/2018 12.00 & Ending time for the time step block
+CN-ELEVC 2 0 0.0000 0.000 0.000 0.000 & Continuity line elevation CN-ELEVC is elevation specified as constant along a line - CCLINE, Elevation ,bc-sal,bc-temp,bc-sed
+CN-QC 1 0 50.00 0.00 0.000 0.000 0.000 & Continuity line inflow CN-QC is standard inflow , CCLINE, Total Flow, Dir, bc-sal,bc-temp,bc-sed FL-QC means subsequent steps from file
+ENDSTEP
+ENDDATA