[R-sig-Geo] Fwd: TRMM data Analysis

[sownal chand]

https://drive.google.com/file/d/1krtBKQUIhupvYKEG1gBV8bgQaQxx96pb/view?usp=drivesdk


Attached is the link for the data


---------- Forwarded message ---------
From: sownal chand <sownalchandfms using gmail.com>
Date: Fri, Sep 2, 2022, 10:10
Subject: TRMM data Analysis
To: <r-sig-geo using r-project.org>


Dear sir/madam

I have a data set in NCDF file and was visualizing the data precipitation
data set but ran into some small errors.
The code is below:
*******************************************************
library(ncdf4)

ncin <- nc_open("c:/Users/Downloads/PreNCDF.nc")

#get longitutes and latitudes

lon <<- ncvar_get(ncin,"lon")

nlon <- dim(lon)

lat <<- ncvar_get(ncin,"lat", verbose=F)

nlat <- dim(lat)

print(ncin)
#get the time variable and attributes

tori <<- ncvar_get(ncin,"TRMM_3B42RT_7_precipitation") #for time original

tunits <<- ncatt_get(ncin,"TRMM_3B42RT_7_precipitation","units")

nt <- dim(tori)


View(nt)
#get the variable and its attributes

TRMM.array <<- ncvar_get(ncin,"TRMM_3B42RT_7_precipitation") #dimensions
(lon,lat,time)
print(TRMM.array)

dim(TRMM.array)

#m is the time slice you want to plot

tmp.slice <- TRMM.array[, , 1]

image(lon, lat, tmp.slice, col = rev(brewer.pal(10, "RdBu")))

#betterway to plot

TRMM.Precip <- tmp.array [, , m]  #note this is the same as above

mapCDFtemp <- function(lat,lon,precip) #model and perc should be a string

{

  titletext <- "title"

  expand.grid(lon, lat) %>%

    rename(lon = Var1, lat = Var2) %>%

    mutate(lon = ifelse(lon > 180, -(360 - lon), lon),

           precip = as.vector(precip)) %>%

    #mutate(tas = convert_precip(, "k", "c")) %>%



    ggplot() +

    geom_point(aes(x = lon, y = lat, color = precip),

               size = 0.8) +

    borders("world", colour="black", fill=NA) +

    scale_color_viridis(na.value="white",name = "precip") +

    theme(legend.direction="vertical", legend.position="right",
legend.key.width=unit(0.4,"cm"), legend.key.heigh=unit(2,"cm")) +

    coord_quickmap() +

    ggtitle(titletext)

}


*******************************************
Would really appreciate if it could be corrected and really appreciate your
assistance

thanking you in advance
sownalc

	[[alternative HTML version deleted]]