remo data analysis and comparison to observations

larsbuntemeyer commented 3 years ago

How to compare to gridded observations

How to handle masks?

How do we handle masks. Do we use the original REMO land sea mask before regridding? How to avoid regridding artifacts on the coastlines when REMO has no sea topography?

When to do the regridding?

Right now, we did the regridding before any further computation. Can the regridding be done in the end, e.g., on the the seasonal means?

Height correction

How does the heigh correction work correctly, e.g., what is height1 and height2 exactly? is my documentation here correct?

def height_correction(height1, height2):
    """Returns height correction in [K].

    Citation: trocken - adiabatischer Gradient, ludwig fragen.... link?

    Parameters
    ------------

    height1: array like
        topography of model data [m]
    height2: array like
        topography of *true* data [m], e.g., observations

    Returns
    --------

    height correction: array like
        temperature height correction due to different topographies
    """
    return (height2 - height1) * 0.0065

see ecmwf report page 11

larsbuntemeyer commented 3 years ago

@KatharinaBuelow just to let you know about my questions here.

larsbuntemeyer commented 3 years ago

check soil spin upd on different regions, srex regions....
erst am ende remappen und höhenkorrektur am ende machen...
bei precipitation brauchen wir conservative remapping.
wir nutzen als gemeinsame maske die remo maske remapped auf die cru daten...
zusätzlich seasonal cycle vergleichen...
sst ansehen (code 54) in den forcing daten... keine missing values einlesen.
zusätlich mit EOBS Daten vergleichen....

larsbuntemeyer commented 3 years ago

we need bounds in remo output for conservative remapping, where do we get those?

larsbuntemeyer commented 3 years ago

prototype notebook here: https://nbviewer.jupyter.org/github/remo-rcm/pyremo/blob/master/notebooks/remo-dataset.ipynb

larsbuntemeyer commented 3 years ago

ludwigs prototype: /pf/g/g300100/remo-dataset.ipynb

ludwiglierhammer commented 3 years ago

Ich hab folgende Änderungen vorgenommen:

klimatologische Jahresgänge
Vergleich mit E-OBS und ERA5
Höhenkorrektur: (REMO - OBS) * 0.0065
Vergleiche: REMO - OBS (war vorher umgekehrt)
Reihenfolge: Mitteln, Remappen, Höhenkorrektur
REMO Landseemaske zum Maskieren

remo-rcm / pyremo