Welcome, Guest
Username: Password: Remember me

TOPIC: drogue file output seems very short

drogue file output seems very short 4 years 11 months ago #35085

  • shivanesh.rao
  • shivanesh.rao's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 37
  • Thank you received: 1
Hi all,

I'm stick at an odd place with my telemac2d & tomawac coupled model.

I have a fortran file in tomwac.cas that does time-varying (spatial uniform) wave heights along the boundary. I have another fortran file in telemac2d.cas that adds drogues.

When I run the simulation, the drogue output file has only a 3 line output (attached). But I specified:
# of drogues =200
printout period for drogues=10
drogue file="drogues_output'
stochastic diffusion model=1

The flot.f subroutine in my fortran file has
IF (LT.EQ.1) THEN
CALL ADD_PARTICLE(355729.125D0,6298045.500D0,0.D0,LT,NFLOT,...)
ENDIF

!!NO CALL DEL_PARTICLE(...)

I added WRITE(*,*) "IN LOOP" into the IF loop and that shows up in log, so it is going through the loop.

Should I be expecting 200 lines in the drogue output file (even if the X,Y is fixed)? I'm a bit confused on what to expect from the tec file?

Thanks,
Shiv
PS: Unfortunately, I can't plot this atm because the tecplot eval product has network errors to activate product...maybe the company firewall.
The administrator has disabled public write access.

drogue file output seems very short 4 years 11 months ago #35094

  • shivanesh.rao
  • shivanesh.rao's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 37
  • Thank you received: 1
A bit more information...

1. Drogue seems to jump to the boundary of domain after 1st time-step (~10-15km)
2. The subsequent drogue position (LT>=2) don't appear to change with time, i.e., XFLOT and YFLOT
3. The drogue output tecplot file has only 1 line of data.

The telemac2d.cas file introduces the drogues via:
NUMBER OF DROGUES=1
PRINTOUT PERIOD FOR DROGUES=1
DROGUES FILE='drogues_output.tec'
STOCHASTIC DIFFUSION MODEL=1

The fortran file (for telemac2d) has the flot.f (from same version of telemac2d) subroutine that adds the particle via:
      IF(LT.EQ.1) THEN
       WRITE(*,*) 'IN THE ADD LOOP'
        CALL ADD_PARTICLE(355729.125D0,6298045.5D0,0.0D0,1,NFLOT,
     &                    NFLOT_MAX,XFLOT,YFLOT,YFLOT,TAGFLO,
     &                    SHPFLO,SHPFLO,ELTFLO,ELTFLO,MESH,1,
     &                    0.D0,0.D0,0.D0,0.D0,0,0)      
      ENDIF
      WRITE(*,*) 'DROGUE COORDINATES',XFLOT,YFLOT
!

Then the model runs...but the XFLOT and YFLOT values do not change in the WRITE to the log file. In fact, the first XFLOT,YFLOT is the value I assign, but the next time-step, the drogue position has moved to the other end of the domain and remains at that position for the rest of the simulation.
IN THE ADD LOOP
DROGUE COORDINATES 355729.12500000000 6298045.5000000000
NUMBER OF SUB-ITERATIONS : 29
...
...
NUMBER OF SUB-ITERATIONS : 20
WIND AT T= 20.000000000000000 UAIR= 0.13333333333333333 VAIR= -0.10000000000000001
DROGUE COORDINATES 339392.97361245431 6316790.5233478788
NUMBER OF SUB-ITERATIONS : 29
...
...
NUMBER OF SUB-ITERATIONS : 20
WIND AT T= 30.000000000000000 UAIR= 0.20000000000000001 VAIR= -0.15000000000000002
DROGUE COORDINATES 339392.97361245431 6316790.5233478788
NUMBER OF SUB-ITERATIONS : 29
...
...
NUMBER OF SUB-ITERATIONS : 27
Attachments:
The administrator has disabled public write access.
Moderators: pham

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