--- Welcome to the official ADCIRCWiki site! The site is currently under construction, with limited information. ---
For general information, see the ADCIRC site, ADCIRC FAQ, or the ADCIRC Wikipedia page. For model documentation not yet available on the wiki, see the ADCIRC site. New content is being continuously added to the ADCIRCWiki, and material from the main ADCIRC site will be gradually transitioned over to the wiki.
Talk:TAU0
Notes from code dive
I dug through the code, came up with some inner workings that aren't documented. The following info needs to be added:
- spatially-averaging is done to tau0 in: SUBROUTINE CalculateTimeVaryingTau0(TK, NNeigh, NeiTab, NP)
- I see Casey's fingerprints
- tau0=-3, -5, -6, -7: looks to average tau0 across all immediate (?) nodal neighbors
- Rick said tau0=-5 is ineffective
- tau0=-3: the actual tau0 is capped at 0.2
- tau0=-6 is a time-averaging:
C jgf47.33 Perform time averaging of tau0 if requested.
IF (TimeAveragedTau0) THEN DO I=1, NP TAU0VAR(I) = 0.5d0*TAU0VAR(I) + 0.5d0*LastTau0(I) LastTau0(I) = TAU0VAR(I) ENDDO ENDIF
- tau0=-7 is a different time-averaging:
C jgf48.42 Perform backloaded time averaging of tau0 if requested.
IF (BackLoadedTimeAveragedTau0) THEN DO I=1, NP TAU0VAR(I) = AlphaTau0*TAU0VAR(I) & + (1.d0-AlphaTau0)*LastTau0(I) LastTau0(I) = TAU0VAR(I) ENDDO ENDIF
- AlphaTau0 is hard-coded to 0.25, so 75% last tau0, 25% current one
T$Taylorgasher (talk) 00:10, 16 February 2019 (UTC)
- OK, done.Taylorgasher (talk) 13:49, 19 February 2019 (UTC)