Download W2 Version 3 Release Notes - Civil & Environmental Engineering

Transcript
#
Code: W2
or PREW2
or GUI
Fix
or
Enhancement
Type
Description of Bug/Enhancement
Date
Bug
Fixed
or
Enhancemen
t added
25
W2
Initial
WL
Calculation
Changed SLOPE to SLOPEC in init—u-elws.f90 routine
since the normal depth should be based on SLOPEC.
4/5/2013
END IF
FUNCVALUE=FLOWXAREA*HRAD**0.6667*SLOPEC(JB)**0.5/FMANN
! SW 4/5/2013
RETURN
END SUBROUTINE MANNINGS_EQN
Also changed KB(I)-1 to KB(I)+1 for ELWS:
IF(ABS(DX).LT.XACC .OR. FMID.EQ.0.)THEN
ELWS(I)=RTBIS+EL(KB(I)+1,I)
! SW 4/5/13
RETURN
Also changed KTTOP from REAL to an INTEGER:
REAL :: XAREA, WSURF
INTEGER :: KTTOP
! 4/5/13 SW
! 4/5/13 SW
25
W2
Output
for
pumps,
spillways,
gates
If the LAT option was chosen, the output files index for
JWD was incorrect. This may have affected output
temperatures and concentrations.
5/17/2013
26
PRE-W2
Mass loading
calculation
6/21./2013
27
W2
Assorted code
updates
There were cases where the preprocessor bombed
while calculating the mass loading for output to the
pre.opt file. This error has been fixed.
Minor format errors (that were ignored by compiler),
update to code comments, and faster code
initializations to speed up model performance were
performed in several subroutines: input_PAR.f90,
temperature_PAR.f90, transport_PAR.f90, update.f90,
and w2_37_win.f90.
6/21/2013
An example of an initialization code speed up from
temperature_PAR.f90:
New code:
DO K=KT,KB(I)
AT(K,I) = 0.0D0; CT(K,I) = 0.0D0; VT(K,I) = 0.0D0
! SW CODE SPEEDUP 6/15/13
ENDDO
Old code
AT(:,I) = 0.0D0; CT(:,I) = 0.0D0; VT(:,I) = 0.0D0
W2 Version 3.71 Release Notes
11