[ERROR] sign error in rsutaf and rsutafcs

[DirkOlivie]

Describe the error
There is an error in the fields rsutaf and rsutafcs. They should have the inverse sign from what they currently have.

Way to correct
Is there a place where we can share this information for people using NorESM2 data from ESGF? (saying these two fields should be multiplied by -1 to have the correct physical meaning)

At a later state, we can maybe recmorize these two fields, but that can/should possibly wait.

[YanchunHe]

Hi Dirk, does this apply for all datasets for rsutaf and rsutafcs, right?

@monsieuralok should be able to report an errtum at ESGF doc.

[DirkOlivie]

Hi Yanchun, yes, it applies to all datasets.

[adagj]

Hi,
I made a new issue, but now I see that this issue already exists so I'll add here and close the other one.
I don't agree that -1 should be multiplied to the fields, I think the problem rater is related to that the cmorized AERmon variables rsutaf and rsutcsaf are the difference from rsdt and not the actual variables. E.g. rsutaf (cmorized in noresm) => rsdt - rsutaf

This is different from what is outputted in other models and not described in the netcdf file

rsutaf:standard_name = "toa_outgoing_shortwave_flux" ;
rsutaf:long_name = "TOA Outgoing Aerosol-Free Shortwave Radiation" ;
rsutaf:comment = "Flux corresponding to rsut resulting from aerosol-free call to radiation, following Ghan (ACP, 2013)" ;

Way to correct:
To cmorize the actual field and not the difference from rsdt

[YanchunHe]

Hi @adagj, what is the actually fields name for these two variables rsutaf and rsutcsaf?

Are there some conclusion for the solution?

[adagj]

Hi @YanchunHe ,
I think the best solution is to retract the files and publish new files. The cmorized rsutaf should just be the outgoing shortwave flux for aerosol free air, and not the net radiation minus the outgoing shortwave flux for aerosol free air.

From the NorESM2 rawdata:

• The cmorized rsutaf says
  rsutaf:original_name = "FSNT_DRF" ;
  rsutaf:original_units = "W/m^2" ;
  rsutaf:standard_name = "toa_outgoing_shortwave_flux" ;
  rsutaf:long_name = "TOA Outgoing Aerosol-Free Shortwave Radiation" ;
  rsutaf:comment = "Flux corresponding to rsut resulting from aerosol-free call to radiation, following Ghan (ACP, 2013)" ;

• The FSNT_DRF from NorESM2 raw data
  FSNT_DRF:units = "W/m^2" ;
  FSNT_DRF:long_name = "Total column absorbed solar flux (DIRind)" ;

And the total column absorbed solar flux is not the same as TOA outgoing aerosol-free shortwave radiation.
I think what should be cmorized is the difference from the net, so you should cmorize FSNT_DRF - FSNT, and the same with rsutafcs -> cmorize FSNTCDRF - FSNTC

But I think we need @DirkOlivie or @oyvindseland to approve that is the correct method :-)

[adagj]

Hi @YanchunHe .
The calculation turned out a bit more complicated than first assumed, at least by me. The cmorized variables should be top-of-atmosphere fluxes, and not top-of-model fluxes which I think FSNT_DRF is, even if the long name suggests differently. So the question on how to deal with these variables remain unanswered and I don't know what is the best method. I don't know how to calculate the TOA fluxes from the TOM output we have available at the moment, but maybe @DirkOlivie or @oyvindseland have a good suggestion...

You can also see a relevant discussion started by Dirk here: NorESMhub/NorESM#418

[adagj]

@DirkOlivie ,
I checked what variables CESM2 uses to calculate TOA fluxes, but unfortunately the original name is not provided. However, CMCC uses CAM, as far as I understand, and uses different (cmor) calculations than us (however I can not find any aerosol free variables, but I checked clear sky)

• rsutcs (TOA Outgoing Clear-Sky Shortwave Radiation)
  CMCC-CM2-SR5: rsutcs:original_name = "FSUTOA+SWCF"
  NorESM2-LM: rsutcs:original_name = "SOLIN-FSNTOAC"

• rlutcs (TOA Outgoing Clear-Sky Longwave Radiation)
  CMCC-CM2-SR5: rlutcs:original_name = "FLUTC"
  NorESM2-LM: rlutcs:original_name = "FSNTOAC-FSNTC+FLNTC"

[YanchunHe]

Any conclusion here?

[YanchunHe]

Could any one summerize what to do for this?

Retract the data?
any fix for the cmorization?

[YanchunHe]

So, are the following expression correct?

rsutaf  = FSNT_DRF-FSNT
rsutcsaf = FSNTCDRF-FSNTC

[YanchunHe]

Please retract all rsutaf and rsutcsaf, @monsieuralok

[monsieuralok]

@YanchunHe should we remove rsutaf and rsutcsaf from all experiments, also from LM and MM and all versions?

[YanchunHe]

@YanchunHe should we remove rsutaf and rsutcsaf from all experiments, also from LM and MM and all versions?

I think so, thanks!

[adagj]

@YanchunHe @monsieuralok @DirkOlivie

Hi,
I had a brief discussion with Dirk and we suggest not to retract the data. The data is not wrong, if used correct. But we probably should add a description or an errata and explain how the data should be used (multiplied with -1 and only look at differences from the control simulation). Do you know where we can add such a description?

Best,
Ada

[monsieuralok]

@adagj We need to add description in errata and provide list of all affected datasets. @adagj could you add exact description what should we mention? if I am correct you want to say rsutaf = rsutaf -1 and rsutcsaf=rsutcsaf-1 . For control simulations data are correct or ..? or comparing control simulations from other data providers?

[adagj]

Hi,
I think it is best if @DirkOlivie write such a description since he knows these data best. I can provide a list of all affected datasets

[adagj]

The cases are:

NorESM2-MM :

• 1pctCO2
• abrupt-4xCO2
• historical (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• piControl
• piClim-4xCO2
• piClim-aer
• piClim-control
• ssp126
• ssp245 (r1i1p1f1, r2i1p1f1)
• ssp370
• ssp585

and for NorESM2-LM :

• hist-piAer (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• hist-piNTCF (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• histSST-piAer (r1i1p1f1, r1i1p2f1)
• histSST-piNTCF (r1i1p1f1, r1i1p2f1)
• histSST (r1i1p1f1, r1i1p2f1)
• piClim-2xDMS (r1i1p1f1, r1i1p2f1)
• piClim-2xdust (r1i1p1f1, r1i1p2f1)
• piClim-2xss (r1i1p1f1, r1i1p2f1)
• piClim-2xVOC (r1i1p1f1, r1i1p2f1)
• piClim-BC (r1i1p1f1, r1i1p2f1)
• piClim-CH4 (r1i1p1f1, r1i1p2f1)
• piClim-N2O (r1i1p1f1, r1i1p2f1)
• piClim-NTCF (r1i1p2f1)
• piClim-O3 ( r1i1p1f1, r1i1p2f1)
• piClim-OC (r1i1p1f1, r1i1p2f1)
• piClim-SO2 (r1i1p1f1, r1i1p2f1)
• ssp370-lowNTCF (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• ssp370SST-lowAer (r1i1p1f1)
• ssp370SST (r1i1p1f1)
• 1pctCO2-bgc (r1i1p1f1)
• 1pctCO2-rad (r1i1p1f1)
• esm-1pct-brch-1000PgC (r1i1p1f1)
• hist-bgc (r1i1p1f1)
• 1pctCO2-cdr (r1i1p1f1)
• esm-pi-cdr-pulse (r1i1p1f1)
• esm-pi-CO2pulse (r1i1p1f1)
• amip-4xCO2 (r1i1p2f1)
• amip-p4K (r1i1p2f1)
• 1pctCO2 (r1i1p1f1)
• abrupt-4xCO2 (r1i1p1f1)
• amip (r1i1p1f1, r1i1p2f1)
• esm-hist (r2i1p1f1)
• esm-piControl (r1i1p1f1)
• historical (r1i1p1f1, r1i1p4f1, r2i1p1f1, r3i1p1f1)
• piControl (r1i1p1f1)
• hist-aer (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• hist-GHG (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• hist-nat (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• ssp245-aer (r1i1p1f1)
• ssp245-GHG (r1i1p1f1)
• ssp245-nat (r1i1p1f1)
• hist-spAer-all (r1i1p1f1, r2i1p1f1)
• piClim-4xCO2 (r1i1p1f1, r1i1p2f1)
• piClim-aer (r1i1p1f1, r1i1p2f1)
• piClim-anthro (r1i1p1f1, r1i1p2f1)
• piClim-control (r1i1p1f1, r1i1p2f1)
• piClim-ghg (r1i1p1f1, r1i1p2f1)
• piClim-histaer (r1i1p1f1, r1i1p2f1, r2i1p1f1, r2i1p2f1, r3i1p1f1, r3i1p2f1)
• piClim-histall (r1i1p1f1, r1i1p2f1, r2i1p1f1, r2i1p2f1, r3i1p1f1, r3i1p2f1)
• piClim-histghg (r1i1p1f1, r1i1p2f1, r2i1p1f1, r2i1p2f1, r3i1p1f1, r3i1p2f1)
• piClim-histnat (r1i1p1f1, r1i1p2f1, r2i1p1f1, r2i1p2f1, r3i1p1f1, r3i1p2f1)
• piClim-lu (r1i1p1f1, r1i1p2f1)
• piClim-spAer-aer (r1i1p1f1, r1i1p2f1)
• piClim-spAer-anthro (r1i1p1f1, r1i1p2f1)
• piClim-spAer-histaer (r1i1p2f1)
• ssp126 (r1i1p1f1)
• ssp245 (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• ssp370 (r1i1p1f1, r2i1p1f1, r3i1p1f1)
• ssp585 (r1i1p1f1)

[monsieuralok]

@DirkOlivie @YanchunHe Please could you provide description for these datasets only mentioning
calculate rsutaf and rsutcsaf like below would be sufficient
rsutaf = FSNT_DRF-FSNT
rsutcsaf = FSNTCDRF-FSNTC

or it should be rsutaf = (rsutaf)x-1 and rsutcsaf=(rsutcsaf)x-1