IMD High resolution gridded daily rainfall data (1951-2007) - version 3 This data product is a the third version of IMD gridded daily rainfall data developed in 2007. We have updated the data upto 2007. This data is arranged in 35x33 grid points. All the gridded data are in the directory “data” in the CDROM. The "data" directory contains 3 subdirectories viz. "binary" ,"ascii" and "stn_grd". The "binary" folder contains 57 binary data files one each for 57 years (1951-2007). The data file is named as rf_1951.grd, rf_1952.grd etc. The relevant Grads control file, rf.ctl is also included in the data directory. For leap years, data for 366 days are included. Sample Fortran and C programs to read the data are given below. There are 57 data files in ascii format in the subdirecory "ascii", one each for 57 years (1951-2007). The data file is named as drf_1951.prt, drf_1952.prt etc. A Sample Fortran program to read the data is also given in the subdirectory. Number of stations per grid are stored in the subdirectory "stn_grd" There are 57 binary files in the subdirecory , one each for 57 years (1951-2007). The data file is named as stn1951.grd, stn1952.grd etc. The relevant Grads control file, rfstn.ctl is also included in the data directory The unit of rainfall is in mm. ======================================================== PROGRAM READ PARAMETER(ISIZ=35,JSIZ=33) DIMENSION RF(366,ISIZ,JSIZ) OPEN(7,FILE='E:\DATA\RF_1951.GRD', + FORM='UNFORMATTED',ACCESS='DIRECT',RECL=ISIZ*JSIZ*4, + STATUS='OLD') C TAKE NDAY=366 FOR LEAP YEARS NDAY=365 DO IDAY=1,NDAY READ(7,REC=IDAY)((RF(IDAY,I,J),I=1,ISIZ),J=1,JSIZ) ENDDO STOP END =========================================================== Sample C program to read the data is given below: ==================================================== /* This program reads binary data for 365/366 days and writes in ascii file. */ #include main() { float rf[35][33]; int i,j ,k; FILE *fptr1,*fptr2; fptr1=fopen("d:\\jb\\rf_1951.grd","rb"); // input file fptr2=fopen("d:\\jb\\rf_1951.prt","w"); if(fptr1==NULL) { printf("Can't open file"); return 0; } if(fptr2==NULL) { printf("Can't open file"); return 0; } k=1; while(fread(&rf,sizeof(rf),1,fptr1)==1) { fprintf(fptr2,"\nrecord = %d \n",k); k++; for(i=0 ; i < 35 ; i++) { fprintf(fptr2,"\n") ; for(j=0 ; j < 33 ; j++) fprintf(fptr2," %6.1f",rf[i][j] ); } } fclose(fptr1); fclose(fptr2); return 0; } /* end of main */ ===================================================== The data on the grid is from South to North. i.e. J=1 corresponds to lat. 6.5 N, J=2 corresponds to 7.5 N and so on. The first data in the record is at 6.5N and 66.5E, and the second is at 6.5N, 67.5E etc The reference paper is given in the directory “DOC”. The relevant reference research paper is : 1. Rajeevan, M., J. Bhate, J.D. Kale and B. Lal, 2006 : A high resolution daily gridded rainfall for the Indian region : analysis of break and active monsoon spells, 2006, Current Science, Vol.91, No.3 pp 296-306. These are also included in the CD ROM as PDF file “ref_paper.pdf” in the doc directory. We would be grateful if you cite this paper in any publiction you may produce using these data and also send us a reprint of your paper after it is published. Disclaimer Although every care has been taken in preparing and testing the data set, India Meteorological Department cannot guarantee that the data are correct in all circumstances. India Meteorological Department also does not accept any liability whatsoever for any error or omission in the data, or for any loss or damage arising from its use. For comments and questions on the data or if you encounter any problems, please contact: National Climate Centre India Meteorological Department Pune, INDIA. 411 005. Phone: 091-20-2553 5211 / 2553 5877 FAx: 091-20-2553 5435. E Mail: ncc@imdpune.gov.in