Welcome, Guest
Username: Password: Remember me

TOPIC: How to iterate the boundary nodes for HN without using NBOR

How to iterate the boundary nodes for HN without using NBOR 11 years 8 months ago #8010

  • gourish
  • gourish's Avatar
Hi,

I am trying to iterate through the boundary nodes in the telemac2d function Q for the variable HN, to check for minimum water depth at previous time step.

However it is required that the array NBOR is passed to the function. Currently Q has only IFRLIQ variable is passed to it.

How can I iterate through the boundary HN values without using NBOR.
I thought of doing the following way, which of course requires NBOR.
      HMIN = -1;
      DO K=1,NPTFR
        IF(NUMLIQ(K).EQ.I) THEN
            if( HN%R(NBOR(K)) .LT. HTRESH(I) ) THEN
                  IF( HN%R(NBOR(K)) < HMIN) then
                        HMIN = HN%R(NBOR(K))
                  ENDIF
            ENDIF
        ENDIF
      ENDDO

Is there any other method.

Thanks and regards,
Gourish

File Attachment:

File Name: qmodf.f
File Size: 5 KB
The administrator has disabled public write access.

How to iterate the boundary nodes for HN without using NBOR 11 years 8 months ago #8012

  • riadh
  • riadh's Avatar
Hello Gourish

You can add (before IMPLICIT NONE):
USE BIEF_DEF, ONLY:NBOR

And the code you've written should work.

With my best regards

Riadh
The administrator has disabled public write access.

How to iterate the boundary nodes for HN without using NBOR 11 years 8 months ago #8018

  • jmhervouet
  • jmhervouet's Avatar
Hello,

Actually I am not sure that NBOR has been given an alias, you may reach it with:

USE DECLARATIONS_TELEMAC2D, ONLY : MESH

and then your NBOR(K) would be : MESH%NBOR%I(K)

With best regards,

Jean-Michel Hervouet
The administrator has disabled public write access.

How to iterate the boundary nodes for HN without using NBOR 11 years 8 months ago #8023

  • gourish
  • gourish's Avatar
Hi Jean-Michel Hervouet,

You are correct, I used the method as mentioned by you and it worked perfectly.

Thanks and Best regards,
Gourish
The administrator has disabled public write access.
Moderators: pham

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