读TRMM的3B43和3B42可以使用相同的程序

《读TRMM的3B43和3B42可以使用相同的程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、读TRMM的3B43和3B42可以使用相同的程序，进行时间统计计算时注意它们是3小时和日值数据就行。附件3B43文档中有以IDL语言编写的读3B43的HDF文件的程序，用到3B42的HDF文件是没有问题的。建议先装IDL语言（软件），或ENVI（基于IDL平台，因此有IDL环境），然后把数据输出到TXT文件中（ASCII码格式），再用GRADS画图。或用IDL进行数据处理、画图，均可。SATE_TRMM_3B43_CASE.pro程序文件proSATE_TRMM_3B43_CASE;1*1440*400FILENAME='

2、SATE_L3_TRM_MUTDS_MWB_3B43_GLB_V6-199807.HDF'print,'*',FILENAMESD_ID=HDF_SD_START(FILENAME,/READ)P='precipitation';range:0-10,unit:mmhr-1SDS_INDEX_P=HDF_SD_NAMETOINDEX(SD_ID,P)SDS_ID_P=HDF_SD_SELECT(SD_ID,SDS_INDEX_P)HDF_SD_GETDATA,SDS_ID_P,DATA_PHDF_SD_ENDACCESS,

3、SDS_ID_PHDF_SD_END,SD_IDPRINT,DATA_P(0,0)print,size(DATA_P)OPENW,LUN,'0_50N_0_180E_199807_output.txt',/GET_LUN;0-50N,0-180EeastasiaFORJ=0,399DOBEGINFORI=0,1439DOBEGINIFJgt200THENBEGIN;0-50N(J>200)IFIgt720THENBEGIN;0-180E(1440

4、4.0-50.125,(I+1)/4.0-180.125,DATA_P(J,I),$format='(2(f7.3,2x),f7.3,2x)';lat,lon,precipitationENDIFENDIF;ELSEBEGIN;180-360E-->-180~0,,AmericaRINTF,LUN,(J+1)/4.0-50.125,(I+1)/4.0-180.125,DATA_P(J,I,0),$;format='(2(f6.2,2x),f6.4,2x)';ENDELSE;1-1440,I,-->0.125-359.875

5、,Longitude;1-720->(Europea,America),0.125-179.875W,Lon=I/4.0-180.125<0;721-1440->(Asia,Australia)0.125-179.875E,Lon=I/4.0-180.125<180;;1-400,J,-->49.875N-49.875S,Latitude;1-200->Lat=J/4.0-50.125<0,49.875S-0.125S;201-400->Lat=J/4.0-50.125>0,0.125N-49.875NENDFORENDF

6、ORCLOSE,LUNprint,'end'END附件中文档所带的FORTRAN程序未用过。但下面的程序读BIN格式的3B4X数据应该没问题。PROGRAMTRMMreal*4data1(1440,400)realdata(1440,400)reallon(1440)reallat(400)open(10,file='3B42.000101.00z.6.precipitation.bin',+form='unformatted',+access='DIRECT',status='OLD',recl=1440*400*4)o

7、pen(12,file='rain.out',form='unformatted',+access='DIRECT',status='unknown',recl=1440*400*4)cread(10,rec=1)((data(ii,jj),ii=1,1440),jj=1,400)cdo10ii=1,1440do10jj=1,400lon(ii)=-179.875+0.25*(ii-1)lat(jj)=-49.875+0.25*(jj-1)Cwrite(*,*)lon(ii),lat(jj),data(ii,jj)data

8、1(ii,jj)=data(ii,jj)10continuewrite(12,rec=1)((data1(ii,jj),ii=1,1440),jj=1,400)cclose(10)cend==========================================================