^{1}

^{1}

^{1}

^{1}

^{1}

^{1}

^{1}

^{1}

^{2}

^{3}

The version 2.1 of PETROMAR-3D model, created in the Center for Marine Meteorology of the Meteorology Institute of Cuba, is presented. This Lagrangian model has been designed to describe the physical processes of marine oil spills in the face of multiple scenarios of the marine environment. Although it is applicable to any part of the world, it is mainly designed for inter-American seas. The novelty has been to integrate the processes of drift and weathering into a model, with updated methods that incorporate 3D phenomena, a very favorable situation to achieve an operating system in Cuba and the region for the immediate and medium term. Python was chosen as the programming language because it has advanced libraries for numerical modeling, automation work and other useful tools for pre-and post-processing. By means of adapters, an important number of atmospheric, hydrodynamic and wave models have been considered to create the scenarios efficiently. The modular distribution in which the code has been created facilitates its use for other dispersion analysis and biophysical applications. Finally, a set of simple images are presented, aimed at informing decision-makers in order to mitigate the effects of the spill on the environment.

The presence of Exclusive Economic Zone (ZEE as abbreviated in Spanish) of the countries that make up the Gulf of Mexico, together with the transit through its territorial waters of millions of barrels of oil daily, has promoted scientific research on the processes involved during a spill. In this regard, one of the efforts made is the creation of the PETROMAR-3D model, whose system is ready to be available in the region of the inter-American seas.

This system has been used in Cuba for the main practical exercises in the face of oil prospecting by some international companies. Additionally, it was used in the event of contingencies such as the Macondo well, an event that not only made profound damages to environment in several ways, but also generated countless important investigations in the region [

The main antecedents of this work are a national scientific project completed in 2009 and the experience gained from the Gulf Spill in 2010, with own evaluated daily forecasts made for Ministry of Technology and Science functionaries and other national decision-makers. Version 1.0, 1.1 and 1.2 of a drift model, designed and implemented by the authors named PETROMAR [

The main outputs of the last project, performed by the team and finished in November 2017, are presented in this article. It was initially proposed to update and improve the algorithms and methods of previous versions of the model, with the introduction of the weathering processes and description of the blowout, to obtain more real and objective scenarios.

The main update made is related with the three-dimensional conception of this version of the model, expressed in its modules that describe the ascent of the pollutant from the depths to the surface and weathering processes. Other improvements were made to the drift module, so the three modules work jointly and in harmony.

In the calculation of drift and weathering of the slick there are, on the one hand, many empirical functions that are still used with well-defined parameters and, on the other hand, numerical models that apply equations of fluid dynamics. Until now, both tendencies have elements of significance and that is why neither of them should be discarded as a possible solution [

- If the spill occurs on the sea surface, the runoff is calculated according to Fay’s equations [

- At the same time, Evaporation processes (with the PC method), vertical dispersion (with Delvigne and Sweeney’s proposal) and emulsification (through Eley’s law of change) are calculated (from the beginning). All this combination of processes will favor, for the first time, that time based parameters such as density, viscosity, thickness and concentration of the oil slick could be calculated.

- If the spill takes place in the sea depths, everything will be executed in a similar way, except for two main aspects: the spreading of the oil on the surface will be the result of the value of the radius of the plume around its axis. The parameters of the slick’s position and dimension are described in the blowout process. Other aspect is that the values of density and viscosity of the pollutant in the surface movement will be linked to the upward movement, introduced by the method proposed by Zheng and Yapa [

The system will do the above in four stages:

1) Data selection and interpretation.

2) Creation of the scenarios (framed by dates, zones and supplier models, starting with a configuration file that allows the entry of initial parameters, see Annex).

3) Calculation of the trajectory considering the initial coordinates (X, Y, Z), the volume spilled or the rate of spillage, the type of oil and other aerial or satellite information are also initialized within the configuration file (see Annex).

4) Graphic representations and storage of the numerical data of the results are tasks carried out through the Matplotlib and Basemap libraries.

The Lagrangian method to describe the processes in the slick causes the concentration to be calculated indirectly. This process is done considering the amount of existing virtual particles (with each of its macro-scalar features) in a pre-determined control volume (hereafter CV) around each node of the mesh, as well as its macro-scalar properties that are modified by the environment. Most of the models developed in the first decade of the 21^{st} century and even during the first years of the second, present similar ideas. Despite this, there is a novel idea to carry out the reconstruction of the oil concentration field from the distribution process of advection-diffusion and transformation processes of the lagrangian particles [

This model is developed based on the Lagrangian concept proposed by Lee and Cheung [

1) The plume cross-section is circular and perpendicular to the trajectory.

2) The control volume has the shape of an inclined cone.

3) The effects of oil viscosity are not included in the calculations due to the turbulent plume development.

4) Forced penetration of ambient fluid into the plume occurs from its front.

5) The model is of the quasi-stationary type.

1) Governing plume equations

In this work, a Lagrangian element is defined as a control volume moving through the central plume line. Element moves with the speed V of its central line. Thickness and mass are defined by h = V d t and m = ρ π b 2 h , been ρ and b the density and radius of the volume respectively. The following governing equations are applied to the control volume.

a) Conservation of mass equation for a plume

d m d t = ρ a Q e − ∑ i n ( d m i d t − d m d d t ) (1)

In the Equation (1), ρ a is the density of the ambient fluid, Q e = Q f + Q s the volume flow entering the control volume, due to forced penetration Q f and induced shear penetration Q s . The second term is the summation of the differences of the loss of mass due to dissolution and turbulent diffusion.

In practice, it is convenient to define the term Q e as a result of a constant of proportionality E multiplied by an environmental function.

Q e = E ⋅ f ( ρ a , ρ p , v a , v p ) (2)

This constant is empirical and must be adjusted through experiments. The aforementioned function f in (2) is proportional to ambient density ( ρ a ), CV density ( ρ p ), ambient speed ( v a ) and CV speed ( v p ).

b) Momentum conservation

d ( m V ) d t = V a d m d t + m Δ ρ ρ g k − 2 ρ b h C D ( | V | − V ′ a ) 2 | V | | V ′ a | (3)

In the above Equation (3), V a is the average vector of the flow velocity over the exposed plume surface while V ′ a is a projection of V a in the direction of V and k is a unitary vector in the vertical direction. On the other hand, the deficiency of density between oil and the ambient fluid Δ p is presented as the difference Δ p = p a − p and C D is a drag coefficient.

The first term on the right in the equation represents the momentum of the mass penetrating the plume, the last two terms represent the net forces acting on a control volume.

c) Conservation of scalar magnitudes

d ( m I ) d t = I a d m d t − 2 a π ρ b h k I − I a b (4)

In this case, I is a symbol representing any of the scalar magnitudes that characterize the plume, for heat I = c p T , for salinity I = S and for oil mass concentration I = C . Diffusive constants, such as heat and salinity diffusivity or oil concentration, are represented by k.

In the case of concentration C, an additional term is necessary to take into account the dissolution rate

∑ i n d m i d t . (5)

This continuity equation states that scalar changes in the control volume are equal to the scalar magnitude that mainly penetrates through the mass entering the control volume, and the loss of this magnitude due to turbulent diffusion.

d) Equation of state

The set of governing equations needed to model the behavior of a submerged plume can be completed with an equation of state. Its purpose is to describe the changes in density due to salinity, temperature and concentration taking the general form

ρ = ρ ( T , S , C ) . (6)

Exact shape of the above equation depends on the type of fluid in the plume, for example, the density of hot water released into seawater depends on the temperature and salinity of the mixed water.

Model selected includes the main processes affecting plume development and correctly represents the physical behavior of the oil that ascends in a broad spectrum of initial and environmental conditions. Because of its formulation, it ceases to be valid once the vertical velocity is zero or very close to it, known as the terminal level (which may be found below the surface). In these cases, it is necessary to take into account other processes that affect oil, such as horizontal spreading.

e) Implementation details

Two classes were implemented to programmatically encapsulate the methods used in this model. One for the numerical processing of the variables of interest on the control volume; the second was implemented for the visualization of the results providing a wide variety of graphical outputs to facilitate the analysis and verification of the results.

The model was implemented in the Python programming language version 2.7, using the integrated development environment Spyder 3.0 (Scientific Python Development Environment). Python Seawater module [

In order to plot the results, a Python module was implemented where MATPLOTLIB library [

The drift processes are governed by the equations

d x i d t = u a ( x i , t ) + u d ( x i , t ) (7)

and

d y i d t = v a ( y i , t ) + v d ( y i , t ) , (8)

considering the motion of a particle i in each axis. The u a and v a are the components of the speed of the particle, consequence of the advection and in the other hand, u d and v d are those referring to diffusion. Thus, putting equations (6) and (7) in perspective, you would obtain that the advection part could be modelled like

u a ( x i , t ) = c c o r u c o r ( x i , t ) + c v t o u v t o ( x i , t ) (9)

and

v a ( y i , t ) = c c o r v c o r ( y i , t ) + c v t o v v t o ( y i , t ) , (10)

in which u c o r and v c o r are the current components and u v t o and v v t o the wind ones. Coefficients c c o r and c v t o represent the interaction between the current and the wind with the oil slick. They were selected, taking into account the values obtained in tropical laboratories exercises according to Betancourt [

Furthermore, if we incorporate a random walk method, such as the Brownian motion, into these formulations, the position of the particles at any time would be now described by the stochastic differential equations

x i n + 1 = x i n + ( c c o r u c o r ( x i n , t ) + c v t o u v t o ( x i n , t ) ) Δ t + 12 D h Δ t [ R 1 ] 0 1 cos ( 2 π [ R 2 ] 0 1 ) (11)

y i n + 1 = y i n + ( c c o r v c o r ( y i n , t ) + c v t o v v t o ( y i n , t ) ) Δ t + 12 D h Δ t [ R 1 ] 0 1 sin ( 2 π [ R 2 ] 0 1 ) (12)

In these formulations, D h represents the horizontal diffusion coefficient, choosing 10^{5} m^{2}/s for the calculations [_{1} and R_{2} are independent random variables within the [ 0 , 1 ] interval.

The spills described by the system could be instantaneous or continuous (considering the type of source), on the surface or subsurface (taking into account the depth) and in a fixed location.

1) Implementation details

This module has many details from its previous versions. The programming language selected for this version to re-elaborate the codes of the algorithms was Python v2.7.6. This language covers the majority of the entire code and serves as a coordinating element between different scientific libraries used, mostly written in FORTRAN and C, languages designed for high performance. Numpy, Scipy, Matplotlib, Trapezoid Map Finder, Basemap, Cmocean and Seawater are just examples of them.

Precisely, the ability of Python to interact with other languages to optimally solve specific problems, is one of the main arguments that support its choice as the main language. This feature makes the drift module have an efficient implementation.

On the other hand, the fact that Python and the mentioned scientific libraries are supported by multiple platforms, makes this implementation highly subjected to be implemented in them. The current version was implemented for Linux operating system (Ubuntu 14.04 or superior to be more exact) to facilitate the automation processes.

Three mass transferring processes will be explained below. Two of them (evaporation and vertical dispersion) reduce the initial volume of the pollutant; the first transfers it to the atmosphere and the second to the water column underlying the spill. The third, emulsification, occurs when to water column does not enter any more oil molecules (saturation state) resulting in the mixture. Therefore, the latter increases the initial volume and makes the modeling more complex by giving rise to a new substance with very complex characteristics. Details are given below.

1) Evaporation

Evaporation is the most important mass removal process in an oil spill. It is controlled by the mass transfer coefficient, wind speed, oil diffusivity (Schmidt number) and vapor pressure; although other authors also include water temperature.

The authors of this paper choose the pseudo-component (hereafter PC) model to describe this process. Therefore, the hydrocarbon mixture is divided into several fractions from boiling point data. Gros et al., Kotzakoulakis and George, and Spanoudaki et al. [

is constructed from the first-order approximation of d T d f using empirical correlation based on its API.

T 0 = 457.16 − 3.3447 API (13)

d T d f = 1356.7 − 247.36 ln API (14)

First five PCs of equal volume are created and their boiling points are determined by

B P i = T 0 + d T d f i − 1 2 5 . (15)

Vapor pressure of each PC is based on the Antoine equation [

ln P i p 0 = Δ S i ( B P i − C 2 , 1 ) 2 Δ Z b R B P i [ 1 B P i − C 2 , i − 1 T − C 2 , i ] (15)

in which

C 2 , i = 0.19 B P i − 18 (16)

and

Δ S i = 8.75 + 1.987 log B P i . (17)

The relative molar volume V ¯ i of each PC is obtained from a correlation between the molar volume and the boiling point for a series of normal alkanes (C3 to C20).

V ¯ i = 7.0 e − 5 − ( 2.102 e − 7 B P i ) + ( 9.494 e − 7 ( B P i ) 2 ) (18)

Effective molecular weight of each PC is obtained in the same way.

M W i = 0.04132 − ( 1.985 e − 4 B P i ) + ( 9.494 e − 7 ( B P i ) 2 ) (19)

Reference mass transfer coefficient K i 0 is defined at a wind speed of 1 m/s, and a pool diameter of 1 m. A single Schmidt number value is used for all components and is based on the hydrocarbon moles weighted average.

K i 0 = 0.0048 ( 1.3676 ( 0.018 M W a v g ) 1 2 ) (20)

M W a v g = ∑ i = 0 n c p M W i V i ( 0 ) V ∑ i = 0 n c p V i ( 0 ) V (21)

In addition, the exposure evaporation matrix is constructed consisting of four columns: exposure evaporation, volatility factor Z, evaporated fraction F 2 and vapor pressure.

Volatility factor is a measure of how volatile the PC is and is defined by

Z ( t ) = ∑ i = 0 n c p P i V t R T N ( t ) . (22)

The evaporation exposure values are selected in such a way that the separation between them increases with the index of each PC. The minimum is set to zero, and the maximum was obtained using the maximum time and maximum wind speed.

ϕ max = t max K 0 U max 7 9 (23)

The incremental change in evaporation due to exposure also increases with the index of each PC as follows.

Δ ϕ ( j ) = 2 ϕ max j j max 2 (24)

The differential equations describing evaporation are then solved in these preset values of evaporation due to exposure.

d V i d ϕ = − V ¯ ( ϕ ) ( P i R T ) V i ( ϕ ) (25)

A numerical approximation in terms of ln V is used:

V i ( ϕ i + 1 ) = V i ( ϕ i ) e − V ¯ ( ϕ i ) P i R T ( ϕ i + 1 − ϕ i ) . (26)

Once the PC volumes are at certain discrete values of evaporation due to exposure, the volatility factor, the evaporated fraction and the vapor pressure of each PC are easily calculated by the following ratios.

● Volatility factor:

Z ( ϕ ) = ∑ i = 0 n c p P i V i ( ϕ ) R T N ( ϕ ) . (27)

● Evaporative fraction:

F z ( ϕ ) = ∑ i = 0 n c p V i V 0 . (28)

● Vapor pressure is the sum of the partial pressure of the individual components such that

P T = ∑ i = 0 n c p n i P i ∑ i = 0 n c p n i . (29)

The last one is calculated from the amount of hydrocarbon that evaporates over time, using the previously created exposure evaporation so

d V V d t = − { K 0 ( 1 − Y ( t ) ) U 7 9 h ( t ) } Z ( ϕ ) . (30)

The integration time was carried out as follows:

V ( t n + 1 ) = V ( t n ) e [ − { K 0 ( 1 − Y ( t n ) ) U J 9 h ( t n ) } Z ( ϕ ( t n ) ) ( t n + 1 − t n ) ] . (31)

2) Vertical dispersion

Natural dispersion of oil in a spill is understood as the process of fragmentation of oil into droplets due to the breaking of the waves on the slick. According to Delvigne and Sweeney [

Q = C d i s p D e 0.57 f b w V d i s p , (32)

where C d i s p is an experimentally determined parameter, f b w is the fraction of breaking wave per wave period, and V d i s p is the volume of oil introduced into the column per unit volume of water.

The energy dissipation of the wave per unit area is defined by

D e = 0.0034 ρ w g H r m s 2 , (33)

in which H r m s 2 is the mean root of the wave height in meters, assuming the significant spectral height of wave h 0 as H r m s = 0.707 H 0 .

V d i s p is obtained by integrating the product of the droplet volume and the frequency distribution of the droplets on the oil volume. In practice, the integration is carried out between the minimum and maximum size that the droplet could reach, considering the experimental analyses.

V d i s p ∝ ∫ d min d max N ( δ ) δ 3 d δ (34)

Larger droplets, whose imposed maximum allowable size is 70 microns, will refloat in less time than necessary for the stain to pass through the area covered by the dispersed oil, thus reintegrating into the oil slick (forming the tail). Droplets 70 microns or less in size remain in suspension for a time determined by the buoyancy and depth of penetration into the water column. The minimum droplet size is generally considered to be 5 microns, from which the droplets are considered to be precipitated to the seabed [

N ( δ ) is the number of droplets of oil per unit volume of water per droplet size measured by

N ( δ ) = N 0 ( δ 0 δ ) 2 3 (35)

where N 0 and δ 0 are experimental reference values.

3) Emulsification

Wave energy spontaneously mixes hydrocarbon with seawater, forming a “water-in-oil” emulsion. In a series of laboratory experiments with crude oils [

d S d t = k s ( 1 − S S max ) , (36)

where the interfacial parameter k s is sensitive to wave energy; S and S max are the oil-water interfacial area and the maximum interfacial area, respectively. The water fraction contained in the oil slick is related to the interfacial area and to the mean diameter of the water droplets d w such that

Y = S d w 6 + S d w . (37)

4) Change of physical-chemical properties

Evaporation and emulsification increase the density and viscosity of the slick, even when the oil is freshly spilled. Most hydrocarbons and their by-products are more viscous than water. Emulsified hydrocarbons may have a viscosity of several orders of magnitude greater than the corresponding fresh hydrocarbon and often behave like a non-Newtonian fluid. Density will also increase due to the weather, but not to the point of losing the buoyancy of the oil.

5) Density calculation

The two circumstances in which slick sinking has been recorded are when the slick mixes with the mineral sediment in the water column or when it burns in-situ, generating a residue that may be of high density. The weathering algorithms used will not identify any hydrocarbon that is denser than the surrounding water, as these will simply not be modeled.

A common way to present the density of a liquid is to compare it with the density of fresh water, giving a specific gravity of a fluid in relation to the water. However, in the oil industry, the common system of units is the API. The two references are related by the formula:

G e = 141.5 131.5 + API . (38)

Three factors affecting the slick density are considered: water temperature, degree of emulsification, and evaporated hydrocarbon fraction. The formula used that estimates this parameter in comparison with the density of fresh oil is

ρ = Y ρ w + ( 1 − Y ) ρ r e f [ 1 − c 1 ( T − T r e f ) ] ( 1 + c 2 f e v a p ) , (39)

where c 1 and c 2 are empirical constants, ρ is the hydrocarbon density, ρ w the water density, ρ r e f the reference hydrocarbon density, Y the emulsion fraction, T the water temperature, T r e f the reference temperature and f e v a p the evaporated hydrocarbon fraction.

6) Calculation of viscosity

Viscosity is the property that determines how freely a fluid flows and may be physically regarded as the internal friction of the fluid. Viscosity is significantly affected by temperature and increases with the time the oil has been outdoors. For a Newtonian fluid it is defined by the relationship between stress and shear rate. For a dimensional flow in the X direction, the τ stress is proportional to the velocity gradient.

τ = − η d v x d y (40)

Taking in account the mentioned temperature sensibility of viscosity, Andrade correlation is generally valid for small temperature variations [

ln η 2 η 1 = C T ( 1 T 2 − 1 T 1 ) (41)

The proportionality constant C T is taken with the value 5000 K.

As oil evaporates, the relative concentrations of the various components change. Perry’s rule [

ln η max = ∑ i = 0 n c p w i ln η i (42)

The following relationship given by Mackay et al. [

ln η η 0 = C B F z (43)

The empirical proportionality constant C B = 5 is used for this research.

Mooney’s equation [

η η 0 = e ( 2.5 Y 1 − C Y Y ) (44)

Nguyen [

C η ( t ) = d min Y max d ( t ) ( 1 − 2.5 Y ln R η ) (45)

where d ( t ) is the diameter of droplets based on time, d min the minimum diameter of droplets, F z ( t ) the fraction of evaporated oil based on time, R n the maximum emulsion radius, T the temperature (K), v the fluid velocity, Y ( t ) the fraction of water content based on time, Y max the fraction of maximum water content, n max the maximum dynamic viscosity and n ( t ) the dynamic viscosity based on time.

7) Module implementation

The libraries “matplotlib”, “basemap” and “cmocean” were used, mainly in the preparation of the new graphic outputs.

Weather data used in the creation of the scenarios come from the following atmospheric models.

● WRF: Weather Research and Forecasting [

● MM5: Mesoscale Model 5 [

● NAM-CARIBE: North American Mesoscale-Caribbean Forecast System, downloaded from Internet. The domain considers the entire Greater Caribbean and presents a spatial resolution of 12 km and a time of 3 hours [

Among the models and projects that provide the variables of marine circulation, those mentioned below were chosen.

● ROMS: Regional Oceanic Modeling System [

● HYCOM: Hybrid Coordinate Ocean Model, downloaded from Internet. The horizontal grid has an average space between nodes of 7 km and 32 vertical layers [

● NCOM: Navy Coastal Ocean Model, downloaded from the Internet, is the U.S. Navy’s worldwide operating ocean model that includes several domains within which the Gulf of Mexico and the Caribbean Sea [

● CW: COSTWATCH, downloaded from the Internet) this project does not yields predictive scenarios, just diagnostic. However, its values have been corrected for their drifter service, improving the quality of the data [

In order to consider the parameters of the waves, the outputs of the following models were incorporated.

● WW III: WAVEWATCH III [

● SWAN: Simulating Waves Nearshores [

The bathymetric information for the Gulf of Mexico, Caribbean Sea and Near Atlantic, referring to bathymetric and coastline data, was obtained by integrating two sources. In first place The GEBCO project [

Characteristics of the basemap were conformed creating the regions and other components in the following way: scale 1:250,000, WGS84 projection, latitude range from 7.20 degrees to 33 degrees latitude north and longitude range from 61 degrees to 101 degrees longitude west. The model has twelve pre-conceived zones related to the Interamerican Seas, spatially designed by finite elements [

Simple ways to illustrate the graphical outputs results is one of the advantages of PETROMAR-3D model. The purpose is to give as accurate as possible forecasts to facilitate decision-makers quick and effective analyses and responses.

Images describing the upward movement are illustrated in

Outputs were obtained taking the value of the constant E equal to 0.01, as the best fit for the northwestern area of Cuban territorial waters.

Horizontal displacement is one of the most important variables for those responsible for taking action in the event of disasters. Its variations for different values of parameter E (

In dotted lines are shown the results of the ASCENSOMAR model, proposed by Pérez [

Additionally, the new curves reflect, to a better extent, the behavior of the particle in its ascent, mainly the fact of having the concavity downwards (the slope to the curve decreases with height) which is theoretically more accepted. The displacement in the near spill is very limited because the exit speed is imposed and then gradually reduced to a limit speed.

Once the calculations of each weathering process are attributed to the virtual

particles of the Lagrangian method, the oil concentration values will be obtained at all times considered in the scenario as a sample of the existing intermodular relationship.

Numerical and graphical outputs are yielded every one hour giving information of state and position of all particles conforming the modelled oil slick which are been transported. The first ones are contained in a file with the variables of each run: initial time, final time and time of greatest arrival of the pollutant to a coastal sector, as well as location of virtual particles with their macroscale variables every one hour. This kind of output becomes very convenient because it can be used as input for other analysis and post-processing tools like geographic information systems to obtain, for example, the affected coastal sector in kilometers. In addition, the model can suggest the type of optimal response technique to be used over time to mitigate the environmental impact of the spill.

Graphical outputs are released to illustrate not only all the numerical releases but also other useful information like cartographical maps and meshes configured on PETROMAR-3D.

Figures below are results of a 5000 bb spill from an average oil of 30 API. They are very important for users because of their integration of useful information in a simple way. This feature allows them to take quick actions with respect to the incident.

The version 2.1 of PETROMAR-3D model has a module that simulates a blowout in the sea depths. The enhancement of this process includes the characteristics of the type of oil, marine environment and the spill itself. With the introduction of weathering, the spill is better represented and the rest of the processes are simulated in a more realistic way. With these conditions it is possible to represent the concentration of the slick based on time, favoring a better analysis of the scenario.

Given the choice of the mentioned theory to describe the physical processes, the mathematical methods and the numerical schemes, the proposed model is considered to have advantages to be used in the simulation of spills in the inter-American seas. So the novelty of the work is directed to three fundamental aspects. First, the design of the region’s areas has been refined with great geographic details and the outputs of the main existing operating models have been considered. The drift, weathering and blowout processes have been unified in the same model, achieving similarity with others of advanced performance. Finally, the graphical and numerical outputs have been adjusted to the demands of decision-makers present in the region and there is still a list of suggestions that authors will consider.

It is recommended to include the influence of waves on surface drift through Stokes’s law [

The authors declare no conflicts of interest regarding the publication of this paper.

Calzada, A.E., Delgado, I., Ramos, C., Pérez, F., Reyes, D., Carracedo, D., Rodríguez, A., Chang, D., Cabrales, J. and Lobaina, A. (2021) Lagrangian Model PETROMAR-3D to Describe Complex Processes in Marine Oil Spills. Open Journal of Marine Science, 11, 17-40. https://doi.org/10.4236/ojms.2021.111002

– Configuration file of PETROMAR-3D

### AUTHOR: CENTRE FOR MARINE METEOROLOGY ###

General:

model_data: /home/petromar/PetroData

output_data: /home/petromar/PetrOut

#############################

### SCENARIO CREATION ###

### INPUT VARIABLES ###

#############################

Scenario:

name_scenario: Macondo_z16_ROMS_Py

initial_date: 22/04/2010

number_days: 30

### DESIGNED ZONES ###

###Golfo_zon; zona01; zona02; zona03; zona04; zona05; zona06; zona07; zona08

###zona09; golfo_zona01; golfo_zona01_06; zona01_02_03; zona01_02_03_06_11; zona01_06; zona03_04; zona04_05_06

###zona07_08; zona08_09

zona_s: zona01_06

### MODELS TO USE ###

### MARINE CURRENT MODELS: Curr ###

### hycom1; hycom2; coastwatch; roms; ncom

curr: roms

### ATMOSPHERE WIND MODELS: wd ###

### mm5; gfs; wrf

wd: 'wrf'

### WAVE MODELS: wv ###

### ww3; swan

wv: swan

### TYPE OF RUNS (tr): FORWARD (1) or BACKWARD (−1)

tr: 1

##############################

### TRAJECTORY CALCULATION ###

### INPUT VARIABLES ###

##############################

Model:

name_scenario: Macondo_z16_ROMS_PY

name_track: t1_rk4

### Runge-Kutta numerical schemes ###

### rk1; rk2; rk4

numeric_schema: rk4

### longitude(geog. degrees): x

### latitude(geog. degrees): y

x: −83.49

y: 23.31

### Kind of the spill (ks): “instantaneous” or “continuous”

ks: instantaneous

initial_time: 10

### Volume magnitude (mag_vol): “m3”, “liters”, “barrels”

mag_vol: barrels

Volume: 500

### Density (Ro): “density”, “specific gravity: ge”,

### “American Petroleum Institute: api”

Ro: api

Value: 20

Oilvis (cSt): 259.0

Oilvisreftem (C): 25.0

#Sea water density (kg/m^{3}): Roa

Roa: 1027

#Water cinematic viscosity (m^{2}/s): visc

visc: 140.267

#Superficial tension coeficient: coef_tension

coef_tension: 27

###############################

### Subsurface spills ###

###############################

salinity_p: 27

temperature_p: 35

depth_p: 100.0

density_p: 893

density15_p: 892

radio_v: 0.4016

### Release speed of blowout ###

v_p(m/s): 1

u_p(m/s): 1

w_p(m/s): 2.1

final_time: 700

constante_E: 0.5

######################## PETROMAR-3D END #######################