Zac2, R - rextester

Zac2

Language: Layout:

data_at_time_1 <- read.table(text = 'Date lon lat uwnd vwnd
2008/12/31 -67.25 -55.75 8.2 -2.6
2009/01/01 -69.75 -55.75 -7.9 1.2
2009/01/02 -69.25 -55.75 1.4 4.9
2009/01/03 -68.75 -55.75 -7.8 2.3
2009/01/04 -68.25 -55.25 8.9 -1.2
2009/01/05 -67.75 -55.25 -9.1 0.8
2009/01/06 -67.25 -55.25 3.0 1.3
2009/01/07 -71.75 -54.75 -1.8 -1.6
2009/01/08 -70.75 -54.75 8.7 0.2
2009/01/09 -70.25 -54.75 -3.9 1.4
2009/01/10 -70.75 -54.75 0.0 0.0
2009/01/11 -69.75 -54.75 -2.8 3.4',
                             header = TRUE)

data_at_time_1

X = data_at_time_1[4]
Y= data_at_time_1[5]
head(X)
head(Y)

m <- 1
for (i in 1:length(X)){

p <- X[i]
q <- Y[i]

temporary_function <- function(p, q)
{
  if ((p > 0) & (q > 0))
  {
     (270 - ((180/pi)*(atan(q/p))))
  } else if ((p > 0) & (q < 0))
  {
    (270 - ((180/pi)*(atan(q/p))))
  } else if ((p < 0) & (q > 0))
  {
    (90 - ((180/pi)*(atan(q/p))))
  } else if ((p < 0) & (q <0))
  {
    (90 - ((180/pi)*(atan(q/p))))
  } else if ((p == 0) || (q == 0))
  {
    90
  } else
  {
    0
  }
}

mapply(FUN = temporary_function, X, Y)
JJJ = ((mapply(FUN = temporary_function, X, Y)))
}
print(JJJ)
i

KKK <- matrix(JJJ, nrow = length(X), byrow = TRUE) 
# KKK <- matrix(JJJ, nrow = n, byrow = TRUE)  # to get results for desired n rows
print(KKK)

[ + ] Show input
Absolute running time: 0.49 sec, cpu time: 0.87 sec, memory peak: 31 Mb, absolute service time: 0,5 sec

Error(s), warning(s):

Warning messages:
1: In if ((p > 0) & (q > 0)) { :
  the condition has length > 1 and only the first element will be used
2: In if ((p > 0) & (q < 0)) { :
  the condition has length > 1 and only the first element will be used
3: In if ((p > 0) & (q > 0)) { :
  the condition has length > 1 and only the first element will be used
4: In if ((p > 0) & (q < 0)) { :
  the condition has length > 1 and only the first element will be used

Date lon lat uwnd vwnd 1 2008/12/31 -67.25 -55.75 8.2 -2.6 2 2009/01/01 -69.75 -55.75 -7.9 1.2 3 2009/01/02 -69.25 -55.75 1.4 4.9 4 2009/01/03 -68.75 -55.75 -7.8 2.3 5 2009/01/04 -68.25 -55.25 8.9 -1.2 6 2009/01/05 -67.75 -55.25 -9.1 0.8 7 2009/01/06 -67.25 -55.25 3.0 1.3 8 2009/01/07 -71.75 -54.75 -1.8 -1.6 9 2009/01/08 -70.75 -54.75 8.7 0.2 10 2009/01/09 -70.25 -54.75 -3.9 1.4 11 2009/01/10 -70.75 -54.75 0.0 0.0 12 2009/01/11 -69.75 -54.75 -2.8 3.4 uwnd 1 8.2 2 -7.9 3 1.4 4 -7.8 5 8.9 6 -9.1 vwnd 1 -2.6 2 1.2 3 4.9 4 2.3 5 -1.2 6 0.8 uwnd [1,] 287.5924 [2,] 278.6371 [3,] 195.9454 [4,] 286.4293 [5,] 277.6790 [6,] 275.0241 [7,] 246.5713 [8,] 228.3665 [9,] 268.6831 [10,] 289.7468 [11,] NaN [12,] 320.5275 [1] 1 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 287.5924 278.6371 195.9454 286.4293 277.679 275.0241 246.5713 228.3665 [,9] [,10] [,11] [,12] [1,] 268.6831 289.7468 NaN 320.5275