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

下载文档 降价啦

15
0
约1.51万字
约 12页
2017-09-01 发布于湖北
举报
版权申诉
保障服务

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

1、本文档共12页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

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

读TRMM的3B43和3B42可以使用相同的程序，进行时间统计计算时注意它们是3小时和日值数据就行。附件3B43文档中有以IDL语言编写的读3B43的HDF文件的程序，用到3B42的HDF文件是没有问题的。建议先装IDL语言（软件），或ENVI（基于IDL平台，因此有IDL环境），然后把数据输出到TXT文件中（ASCII码格式），再用GRADS画图。或用IDL进行数据处理、画图，均可。SATE_TRMM_3B43_CASE.pro程序文件pro SATE_TRMM_3B43_CASE; 1*1440*400FILENAME=SATE_L3_TRM_MUTDS_MWB_3B43_GLB_V6-199807.HDFprint,*,FILENAMESD_ID=HDF_SD_START(FILENAME,/READ)P=precipitation ; range:0-10, unit:mm hr-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,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-180E east asiaFOR J=0,399 DO BEGINFOR I=0,1439 DO BEGINIF J gt 200 THEN BEGIN ;0-50N (J200) IF I gt 720 THEN BEGIN ;0-180E (1440I720), asia, Australia PRINTF,LUN,(J+1)/4.0-50.125,(I+1)/4.0-180.125,DATA_P(J,I), $ format=(2(f7.3,2x),f7.3,2x) ; lat, lon, precipitation ENDIFENDIF ;ELSE BEGIN ; 180-360E -- -180~0,, America RINTF,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, Longitude;1-720 - (Europea, America), 0.125-179.875W, Lon=I/4.0-180.1250; 721-1440 - (Asia, Australia)0.125-179.875E, Lon=I/4.0-180.125180;;1-400, J, --49.875N-49.875S, Latitude;1-200 -Lat=J/4.0-50.1250, 49.875S-0.125S; 201-400 - Lat=J/4.0-50.1250, 0.125N-49.875NENDFORENDFORCLOSE,LUNprint,endEND附件中文档所带的FORTRAN程序未用过。但下面的程序读BIN格式的3B4X数据应该没问题。PROGRAM TRMMreal*4 data1(1440,400)real data(1440,400)real lon(1440)real lat(400)open(10,file=3B42.000101.00z.6.precipitation.bin,+ form=unformatted,+ access=DIRECT,status=OLD,recl=1440*400*4)open(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)cdo 10 ii=1,1440do 10 jj=1,400lon(ii) = -179.875+0.25*(ii-1)lat(jj) = -49.875+0.25*(jj-1)C write(*