Recipe: Matlab program converts cubed-sphere to latlon

function [lons lats T]= CSnative(fname,varname,time,level)
   ncid=netcdf.open(fname,'NOWRITE');
   vinfo=ncinfo(fname,varname);
   S=size(vinfo.Size);
   ydim=vinfo.Size(1);
   xdim=vinfo.Size(2);
   nf=vinfo.Size(3);

   lats=ncread(fname,'lats');
   lons=ncread(fname,'lons');
   if( S(2)==5)
      T=ncread(name,varname,[1 1 1 level time],[ydim xdim nf level time]);
   else
      T=ncread(name,varname,[1 1 1 time],[ydim xdim nf time]);
   end
   netcdf.close(ncid);