Temperature Layer

The Temperature Layer code processes the CR2 temperature product to obtain the mean temperature layer between 1997 to 2016. The initial_year and final_year defines the initial and last year of execution respectively. The initial_year_temp and the final_year_temp as the initial and last year of execution of the CR2 product respectively.

The layer_temp is defined as the product of the average mean temperature defined by the CR2 climate research centre. After the definition of the Temperature raster folder for the input_temp, we read the first layer of the netCDF file. Each netCDF band is a month and the band corresponding to January of the first year of the simulation is calculated.

temp <- raster(paste0(input_temp, "/", layer_temp), band = 1)
dif_year <- initial_year - initial_year_temp
initial_band <- (dif_year * 12) + 1
final_band <- ((final_year_temp - initial_year_temp) * 12) + 12

Later we calculate the mean temperature by adding all the bands and dividing them by the number of bands.

for (i in initial_band:final_band) {
  temp <- raster(paste0(input_temp, "/", layer_temp), band = i)
  if (i == initial_band) {

    sum_temp <- temp
    count <- 0
  } else {
    sum_temp <- sum_temp + temp
    count <- count + 1
  }
  setTxtProgressBar(pb, i)
}

mean_temp <- sum_temp / count

The maximum and minimum temperature is calculated and the output is written to a raster file.