Hi, all
Thanks you for your quickly respond,
Yes, I fixed it and the result is good.
I suggest that we should modify subroutine "TIDAL_MODEL_T2D" to be make sure water level and velocity are imposed on boundary points responding to code number (566).
Sample:
DO K=1,NPTFR
IPTFR = BOUNDARY_COLOUR%I(K) !code number (566)
IF((IPTFR>=1.AND.IPTFR<=44).OR.(IPTFR>=2680.AND.IPTFR<=2691))THEN
IF(NUMTIDE%I(K).GT.0) THEN
! POSSIBLE SMOOTHING AT THE BEGINNING
IF(AT.LT.1800.D0) THEN
UBTIDE%R(K) = UBTIDE%R(K)*(AT/1800.D0)
VBTIDE%R(K) = VBTIDE%R(K)*(AT/1800.D0)
ENDIF
IF(LIUBOR%I(K).EQ.KENTU) THEN
UBOR%R(K) = UBTIDE%R(K)
VBOR%R(K) = VBTIDE%R(K)
U%R(MESH%NBOR%I(K)) = UBOR%R(K)
V%R(MESH%NBOR%I(K)) = VBOR%R(K)
ENDIF
IF(LIHBOR%I(K).EQ.KENT) THEN
HBOR%R(K) = HBTIDE%R(K)
H%R(MESH%NBOR%I(K)) = HBOR%R(K)
ENDIF
ENDIF
ENDIF
ENDDO
One more time, thanks.
Best regards,
NgocAnh.