open TELEMAC-MASCARET The mathematically superior suite of solvers

The model crushed when I tried to use Velocity Profile Option 5.

The erro message is "Fortran runtime error: Array reference out of bounds for array 'zmin', lower bound of dimension 1 exceeded (0 < 1)"

The crush pointed to Line 141 in bord.f

DO K=1,NPTFR
IFR=NUMLIQ(K)
ZMIN(IFR)=MIN(ZMIN(IFR),ZF(NBOR(K))+H%R(NBOR(K)))
ENDDO

I had a look at bord3d.f and it has a check statement for the array bound.

DO K=1,NPTFR2
IFRLIQ=NUMLIQ%I(K)
IPOIN2=NBOR2%I(K)
IF(IFRLIQ.NE.0) THEN
ZMIN(IFRLIQ)=MIN(ZMIN(IFRLIQ),ZF%R(IPOIN2)+H%R(IPOIN2))
ENDIF
ENDDO

should this check statement also be incorporated in bord.f? Thanks.

Edwin

Hello,

You are perfectly right (and we appreciate your remarks...). It will be done in version 6.3 (release in a few days). We just missed a test case to see that because we compile here with bound checking on Intel compiler. Note that this is very useful for development but maybe not efficient for applications. We sometimes observe a factor 3 between Intel with all checking (bounds, traceback) and Intel without.

Answer on radiat.f in Tomawac still pending...

With best regards,

Jean-Michel Hervouet