It is also not surprising that the Tibetan plateau
specific humidity at this pressure level is considerably higher than the zonal mean (3.5 - 4
times larger during the summer).
library (ncdf) ncep500 = open.ncdf («NCEP.nc») ncep500 [1] «file NCEP.nc has 4 dimensions:» [1] «lat Size: 73» [1] «level Size: 1» [1] «lon Size: 144» [1] «
time Size: 361» [1] «----------------» [1] «file NCEP.nc has 1 variables:» [1] «short shum [lon, lat, level,
time] Longname: Monthly Mean of
Specific Humidity Missval: 32766 ″ nc500 = get.var.ncdf (ncep500) dim (nc500)[1] 144 73 361 nc25SN = nc500 -LSB-, mean (c (27:47)-RRB-,] dim (nc25SN)[1] 144 361 SN25 = colMeans (nc25SN) x = 1:361 lmSN25 = lm (SN25 ~ x) summary (lmSN25) acf (residuals (lmSN25)-RRB- $ acf [2][1] 0.5715495 plot (x, SN25, type =» b», main =»
Specific Humidity 25S to 25N from 01/1979 to 01/2009 ″, xlab = «Months Starting at 01/1979 ″, ylab =» Secific
Humidity», col =» dark red»)