Welcome, Guest
Username: Password: Remember me

TOPIC: Liquid boundary file and TPXO database file

Liquid boundary file and TPXO database file 10 years 2 months ago #14102

  • anh_siwrp
  • anh_siwrp's Avatar
Hi, all.

In case 1: only one liquid boundary (566), i use Telemac2d to read from TPXO database file, It is ok.
But in case 2: there are 2 liquid boundary in my domain
(1): 566 (sea) and i want to read from TPXO database file ;
(2): 455 (up stream) Q is read from Liquid boundary file (.qsl).

I don't know how to set up liquid boundary file to do that? Please, give me advices. Thanks.

Regards,
NgocAnh.
The administrator has disabled public write access.

Liquid boundary file and TPXO database file 10 years 2 months ago #14104

  • gourish
  • gourish's Avatar
Hi,

As per you case 2,you may try putting following entries in the steering file
LIQUID BOUNDARIES FILE = 'QSL.DAT'

PRESCRIBED ELEVATIONS = 0;0
PRESCRIBED FLOWRATES = 2.;2.0
PRESCRIBED VELOCITIES = 0;0

OPTION FOR LIQUID BOUNDARIES = 2;1

and in the "QSL.DAT" file, insert the discharge as
T      Q(2)
t      m3/s
0      2.5
1800   2.6
3600   2.6

hope it helps

With Best regards,
Gourish
The administrator has disabled public write access.

Liquid boundary file and TPXO database file 10 years 2 months ago #14119

  • anh_siwrp
  • anh_siwrp's Avatar
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.
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.