Dimension size mismatch

Dear NCL, I am computing 50 -99.9 percentile of precipitation using historical and rcp85 data from CMIP5 projects and found following fatal error.  You can find the error line is highlighted in the script  Appreciated for your help.   Please note I have gone through the R Manual of NCL and I have still problems about dimension. Any recommendation for further reading is really helpful. RegardsAzizul ***********fatal:Dimension size mismatch, dimension (0) of left hand side reference does not have the same size as the right hand side reference after subscripting. fatal:["Execute.c":8567]:Execute: Error occurred at or near line 63 in file mmp1.ncl************ Error part of the scripts infile1=addfile(inname1,"r") infile2=addfile(inname2,"r") infile3=addfile(inname3,"r") infile4=addfile(inname4,"r") infile5=addfile(inname5,"r") infile6=addfile(inname6,"r") infile7=addfile(inname7,"r") infile8=addfile(inname8,"r") infile9=addfile(inname9,"r") infile10=addfile(inname10,"r") infile11=addfile(inname11,"r") infile12=addfile(inname12,"r") infile13=addfile(inname13,"r") infile14=addfile(inname14,"r") infile15=addfile(inname15,"r") infile16=addfile(inname16,"r") infile17=addfile(inname17,"r") infile18=addfile(inname18,"r") pct=fspan(50.,100.,2) prsa=new((/9,dimsizes(pct)/),"float") diffsa=prsa pr20in1=infile1->pr (:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr21in10=infile10->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr20avg1=dim_avg_Wrap(dim_avg_Wrap(pr20in1)) pr21avg1=dim_avg_Wrap(dim_avg_Wrap(pr21in10)) printVarSummary(prsa(0,:)) printVarSummary(pr20avg1(:)) prsa(0,:)=pr20avg1(:) diffsa(0,:)=pr21avg1(:) pr20in2=infile2->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr21in11=infile11->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr20avg2=dim_avg_Wrap(dim_avg_Wrap(pr20in2)) pr21avg2=dim_avg_Wrap(dim_avg_Wrap(pr21in11)) printVarSummary(prsa(1,:)) printVarSummary(pr20avg2(:)) prsa(1,:)=pr20avg2(:) diffsa(1,:)=pr21avg2(:)-pr20avg2(:)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu May 22 14:13:43 2014