ftp://ftp.cpc.ncep.noaa.gov/precip/CMORPH_V1.0/RAW/8km-30min/

Hello Bipasha,

Which software or program language do you normaly use for handling CMORPH data?

I  wrote a code for handling CMORPH data in MATLAB, though  a fortran read code is provided with the dataset

Hello Bipasha

can you share the matlab code for read cmorph 8km-30min. Regards 

pfa

Thanks

Related to this topic I was hoping somebody could quickly provide some information about translating coordinates provided by NOAA to the gridded output. Thanks to the suggestion here I was able to read the CMORPH data, now i'm wondering how to access the correct row and col index that corresponds to the correct latitude, longitude coordinates (e.g. from google maps). 

The rainfall data results in a matrix of 4948 x 1649 (lon x lat dimension), I had the following transformation in mind:

Latitude = 31.98;

Longitude = 35,92;

lat_index = round(abs(59.963614 - Latitude) / gridspacing_lat);

lon_index = round(abs(0.036378335 - Longitude) / gridspacing_lon);

Here i end up with longitude&latitude index 481, 385. If i check the precipitation with known stationdata it doesn't even closely correspond. Any idea where i'm going wrong?

Any help is much appreciated! 

pl. check whether it should be    -59.963614 - Latitude (latitude starting point with sign)

Thank you for your quick answer Bipasha.

I will try this, but because the CMORPH description states that: "The arrays are oriented from North to South, beginning from latitude 59.963614N and from West to EAST from

longitude 0.036378335E" 

So I assumed the upper left colum starts at 60 and the first row starts at 0.03, not sure why it would be otherwise...

Regards

Job

Hello All,

I tried Bipasha's code in Matlab to load 8 km-30 min CMORPH data but I am getting error message. Job: I also looked at your code and how were you able to read and import CMORPH data into Matlab? (i.e. the file you named All_files in previous Matlab News group).  Would you mind sharing that code as well?  I am new to this and greatly appreciate any suggestions.

@Liya : I am using the code for  reading the CMORPH data (2013, 2014 data). It works alright. The CMORPH file is in archived as bz2, which we have to extract first , using any unzip software bunzip / Winrar etc. Then the above  matlab code has to be run. 

@Liya, Great to hear others are working on it! I'm trying to work with the half-hourly data and extract per location (in Lat Lon extracted from google maps). For some reason the data (for example 2010) is not comperable with station data at different places in the world. And then i'm not talking about 10 - 40% wrong but not even near seasonal changes or monthly average values. I build a little bit further on @Bipasha her approach, hopefully this makes sense (see attached). 

Love to hear more experiences about extracting half hourly data per location. I'm still confused wether the correction from lat/lon to row/col index is correct! 

Good luck!

Thank you both! I am trying to improve on the codes shared here. I am having difficulty locating the most up to date metadata for the CMORPH data. I was able to find the data description from 2010 (ftp://ftp.cpc.ncep.noaa.gov/precip/global_CMORPH/README.cmorph.8km_30minute). but not sure if this still holds true. I am particularly worried about the precipitation values being multiplied by "0.2" and the difference b/n CRT and RAW data.

Hi Liya,

I'd be interested to see what you improved so far. Location access gives you estimates that coincides with for example station data? 

Unfortunately i'm not able to tell you if the description is up to date, i'm sure hoping it is ;)

I'm fairly confident about the scaling factor of 0.2 but agree about the binary data format versus the raw format. Even in FORTRAN it is not straightforward to access. 

Good luck!

Hi.

You can download the monthly data here...

https://climexp.knmi.nl/select.cgi?id=someone@somewhere&field=cmorph_monthly

regards

rodrigo

Great !!!!

Do you know the scaling factor to apply to get the data in mm ?

No. Did you find it ?

r