hdf格式卫星数据（modis数据）处理程序出错，希望哪位能帮忙看看哪儿出错了

最初的向往

用matlab处理modis卫星数据时，命令窗口显示程序第55行出错，55行的程序是：z=0.39log(COD)-0.01  我看了下下真不知道错哪儿了？求哪位大咖给看看指点一下。
下面是程序
clear all;close all;clc;
filename = 'MOD06_L2.A2017098.0435.006.2017098135158.hdf';


file_info = hdfinfo(filename);

lon = hdfread(filename,'Longitude');                            %经度
lat = hdfread(filename,'Latitude');                             %纬度
%----------------经纬度坐标插值--------------------
y=1:406;x=1:270;
yi=linspace(min(y),max(y),2030);
xi=linspace(min(x),max(x),1354);
[X,Y]=meshgrid(x,y);
[Xi,Yi]=meshgrid(xi,yi);
Z1=lon;
Longitude=interp2(X,Y,Z1,Xi,Yi,'spline');
Z2=lat;
Latitude =interp2(X,Y,Z2,Xi,Yi,'spline');
%---云顶温度,k----------------------------------------------
%temperature = hdfread(filename,'Cloud_Top_Temperature');        %云顶温度
temperature = hdfread(filename,'cloud_top_temperature_1km');    %云顶温度
scale_factor = 0.009999999776482582;
add_offset = -15000.0;
temp = (temperature - add_offset)*scale_factor;                 %温度数据转换,单位：K
% ---云顶高度,m-----------------------------------------------
% height=hdfread(filename,'Cloud_Top_Height');                  %云顶高度
height=hdfread(filename,'cloud_top_height_1km');  
scale_factor = 1.0;
add_offset = 0.0;
zt = (height - add_offset)*scale_factor;                          %高度转换，单位：m

%----有效粒子半径，单位：micron -------------------------------
radius=hdfread(filename,'Cloud_Effective_Radius')
scale_factor = 0.009999999776482582;
add_offset = 0.0;
re = (radius - add_offset)*scale_factor;                          %1*1，有效粒子半径，单位：micron
%-----------------------------------------------------
thickness=hdfread(filename,'Cloud_Optical_Thickness')
scale_factor = 0.009999999776482582;
add_offset = 0.0;
COD = (thickness - add_offset)*scale_factor;                      %1*1
%-freezing level,m-----------------------------------------------------
zfr=zt+(temp-273.15)/6.5
%-几何高度,km?----------------------------------------
z=0.39log(COD)-0.01
%-cloud base altitude-----------------------------------------
zb=zt-z

%-LWP,1*1,g/m^2------------------------------------------------
path=hdfread(filename,'Cloud_Water_Path');  
scale_factor = 1.0;
add_offset = 0.0;
lwp=(path - add_offset)*scale_factor;
if zb<zfr
    slwp=lwp(zt-zfr)/z
else
    slwp=lwp
end
% --------------------------------------------------
contourf(Longitude,Latitude,slwp);
%contourf(lon,lat,ht);
colorbar
axis([120 125 30 34 ])

最初的向往

自己顶一个

river

z=0.39log(COD)-0.01是不是少了一个 * ？
z=0.39*log(COD)-0.01
我不会matlab，我猜的，哈哈

愛、星晴

river 发表于 2017-10-16 18:03
z=0.39log(COD)-0.01是不是少了一个 * ？
z=0.39*log(COD)-0.01
我不会matlab，我猜的，哈哈

我觉得也是少了一个乘号

roby

看看，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，

゛绍〤z、

少了个乘号啊。。

最初的向往

愛、星晴 发表于 2017-10-17 05:02
我觉得也是少了一个乘号

我加了乘号，但是还是不行！还是显示”Error in modis_read (line 58)  z=0.39*log(COD)-0.01“。你看看还有其他错误吗？

最初的向往

゛绍〤z、 发表于 2017-10-17 07:50
少了个乘号啊。。

我加了乘号，但是还是不行！还是显示”Error in modis_read (line 58)  z=0.39*log(COD)-0.01“。要不你再给看看？

最初的向往

roby 发表于 2017-10-17 07:12
看看，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，

看出来啥毛病了吗？

最初的向往

快来人给看看吧！我一直在这里等你们