#==============Import COSMO Output Data=================

from qgis.core import *
import os.path, glob

#select variable
var="CO2_RA" #varibale name in NetCDF
path = 'Z:/burc/COSMO/20200627_04/*.nc' # Change this base path

#load all files in folder
for file in glob.glob(path):
    #string for NetCDF:
    uri = 'NETCDF:"' + file +'":'+ var
    #filename with date string
    fname = str(int(''.join(filter(str.isdigit, file))))
    if len(fname) > 10:
            fname = fname[10:]
    #load raster layer
    rlayer = QgsRasterLayer(uri,fname)
    if not rlayer.isValid():
        print("Layer failed to load!")
    else:
        QgsProject.instance().addMapLayer(rlayer)

#set coordinate reference system
    crs = QgsCoordinateReferenceSystem()
    #proj4string to define crs:
    crs.createFromProj("+proj=ob_tran +o_proj=longlat +o_lon_p=0 +o_lat_p=43 +lon_0=10")
    #apply crs on raster layer
    if not crs.isValid():
        print("CRS is not valid!")
    else:
        rlayer.setCrs(crs)
    QgsProject.instance().addMapLayer(rlayer)

#load style for variable
    style_uri = "Z:/burc/GIS/COSMO/"+var+"_style.qml"
    rlayer.loadNamedStyle(style_uri)