[ncl-talk] HOW TO READ TRMM 3B42 data?

[闯关东]

Dear Sir,

I am reading the TRMM 3B42RT data with ncl, the data is writen in a direct format with only one record. But in this record, there are a header, 3 short integer and 1 byte arrays.


I am trying to read it with fbindirread, but the data is not OK.


The fortran code to read the file is as following:


            READ ( iunit, REC=1, ERR=100, IOSTAT=iostat ) head (1:2880),
     +          ( ( i2p (i, j), i = 1, LN ), j = 1, LT ),
     +          ( ( i2e (i, j), i = 1, LN ), j = 1, LT ),
     +          ( ( np1 (i, j), i = 1, LN ), j = 1, LT ),
     +          ( ( i2u (i, j), i = 1, LN ), j = 1, LT )



The NCL code i want to read the file is as followed,


head = fbindirread("3B42RT.2013071318.7.bin",0,(/2880/),"Character")
trmm1 = fbindirread("3B42RT.2013071318.7.bin",1,(/480,1440/),"short")
trmm2 = fbindirread("3B42RT.2013071318.7.bin",2,(/480,1440/),"short")

trmm3 = fbindirread("3B42RT.2013071318.7.bin",3,(/480,1440/),"byte")
trmm4 = fbindirread("3B42RT.2013071318.7.bin",4,(/480,1440/),"short")


If i set the parameter "rec_num" 0,1,2,3,4 all to 0, the results is also wrong.







How to read this file? Attached are the fortran code and a sample data.

I am waiting for your quick reply. Thanks a lot.

--
yours,
Best regards,

Yuchun
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150710/33edf001/attachment.html