TRMM数据校正代码matlab

clearclc
for year=1998:2016
    for  month=1:1
        if month<10
            monthstr=cat（2‘0‘num2str（month））;
        else
            monthstr=num2str（month）;
        end
        yearstr=num2str（year）;
        filename=cat（2‘3B43.‘yearstrmonthstr‘01.7.HDF‘）;
        precipitation = hdfread（filename‘/Grid/precipitation‘）;
        rain_90rot = rot90（precipitation）;
        [hw]=size（precipitation）;
        lat=[hw];
        lon=[hw];
        
        for i_lat = 1:w
            for j_lo