##–LLAMADO DE LAS LIBRERIAS–##

library(DT)

##–LAMADO DE LA BASE DE DATOS–##

Aqui se llama a la base de datos y se atribuyen valores

setwd("C:/Users/isabe/OneDrive/Desktop/SEMESTRE 6/DISEÑO Y MODELACIÓN DE CULTIVOS/MODELACIÓN/MODULO 1")
tolima=read.table("13_Tolima_Espinal.csv", header = TRUE, sep = ",")
nueva.tolima=tolima[,1:5]
DT:: datatable(nueva.tolima)

Obtencion de variables para el calculo de Area foliar

DTU

nueva.tolima$DTU=((nueva.tolima$Tem..Maxima..C.+ nueva.tolima$Tem..Minima..C.)/2)-10

INODE

nueva.tolima$INODE=nueva.tolima$DTU/100

Sumatoria del INODE, denominado MSNN

INODE=length(nueva.tolima$INODE)
INODE
## [1] 365
nueva.tolima$MSNN[1]=nueva.tolima$INODE[1]

for (i in 2:365){
  nueva.tolima$MSNN[i]=nueva.tolima$MSNN[i-1]+nueva.tolima$INODE[i]
}

###—CREACION DE UNA OPERACION PARA AUTOMATIZAR LA SEGMENTACION DE LOS DATOS DE INTERES—###

###---NUMERO DE NUDOS---###
nudos=16
df=nueva.tolima$MSNN
for (i in INODE){
  icut=which(df<=nudos)
}
icut1=length(icut)
data.tolima=nueva.tolima[0:icut1+1,]
names(data.tolima)[3:5]=c("R.Solar", "Tmax", "Tmin") 
fecha=data.tolima[icut1,1]
as.character(fecha)
## [1] "3/21/2010"
print(c("los", nudos,"nudos se alcanzan el", fecha))
## [1] "los"                  "16"                   "nudos se alcanzan el"
## [4] "3/21/2010"

###—CALCULO AREA FOLIAR—###

#PARAETROS
#PLA1
paramA1=0.0235
paramB1=0.166
#PLA2
paramA2=0.08
paramB2=0.166

data.tolima$PLA1=paramA1*exp(paramB1*data.tolima$MSNN)
data.tolima$PLA2=paramA2*exp(paramB2*data.tolima$MSNN)

data.tolima=data.tolima[,c(1,3,4,5,6,7,8,9,10)]
PDEN=8
data.tolima$IAF1=data.tolima$PLA1*PDEN
data.tolima$IAF2=data.tolima$PLA2*PDEN
DT::datatable(data.tolima)

###—Graficas—###

plot(data.tolima$MSNN, type="l",xlab = "Dias",ylab = "Nudos acumulados")

plot(data.tolima$PLA1,type = "l", ylim= c(0,max(data.tolima$PLA2)),xlab = "Dias",ylab = "Indice de Area Foliar")
lines(data.tolima$PLA2, type = "l", col= "red")

CONCLUSIONES

La oferta agro-climatológica (temperatura y luminosidad) de las dos localizaciones es muy competitiva debido a que el desarrollo de 16 nudos en las dos se obtiene al día 82 en Neiva y en Espinal al día 81, teniendo como parámetros a=0.0235 y b=0.166 para el desarrollo del área foliar de la planta

El índice de área foliar (LAI) de NEIVA-HUILA converge en la parametrización a y b en el nudo número 15. Esto varía debido a la parametrización que se tuvo en cuenta. Se alcanzaron áreas foliares que no superaron los 1.2m2/planta en ambas localidades, lo que sigue demostrando la similitud de temperaturas en ambas localidades.

Los DTU obtenidos fueron 18.25°, logrando un desarrollo diario de nudos de un 18.25%; lo que indica que se requieren 5 días y medio aproximadamente para alcanzar el PHYL de 100°C, en ambas localizaciones