%Hallo Fellas
%1) Ersetzt mit 'Suchen und ersetzen' Filelocation+Filename mit dem 
%   Dateipfad und Dateinamen der zutreffenden NetCDF-Datei etwa sowas
%   " 'M:\zzz_Arbeit RLI\__Winlass\___DAten\wspd.50m.2010.nc' " , 
%2) Ersetz mit"Suchen und Ersetzen" TimeLine durch den Wert den ihr in der
%   time-Variablen findet (z-Achse)
%3) und mit "Suchen und Ersetzen" auch 'Laufvariable' duch die, die ihr
%   auslesen wollt,... siehe --> 'ncdisp('Filelocation+Filename','/')'

%nc Open
ncid=netcdf.open( 'Filelocation+Filename','NOWRITE')%Nowrite könnt ihr lassen.

%Visually inspect a NetCDF file, hiding the attributes.
ncdisp('Filelocation+Filename','/')% Vorwärts-Shlash steht für 'alle', könnt ihr lassen

%lesen der Daten
data=ncread ('Filelocation+Filename','Laufvariable'); %hier x,y,time

%dann um die Zeitachse zusammenzuschieben
ave=sum(data,3)/TimeLine % data,3 steht für bei data die 3 Achse


%  wspd          
%          Size:       187x148x52560
%          Dimensions: x,y,time
%          Datatype:   int16
%          Attributes:
%                      long_name     = '4xDaily wind speed'
%                      units         = 'm/s'
%                      valid_range   = [0.00e+000 2.00e+002]
%                      var_desc      = 'Wind Speed\n'
%                      scale_factor  = 0.01
%                      add_offset    = 100
%                      missing_value = 3.28e+004
%                      actual_range  = [0.00e+000 5.77e+001]

ave=ave*0.01+100  %also geteilt durch den Scale_factor + add_offset

%nun kann man ave in ein ASCII file übertragen und fertig
% der ASCII HEDAR solle so Ausssehen (ohne Leerzeichen, am Anfang)
%    ncols         Anzahl Reihen
%    nrows         Anzahl Spalten (ohne Hedar)
%    xllcorner     Startpunkt x-Achse
%    yllcorner     Startpunkt y-Achse
%    cellsize      0,00083333333333333   (das 1km, bei 20km halt 0,17142857)
%    NODATA_value  -9999
%dann eure daten ins ASCII-File schupsen und ab damit in ArcGIS, umwandeln in Rasterdateien etc.

%Als Beispiel:
%
%%nc Open
%ncid = netcdf.open( 'M:\zzz_Arbeit RLI\__Winlass\___DAten\wspd.50m.2010.nc','NOWRITE')
%
%%Visually inspect a NetCDF file, hiding the attributes.
%ncdisp('M:\zzz_Arbeit RLI\__Winlass\___DAten\Windspeed\wspd.50m.2010.nc','/')
%
%%lesen der Daten
%data=ncread ('M:\zzz_Arbeit RLI\__Winlass\___DAten\wspd.50m.2010.nc','wspd');
%
%ave=sum(data,3)/52560;
%ave=ave*0.01+100;
%
%das Simikolon nicht vergessen, sonst dauerst zu lang ;)
% und so habt ihr die NetCDF_Datein sauber ausgeleden.