NorESMhub / noresm2cmor

A command line tool for cmorizing NorESM output
http://noresmhub.github.io/noresm2cmor/
5 stars 16 forks source link

[ERROR] in wap ssp585 #316

Closed YanchunHe closed 1 year ago

YanchunHe commented 2 years ago

Describe the error Quote from Birthe

the wap is given on 4 Plevels, but all the values in the vertical are equal.

Way to correct

Is it possible to get wap data for this scenario that is more reasonable? If not,, is there any reasonable wap fields for another model simulation on 6 hour time resolution?

YanchunHe commented 2 years ago

I can't track where is the cause of the problem now. but I redo the cmorization and it seems OK.

The fixed dataset is here: /projects/NS9034K/CMIP6/ScenarioMIP/NCC/NorESM2-MM/ssp585/r1i1p1f1/6hrPlev/wap/gn/v20220505

I also noticed that all (just check randomly) the 3D 6hrPlev data are wrong with identicial data on plevs.

This applies for wap and hus.

This should apply for the experiments with 6hrPlev dataset available, including:

NorESM2-LM:

NorESM2-MM:

(a late update on this, 20221130: weired that only the cmorized NorESM2-MM ssp585 has different values in Plev, but others still have the same value on Plevs).

YanchunHe commented 1 year ago

Just find that 'PS' are missing in the 6 hourly averaged files, i.e., the cam.h2 streams.

It is only included in the SSP585 with high frequency output run, NSSP585frc2_f09_tn14_20200919, for FLEXPART by Dirk. (not in the previous NSSP585frc2_f09_tn14_20191105), see https://github.com/NorESMhub/noresm2cmor/issues/151

All other NorESM2-MM, LM runs miss PS. (not check all, but some of them)

That is why all the PS has been setting as '0' and gives wrong interpolation from hybrid to pressure levels.

However, PS does available in the instant 6 hourly data (cam.h3).

We may use PS from 6 hourly instant frequency as a workaround?

How do you think? @oyvindseland @DirkOlivie @matsbn

monsieuralok commented 1 year ago

@YanchunHe I have published errata issue on that : -https://errata.es-doc.org/static/view.html?uid=590fbf04-3fca-8d3d-a286-d7876f4b2970 datasets are not removed yet as again I am having problem with ESGF node.

monsieuralok commented 1 year ago

@YanchunHe retracted

YanchunHe commented 1 year ago

Thanks, @monsieuralok !

I keep this open, since the new data has not been generated and republished yet, so it is not closed.

oyvindseland commented 1 year ago

Why do you want to calculate additional fields than we have published though. The 3 3D 6 hour average fields we have output Q,Omega and cloud volume do vary with height in the output (checked 3 random files) Any additional fields should be calculated by the user I think

YanchunHe commented 1 year ago

The 3D 6 hourly output, Q, Omega are all right in the original output on sigma levels, but the CMIP6 requires them on interpolated 4 Plevs, corresponding to 'hus' and 'wap'. To calculate the latter with Q and Omega on original sigma levels, we need PS data, which is missing.

Only 'hus' is published on simga levels, 'wap' is not required and thus not cmorized and published.

Fine of me if you don't think it is necessary to calculate these additional fields.

Then please close this issue.

YanchunHe commented 1 year ago

OK, according to the discussion. I close this issue now. reopen when needed.