1) How is Vertical viscosity and diffusivity calculated in EFDC model ?
In EFDC model, subroutine CALAV calculates vertical viscosity and diffusivity. There are three options available to calculate the vertical viscosity and diffusivity. The three methods are :
a) Galerpin
IF(ISTOPT(0) == 1)THEN
SFAV0= 0.392010
SFAV1= 7.760050
SFAV2=34.676440
SFAV3= 6.127200
SFAB0= 0.493928
SFAB1=34.676440
RIQMIN=-0.999/SFAB1
ENDIF
b) Kantha and Clayson (1994)
IF(ISTOPT(0) == 2)THEN
SFAV0= 0.392010
SFAV1= 8.679790
SFAV2=30.192000
SFAV3= 6.127200
SFAB0= 0.493928
SFAB1=30.192000
RIQMIN=-0.999/SFAB1
ENDIF
c) Kantha (2003)
IF(ISTOPT(0) == 3)THEN
SFAV0= 0.392010
SFAV1=14.509100
SFAV2=24.388300
SFAV3= 3.236400
SFAB0= 0.490025
SFAB1=24.388300
RIQMIN=-0.999/SFAB1
ENDIF
a) Galerpin
IF(ISTOPT(0) == 1)THEN
SFAV0= 0.392010
SFAV1= 7.760050
SFAV2=34.676440
SFAV3= 6.127200
SFAB0= 0.493928
SFAB1=34.676440
RIQMIN=-0.999/SFAB1
ENDIF
b) Kantha and Clayson (1994)
IF(ISTOPT(0) == 2)THEN
SFAV0= 0.392010
SFAV1= 8.679790
SFAV2=30.192000
SFAV3= 6.127200
SFAB0= 0.493928
SFAB1=30.192000
RIQMIN=-0.999/SFAB1
ENDIF
c) Kantha (2003)
IF(ISTOPT(0) == 3)THEN
SFAV0= 0.392010
SFAV1=14.509100
SFAV2=24.388300
SFAV3= 3.236400
SFAB0= 0.490025
SFAB1=24.388300
RIQMIN=-0.999/SFAB1
ENDIF