remo data analysis and comparison to observations

[larsbuntemeyer]

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]

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

[larsbuntemeyer]

• 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]

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

[larsbuntemeyer]

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

[larsbuntemeyer]

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

[ludwiglierhammer]

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