MATLAB处理探空数据

zzw9701

请问有人处理过怀俄明大学的sounding数据吗？他的sounding数据文件大小不一致，有头文件和后面的文字描述文件，很不好批量处理。。。附件里的TXT只是一次数据。。。有没有人能讲一下怎么处理的。

croton

for isite = 1:length(site)
    disp(['===> processing ',site{isite}]);
    file = [input_dir,site{isite},'\',num2str(target_year,'%.4i'),'-',num2str(target_month,'%.2i'),'.txt'];
    fileID    = fopen(file,'rt');
    fseek(fileID,0,'eof');
    filelength = ftell(fileID);
    fclose(fileID);
   
    fileID    = fopen(file,'rt');
    clear data;
    for ifileline=1:filelength
        data{ifileline,1} = fgetl(fileID);
        if strfind(data{ifileline,1},[site_name{isite},' Observations'])
            data_beg = ifileline;
            temp     = strfind(data{ifileline,1},[site_name{isite},' Observations']);
            hour = data{ifileline,1}(temp+length([site_name{isite},' Observations'])+4:temp+length([site_name{isite},' Observations'])+5);
            day  = data{ifileline,1}(temp+length([site_name{isite},' Observations'])+8:temp+length([site_name{isite},' Observations'])+9);
        end
        if strfind(data{ifileline,1},'Station information and sounding indices')
            data_end = ifileline;
            fid=fopen([output_dir,site{isite},'\',num2str(target_year,'%.4i'),num2str(target_month,'%.2i'),day,hour,'.txt'],'wt');
            for idataline = data_beg+2:data_end-1
                temp = data{idataline,1};
                fprintf(fid,'%s',temp);
                fprintf(fid,'\n');
            end
            fclose(fid);
        end;
    end
   
%     while ~feof(fileID)
%         data{i,1} = fgetl(fileID);
%     end
    fclose(fileID);
end


之前自己处理时写的一段程序。
大致思路就是虽然数据里面有头尾，但是都还是规律的，可以按照站点编号（即程序中的site）来处理。
楼主可以参考一下

zzw9701

croton 发表于 2019-5-20 09:57
for isite = 1:length(site)
    disp(['===> processing ',site{isite}]);
    file = ;

太谢谢你了！！！

一颗醉甜的草莓

请问处理探空数据时 前几层超过1000hpa 是舍弃还是怎么处理呢

zzw9701

一颗醉甜的草莓 发表于 2020-4-8 09:02
请问处理探空数据时 前几层超过1000hpa 是舍弃还是怎么处理呢

我是继续保留的，如果你只需要1000hPa及以上的数据，可以不读

一颗醉甜的草莓

zzw9701 发表于 2020-4-8 10:47
我是继续保留的，如果你只需要1000hPa及以上的数据，可以不读

那一般计算对流参数 是直接从第一层算呢 还是从1000hpa算 比如cape

zzw9701

一颗醉甜的草莓 发表于 2020-4-8 14:39
那一般计算对流参数 是直接从第一层算呢 还是从1000hpa算 比如cape

我没算过CAPE值。。。

wwggiser

你好，我请假个问题，感谢了。我看怀俄明大学探空数据内显示：Precipitable water [mm] for entire sounding: 10.98。这个结果做水汽含量的验证就可以直接用了吗？

zzw9701

wwggiser 发表于 2021-4-10 16:07
你好，我请假个问题，感谢了。我看怀俄明大学探空数据内显示：Precipitable water [mm] for entire soundin ...

抱歉，我没用过这个值。。。

wwggiser

zzw9701 发表于 2021-4-11 20:55
抱歉，我没用过这个值。。。

感谢你及时回复，辛苦了。我再找找其他资料。