matlab .nc数据格式转tif格式

close all;
clear all;
clc

% ff = ‘R:\study\hhl\lst\lst\cruncepv5_tair_1980.nc‘;
% ncid = netcdf.open（ff‘NOWRITE‘）;
% [numdims numvars numglobalatts unlimdimID] = netcdf.inq（ncid）;
% for i = 0:numvars-1
%     [varname xtype dimids numatts] = netcdf.inqVar（ncidi）;
%     disp（[‘--------------------< ‘ varname ‘ >---------------------‘]）
%     flag = 0;
%     for j = 0:numatts - 1
%         attname1 = netcdf.inqAttName（ncidij）;
%         attname2 = netcdf.getAtt（ncidiattname1）;
%         disp（[attname1 ‘:  ‘ num2str（attname2）]）
%         if strcmp（‘add_offset‘attname1）
%             offset = attname2;
%             pause;
%         end
%         if strcmp（‘scale_factor‘attname1）
%             scale = attname2;
%             flag = 1;
%         end        
%     end
%     disp（‘ ‘）
% 
%     if flag
%         eval（[varname ‘= double（double（netcdf.getVar（ncidi））*scale + offset）;‘]）
%     else
%         eval（[varname ‘= double（netcdf.getVar（ncidi））;‘]）   
%     end
% end

% % varname2 = netcdf.getVar（ncid2）;
% varname3 = netcdf.getVar（ncid2）;




% inway = ‘H:\Hehuili\ERA\ERA\Tavg‘;
inway = ‘R:\study\hhl\precipitation\pre‘;

%out_way = ‘J:\HHL\CRU_NCEP\TIFF‘;
%mkdir（out_way）;
out_way = ‘R:\study\hhl\tif\precipitation‘;
mkdir（out_way）;

yr1 = 1980;
yr2 = 2002;


% coor_way = ‘H:\Hehuili\ERA\HDR‘;
% hdr = sprintf（‘%s\\ERA_T.hdr‘coor_way）;
% prj = sprintf（‘%s\\ERA_T.prj‘coor_way）;
month1 = [31 28 31 30 31 30 31 31 30 31 30 31];
% month = {‘01‘‘02‘‘03‘‘04‘‘05‘‘06‘‘07‘‘08‘‘09‘‘10‘‘11‘‘12‘‘13‘‘14‘‘15‘‘16‘
%     ‘17‘‘18‘‘18‘‘19‘‘20‘‘21‘‘22‘‘23‘‘24‘‘25‘‘26‘‘27‘‘28‘‘29‘‘30‘‘31‘‘32‘‘33‘
%     ‘34‘‘35‘‘36‘‘37‘‘38‘‘39‘‘40‘.‘41‘‘42‘‘43‘‘44‘‘45‘};
%days={}
% month=cell（112）;
% 
% for ii =1:12
%     month（ii）=num2cell（ii）;
%     
% end
month1=[4];
for year = yr1:yr2
%     for mm = 1:12
%         mm1 = char（month（mm））;
%         ff = sprintf（‘%s\\tavg_%d%s.nc‘inwayyear