r - error when specifying latitute dimension when reading Netcdf files -


i want read 2 netcdf files (climate data 2 years) , combine data of 1 single grid (lon [15] , lat[15]) array (named prec_year), script below.

i got error:

"error in r_nc_get_vara_double: netcdf: index exceeds dimension bound var: precip  ndims: 3   start: 0,4294967284,59count: 366,1,1error in get.var.ncdf   (ex.nc, "precip", start = c(lon[1], lat[15], 1),  :

c function r_nc_get_vara_double returned error"

but when dont specify lat[], example "start = c(lon[1], 1,1)", works fine. please (1) why got error, , (2) how fix it?

many thanks! --------------------

setwd ("d://historic climate data/aphrodyte/prec_025/join 2 years") prec_origin <- c(1:2) filenames <- list.files(pattern = null, all.files = false)  for(i in filenames){ filepath <- file.path("d://historic climate data/aphrodyte/prec_025/join 2 years/",paste    (i,sep=""))     ex.nc = open.ncdf (filepath, readunlim=false)     lon = get.var.ncdf( ex.nc, "longitude")          # coordinate variable   lat = get.var.ncdf( ex.nc, "latitude")          # coordinate variable   day =  get.var.ncdf( ex.nc, "time")            # time variable   #change lon[x] , lat[y] select different files/locations comparison   prec_year = get.var.ncdf(ex.nc, "precip", start = c(lon[1],lat[15],1), count = c(1,1,-1))     prec_origin <- c(prec_origin, prec_year)  #combine data 1 array   close.ncdf (ex.nc)  }  prec_origin <- prec_origin[-(1:2)] # delete first 2 rows- -------------- [1] "file d://historic climate data/aphrodyte/prec_025/join 2     years//aphro_ma_025deg_v1003r1.1952.nc has 3 dimensions:" [1] "time   size: 366" [1] "longitude   size: 360" [1] "latitude   size: 280" [1] "------------------------" [1] "file d://historic climate data/aphrodyte/prec_025/join 2 years//aphro_ma_025deg_v1003r1.1952.nc has 2 variables:" [1] "float precip[longitude,latitude,time]  longname:daily precipitation analysis interpolated onto 0.25deg grids [mm/day] missval:-99.9000015258789" [1] "float rstn[longitude,latitude,time]  longname:ratio of 0.05deg-grids station [%]     missval:-99.9000015258789"

in "start" function need define index number, , not actual value wish start.

i believe looking for:

prec_year = get.var.ncdf(ex.nc, "precip", start = c(1,15,1), count = c(1,1,-1))

lon[1] doesn't produce error because gives coordinate can index (for longitude; 1:360). however, output not want (if understood you're question).

as information, if want select region based on pre-selected coordinates, follows (maybe not shortest code):

lon_select=which(x<=10 & x>=20) lat_select=which(y<=40 & y>=50)  precip=get.var.ncdf(ex.nc,"precip",start=c(lon_select[1],lat_select[1],1), count=c(length(lon_select),length(lat_select),-1))

Comments

Post a Comment

Popular posts from this blog

curl - PHP fsockopen help required -

i trying work on cloudflare api , have found in documentation. unfortunately have no idea how can use in fsockopen or through curl. here example of post request. need know how can make request below data post using curl or fsockopen curl https://www.cloudflare.com/api_json.html -d 'a=cache_lvl' \ -d 'tkn=8afbe6dea02407989af4dd4c97bb6e25' \ -d 'email=sample@example.com' \ -d 'z=example.com' \ -d 'v=agg' here link cloudflare https://www.cloudflare.com/docs/client-api.html#s4.2 try one: $ch = curl_init("https://www.cloudflare.com/api_json.html"); curl_setopt($ch, curlopt_returntransfer,1); curl_setopt($ch, curlopt_followlocation, 1); curl_setopt($ch, curlopt_ssl_verifypeer,false); // handling of -d curl parameter here. $param = array( 'a' => 'cache_lvl', 'tkn' => '8afbe6dea02407989af4dd4c97bb6e25', 'email' => 'sample@example.com', 'z' =...
Read more

HTTP/1.0 407 Proxy Authentication Required PHP -

i working in restricted network/limited access internet. since have use proxy make web service call. here php code calling web service , response . $client = new soapclient("http: //www.xyz.com/abc/wsdls/login.wsdl", array("trace" => 1, "exceptions" => 1, 'proxy_host' => $proxy_host, 'proxy_port' => $proxy_port, 'proxy_login' => $proxy_login, 'proxy_password' => $proxy_password)); response is fatal error: uncaught soapfault exception: [wsdl] soap-error: parsing wsdl: couldn't load ' http://www.xyz.com/abc/wsdls/login.wsdl ' : failed load external entity if copy url error message , paste in browser(with same proxy settings) works. 2nd scenario is, if pass values instead of variables giving me diffarant response $client = new soapclient("http: //www.xyz.com/abc/wsdls/login.wsdl", array("trace" => 1, "exceptions" => 1, 'proxy_host' ...
Read more

java - More than one row with the given identifier was found: 1, for class: com.model.Diagnosis -

below model, service implementation , error trace of code. model class: @entity @table(name = "ph_diagnosis_history") @jsonidentityinfo(generator = objectidgenerators.propertygenerator.class, property = "id",scope=diagnosis.class) public class diagnosis extends abstractentity { @column(name = "diagnosis_notes", length=255) @length(max = 255, message = iprprop.max_len_255_msg) private string diagnosisnote; @enumerated(enumtype.string) private diagnosistypeenum type; @temporal(temporaltype.date) @column(name = "diagnosis_date") @jsonserialize(using = jsondateserializer.class) private date diagnosisdate; @manytoone @joincolumn(name = "patient_chart_id") private patientchart patientchart; public string getdiagnosisnote() { return diagnosisnote; } public void setdiagnosisnote(string diagnosisnote) { this.diagnosisnote = diagnosisnote; } public d...
Read more