#Diretório  raiz
rootDir <- 'C:\\Users\\user\\Desktop\\FRANCISCO\\'
#Diretório da base de dados observados a ser lida
dataDir <- paste(rootDir,'entrada\\',sep='')
#Diretório onde escreverá os arquivos (**esse diretório criei antes de rodar o script)
outputloc <- paste(rootDir,'saida\\',sep='')

cidades = c("abelardoluz","aguaclara")

lat = c(-30, 0)

for (i in cidades){
  
  estacao = i
  a = read.table(paste(dataDir,estacao,'.txt',sep=''), sep = ";", header = T); str(a)

  #--------------------Separar colunas de interesse--------------------#
  a$DATE = as.Date(levels(a$DATE))[a$DATE]; str(a)
  a$DiaJuliano = format(a$DATE, trim = T, '%j') ; a$DiaJuliano
  a$DiaJuliano = as.numeric(format(a$DATE, trim = T, '%j')); str(a$DiaJuliano) 

  #--------------------Calculo das Variaves Astronomicas--------------------#
  head(a)
  corr = pi/180
    
    a$decl = 23.45*sin(corr*((a$DiaJuliano-80)*360/365)); a$decl
    a$hn = 1/corr*acos(-tan(corr*lat[(as.numeric(length(i)))])*tan(corr*a$decl))
    a$N = 2*a$hn/15
    head(a)
  
  write.table(a,file=paste(outputloc,i,'.txt',sep=''),col.names=T, row.names=F, sep="\t")
 }