update_Track {TCHazaRds}

Calculate Additional TC Parameters, and temporally Interpolate Along a Tropical Cyclone Track

Description

Calculate Additional TC Parameters, and temporally Interpolate Along a Tropical Cyclone Track

Usage

update_Track(
  outdate = NULL,
  indate,
  TClons,
  TClats,
  vFms,
  thetaFms,
  cPs,
  rMaxModel,
  vMaxModel,
  betaModel,
  eP,
  rho = NULL
)

Arguments

outdate

POSIX times to be interpolated to

indate

POSIX input times

TClons

input central TC longitude

TClats

input central TC latitude

vFms

input forward velocity of TC

thetaFms

input forward direction

cPs

central pressure

rMaxModel

empirical model for radius of maximum wind calculation (rMax in km)

vMaxModel

empirical model for maximum wind velocity calculation (vMax in m/s)

betaModel

empirical model for TC shape parameter beta (dimensionless Beta)

eP

background environmental pressure (hPa)

rho

air density

Value

list of track data inclining the rMax vMax and Beta.

Examples

paramsTable <- read.csv(system.file("extdata/tuningParams/defult_params.csv",package = "TCHazaRds"))
params <- array(paramsTable$value,dim = c(1,length(paramsTable$value)))
colnames(params) <- paramsTable$param
params <- data.frame(params)
require(terra)
TCi <- vect(system.file("extdata/YASI/YASI.shp", package="TCHazaRds"))
TCi$PRES <- TCi$BOM_PRES
t1 <- strptime("2011-02-01 09:00:00","%Y-%m-%d %H:%M:%S", tz = "UTC") #first date in POSIX format
t2 <- strptime(rev(TCi$ISO_TIME)[1],"%Y-%m-%d %H:%M:%S", tz = "UTC") #last date in POSIX format
outdate <- seq(t1,t2,"hour") #array sequence from t1 to t2 stepping by “hour”
TCil = update_Track(outdate = outdate,
                   indate = strptime(TCi$ISO_TIME,"%Y-%m-%d %H:%M:%S", tz = "UTC"),
                   TClons = TCi$LON,
                   TClats = TCi$LAT,
                   vFms=TCi$STORM_SPD,
                  thetaFms=TCi$thetaFm,
                   cPs=TCi$PRES,
                  rMaxModel=params$rMaxModel,
                   vMaxModel=params$vMaxModel,
                   betaModel=params$betaModel,
                   eP = params$eP,
                   rho = params$rhoa)