(no subject)

    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. Regards Azizul   ********fatal:["Execute.c":7940]:Number of subscripts on left-hand-side do not match                                     number of dimensions of variable: (3),  Subscripts used: (2)   fatal:["Execute.c":8567]:Execute: Error occurred at or near line 63 in file mmp1.ncl*******       oad "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"   begin      lat1     =  18.0    ; Region    lat2     =  25.0    lon1     =  75.0    lon2     =  90.0   inname1="/Projects/CMIP5A/Daily/historical/CCSM4/r2i1p1/pr_day_CCSM4_historical_r2i1p1_18500101-20051231.nc" inname2="/Projects/CMIP5A/Daily/historical/CNRM-CM5/r1i1p1/pr_day_CNRM-CM5_historical_r1i1p1_18500101-20051231.nc" inname3="/Projects/CMIP5A/Daily/historical/CESM1-BGC/r1i1p1/pr_day_CESM1-BGC_historical_r1i1p1_18500101-20051231.nc" inname4="/Projects/CMIP5A/Daily/historical/GFDL-ESM2M/r1i1p1/pr_day_GFDL-ESM2M_historical_r1i1p1_18610101-200512.nc" inname5="/Projects/CMIP5A/Daily/historical/NorESM1-M/r1i1p1/pr_day_NorESM1-M_historical_r1i1p1_18500101-20051231.nc" inname6="/Projects/CMIP5A/Daily/historical/MIROC-ESM/r1i1p1/pr_day_MIROC-ESM_historical_r1i1p1_18500101-20051231.nc" inname7="/Projects/CMIP5A/Daily/historical/FGOALS-g2/r1i1p1/pr_day_FGOALS-g2_historical_r1i1p1_19500101-20051231.nc" inname8="/Projects/CMIP5A/Daily/historical/HadGEM2-ES/r1i1p1/pr_day_HadGEM2-ES_historical_r1i1p1_18591201-200511.nc" inname9="/Projects/CMIP5A/Daily/historical/bcc-csm1-1-m/r1i1p1/pr_day_bcc-csm1-1-m_historical_r1i1p1_18500101-20121231.nc"     inname10="/Projects/CMIP5A/Daily/rcp85/CCSM4/r2i1p1/pr_day_CCSM4_rcp85_r2i1p1_20060101-21001231.nc" inname11="/Projects/CMIP5A/Daily/rcp85/CNRM-CM5/r1i1p1/pr_day_CNRM-CM5_rcp85_r1i1p1_20060101-21001231.nc" inname12="/Projects/CMIP5A/Daily/rcp85/CESM1-BGC/r1i1p1/pr_day_CESM1-BGC_rcp85_r1i1p1_20060101-21001231.nc" inname13="/Projects/CMIP5A/Daily/rcp85/GFDL-ESM2M/r1i1p1/pr_day_GFDL-ESM2M_rcp85_r1i1p1_20060101-21001231.nc" inname14="/Projects/CMIP5A/Daily/rcp85/NorESM1-M/r1i1p1/pr_day_NorESM1-M_rcp85_r1i1p1_20060101-21001231.nc" inname15="/Projects/CMIP5A/Daily/rcp85/MIROC-ESM/r1i1p1/pr_day_MIROC-ESM_rcp85_r1i1p1_20060101-21001231.nc" inname16="/Projects/CMIP5A/Daily/rcp85/FGOALS-g2/r1i1p1/pr_day_FGOALS-g2_rcp85_r1i1p1_20060101-21011231.nc" inname17="/Projects/CMIP5A/Daily/rcp85/HadGEM2-ES/r1i1p1/pr_day_HadGEM2-ES_rcp85_r1i1p1_20051201-22991230.nc" inname18="/Projects/CMIP5A/Daily/rcp85/bcc-csm1-1-m/r1i1p1/pr_day_bcc-csm1-1-m_rcp85_r1i1p1_20060101-21001231.nc"   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.,102) 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(:)     pr20in3=infile3->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr21in12=infile12->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr20avg3=dim_avg_Wrap(dim_avg_Wrap(pr20in3)) pr21avg3=dim_avg_Wrap(dim_avg_Wrap(pr21in12)) printVarSummary(prsa(2,:)) printVarSummary(pr20avg3(:)) prsa(2,:)=pr20avg3(:) diffsa(2,:)=pr21avg3(:)-pr20avg3(:)     pr20in4=infile4->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr21in13=infile13->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr20avg4=dim_avg_Wrap(dim_avg_Wrap(pr20in4)) pr21avg4=dim_avg_Wrap(dim_avg_Wrap(pr21in13)) printVarSummary(prsa(3,:)) printVarSummary(pr20avg4(:)) prsa(3,:)=pr20avg4(:) diffsa(3,:)=pr21avg4(:)-pr20avg4(:)     pr20in5=infile5->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr21in14=infile14->pr(:,{lat1:lat2},{lon1:lon2}) ; (time,lat,lon)  (0,1,2) pr20avg5=dim_avg_Wrap(dim_avg_Wrap(pr20in5)) pr21avg5=dim_avg_Wrap(dim_avg_Wrap(pr21in14)) printVarSummary(prsa(4,:)) printVarSummary(pr20avg5(:)) prsa(4,:)=pr20avg5(:) diffsa(4,:)=pr21avg5(:)-pr20avg5(:)

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