Welcome, Guest
Username: Password: Remember me

TOPIC: U and V at velocity boundary

U and V at velocity boundary 13 years 1 week ago #3017

  • olslewfoot
  • olslewfoot's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 132
  • Thank you received: 3
Dear all

I have been following the method decribed in post #745 and the following trail, to supply values for U and V separately at a liquid boundary (rather than velocity normal to the bounday).

I thought I had it right but I still have compilation errors.

I created a liquid boundary file with required VTU(#) and VTV(#), created functions VTU and VTV (copied from vit) amenduimg FCT as suggested, and included them in my fortran file. I also ammended BORD to read VTU and VTV with:

IF(PROVEL(NUMLIQ(K)).EQ.1) THEN
UBOR(K,1) = VTU(NUMLIQ(K),NBOR(K))
VBOR(K,1) = VTV(NUMLIQ(K),NBOR(K))

There are other declarations required for my functions vtu and vtv (interface_telemac2d.f ??). I would be very grateful for any help with where the declarations should be made or if a different way of achieving the time varying unput of U and V at a liquid boundary is possible.

Thanks
John
The administrator has disabled public write access.

Re: U and V at velocity boundary 13 years 1 week ago #3018

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

you do not need interfaces but since your VTU and VTV are FUNCTION (presumably) returning DOUBLE PRECISION, you need to add them to the declaration of BORD, in addition to the EXTERNAL statement:

DOUBLE PRECISION VTU,VTV
EXTERNAL VTU,VTV

Hope this helps.
If not, please provide the compilation error and your PRINCI file.

Sébastien.
The administrator has disabled public write access.

Re: U and V at velocity boundary 13 years 1 week ago #3020

  • olslewfoot
  • olslewfoot's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 132
  • Thank you received: 3
Thanks Sebastien
They are FUNCTIONS - I have declared them in Bord as you suggest.

That is much better - but I have one compilation error as it looks for VIT = VITES(I) in t2dfort.f

t2dfort.f(471): error #6404: This name does not have a type, and must have an explicit type. [VIT]
VIT = VITES(I)
^
compilation aborted for t2dfort.f (code 1)
ERROR : COMPILATION

Can you suggest how to overcome this error please??
John
Attachments:
The administrator has disabled public write access.

Re: U and V at velocity boundary 13 years 1 week ago #3023

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
You may still have the line VIT = VITES(I) in your VTU or VTV functions, which you can replace by VTU = VITES(I), and/or VTV = VITES(I), appropriately.

Also, make sure you have deleted the follwoing line from your new functions:
USE INTERFACE_TELEMAC2D, EX_VIT => VIT

Note that it would be easier for me to get your PRINCI file -- if you can share this with us.

Hope this helps.
Sébastien.
The administrator has disabled public write access.

Re: U and V at velocity boundary 13 years 1 week ago #3026

  • olslewfoot
  • olslewfoot's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 132
  • Thank you received: 3
Hi Sebastien

Thank-you - those modifications have worked. The file compiles.

Maybe my fortran file did not upload on the last post - attached again for completeness.


File Attachment:

File Name: scalar_fcoef_2011-11-16.f
File Size: 19 KB


Thanks for your help.
John
The administrator has disabled public write access.

Re: U and V at velocity boundary 13 years 1 week ago #3028

  • sebourban
  • sebourban's Avatar
  • OFFLINE
  • Administrator
  • Principal Scientist
  • Posts: 814
  • Thank you received: 219
Thank you for the PRINCI John. It worked in your previous attempt, it is just that I did not look properly.

Glade it works.
Sébastien.
The administrator has disabled public write access.
Moderators: pham

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