Welcome, Guest
Username: Password: Remember me
  • Page:
  • 1
  • 2

TOPIC: More than 99 liquid boundaries

More than 99 liquid boundaries 2 years 7 months ago #40190

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
AS you run in //, If the boundary is not in the master node, it'doesn't appear in the listing file I guess...
For the verification, try to run in sequential

Regards
Christophe
The administrator has disabled public write access.

More than 99 liquid boundaries 2 years 7 months ago #40191

  • TomS
  • TomS's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 71
  • Thank you received: 8
Unfortunately, running in serial did not resolve the issue.

I have a potential workaround that seems to be successful. As it is only my inlet flowrates that vary with time and the scalars (salinity/temperature) are constant I have set only the flowrates in the liquid boundaries file. I then include the inlet salinity and temperature values in the actual case file as prescribed values.

This appears to work fine for 100 boundaries so I will test with more.

Thanks,

Tom
The administrator has disabled public write access.

More than 99 liquid boundaries 2 years 7 months ago #40200

  • pham
  • pham's Avatar
  • OFFLINE
  • Administrator
  • Posts: 1559
  • Thank you received: 602
Hello Tom,

In your log_I3_IRANK3.txt, it seems that the error occurs in TR3 function at line 87. I cannot see your whole tr3.f file but does it correspond to
WRITE(FCT(IRANK:IRANK+1),FMT='(I2)') I
just after line
ELSEIF(I.LT.1000) THEN ?

Can you try the enclosed file please and tell if your issue still occurs?

File Attachment:

File Name: tr3_XXX_liq_bnd.f
File Size: 4 KB


I have changed the line by:
WRITE(FCT(IRANK:IRANK+2),FMT='(I3)') I
in addition to
IRANK=IRANK+3

Be careful if you have more than 100 liquid boundaries and more than 10 tracers, you may have to change the length of FCT from 9 to 10 and also change other lengthes in other subroutines/functions.

Hope this helps,

Chi-Tuan

PS: you should avoid using tabulations in your ASCII files. You may encounter problems with TELEMAC.
The administrator has disabled public write access.

More than 99 liquid boundaries 2 years 7 months ago #40201

  • TomS
  • TomS's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 71
  • Thank you received: 8
I appreciate you doing this, Chi-Tuan.

I am still having issues when using the modified tr3.f code that you sent. Examples of the fluxes produced are:

FLUX BOUNDARY 81 : 4.689003
FLUX BOUNDARY 82 : 3.906000
FLUX BOUNDARY 83 : 4.427991
FLUX BOUNDARY 84 : 2.672994
FLUX BOUNDARY 85 : 2.466003
FLUX BOUNDARY 86 : 4761827.
FLUX BOUNDARY 87 : -0.000000
FLUX BOUNDARY 88 : -0.000000
FLUX BOUNDARY 89 : -0.000000
FLUX BOUNDARY 90 : -0.000000
FLUX BOUNDARY 91 : -0.000000
FLUX BOUNDARY 92 : -0.000000
FLUX BOUNDARY 93 : -0.000000
FLUX BOUNDARY 94 : -0.000000
FLUX BOUNDARY 95 : -0.000000
FLUX BOUNDARY 96 : -0.000000
FLUX BOUNDARY 97 : -0.000000
FLUX BOUNDARY 98 : -0.000000
FLUX BOUNDARY 99 : -0.000000
FLUX BOUNDARY 100 : -0.000000
FLUX BOUNDARY 101 : -0.000000
FLUX BOUNDARY 102 : -0.000000
FLUX BOUNDARY 103 : -0.000000
FLUX BOUNDARY 104 : -0.000000
FLUX BOUNDARY 105 : -0.000000
FLUX BOUNDARY 106 : -0.000000
FLUX BOUNDARY 107 : -0.000000
FLUX BOUNDARY 108 : -0.000000
FLUX BOUNDARY 109 : -0.000000
FLUX BOUNDARY 110 : -0.000000
FLUX BOUNDARY 111 : -0.000000
FLUX BOUNDARY 112 : -0.000000
FLUX BOUNDARY 113 : -0.000000
FLUX BOUNDARY 114 : -0.000000
FLUX BOUNDARY 115 : -0.000000

Whereas, using the workaround with only flowrates in the liquid boundaries files and prescribed scalars in the case file produces the correct fluxes:

FLUX BOUNDARY 80 : 7.018711
FLUX BOUNDARY 81 : 6.195022
FLUX BOUNDARY 82 : 6.195022
FLUX BOUNDARY 83 : 59.34480
FLUX BOUNDARY 84 : 2.773822
FLUX BOUNDARY 85 : 19.08324
FLUX BOUNDARY 86 : 0.2828625E+08
FLUX BOUNDARY 87 : 4.976000
FLUX BOUNDARY 88 : 7.662133
FLUX BOUNDARY 89 : 4.622444
FLUX BOUNDARY 90 : 3.942133
FLUX BOUNDARY 91 : 11.35751
FLUX BOUNDARY 92 : 2.510756
FLUX BOUNDARY 93 : 3.471067
FLUX BOUNDARY 94 : 16.86276
FLUX BOUNDARY 95 : 9.140267
FLUX BOUNDARY 96 : 2.823067
FLUX BOUNDARY 97 : 4.351067
FLUX BOUNDARY 98 : 7.628889
FLUX BOUNDARY 99 : 14.49413
FLUX BOUNDARY 100 : 4.638444
FLUX BOUNDARY 101 : 19.80791
FLUX BOUNDARY 102 : 11.23533
FLUX BOUNDARY 103 : 5.077822
FLUX BOUNDARY 104 : 18.47409
FLUX BOUNDARY 105 : 5.374756
FLUX BOUNDARY 106 : 5.919067
FLUX BOUNDARY 107 : 4.168311
FLUX BOUNDARY 108 : 3.472000
FLUX BOUNDARY 109 : 3.935067
FLUX BOUNDARY 110 : 2.375378
FLUX BOUNDARY 111 : 2.192311
FLUX BOUNDARY 112 : 4.976000
FLUX BOUNDARY 113 : 7.662133
FLUX BOUNDARY 114 : 4.622444
FLUX BOUNDARY 115 : 3.942133

I will persist with this method and, if necessary, revisit the tr3.f issues in the future.

Thanks,

Tom
The administrator has disabled public write access.
  • Page:
  • 1
  • 2
Moderators: pham

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