Welcome, Guest
Username: Password: Remember me

TOPIC: Spatially varying salinity and current forcing at tidal boundary

Spatially varying salinity and current forcing at tidal boundary 1 year 11 months ago #41787

Hello,

I would like to implement salinity and current forcing at several locations along the tidal boundary in my model. I have a time series of salinity and current data at 5 locations along the tidal boundary. I understand that if I had a time series at just one location along the boundary, I could include it in the liquid boundaries file, but I am not sure how to include temporally varying salinity and current forcing at numerous locations on the same boundary. Based on the manual, it looks like the solution is to modify BORD3D.f. Is this correct, and if so, could you please advise on how to proceed? I am not sure where to begin in terms of modifying BORD3D for my purposes.

Thank you,
Taylor
The administrator has disabled public write access.

Spatially varying salinity and current forcing at tidal boundary 1 year 11 months ago #41795

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
Hello,

That is right - this is a complex and custom setting of your 3D boundary, so the best way is to program your spatial and temporal interpolation through BORD3D.

Please note that you can use the liquid boundary file to load your 5 time series - you can call them MYV(1), MYV(2), .. MYV(5) in the header names, and then make a copy of the VIT3 subroutine, say MYVIT3, changing the test on FCT='MYV( ' ...
You can then call MYVIT3 instead of VIT3 in BORD3D and do your spatial interpolation, knowing that TELEMAC would have done your temporal interpolation.

Good luck,
Sébastien.
The administrator has disabled public write access.
The following user(s) said Thank You: taylor.l.bailey

Spatially varying salinity and current forcing at tidal boundary 1 year 11 months ago #41796

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
Further to my previous note, what I said would alos work in parallel as the liquid file is shared across all processors.

Do post what you have done so we can help.

Sébastien.
The administrator has disabled public write access.
The following user(s) said Thank You: taylor.l.bailey

Spatially varying salinity and current forcing at tidal boundary 1 year 10 months ago #41863

Hello,

Thank you for your reply. I am a bit lost on how to begin to modify the BORD3D subroutine. Could you please provide more details? I see:
UBORL%R(IBORD)=-MESH2D%XNEBOR%R(K)*VIT3(IVIT,AT,N,INFOGR)
VBORL%R(IBORD)=-MESH2D%YNEBOR%R(K)*VIT3(IVIT,AT,N,INFOGR)

which I believe I would change to:
UBORL%R(IBORD)=-MESH2D%XNEBOR%R(K)*MYVIT3(IVIT,AT,N,INFOGR)
VBORL%R(IBORD)=-MESH2D%YNEBOR%R(K)*MYVIT3(IVIT,AT,N,INFOGR)
but I am not sure if that is the only place I need to make a change in BORD3D.

Thank you,
Taylor
The administrator has disabled public write access.
Moderators: pham

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