Index of /~longmtm/funwave/funwave2d2.0
Name Last modified Size Description
![[DIR]](/icons/back.gif)
Parent Directory 13-Jul-2007 14:28 -
![[ ]](/icons/compressed.gif)
Funwave2D248.tar.gz 09-Jun-2005 16:01 15.3M
Funwave2D249.zip 08-Aug-2005 21:54 14.3M
Funwave2D250.zip 18-Apr-2006 21:35 25.6M
![[DIR]](/icons/folder.gif)
Tsunami/ 16-Jul-2005 12:42 -
funwavenotes/ 23-Aug-2006 00:01 -
gallery/ 30-Oct-2005 23:39 -
//-----------
Funwave2D250.zip (windows zip)
//-----------
Sub version 50 of funwave2d2. (04/18/2006) It's based on version 2.49 and the
modifications are:
1)With bugs fixed in directional Wave maker
2)Sponge layer get_FspGsp
3)Inclusion of variables Umean_acc and Vmean_acc for
better calculation of mean depth averaged velocity.
For details, please see the readme.txt file after unzipping the package.
Wen Long, CACR, 04/18/2006
//-----------
Funwave2D249.zip (windows zip)
//-----------
This is the sub version 49 of funwave2d2, i.e, funwave2d2.49 (08/08/2005)
It's created based on funwave2d2.48, and major modifications compared to version 2.48 are:
1)The slots method is eliminated completely, common /slots/ is deleted
2)Wet-dry check algorithm is implemented based on Wen Long's research (April-July,2005)
3)Total water depth HH(i,j) ... are included for the purpose of checking wet and dry
4)filter2D() subroutine changed to be a much simpler procedure
5)subroutine WetDryCheck() added
6)namelist about SlotConstants canceled, new namelist WetDryParam added, with Roughness = film thickness (m) for
dry points of swash zone. Suggested value is 1/1000 wave height
Comments: this version is not well tested. The inputs in the examples may be inconsistent
//---------
Funwave2D248.tar.gz
//---------
This package contains the Funwave2D Version 2, sub-version 48, as a beta release,
2005, May 31st.
New features compared to Funwave1.0 are:
1)Staggered grid system
2)Curvilinear coordinate system
4)Can handle arbitrary shape structures in domain
5)Moving velocity reference level
6)Multiple reference levels
7)Vertical vorticity retained
8)Symmetric boundary conditions
9)One-way wave maker
10)Periodic longshore boundary conditions implemented
To install in linux/unix:
1)
gzip -d Funwave2D248.tar.gz
2)
tar -vxf Funwave2D248.tar
3)
make
The following bugs are fixed as of 05/31/2005
1)bug found in derivatives calculation around boundaries,
in the following subroutines: der_xi_at_v
der_xi_at_v_symmetry
der_eta_at_u
"NO.eq." ---> "Np.eq."
2)bug found in subroutine get_gxx:
if(ir.lt.1)then
il=1
coefx=1.0
endif
-->
if(il.lt.1)then
il=1
coefx=1.0
endif
3)bug found in subroutine get_E_M_xi_M_eta
"*mu_scale" --> "*mu_scale2"
this bug affects stability when grids are not cartesian
around complex boundaries, will lead to noise
4)bug found in expression of pressure calculation, subroutine
B_tmp=0.5*A_tmp**2-h(i,j)*A_tmp
-->
B_tmp=(zeta(i,j)+z_tmp)*A_tmp/2.d0
5)filtering added for F1,F2, G1,G2 for these cross derivatives
are major reason of instability around complex islands/structures
call Filter2D(F1,2,0)
call Filter2D(F2,2,0)
call Filter2D(G1,3,0)
call Filter2D(G2,3,0)
are added to where subroutines get_F1, get_F2, get_G1, get_G2 are called
6)bug found in subroutine der_U_eta_at_u_symmetry_high
elseif(Nphysi.eq.1.OR.Nphysi.eq.4.OR.Nphsyi.eq.31)then
-->
elseif(Nphysi.eq.1.OR.Nphysi.eq.4.OR.Nphysi.eq.31)then
7)in gridtype.f, bug found in index:
"ii=2,ny+1" ---> "ii=2,mx+1"
8)in gridtype.f, the depth is checked to make sure no water domains
less than or equal two grids wide in either x or y direction exist
Also the depth is checked to make sure no water depth is between 0 to 2 mm
The user is responsible for making sure no land domains are less than two
grids wide in either x or y direction. Otherwise, ambiguity of gridtype will
occure.
The following bugs are fixed as of 04/14/2005 (sub version 46):
1) in subroutine der_xi_at_z_symmetry
"trd1(i,j)=(tr1(i,j+1))-(tr1(i,j))" should be
"trd1(i,j)=(tr1(i+1,j))-(tr1(i,j))"
2) in subroutine der_V_xi_at_v_symmetry_high
"i2=j+2" should be "i2=i+2"
3) "ibeg(i,1)" should be "jbeg(i,1)" somewhere
4) in subroutine wave_maker_preparation
"elseif(FreqNum.eq.1.And.DirectNum.ne.1)then !Only frequency spectrum"
should be:
"elseif(FreqNum.ne.1.And.DirectNum.eq.1)then !Only frequency spectrum"
5) in subroutine print_animation
"300 format(I10,I10,1x,f20.10)"
should be
"300 format(f20.10,1x,f20.10,1x,f20.10)"
6)subroutine damping_for_structures is modified and a variable
"DampingWeight(iq,jq)" is added to file "paralong.fh"
7)some typo of comments and texts are corrected such as "fuwave" should
be "funwave",
Comments:
1)If you want to run it in windows, you need to change function call "rand()"
to "random()" in the funwave2d248.f, please search "windows" in the file.
2)For generalized curvilinear grid generation, please download CoastGrid in matlab language
by Dr. Fengyan Shi, go to http://coastal.udel.edu/~fyshi/coastgrid/coastgrid.html
3)For very complicated geometry and bathymetry, you may try to un-comment the
line "call damping_for_structures" in main program of funwave2d248.f
Reason is the code may still have hidden bugs on boundary treatments as
of May 31st, 2005
4)The program is still under extensive tests. Please send your comments to
Prof. James T. Kirby <kirby@coastal.udel.edu>
or Wen Long <longmtm@coastal.udel.edu>
or Dr. Fengyan Shi <fyshi@coastal.udel.edu>
@Copy Right, CACR (Center for Applied Coastal Research), 2005, May 31st