Hi, everyone
I modified the soubroutine NOMVAR-TELEMAC2D and PRERES_TELEMAC2D respectively to give the names of the new variables and how to caculate them.
DO N=1,NPOIN
if(ABS(H%R(I)).LT.0.5D0) then
PRIVE%ADR(1)%P%R(I)=0.D0
else
PRIVE%ADR(1)%P%R(I)=0.0273D0*2650.D0*(U%R(I)**2+V%R(I)**2)/
& (GRAV*H%R(I)) !挟沙力公式S*=0.0681*U**3/(ghw)+0.0295
endif
ENDDO
DO N=1,NPOIN
if(PRIVE%ADR(1)%P%R(N)>CS%ADR(1)%P%R(N)) then
PRIVE%ADR(2)%P%R(N)=1.D0
else
PRIVE%ADR(2)%P%R(N)=0.25D0
endif
ENDDO
DO N=1,NPOIN
PRIVE%ADR(3)%P%R(N)=PRIVE%ADR(1)%P%R(N)-CS%ADR(1)%P%R(N)
ENDDO
There is three new variables needed to be added. Is it right? And whether I can use them directly in other soubroutine, such as difsou.f by added PRIVE%ADR(1)%P%R(N),PRIVE%ADR(2)%P%R(N),PRIVE%ADR(3)%P%R(N).