From e2cf7fcceab7108dd729c07df5bfdd1c71cf0374 Mon Sep 17 00:00:00 2001 From: Ian P King Date: Sun, 18 Aug 2019 11:03:22 +1000 Subject: [PATCH] updates and reorganisation of input --- LINELABELS.xlsx | Bin 0 -> 13926 bytes RMA10PRP/DUMMY.R10 | 124 ++-- RMA10PRP/New Text Document.txt | 4 + RMA10PRP/RMA10PRP.u2d | Bin 0 -> 38 bytes RMA10PRP/RMA10PRP.vfproj | 2 +- RMA10PRP/SRC/BLKRM10.F90 | 4 +- RMA10PRP/SRC/CONVERT.F90 | 44 +- RMA10PRP/SRC/DEFAULTS.F90 | 26 +- RMA10PRP/SRC/FORMFILE - Copy.F90 | 195 +++++++ RMA10PRP/SRC/FORMFILE.F90 | 93 +-- RMA10PRP/SRC/FORMFILEV2.F90 | 937 +++++++++++++++++++++++++++++++ RMA10PRP/SRC/READRM10.F90 | 23 +- RMA10PRP/SRC/RESOURCE10.res | Bin 40988 -> 41032 bytes RMA10PRP/SRC/SHOWCNTRL.F90 | 115 ++-- RMA10PRP/SRC/SHOWFILES.f90 | 73 +-- RMA10PRP/SRC/resid.f90 | 8 +- RMA10PRP/SRC/resource10.rc | 63 ++- RMA10PRP/SUMMYSH.R10 | 52 ++ defaults.R10 | 94 ++++ 19 files changed, 1616 insertions(+), 241 deletions(-) create mode 100644 LINELABELS.xlsx create mode 100644 RMA10PRP/New Text Document.txt create mode 100644 RMA10PRP/SRC/FORMFILE - Copy.F90 create mode 100644 RMA10PRP/SRC/FORMFILEV2.F90 create mode 100644 RMA10PRP/SUMMYSH.R10 create mode 100644 defaults.R10 diff --git a/LINELABELS.xlsx b/LINELABELS.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..a2d9fca078e322dc2d4c6290af70a2afe81011bc GIT binary patch literal 13926 zcmeIZbypov_C0)Y*WeP|xma)wZVB#g7k4MPySqbhcMBGRy9aj&?(RP_@_c48^SpxZ z^dGm^UDf-n>Z-HroIbTnP8tFd0{{(x1pojffGm1T8w)T102K-Vcn5$5*A%v~b_7{F z>MFU}f*iCNT&*mLze0jj=K#Rpp8wzVe|QAS6S`!&nGl0cCEA1&>CEfwKEhTk>3>I~ zWPVRfeu4-5{FM`bX=%&YP^h?E=&(Q@|M9i6*_pCV+prc0IWv}IsTH9Ht-@V!TxxVW zy!(^|!U6VlEQbqcj-&+BbYRb;#;9Q7;q~ljTV|`G1&@|MKePai3-X&dsObf!o=Y7*ru?7eR?u5+yGm$z_y=$UJiVl};)GR3$ur z2yyRrucx7<72e2$0g{_d*79(43?8yZmx`d|2RkPiT1xvQF}w1O9yI6K>)D%BaVa-y z=k{pYvgYDk>A^KJvFUS>8q{$HRXn(N1q4A@e5t+~eV^4=4elzzW`z`wD}t({ z0hGxYP#|SRzhfB}@2AMQ_SUFpJs-|=>*0XwE*m%-^b02{UVRH!|49;~5t&JO2mqiH z1^_^L8yQziMi(1<3j-S)i{C-4Tv^*DgA>D3*X$KyGgI#pX(%5LdO}P}HP39ZBSMKw z&kyr^Idf|0*O%M&8^rkhec?3cp7}exsPS{YxeV`lmkK-C-T~K6rgK~tq{ZYmgXiNP zPafK8v1So*SD0pQve>^M2uD${yD0FgE*|2Z}~v%cYyymnq6*hO@>tNtLQt%#qA zHok`;&FRo(+baN6j&jh2pnsZ}DrIsLANF=3)PnmFzB`yIz&j$x?pLDT!K?xy9Tbrp za2&A6ALCU(M25eGty%5e4r3#HVJUc;Is3mg@A!K_jQ2inoj^QJp}TcJOAojy>ANjfe|XmSufE zd2Ukj5^S9HSP5Ase1eeIdOk>QZhF*N+p)=`_FeI(S4!}Fs9y7QyskP0KE)KAOpP~n ztCf%4UJ1^OA88{K(PhD5&pu0oEsjW1u&sHozLJ`{ufw)6=S8O6yV`8#sUjMt>GsJ( z_`7@=HN^hLgX{O{K@u}~o3&HMO}leaQ%UQtv->b)5gBTsh3;-_QY-!oI1}CPHG-&b z=o4Byc0_9oMY9go=olUB1xD(4ua&>d%n{$eGx>Vs!QyRNDj0wY)A2<>C=V!SMdBSEho3{SYaDvs3@wk-?YtBTvwUW!u=g`p4lL~t zH4itu!YD`aHO$gU(5W4Cbb7vW(i~Q-kO||n8R~s~R&nd&xi%~mQB7r`+xf?j~h1cVL}wW4SeyR?670;lnb&mo>WM&O84@yP2l`w z0Hf)-eQKe=J@`Gwx|Uhu$L;&j)+KM3iaT;5JQGYKNqE0Qn*zi|xnL4Iz6H*bjMEsyeT#)WeuK~U@r5W; z<7~tN(-^@Dz5)qN#9_J7W3LvirRJ=BJkOJ{I|x`TXA^x;BEB7zDrcHhv1r?%$I-MK zf^nhvN425z(IlK&Hr$R30chc#cl{-p+nnkFExPZ-49l_+*`58k8j=|7o6M;ZztFDD z_m7O%G26U$4b;u=3j9|enmG|64L|?@dbt1qoVO?bv4T68fl8|l34xetKs}C8OE(zKnP?dge6wvL z%0w#Ihf1+WkJ%kWH+j)Fj;l`n!8PA9qr{Kxruj~z7*V!5Uylz*EqS?bN6%XpT>5YH zpYr&ky)~X*TyC>h_LuG(tRB(SIJ1cQsyttvIZ`t?s*ZEH)|_pQ(Zm{8_9w2+w)R2x zUOWw>b@JRtR11z8pOHsTI2*A}Waj4g8kbfIwO!J0h24)TVij9$zI^xgT4~Q@{E3`y zF+!L&-E`ZO{@KM9;+<^kbLFdkM3i{&=*l=_|H%yRvHt3s@}#BMn(Ubk^y^M?f6AJ> zEcX$6(aEyaqh>X4P0I8BIMt=?EoIS-VUaFz+K?AMxiWGi3v1G=nBlD6 zr`AKq+YWaxoL8#hU%d|&H#^4c+NJ4VwZ?3WhmXz+vr4sQbq^L)8y-$w?OxmbIm%7g zxEq#&&kdcMvOf36_qN03cpoi}f4F}vw#o9cx;-wm*&s3JJ!oyxU9{>nzjc>mSm$^> z`PN;PF4j59C31bXmF*Pb7xPjObx>6KSUCEM7*q(@wj41h9xYt+tjCUbP{Vz#oeyh) zX2*@0?9`R>Vx!w&EShv^wOYM>)@>aZXlzZMd%up=MByqCJil@o?c~1t(XzYS=2Bs4 z#nGJ!{FR8%xmsdllpyW&*ZBB~!z=s)Q|sts%PI1!W2ddT3Rs&B#`<{pkEh0=y40cd znvTsL{)rLem8;Rsml@9}V&09x-6kP*tKX2_=63uY? z)C<2*a-htSF;Lstm?ffiO|@hnz`oncWz040m&+=Ml;wMDdeB!7NEHY$n^R zy!v)8U9n>&4Ez-Pc#gCvgSxptgpxL(Kr{atE2O$mWvICGxVCTR@&swl%nj)K8?Oi+u2;Ya>Ru+OAq zJM+aqczq-%xDpk@VMh}uQRo>;KSX2T(}qTUTr31bRQu$QOamfhj!c{kZt|v6ObTLl zb&8T7)AdGIFG0}Uc*IsWgCo6DVvc>ox+Wfs*u5t(L)WAcSB2{E-T|qK@#}nPh}bCr z6tYI{%#jZ%MAcn@J0#I_2?i;F8PL)h3c~n?M6oF&n z`to#W?@|7{6Dfssl86?u)f`qR_`aGI+A!a~nqk*?g_aqfgN25}RtjBrZ-J5XC=S@* z6?|SUdK%3r**=>zDKLgTYLq1cELJFak(66y#0W$WNie)f<)KRdwxnl3O<~c6st!Ue zk{To$>Cr=0vrs6|+g3}qL9Nt5g^w0pJZP`J8=S>F$#c|_MI2CD~>nuZ-UmfVJ#g0}R{6xQHuQ%SGw0)dx0!TF-(U0VQCX8b_KALtu&Sxf9i4e`Gk& z7zZ>G4Dv~&8&QMx*~!3G*LPdQ%HFEVfUGDtY7mWpP#^kX8e)+eH$>WCnI)n}Y;_1E zXPO%28dqQ-&8%OcWx@i_0zoE;nyb$<$t>`;qvDvg8GW6_va9=!03oi|M4r0%3|?3i zrbbR4t9BwmFOE*vE}{3Fj19L4o=~T2EmEoGZoEUp8ZVwNB+l|X1|v%%VM;#kQ@!sL}gQ{{~z#pI(u zOwVTg*eulZ1*m^JmsLYxkEl-M=c+|gzeIWq`a2|`qJmrgBJJ-+pxPU*P}K^pYv{_H z9`gHzH84&yygr^zT#EvX&LOceL=RNyH~BquEscgo!MgULK!TJZrZKA*IX5O$N>iko zMIIMnUH!B;^4cPQb=x%Fw9K|%8Fewv4BiM=u*Iy($^Q%)C=B9F~Qu@w?uTxL6yNa7R=^iIaXff}6pt_GXYe!!&zjd23^Ujw2{EYbRQ5T`K2DJ));iG-=ur=pl)w2#C^H3H zC||)@j|MpS)))v->jO>0o+tfoV-MXRr88nH0~(;mlUUG0Ms29S z6O);lf1r(TM$E4;shJ^-IN_Bx3&o&ntYXFzG9ki)UEm2(3Id_0`_>0aFEX9_LGV0Z zSaUXSRcS6zn-W7>=QGhyR~mBi(+%eOjleFFHa1^u!bU*b7%_2^wVLu`#O;+kI;xPV z+2ei>1TXktU(glkcl=`==hJ&%AFG3936LUXEUv_o+INNsEqGg38j8V%w_bneHslp+ zRBUmGcP)41RUwTYOst2*jk5;f8o$jP%X_k5%r`1TCnmojwB%m=NV1PR8dK2&D#6b)G{g+g~v(nClu_esHHB%-GI^eD2YwQLztCh#tD? z62p|d^N^9OSH)e-bEvKJY4s-xmugHQ{yZ`K5)wUhJt(G(w$Vi!3L~h|h8CvN_;r!} zL5GEu=-OjEB-Rme`Bt(=4c8#zUgreh&nn26%s`ZwjJZ%ustUvPK%P~ODBCwNY>b%? z_hm4}Vnb?`UpobN^dRCd!+nV+2!F? z&87F79ZA6smNVkxYKhBgn;o(SvqLrJf#n*5SOi9}#*DJH@-`Ke!=2sA-4cFh(gj?a2b*SlH%3r5#t8RosCrD zTjOf?U~s`zlC1lL0Q44Ul8@ckYUq)($xF}Vh(?hpERf~PPwZIWdiO7;`S+FaDTr?? z@_LOBP#WI~O);4s^3F6L)rYZLfj`>TL-}1tTtj$cNRAQTfH)|wiO6q3wrDji`9^kh zALPOteRC$5_EL`Ew`tV*#Nr(+J3C^dkIC=E=h5qQzmLhEY`LLHdxvVtyXdU#eJ|*3 z$?oCPa!2%EGxsiR)DD@qx=SDN?@ffiW$uoqAS)2#-+up=v>&ODhU0T!bYedUAUe1_ zv28_Dtgeq+#;=eYWhM|bG#<&Td|-)d!+!@(!F8=bonMeFWW$#*^$`N8{VE!cyndf( zl4^9GtYU{e-BL=!87W-+rM;-=Hp}bFrQ`Wjd&(E_6oemfNxB_!7kn8J^C{%hN!}Km z)T|z|5(yD#AUNxL1o|iGrED*IG^;__{&|Rd5zE{pB1x_#TC85!pxEZ`#+<%ulo9&g^(hs*HlYuC~|AI zx@(0(0{0Y72uT(%^s|_G-fy~3-BP2Ydzs!-Jnv#RF=9=sP>Sq8s|#D_yt1@>((_Q8 zbhy~>GM^%b_1ADO(1j@0>c2jG*63`UDgrZQ^-C;loaa7Fw2sk7l6#oVZB;#09>!Xo z(ABPMAuNp!oXb)2t_~tX{I)KNP$fviuhku#XfI`$zhHYD&A5CFQrVG$;wLpyjvbAN6Bd4e54Zxk6H9)gwqB?jzQT>GLC@f)#wWU* zSq@^xEI#*_v;E89PCoCez58z6+Li`dX=@k?Ma7fZXx-Q6XY5to&X==;bB_HPs`DQ` zolp0}$(^s~+}G8c%b1KEo);JU@z`sg*SoUOgo{fAgehkvp*3ff8-}Np;rhPFkLO?l znLSKGJ#in=_MD4Ph6liPG>H#>I8`72SRzpCgKFdLv{?8?blyT9m0BfM1GQ`z;f&+L zygew2-rFO0ly>ID)XO zRXAj4j*n&YQ+4Ro32K#p$TWqpy&i=EL%py;$nXykalqom=&Ds!cd^ARpNA2d{#!qo zyh2KL`i4yzkK_GP+Y0SPa6AVSUS+{d4g57iH)6~D;9Mq+FJfS zp2&)nHs?-;(wEFAqt(}sGk+A9}$zQNeJhW};ylm?F^>kPH zjyoW6JdDa53)49IpN#S4op9uJk%cgH zOagQkQq?9M)h6W*l5;i0>~`RptS~WMcW%~2+joKm7X^ci2@-YgR|K_m*Ag(xXivo9s3;zSRR^?idmY!fXRtHu4hJ)G`2SLhb98BOf3ip|crl3CVw zF8WEh*tNkT3EW!oGDFbJu_^8eJwRQccOnw&B*9lJI#nN3evPfPk90)ocBqa#1FkGx zWM)AtIPzICM7|q}m9{UXy=#Q!8p+KD2g+ab7jIneUe)EMbtqC$Pb{p}A%nSIl% zF=Kb`T{iWGOy}CdrnnNK$hgDhoWG_RW&L0atsI$@r6aacm-maU zH)oN;o}7?PG8y&56#hU0yI%Nmym}<;S{oM?@YRdP5@tne9+-fq--)W^uaJw@SQhAv zqcZ$lt##l;i?Xg4E3bt-efR*Tr9UC+Y*y7j$8fT>0~`%R3lIOKJA&P?@vzyUwcufQ zibkCI*f5l%ZslwqXxzX@4cRV}H*xfnJ5Z-W^AzdJO7$cG?^!z;67u4#wcH^Ge1m*}&QNcM#?z<+aBJS$`1)VBKHPGBEc4&?1{I(I0HnW% zUk67wOOV6wYFWOrmBkVh>I?t8H|oR5kn(9b5JKgW)~vaVR32kL>S0* zoW5glty8*U3E`$0G(zGd$KI{En+8SFxXdDXoM@fm$JoFg$335hR1>LRU;7e+L5|!_ z6PrNCOH6NcFF+iZw%Z6N1;jtNOG7xh;lmQE)mpYmV0Q%im}7r74gF?nZZMN-7`OUJ z8ogaGsWNv0ZyQ=^p;MdexN9V1>Lv!hJx3x62e%cz$a5)`565iT-|oRFYf+IQT40?B z!R%XX)3vF-e_DE$-R2N7r1%JNTnvy6kjzTbF!d*pfXO`M3zK3OSyt$p3wo?hlVZ5> z=edc}ITp^?o^U@U&+6B&u?<16Gr)8QUmsA68E4LepU?4suD2m>_mjledua9+u|ZrB zw`1Y{_z>Xs?PsCiAph;w-jHIKIgAL>n&tC>BSfdgF+p*55oQs7Me$hyA~6-j5KSs- zwAVxDEcU#1q!Jy%{z(~PH{Zo)&NR7!E`qb5$bO}MxTCN3D*dg~Tjt!G~*J(ZOfs9e?W_hoY&y9rR-@-=CVKln2FjqfbI4KQs$9$^nZ zTo$iUaru`}bV^JQ9ID6o)NmNpm190xCzmNR zqimhs%c3+>amGK6KH&lbahjXC5k=xu_NZyz<&RH9$o_Oz&-9LS1h#n|tM<~{1 zBxTZN`6Uxl6uq}>-^eW)V=%=1V#tJ)T%MJeGs~maR*R);R=zAx$dEAva;XzTFiI-s zGF$Kk|G(`eWY&7|ox%YCOr!t+#y{NIRNo$Cr08gGW^M8}ukO)UvR~pv_gpJ|hDjhQ zJd(s1&RUuI8FfrEn2$cHkxXL7mCBJ8Kwj=9FG1b-bnT8uT!J9+;oR7QNTE^lW8S3a zgCj>m_*&Lxi06(v!_@Z0-jVig$T_)$a4bDTv~|1SE|EBWvKRV;eA?M?9@2QUSL?%} zF2jWu%hBhH;k<%`GViRcUtMwiUM~#V=aAVW)ATnKFC$n*_q1DHDU4jnc8*_Gk2?e+ zJ!1)Kh_$KEcRRejdwz!A)&80**nZkOAjMY2z8J#S$rvgNSA39`)~ecd?VikaHRj23 zIlMyI$;;aE%J|eS++wAk-9E1PD0gxV({5qdFh_*c-#4{kdhp3h03%y7FBG&!5GKIc zzA*9VLH#oAOzF8)igtUPp1zRlvFM=PWI|0ddaTou#%tCNYCrSqQe0~Ktn*0(j3 z)kzlku*9io1zl%~uiubEZ-^@OCvADJ%!R}a#@P*vlG1C)&aII2B>2l>c^Y>kiYUL* zyyR_AL4co3eDFP)f!7|-`UVZC6!=!R8tFTcYr@lAW8#yytJ7NNdl{@BGRN^@%$`5p ztM^x0K@4{#n&|G>K?_?^rusT~cfcT$T5n9;g#D`%HPJ%l;c_HWgq+@q)`x{ksx5UB z*5#sr<&VhtM;yN> zEjV#lGx9jtEK9gdX7HnRCv$+$HMkUUb0Wftp<~}e8m*j+SIbQX_793vR;>5^ zE_QQC1n$odRE@nQal%NIHEF;F3~~oRlqYlO4auX=A2)?Hsy}@=&!9w?cEz!lBfs}x z0PlJ?IKMvSefe~G{+zzPyBpu_|FAVyt|)9Q#ozwHdY$X~FtN?J{U+vd@`@GSvj1L= zF&l5!oR0&V%e~#JpM6^V(IzFp7N3hIGfxt3Klti&E;hSnBAI`8F7GovW@Oh{FR)7S#VyFqiHdU;a$(@5}+>CyZ!o|RC_ayjRrP_Gv6 z`5BO>)rbpKsx!{uSKkYG&Zmpu&*9uQ?z!NT(b8E*G^kv3D}X-#ug!*=n~ zYn$yNgjGZWQ>Ze3MLvNOkJ(#U_eF1oQs`GffgT>DUo|y2LDA9nASt?!aCnwr8ko*F zuk->RilBfz9&C#|yyxH(wZcWn*1H^0QC{?#7&VX11uvELLr~u1d_JqbZUR540k~Hs z@*GiU00A0YxmAs4v&2-y*0x*yPjwcP79ud@W&rP9NQUsuB{@|WU4rEMD%W3RD6zf* zA9crsKO%`((N3e9!QRNklme7#FPwsmHi2^vo~T`X)L@Qm*l4N^0bj!mBorg*fkAJ_nZR}X&|F}8t>Ag;Wb-0t8)BZ91}$e@vs}h&jTa>Z6_sz=MpCsG zTu67L$)=GF4~<4mLNmXBkcI$x(v*q^$TEvFmR=EQCflMP1?N0_EGla5w8(jjNF~l8 z`L?HzIeJq-YoD+8=C(L)2S4F zE-tns(d`=LPwA{aRHL#SfNYMI-p29~tNzxDOH1h+9RT|T`C9CZ#sjZ~d5yOtCBic@ zj|1aWNWaHD!ag%Wa2Z0!FzlKEK0}3Aph>4{oAX15-U{Vv6%5yVTY>E}kBN_e{JSLR zGn-1fXI`mUX{B^$RaJ);5!TqBU}u^ns7y#3$fN?~+kD*hUoVV2Hca$;E;H~3t(cbN zs-sSQl`Tsi?%+R3y!e^eDNP!zH<^vs1f2xb7_P7AieD&vR?{>#4H>EqF;qTNdh`#X z*n-DT+quL<>Cu{-0{`x9VeZPt;B28&CY6Z}#49U>v4I&t(zbR))vdi1-IQ zOGu!ePwo8p?bMBXD}j#*XNP`6mYN@&v9Z*~Pq*~RthuCs8fG^kQj9(e1B`&>f}o*A zZ^E8t9vS9%gxM-MvJ7&F;(eg}3ptXEa{LcJq=R|=M(fkF+b^~@v!y?^6J6Puo2bQZ zi_$&S^2|TR3Lp2&y~~x!$t-B+`jsG|JggGFLLsE6K#!CTCkDwk_%-KDPcl94O-f~K#8U4Tuygh`6Wd{^f5F!yJ^Sk9CO<_6w< zbvq^I`9WIRX@{a!*G)nUcaHlQPa#%^tr0_((;Q0FJX2e^4;7AdSisV^mM)*Vs9ryM zQNfY$KCyO$>Ma6wlm}7usONx6@%1n>CWMwN(E#UmdryDsR_aK4Su!=%Hcl<%jKTTt z)oiQ_RUwc%t@icNVrc{Y6VLaLsG-bkz+O|W^-#mBA{45fIjt)@HyLKEgs~sCQRn=H zG|duGuk0$!4+vi^*kT%b57Nn$ouR*_0!whTgsL#j{1(P?1<8d4KXwEgR0ff($8b;i z?^E%lB$T#J=Q(+TYZbv6G*4_{wqszNrXA=D0wP^fc4W|OK? zu3o}zUktf`((4pzDtm19ya5iEl)C{|X=^Kj`XRU6Y6PRd=7|9tZcvkcUV+YE${b@E z8fTg4=G3mFK9dPZehe&{lYS{DD*EM$t~hSye!AU5jnj!cO8@0}wLh^7iQ(vR5Q1)O zv6B{Es(mw;U2#pJKzeKCbOv8JHQ$N5Q?tMedeu`8_^aB!IX9~LC(z1wD$PZ^K#@f^ z>|Tb-z)QKLZgDw^K~kAyv}movCh;r!G4Fxrqs+ClI=bt0qu0HV{3!Od?y%A+ap0D5 zpNsp$=31$2r(-=ek2r}Kq<2>|kSn?SI)z#z*7QdZb~lof@BA@_qztkI)!HoN{({amb%GlA6+{!^kG zGbvEp&eLHnNzpQ4{O&OZcO9jQE%ilT_;*xSKo$be$=9`h59+V7oxv|7n`s>_IzJ~W zyP}|ZH2nE95nNe&Y6FlCxKXz8phG~XxrAp*#HZtSj7?p_%kZ`;KDmC)n0l3aT8r#> zWM|CvR^R)5I>26`pOV*Ms=?S$Yfy)Ts1b7Mkpdq`$_oktcSwx&Iw8}s^vrY2z88X} zLti27PxHV#!5BJz@9I(`IJ3J7gh58;Tvr$>i_bP-NF5u z)TH}1=Km=HLdOaEA@MdVx^G{iXm1iABO60GdmCE^MnfBW(C=w^)2jYY`SUh(5pkn$ zI%?E_-#Y4T_k;ZXs>FP7DcOnuV%VP1+9oTJBeh17K9Bqq=hHS(J{m=yyO@v`F|b z7NXZTsh~w5E&1mduMy_g z%v7!&)YrnmqzuPkr$1ViOz)1aT-dco*;T2mRps#Is>^h-gStA-)a#| z1qsmEOm?$q<~rg8I0&TNZ;&JO8L^UW-HY?xt9lbw0Yj#rUQvtGnJDB%#N?x0q4SwL z9?MEh<_$H#;2f|k47{7 zy%EhF4H^fDS{`P>+^~R{cyvd>=B|!ErpGEq=of?s3OsJ3XgAtM4Ry65-&$oX!81j{ zW#6IUW$JhGVZJI5rWZ9IFVlXy4Tez##xT>d7n>&J1*C%%S^C(CMep~vJ~#sc{HVH= z2)bLG->0A2bHOXvGK@gKlQ7R9`x!*PzUG`q>iMxP<6Si3yVHQWQIF{nuC=>u)_sA~ zhEsq0;2BP_&PeyZN4ZB$(X+;}%dhj!v_kO=wpzC}#3xgfd83qvGA)FxHfgt!9#6sh zyTd!PU!PP8TgS}(vR|5?-s0Ea9T^N9_*Uxp?`LxUIdT8I|CbXwa?<~9;NK4n{1f>5 zp8v*pzzC-6@Z^Z33km>;!~YHZ|Dy=~)6Sndz`rboBmRG# z_zzX!pH}|dll;qy=$p9k?N$EVq5RXppH`U(oosn$~aV zpY^OiE&Q3{{AHn@{%;okE#dhS{qLdpFDw8c#Q*^OM_~RF{_jiazrqKZ{ssPzg;h=( V>diNPpG-jjbiBp>xDUVo{Xb2HAi4km literal 0 HcmV?d00001 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 e69de29bb2d1d6434b8b29ae775ad8c2e48c5391..d45a425983bd033bbb2491c43c3fb8924a021190 100644 GIT binary patch literal 38 ncmZQzU|`?_;vj}ph8!SHWGDubH4NblnGAUh#tcS4IR-`mTm=N+ literal 0 HcmV?d00001 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 931bcb5c077f9cca4cdf6c297adadd97ca5da64b..a46b84027e9cb9ecef257edb92f0191b3ee8a8dc 100644 GIT binary patch delta 116 zcmbPpfa$~mrU@GuB_?jnX4Kp)%UCJK-x$CE1nLai8G;yG7`P{Yl$PIoL|lc1S(qVq z^J{5oR(2x=Fxbqi_@0?@$>hH(>XUgQ1SZ#L$xJ>WA~!itl}`|41SniO%UCJKzcGLT2xJ(JFvKv(GH_4+C@sJFh`0(1vlv6_ w=GW5Ftn7vgAYizeSMfbFW5wjZD(aI10tF_o(2`+{njEOA&8V<>rRw}D055?X>;M1& 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