The experimental work was developed in a greenhouse with 8640 m² of covered area, belonging to a farm dedicated to cut rose production (5º03'22.16''N, 73º55'49.07''W, 2583 masl) located in Cogua Municipality (Cundinamarca). The evaluated greenhouse consisted of 24 spans each with a width of 6.0 m (Figure 1A), the minimum and maximum heights under the channel were 3.2 and 6.4 m, respectively. The longitudinal distance of the greenhouse was 60 m and was oriented in the northeast-southwest direction (NE-SW). Each span had a moveable roof ventilation opening 1.0 m wide and 48 m long and fixed lateral openings 1.9 m wide on each of the four sides of the structure, for a total natural ventilation area of 1927 m².

The exterior meteorological conditions were evaluated for 24 hours a day, within the period between 00:00 hours on March 29 and 23:00 hours on May 2, 2016. The variables were recorded through of a weather station (Vantage Pro2 Plus, Davis Instruments, Hayward CA), which integrated sensors for global radiation, temperature, relative humidity, precipitation, wind speed and direction with a registration frequency of 10 minutes.

The CFD numerical simulation allows solving the governing equations of fluid flow using the finite volume method (Molina et al., 2009). These equations can be represented as convection-diffusion equations of a fluid for three conservation laws, which include the moment, energy and transport equations of a compressible fluid and in a three-dimensional (3D) field and are expressed as follows:

where ρ is the fluid density (kg m^-3), ∇ is the nabla operator, φ represents the concentration variable, $\overrightarrow{v}$ is the velocity vector (m s^-1), Γ is the diffusion coefficient (m² s^-1) and S represents the source term (Piscia et al., 2012). The turbulent nature of the airflow was simulated using the standard turbulence model k-ε, which is based on two main equations, one for k representing kinetic energy and the other for ε representing the rate of dissipation in time and volume units. The transport equations for k and ε can be modeled as:

where μ is the viscosity and ${\mu }_t$ is turbulent viscosity in (kg m^-1s^-1), ${\sigma }_k \text{and} {\sigma }_{ϵ}$ are the turbulent Prandtl numbers of Prandtl for k and ε. G _k is the generation of turbulent kinetic energy due to velocity gradients, G _b is the generation of turbulent kinetic energy due to buoyancy, Y _M is the fluctuating dilation in turbulence due to the global dissipation rate and v is the kinematic viscosity coefficient. C _1ε , C _2ε , C _μ , σ _k and σ _ε are constants with default values. This standard model k-ε has been widely used and validated in studies focused on greenhouses demonstrating an adequate precision (Fatnassi et al., 2006; Katsoulas et al., 2006). All the simulations considered the energy equation, which allowed analyzing the scalar field of temperatures inside the greenhouse. Likewise, the phenomenon of buoyancy was modeled. A phenomenon that drives the movement of air under the condition of external speed was made by means of the Boussinesq approach, which is described as follows:

where ρ _ref is the constant density of the flow; T _ref is the actual temperature (°C) and β is the coefficient of thermal expansion for air (β = 0.00329°C^-1). The Boussinesq approximation is valid if the temperature differences that appear in the computational domain are not too large β(T-T _ref ) < 1, a situation that generally occurs in the microclimate study of greenhouse (Baeza et al., 2006).

The selected radiation model was that of discrete ordinates (DO) with angular discretization. This model allows the analysis of climate under nighttime conditions, simulating and resolving the radiation phenomenon from the greenhouse floor to the external environment. For this purpose, the sky was considered a black body with an equivalent temperature (T _C ) of 0°C for a predominant scenario of wet and clear nights (Iglesias et al., 2009).

The ANSYS ICEM CFD preprocessing software (v. 17.0) was used to generate a large computational domain composed of the greenhouse and its surroundings. In order to guarantee the no affectation of the numerical solution of the flow field outside the greenhouse and to allow an appropriate definition of the atmospheric boundary layer (Rico, 2011) The dimensions of the computational domain were 354, 260 and 70 m for the x, y and z axes, respectively (Figure 1C). This size was determined following the guidelines for calculating the wind environment around the buildings using CFD (Tominaga et al., 2008). The computational domain consisted of an unstructured mesh of 76.651.871 volumes discretized in space, this number of elements was obtained after verifying the independence of the numerical solutions of the airflow at mesh sizes with a higher and lower number of elements according to the procedure reported by He et al. (2017). A fundamental criterion to establish the accuracy of the solutions is to evaluate the quality of the mesh (Figure 1B). The quality parameters evaluated were the cell size and the cell-to-cell size variation, finding that 93,1% of the cells in the mesh were within the high-quality range (0,95-1). The orthogonal quality was also evaluated, where the minimum value obtained was 0.94, results that were classified within the high-quality range (Flores et al., 2015). The convergence criteria of the model were established in 10^-8 for the energy equation and 10^-6 for the continuity, momentum and turbulence equations (Baxevanou et al., 2017).

The upper limit of the domain and the surfaces parallel to the flow were set with boundary conditions of symmetrical properties so as not to generate frictional losses of the airflow in contact with these surfaces. The simulations considered the atmospheric characteristics of Cogua Municipality such as atmospheric pressure of 74992 Pa, air viscosity equal to 1,7E-05 kg m^-1 s^-1 and gravity of 9.81 m s^-2. Other properties for polyethylene and agricultural soil such as specific heat (Cp), thermal conductivity (k) and density (ρ) were established according to Villagrán et al. (2012). The perimeter boundaries of the computational domain were established as limits of air intake or pressure output as the case may be. A uniform profile of wind speed was considered, evaluating speeds with values between 0.001 and 1.21 m s^-1 and average temperature values for each hour of the day. These values were established based on the climatic information collected during the experimental measurement period.

Other input parameters that are necessary to feed the radiation model are the optical properties of the roofing material such as, absorption coefficient (α) of 0,69, transmission coefficient (Ƭ) of 0,19 and reflection coefficient (ρ) of 0,11. The model did not include any crop and a maximum hermetic seali ng of the greenhouse was assumed. These simplifications are valid since they are applied to each of the simulated cases and the errors that may arise will have the same degree of magnitude for each scenario.

The evaluation was carried out assuming the average meteorological conditions for each hour evaluated (Table 1) and the standard ventilation configuration of the greenhouse studied, which, for this case, was lateral ventilation combined with roof ventilation for daytime hours (6-18 hours) and roof ventilation for nighttime conditions (19-5 hours).

The validation of the CFD model was made through the comparison between the results of the simulations and the experimental records for the same temporal frequency of the temperature inside the real greenhouse during the evaluation period. The temperature of the greenhouse air was recorded by means of thirty T-type thermocouples (copper-constantan) connected to an equal number of data recorders (Cox-Tracer Junior, Escort DLS, Edison, NJ) deployed uniformly at 1,5 m from the ground surface and recording data every 10 minutes.

The model validation was made by comparing the temperature values measured experimentally and those obtained through CFD simulations, for the 24 hours of the day. Figure 2 shows the behavior of the temperature for each evaluated hour, showing a good fit between the two methodologies used and the same trend for the measurement period. Quantitatively, when comparing the measured and simulated data, an average absolute percentage error (MAPE) of 4,26% was obtained, an absolute mean error (MAE) of 0,75°C and a mean square error (MSE) of 0,94°C, values that guarantee that the use of the CFD model is adequate for the evaluation and description of the thermal behavior of the greenhouse.

Table 2 shows the time evolution of the interior greenhouse temperature for the daytime period. It is observed that the variable presented a minimum value of 12,17°C for hour 6 and reached a maximum value of 24,31°C at hour 14, this temporary increase is driven by the solar radiation. The thermal difference between the interior and the exterior of the greenhouse (ΔT = T _internal -T _external ) exhibited a minimum value of 0,88°C for the hour 6. This value is consistent since the indoor air comes from an energy loss process during the night hours and additionally, the solar radiation for the hour 6 is generally low in the intertropical zones, therefore, the energy gain inside the greenhouse is also low. The maximum value of ΔT was 5,93°C for hour 14, a value that is generated as a function of the thermal gain inside the greenhouse over the course of 7 hours (7-14 hours).

Temperature is one of the main parameters that affects the quality and growth of the flower stems. In the case of rose cultivation, it is recommended that optimum temperatures during the daytime period should be between 21 and 24°C Yong (2004). For the evaluated case, it is observed that this is fulfilled for the period between 11 and 16 hours; however, for the remaining hours the temperature are lower than recommended; this translates into longer vegetative cycles and lower stem production (Zieslin and Mor, 1990). Additionally, temperatures between 15 and 20°C under high humidity conditions favor the germination of Botrytis cinerea conidia (Restrepo, 2010), a disease that generates great economic losses for rose crops in the Bogota savanna.

Figure 3 shows the thermal distribution patterns in a greenhouse section at 1.5 m from ground level for the hours 8, 12 and 16 (Figure 3 A-C). For the eight hour, a homogeneous thermal behavior can be observed with an average value of 17,52°C (Figure 3A). For the twelfth hour an average temperature of 23,05°C was obtained, under this condition a heterogeneous behavior was observed in 40% of the evaluated area. Higher temperatures were found near the lateral zones with averages of 25,1°C and lower temperatures of 21,1°C in the central zone of the greenhouse, just in the areas of influence of the roof vents (Figure 3B). These heterogeneous conditions are not suitable for the production of flowers since they generate differentiated production volumes of stems in quantity and quality. For the sixteenth hour, the average interior temperature was 21,17°C, observing some areas with temperature values of 22,8°C near the lateral and frontal ventilation openings (Figure 3C). This allows us to conclude that the thermal conditions inside the greenhouse are not thoroughly homogeneous and the thermal behavior for the first hours of the day period are the most unfavorable for the production of roses.

Table 3 shows the average temperature values calculated for the night hours. The maximum temperature value of 16,19°C was registered for the hour 19, moment from which the behavior is decreasing as night passes. This is mainly due to the fact that under passive greenhouses an irradiative cooling is generated from the ground towards the external environment promoted by the confinement of the internal air and the high transmission of the plastic cover of the infrared radiation (Majdoubi et al., 2016).

The greenhouse has an acceptable thermal hermetic sealing, which allowed the generation of a positive ΔT for the hours evaluated, showing maximum and minimum ΔT values of 1,84 and 0,37°C for the first and last hour of the night, respectively. Taking into account that the general conditions were clear and humid it should be mentioned that the greenhouse presented a thermal behavior that limited the occurrence of the thermal inversion phenomenon that under these prevailing meteorological conditions is frequent within Colombian greenhouses.

The general recommendation for the cultivation of rose is to be able to guarantee night temperatures that oscillate between 15 and 17°C (Yong, 2004).The analysis of the data indicated that these recommended values are only obtained for hours 19 and 20. The following hours presented average values below the recommended value which can generate relevant effects on the growth and development of the plant; in productive terms, under these conditions it can be usual to obtain deformed flowers with a number of undesirable petals (Yong, 2004).

The qualitative analysis of thermal distribution patterns in a horizontal plane at 1.5 m above ground level for hours 20, 0 and 4 is presented in Figures 4 A, B and C. In general, a homogeneous thermal behavior is observed for the three hours evaluated in approximately 95% of the evaluated area.

El trabajo experimental se desarrolló en un invernadero de 8640 m² de área cubierta, perteneciente a una finca productora de Rosa (5º03’22.16’’N, 73º55’49.07’’W, 2583 msnm) ubicada en el municipio de Cogua (Cundinamarca).El invernadero evaluado estaba compuesto de 24 naves cada una con una luz de 6,0 m (Figura 1A), las alturas mínimas y máximas bajo canal fueron de 3,2 y 6,4 m respectivamente, la distancia longitudinal del invernadero fue de 60 m y estaba orientado en sentido noreste-suroeste (NE-SW). Cada nave disponía de una apertura de ventilación cenital móvil de 1,0 m de ancho y 48 m de longitud y unas aperturas laterales fijas de 1,9 m de ancho en cada uno de los cuatro costados de la estructura, para un área total de ventilación natural de 1927 m².

Se evaluaron las condiciones meteorológicas exteriores para las 24 horas del día, dentro del periodo comprendido entre las 00:00 horas del 29 de marzo y las 23:00 horas de 02 de mayo del año 2016. El registro de las variables se realizó por medio de una estación meteorológica (Vantage Pro2 Plus, Davis Instruments, Hayward CA), que integraba sensores de radiación global, temperatura, humedad relativa, precipitación, velocidad del viento y dirección y la frecuencia de registro fue cada 10 minutos.

La simulación numérica CFD permite resolver las ecuaciones gobernantes del flujo de fluidos utilizando el método de volumen finito (Molina et al., 2009). Estas ecuaciones pueden ser representadas como ecuaciones de convección-difusión de un fluido para tres leyes de conservación, que incluyen las ecuaciones de momento, energía y transporte de un fluido compresible y en un campo tridimensional (3D) y se expresan así:

donde ρ es la densidad del fluido (kg m^-3), $\nabla$ es el operador nabla, ϕ representa la variable de concentración, $\overrightarrow{v}$ es el vector de velocidad (m s^-1), Γ es el coeficiente de difusión (m² s^-1) y S representa el término fuente (Piscia et al., 2012). La naturaleza turbulenta del flujo de aire se simuló utilizando el modelo de turbulencia estándar k-ε, el cual se basa en dos ecuaciones principales, una para k que representa la energía cinética y otra para ε que representa la tasa de disipación en tiempo y volumen unitarios. Las ecuaciones de transporte para k y ε pueden ser modeladas como:

donde $\mu$ es viscosidad y ${\mu }_t$ es viscosidad turbulenta en (kg m^-1s^-1), ${\sigma }_k y {\sigma }_{ϵ}$ son los números turbulentos de Prandtl para k y $\epsilon$ , G _k es la generación de energía cinética turbulenta debida a los gradientes de velocidad, G _b es generación de energía cinética turbulenta debida a la flotabilidad, Y _M es la dilatación fluctuante en la turbulencia debida a la tasa de disipacion global y $v$ es el coeficiente de viscosidad cinemática. $C_{1ϵ}, C_{2ϵ}, C_{\mu }, {\sigma }_k y {\sigma }_{ϵ}$ son constantes con valores predeterminados. Este modelo estándar k-ε ha sido ampliamente usado y validado en estudios enfocados a invernaderos demostrando una precisión adecuada (Fatnassi et al., 2006; Katsoulas et al., 2006).Todas las simulaciones consideraron la ecuación de energía lo cual permitió analizar el campo escalar de temperaturas al interior del invernadero. Así mismo se modeló el fenómeno de flotabilidad, fenómeno que impulsa el movimiento del aire bajo condición de velocidad exterior se realizó por medio de la aproximación de Boussinesq, que se describe de la siguiente forma:

donde ρ _ref es la densidad constante del flujo; T _ref es la temperatura real (°C) y β es el coeficiente de expansión térmica para el aire (β = 0.00329 °C^-1). La aproximación de Boussinesq es válida si las diferencias de temperaturas que aparecen en el dominio computacional no son demasiado grandes β(T-T _ref ) < 1, situación que generalmente ocurre en el estudio de microclima de invernaderos (Baeza et al., 2006).

El modelo de radiación seleccionado fue el de ordenadas discretas (DO) con discretización angular. Este modelo permite realizar el análisis de clima en condiciones de periodo nocturno, simulando y resolviendo el fenómeno de radiación desde el suelo del invernadero hacia el ambiente exterior. Para tal fin se consideró el cielo como un cuerpo negro con una temperatura equivalente (T_C) de 0 °C para un escenario predominante de noches húmedas y despejadas (Iglesias et al., 2009).

El software de preprocesamiento ANSYS ICEM CFD (v. 17.0) se utilizó para generar un gran dominio computacional compuesto por el invernadero y sus alrededores, esto con el fin de garantizar la no afectación de la solución numérica del campo de flujo fuera del invernadero y para permitir una definición apropiada de la capa límite atmosférica (Rico, 2011). Las dimensiones del dominio computacional fueron de 354, 260 y 70 m para los ejes x, y, z, respectivamente (Figura 1C). Este tamaño se determinó siguiendo las pautas para el cálculo de CFD del entorno eólico alrededor de los edificios (Tominaga et al., 2008). El dominio computacional estuvo compuesto por una malla no estructurada de 76,651,871 volúmenes discretizados en el espacio, número de elementos que se obtuvo luego de verificar la independencia de las soluciones numéricas del flujo de aire a tamaños de malla con un número superior e inferior de elementos de acuerdo al procedimiento reportado por He et al. (2017). Un criterio fundamental para establecer la precisión de las soluciones consiste en evaluar la calidad de la malla (Figura 1B). Los parámetros de calidad evaluados fueron el tamaño de las celdas y la variación del tamaño de celda a celda encontrando que un 93,1% de las celdas de la malla estaban dentro del intervalo de calidad alta (0,95-1). También se evaluó la calidad ortogonal, donde el valor mínimo obtenido fue de 0,94, resultados que se clasifican dentro del rango de alta calidad (Flores et al., 2015). Los criterios de convergencia del modelo fueron establecidos en 10^-8 para la ecuación de energía y en 10^-6 para las ecuaciones de continuidad, momento y turbulencia (Baxevanou et al., 2017).

El límite superior del dominio y las superficies paralelas al flujo fueron fijados con condiciones de frontera de propiedades simétricas para no generar pérdidas de fricción del flujo de aire en contacto con estas superficies. Las simulaciones consideraron las características atmosféricas del municipio de Cogua tales como presión atmosférica de 74992 Pa, viscosidad del aire igual a 1,7E-05 kg m^-1 s^-1 y gravedad de 9,81 m s^-2. Otras propiedades del polietileno y del suelo agrícola como calor especifico (Cp), conductividad térmica (k) y densidad (ρ) fueron establecidas de acuerdo con Villagrán et al. (2012). Los límites perimetrales del dominio computacional se establecieron como límites de entrada de aire o salida de presión según el caso evaluar. Se consideró un perfil uniforme de velocidad del viento evaluando velocidades con valores entre 0,001y 1,21 m s^-1 y valores de temperatura medios para cada hora del día. Dichos valores fueron establecidos a partir de la información climática recopilada en el periodo de medición experimental.

Otros parámetros de entrada que son necesarios para alimentar el modelo de radiación son las propiedades ópticas del material de cubierta tales como, coeficiente absorción (α) de 0,69, coeficiente de transmisión (Ƭ) de 0,19 y coeficiente de reflexión (ρ) de 0,11. El modelo no incluyó cultivo alguno y se asumió una hermeticidad máxima del invernadero. Estas simplificaciones son válidas puesto que son aplicadas a cada una de los casos simulados y los errores que puedan derivarse tendrán el mismo grado de magnitud para cada escenario.

La evaluación se realizó asumiendo las condiciones meteorológicas medias para cada hora evaluada (Tabla 1) y la configuración de ventilación estándar del invernadero estudiado que, para este caso, fue ventilación lateral combinada con ventilación de techo para las horas diurnas (6-18 horas) y ventilación de techo para las condiciones nocturnas (19-5 horas).

La validación del modelo CFD se realizó a través de la comparación entre los resultados de las simulaciones y el registro experimental para la misma frecuencia temporal de la temperatura dentro del invernadero real durante el periodo de evaluación. La temperatura del aire del invernadero se registró por medio de treinta termopares tipo T (cobre-constantan) conectados a un número igual de registradores de datos (Cox-Tracer Junior, Escort DLS, Edison, NJ) desplegados de manera uniforme a 1.5 m de la superficie del suelo y registrando datos cada 10 minutos.

La validez del modelo se realizó mediante la comparación de los valores de temperatura medidos experimentalmente y los obtenidos a través de las simulaciones CFD, para las 24 horas del día. La Figura 2 muestra que el comportamiento de la temperatura para cada hora evaluada presenta un buen ajuste entre las dos metodologías usadas y a su vez una misma tendencia para el periodo de medición. Cuantitativamente al comparar los datos medidos y simulados se obtuvo un error porcentual absoluto medio (MAPE) de 4,26%, un error medio absoluto (MAE) de 0,75 °C y el error cuadrático medio (MSE) de 0,94 °C, valores que permiten garantizar que el uso del modelo CFD es adecuado para la evaluación y descripción del comportamiento térmico del invernadero.

En la Table 2 se presenta la evolución temporal de la temperatura interior del invernadero para el periodo diurno. Se observa que la variable presento un valor mínimo de 12,17 °C para la hora 6 y alcanzó un valor máximo de 24,31 °C en la hora 14, este aumento temporal está influenciado por la radiación solar. El diferencial térmico entre el interior y el exterior del invernadero (∆T= T_{media interior}-T_{media exterior}) exhibió un valor mínimo de 0,88 °C para la hora 6. Este valor es coherente puesto que el aire interior del invernadero viene de un proceso de pérdida energética durante las horas nocturnas y adicionalmente la radiación solar para la hora 6 es generalmente baja en las zonas intertropicales, por lo tanto, la ganancia energética en el interior del invernadero es baja también. El valor máximo de ∆T fue de 5,93 °C para la hora 14, valor que es generado en función de la ganancia térmica en el interior del invernadero en el transcurso de 7 horas (7-14 horas)

La temperatura es uno de los principales parámetros que incide en la calidad y crecimiento de los tallos florales. En el caso del cultivo de rosa es recomendable que las temperaturas óptimas en el periodo diurno se encuentren en valores que oscilen entre 21 y 24 °C según Yong (2004). Para el caso evaluado se observa que esto se cumple para el periodo comprendido entre las 11 y las 16 horas, para las horas restantes el valor de temperatura presenta niveles inferiores al recomendado; esto se traduce en ciclos vegetativos más largos y una menor producción de tallos Zieslin y Mor (1990). Adicionalmente, temperaturas entre 15 y 20°C y en condiciones de alta humedad favorecen la germinación de conidios de Botrytis cinérea Restrepo (2010), enfermedad que genera grandes pérdidas económicas en los cultivos de rosa de la sabana de Bogotá.

La Figura 3 presenta los patrones de distribución térmica en una sección vista de planta a 1,5 m del nivel del suelo para las horas 8, 12 y 16 (Figura 3 A-C). Para la hora 8 se puede observar un comportamiento térmico homogéneo con un valor medio de 17,52 °C (Figura 3A). Para la hora 12 se obtuvo un valor medio de temperatura de 23,05 °C, bajo esta condición se observa un comportamiento heterogéneo en un 40% del área evaluada, encontrando temperaturas más elevadas cerca de los costados laterales con valores medios de 25,1 °C y temperaturas más bajas con valores aproximados de 21,1 °C en la zona central del invernadero justo en las áreas de influencia de las ventilaciones de techo (Figura 3B). Estas condiciones heterogéneas no son adecuadas para la producción de flores ya que generan volúmenes de producción de tallos diferenciados en cantidad y calidad. Para la hora 16 el valor medio de temperatura interior obtenido fue de 21,17 °C, observando algunas zonas con valores de temperatura de 22,8 °C cerca de las aperturas laterales y frontales de ventilación (Figura 3C). Lo anterior permite concluir que las condiciones térmicas en el interior del invernadero no son totalmente homogéneas y el comportamiento térmico para las primeras horas del periodo diurno son las más desfavorables para la producción de rosa.

En la Table 3 se pueden observar los valores medios de temperatura calculados para las horas nocturnas. El máximo valor de temperatura de 16,19 °C se registró para la hora 19, momento a partir del cual el comportamiento es decreciente a medida que trascurre la noche. Lo anterior se debe principalmente a que en los invernaderos pasivos se genera un enfriamiento radiativo desde el suelo hacia el ambiente exterior promovido por el confinamiento del aire interno y la alta transmisión de la cubierta plástica a las radiaciones infrarrojas (Majdoubi et al., 2016).

El invernadero presenta una hermeticidad térmica aceptable, lo que permitió la generación de ∆T positivos para las horas evaluadas, presentando unos valores máximos y mínimos de ∆T de 1,84 y 0,37 °C para la primera y última hora de la noche respectivamente. Teniendo en cuenta que las condiciones generales fueron despejadas y húmedas se debe mencionar que el invernadero presentó un comportamiento térmico que limitó la ocurrencia del fenómeno de inversión térmica que bajo estas condiciones meteorológicas predominantes es frecuente en los invernaderos colombianos.

La recomendación general para el cultivo de rosa es poder garantizar temperaturas nocturnas que oscilen entre los 15 y 17 °C (Yong, 2004). El análisis de los datos indicó que estos valores recomendados solo se obtienen para las horas 19 y 20. Las horas posteriores presentan valores medios por debajo del valor recomendado lo cual puede generar afectaciones relevantes en el crecimiento y desarrollo de la planta, en términos productivos bajo estas condiciones predomina la obtención de flores deformadas con un número de pétalos no deseables (Yong, 2004).

El análisis cualitativo de los patrones de distribución térmica en un plano a 1,5 m sobre el nivel del suelo para las horas 20, 0 y 4 se presenta en las Figuras 4 A, B y C. En general se observa un comportamiento térmico homogéneo para las tres horas evaluadas en aproximadamente un 95% del área evaluada.