Welcome, Guest
Username: Password: Remember me

TOPIC: Oilspill

Oilspill 10 years 5 months ago #13263

  • smk.svks
  • smk.svks's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 93
  • Thank you received: 1
Dear sir,
The examples given for the oil spill modules are for instantaneous oil spills, how can i simulate a continuous oil spill???

-thanks,
Sathish
The administrator has disabled public write access.

Oilspill 10 years 5 months ago #13268

  • gourish
  • gourish's Avatar
Hi,

In the given example, there is a fortran file which modifies the OIL_FLOT subroutine to give the time and location of the float. I think you may have to modify the subroutine to give a range of time-steps where you want to simulate the continuous oils spill, and release the drogues at these particular time steps.

Correspondingly the NUMBER OF DROGUES should be equal to the "total number of drogues spilled per time step" x "Total number of time step for spill".
I have not tried this method, but possibly I think may be a solution of getting a continuous spill.

Check the document www.opentelemac.org/index.php/manuals/su...looms-and-oil-spills , page no. 53 for more details.

Hope this helps,
With Best regards,
Gourish
The administrator has disabled public write access.

Oilspill 10 years 5 months ago #13274

  • goeury
  • goeury's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 16
  • Thank you received: 6
Dear Sathish,

If you want to model a continuous oil spill, you need to modify the subroutine OIL_FLOT.f to make a release of some particles at some time steps. Gourish gives you the good response. So, I give you a little example.
If you want to make a release of 10 particles every 10 times steps until the time step number 50, you need to write in OIL_FLOT.f :

==============================================================================
00156 IF(LT.EQ.10.OR.LT.EQ.20.OR.LT.EQ.30.OR.LT.EQ.40.OR.LT.EQ.50) THEN
00157
00158 NUM_MAX=0
00159 NUM_LOC=0
00160 COORD_X=0.D0
00161 COORD_Y=0.D0
00162 NUM_MAX=10
00163 DO K=0,NUM_MAX-1
00164 DO J=0,NUM_MAX-1
00165 COORD_X=336000.D0+REAL(j)
00166 COORD_Y=371000.D0+REAL(k)
00167 NUM_GLO=NUM_GLO+1
00168 NFLOT_OIL = 0
00169 CALL ADD_PARTICLE(COORD_X,COORD_Y,0.D0,NUM_GLO,NFLOT_OIL,
00170 & 1,XFLOT,YFLOT,YFLOT,TAGFLO,SHPFLO,SHPFLO,ELTFLO,
00171 & ELTFLO,MESH,1,0.D0,0.D0,0.D0,0.D0,0,0)
00172 IF(NFLOT_OIL.EQ.1)THEN
00173 NUM_LOC = NUM_LOC+1
00174 !=========================================================================
00175 !----INITIALIZATION PARAMETERS FOR THE CALCULATION OF PARTICULE MOTION----
00176 !=========================================================================
00177 PARTICULES(NUM_LOC)%XOIL = XFLOT(1)
00178 PARTICULES(NUM_LOC)%YOIL = YFLOT(1)
00179 PARTICULES(NUM_LOC)%ID = TAGFLO(1)
00180 PARTICULES(NUM_LOC)%SHPOIL(1) = SHPFLO(1,1)
00181 PARTICULES(NUM_LOC)%SHPOIL(2) = SHPFLO(2,1)
00182 PARTICULES(NUM_LOC)%SHPOIL(3) = SHPFLO(3,1)
00183 PARTICULES(NUM_LOC)%ELTOIL = ELTFLO(1)
.....
==========================================================================

Best regards,
Cedric
The administrator has disabled public write access.
Moderators: pham

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