Upto this I have done. Please help me.
read_grads('2016081600-rainfall-3da.ctl','rain');
read_grads('imdmsg1.ctl','rf');
%% Box to read%% Latitude1=5;Latitude2=38; Longitude1=65;Longitude2=100;
Bx(1)=Latitude2; Bx(2)=Longitude1; Bx(3)=Latitude1; Bx(4)=Longitude2; % Box for Interpolation
stx=65;stpx=.036;nx=1001-1; %%% Longitude from CTL file sty=5;stpy=.036;ny=1001-1; %%% Latitude from CTL file
[myX myY]=meshgrid(stx:stpx:stx+nx*stpx,sty:stpy:sty+ny*stpy);