Global model data for "Revealing uncertainties in defining the sensitivity of cloud formation to aerosol changes" by Virtanen A., et al. (2024)

The data includes simulated global model data sampled at the measurement stations. ECHAM-HAM, ECHAM-SALSA, and UKESM data is in CSV text format, NorESM data is in NetCDF format. Python scripts cumulative_M7.py and cumulative_SALSA.py calculate cumulative historgrams for activation dry diameters for ECHAM-HAM and ECHAM-M7.

Usage of Python scripts:
cumulative_M7.py
cumulative_SALSA.py

Puijo_nearest.csv includes UKESM data for:
Dry diameter of a mode, column name: ddry_, unit m
Number concentration of a mode, column name: n_, unit m^-3
Cloud droplet number concentration, column name: CDNC, unit m^-3
Activation supersaturation, Act. SupSat, unit %
Activation vertical velocity, column name: Vert. Vel., unit m/s
Cloud fraction, column name: Cloud Frac., unitless
Air temperature, column name: Air Temp, unit K

Puijo_nearest_rCrit_correct_NEW.csv includes UKESM data for
Critical dry radius of a mode, column name: rcrit_, unit m

M7_.csv and SALSA_ include data for for
Number concentration particle in a mode or a size class larger than , column name: CCN , unit cm^-3
Cloud fraction, column name: Cloud Frac., unitless
Activation vertical velocity, column name: Act. udr, unit m/s
Air temperature, column name: Air temp, unit K
Maximum supersaturation of air in cloud activation, column name: Max SS
Dry diameter of a mode or a size class, column name: Ddry , unit m
Number concentration of a mode or a size class, column name: N , unit cm^-3
Critical dry radius of a mode or a size class, column name: rcrit , unit m
Critical supersaturation of a mode or a size class, column name: scrit , unitless

OsloAero_f19_f19_mg17_act_corr_height_2009-01-01-2011-01-01_w_act_diam.nc includes the following data:
dimensions:
time = 5841 ;
station = 5 ;
mode = 11 ;
variables:
int64 time(time) ;
time:units = "hours since 2009-01-01 00:00:00" ;
time:calendar = "proleptic_gregorian" ;
string station(station) ;
float ACTNL(station, time) ;
ACTNL:_FillValue = NaNf ;
ACTNL:units = "#/cm^3" ;
ACTNL:long_name = "Average Cloud Top droplet number" ;
float ACTREI(station, time) ;
ACTREI:_FillValue = NaNf ;
ACTREI:units = "Micron" ;
ACTREI:long_name = "Average Cloud Top ice effective radius" ;
float ACTREL(station, time) ;
ACTREL:_FillValue = NaNf ;
ACTREL:units = "Micron" ;
ACTREL:long_name = "Average Cloud Top droplet effective radius" ;
float AREL(station, time) ;
AREL:_FillValue = NaNf ;
AREL:mdims = 1 ;
AREL:units = "Micron" ;
AREL:long_name = "Average droplet effective radius" ;
float AWNC(station, time) ;
AWNC:_FillValue = NaNf ;
AWNC:mdims = 1 ;
AWNC:units = "#/cm^3" ;
AWNC:long_name = "Average cloud water number conc" ;
float CLDTOT(station, time) ;
CLDTOT:_FillValue = NaNf ;
CLDTOT:units = "fraction" ;
CLDTOT:long_name = "Vertically-integrated total cloud" ;
float CLOUD(station, time) ;
CLOUD:_FillValue = NaNf ;
CLOUD:mdims = 1 ;
CLOUD:units = "fraction" ;
CLOUD:long_name = "Cloud fraction" ;
float CLOUDFRAC_CLUBB(station, time) ;
CLOUDFRAC_CLUBB:_FillValue = NaNf ;
CLOUDFRAC_CLUBB:mdims = 1 ;
CLOUDFRAC_CLUBB:units = "fraction" ;
CLOUDFRAC_CLUBB:long_name = "Cloud Fraction" ;
float FCTI(station, time) ;
FCTI:_FillValue = NaNf ;
FCTI:units = "fraction" ;
FCTI:long_name = "Fractional occurrence of cloud top ice" ;
float FCTL(station, time) ;
FCTL:_FillValue = NaNf ;
FCTL:units = "fraction" ;
FCTL:long_name = "Fractional occurrence of cloud top liquid" ;
float FREQI(station, time) ;
FREQI:_FillValue = NaNf ;
FREQI:mdims = 1 ;
FREQI:units = "fraction" ;
FREQI:long_name = "Fractional occurrence of ice" ;
float FREQL(station, time) ;
FREQL:_FillValue = NaNf ;
FREQL:mdims = 1 ;
FREQL:units = "fraction" ;
FREQL:long_name = "Fractional occurrence of liquid" ;
float cloud_fraction(station, time) ;
cloud_fraction:_FillValue = NaNf ;
cloud_fraction:mdims = 1 ;
cloud_fraction:long_name = "Liquid cloud fraction used in stratus activation" ;
cloud_fraction:units = "1" ;
double n100(station, time) ;
n100:_FillValue = NaN ;
n100:mdims = 1 ;
n100:units = "m-3" ;
n100:long_name = "N$_{100-10000}$" ;
n100:raw_data_path = "/proj/bolinc/users/x_sarbl/noresm_archive" ;
n100:model = "NorESM" ;
n100:model_name = "NorESM" ;
n100:case_name = "OsloAero_f19_f19_mg17_act" ;
n100:case = "OsloAero_f19_f19_mg17_act" ;
n100:case_name_nice = "OsloAero_f19_f19_mg17_act" ;
n100:isSectional = "False" ;
n100:from_time = "2009-01-01" ;
n100:to_time = "2011-01-01" ;
n100:time_resolution = "hour" ;
n100:history_field = ".h1." ;
n100:pressure_coords = "False" ;
double n30(station, time) ;
n30:_FillValue = NaN ;
n30:mdims = 1 ;
n30:units = "m-3" ;
n30:long_name = "N$_{30-10000}$" ;
n30:raw_data_path = "/proj/bolinc/users/x_sarbl/noresm_archive" ;
n30:model = "NorESM" ;
n30:model_name = "NorESM" ;
n30:case_name = "OsloAero_f19_f19_mg17_act" ;
n30:case = "OsloAero_f19_f19_mg17_act" ;
n30:case_name_nice = "OsloAero_f19_f19_mg17_act" ;
n30:isSectional = "False" ;
n30:from_time = "2009-01-01" ;
n30:to_time = "2011-01-01" ;
n30:time_resolution = "hour" ;
n30:history_field = ".h1." ;
n30:pressure_coords = "False" ;
double n50(station, time) ;
n50:_FillValue = NaN ;
n50:mdims = 1 ;
n50:units = "m-3" ;
n50:long_name = "N$_{50-10000}$" ;
n50:raw_data_path = "/proj/bolinc/users/x_sarbl/noresm_archive" ;
n50:model = "NorESM" ;
n50:model_name = "NorESM" ;
n50:case_name = "OsloAero_f19_f19_mg17_act" ;
n50:case = "OsloAero_f19_f19_mg17_act" ;
n50:case_name_nice = "OsloAero_f19_f19_mg17_act" ;
n50:isSectional = "False" ;
n50:from_time = "2009-01-01" ;
n50:to_time = "2011-01-01" ;
n50:time_resolution = "hour" ;
n50:history_field = ".h1." ;
n50:pressure_coords = "False" ;
double n70(station, time) ;
n70:_FillValue = NaN ;
n70:mdims = 1 ;
n70:units = "m-3" ;
n70:long_name = "N$_{70-10000}$" ;
n70:raw_data_path = "/proj/bolinc/users/x_sarbl/noresm_archive" ;
n70:model = "NorESM" ;
n70:model_name = "NorESM" ;
n70:case_name = "OsloAero_f19_f19_mg17_act" ;
n70:case = "OsloAero_f19_f19_mg17_act" ;
n70:case_name_nice = "OsloAero_f19_f19_mg17_act" ;
n70:isSectional = "False" ;
n70:from_time = "2009-01-01" ;
n70:to_time = "2011-01-01" ;
n70:time_resolution = "hour" ;
n70:history_field = ".h1." ;
n70:pressure_coords = "False" ;
float NACT_FRAC01(station, time) ;
NACT_FRAC01:_FillValue = NaNf ;
NACT_FRAC01:mdims = 1 ;
NACT_FRAC01:units = "%" ;
NACT_FRAC01:long_name = "number fraction of aerosols activated 01" ;
float NACT_FRAC02(station, time) ;
NACT_FRAC02:_FillValue = NaNf ;
NACT_FRAC02:mdims = 1 ;
NACT_FRAC02:units = "%" ;
NACT_FRAC02:long_name = "number fraction of aerosols activated 02" ;
float NACT_FRAC03(station, time) ;
NACT_FRAC03:_FillValue = NaNf ;
NACT_FRAC03:mdims = 1 ;
NACT_FRAC03:units = "%" ;
NACT_FRAC03:long_name = "number fraction of aerosols activated 03" ;
float NACT_FRAC04(station, time) ;
NACT_FRAC04:_FillValue = NaNf ;
NACT_FRAC04:mdims = 1 ;
NACT_FRAC04:units = "%" ;
NACT_FRAC04:long_name = "number fraction of aerosols activated 04" ;
float NACT_FRAC05(station, time) ;
NACT_FRAC05:_FillValue = NaNf ;
NACT_FRAC05:mdims = 1 ;
NACT_FRAC05:units = "%" ;
NACT_FRAC05:long_name = "number fraction of aerosols activated 05" ;
float NACT_FRAC06(station, time) ;
NACT_FRAC06:_FillValue = NaNf ;
NACT_FRAC06:mdims = 1 ;
NACT_FRAC06:units = "%" ;
NACT_FRAC06:long_name = "number fraction of aerosols activated 06" ;
float NACT_FRAC07(station, time) ;
NACT_FRAC07:_FillValue = NaNf ;
NACT_FRAC07:mdims = 1 ;
NACT_FRAC07:units = "%" ;
NACT_FRAC07:long_name = "number fraction of aerosols activated 07" ;
float NACT_FRAC08(station, time) ;
NACT_FRAC08:_FillValue = NaNf ;
NACT_FRAC08:mdims = 1 ;
NACT_FRAC08:units = "%" ;
NACT_FRAC08:long_name = "number fraction of aerosols activated 08" ;
float NACT_FRAC09(station, time) ;
NACT_FRAC09:_FillValue = NaNf ;
NACT_FRAC09:mdims = 1 ;
NACT_FRAC09:units = "%" ;
NACT_FRAC09:long_name = "number fraction of aerosols activated 09" ;
float NACT_FRAC10(station, time) ;
NACT_FRAC10:_FillValue = NaNf ;
NACT_FRAC10:mdims = 1 ;
NACT_FRAC10:units = "%" ;
NACT_FRAC10:long_name = "number fraction of aerosols activated 10" ;
float NACT_FRAC11(station, time) ;
NACT_FRAC11:_FillValue = NaNf ;
NACT_FRAC11:mdims = 1 ;
NACT_FRAC11:units = "%" ;
NACT_FRAC11:long_name = "number fraction of aerosols activated 11" ;
float NACT_FRAC12(station, time) ;
NACT_FRAC12:_FillValue = NaNf ;
NACT_FRAC12:mdims = 1 ;
NACT_FRAC12:units = "%" ;
NACT_FRAC12:long_name = "number fraction of aerosols activated 12" ;
float NACT_FRAC13(station, time) ;
NACT_FRAC13:_FillValue = NaNf ;
NACT_FRAC13:mdims = 1 ;
NACT_FRAC13:units = "%" ;
NACT_FRAC13:long_name = "number fraction of aerosols activated 13" ;
float NACT_FRAC14(station, time) ;
NACT_FRAC14:_FillValue = NaNf ;
NACT_FRAC14:mdims = 1 ;
NACT_FRAC14:units = "%" ;
NACT_FRAC14:long_name = "number fraction of aerosols activated 14" ;
float NCONC00(station, time) ;
NCONC00:_FillValue = NaNf ;
NCONC00:mdims = 1 ;
NCONC00:units = "1/cm3" ;
NCONC00:long_name = "Aerosol mode 0 number concentration" ;
float NCONC01(station, time) ;
NCONC01:_FillValue = NaNf ;
NCONC01:mdims = 1 ;
NCONC01:units = "#/cm3" ;
NCONC01:long_name = "number concentration mode 01" ;
float NCONC02(station, time) ;
NCONC02:_FillValue = NaNf ;
NCONC02:mdims = 1 ;
NCONC02:units = "#/cm3" ;
NCONC02:long_name = "number concentration mode 02" ;
float NCONC04(station, time) ;
NCONC04:_FillValue = NaNf ;
NCONC04:mdims = 1 ;
NCONC04:units = "#/cm3" ;
NCONC04:long_name = "number concentration mode 04" ;
float NCONC05(station, time) ;
NCONC05:_FillValue = NaNf ;
NCONC05:mdims = 1 ;
NCONC05:units = "#/cm3" ;
NCONC05:long_name = "number concentration mode 05" ;
float NCONC06(station, time) ;
NCONC06:_FillValue = NaNf ;
NCONC06:mdims = 1 ;
NCONC06:units = "#/cm3" ;
NCONC06:long_name = "number concentration mode 06" ;
float NCONC07(station, time) ;
NCONC07:_FillValue = NaNf ;
NCONC07:mdims = 1 ;
NCONC07:units = "#/cm3" ;
NCONC07:long_name = "number concentration mode 07" ;
float NCONC08(station, time) ;
NCONC08:_FillValue = NaNf ;
NCONC08:mdims = 1 ;
NCONC08:units = "#/cm3" ;
NCONC08:long_name = "number concentration mode 08" ;
float NCONC09(station, time) ;
NCONC09:_FillValue = NaNf ;
NCONC09:mdims = 1 ;
NCONC09:units = "#/cm3" ;
NCONC09:long_name = "number concentration mode 09" ;
float NCONC10(station, time) ;
NCONC10:_FillValue = NaNf ;
NCONC10:mdims = 1 ;
NCONC10:units = "#/cm3" ;
NCONC10:long_name = "number concentration mode 10" ;
float NCONC12(station, time) ;
NCONC12:_FillValue = NaNf ;
NCONC12:mdims = 1 ;
NCONC12:units = "#/cm3" ;
NCONC12:long_name = "number concentration mode 12" ;
float NCONC14(station, time) ;
NCONC14:_FillValue = NaNf ;
NCONC14:mdims = 1 ;
NCONC14:units = "#/cm3" ;
NCONC14:long_name = "number concentration mode 14" ;
float NMR00(station, time) ;
NMR00:_FillValue = NaNf ;
NMR00:mdims = 1 ;
NMR00:units = "nm" ;
NMR00:long_name = "Aerosol mode 0 number concentration" ;
float NMR01(station, time) ;
NMR01:_FillValue = NaNf ;
NMR01:mdims = 1 ;
NMR01:units = "nm" ;
NMR01:long_name = "number median radius mode 01" ;
float NMR02(station, time) ;
NMR02:_FillValue = NaNf ;
NMR02:mdims = 1 ;
NMR02:units = "nm" ;
NMR02:long_name = "number median radius mode 02" ;
float NMR04(station, time) ;
NMR04:_FillValue = NaNf ;
NMR04:mdims = 1 ;
NMR04:units = "nm" ;
NMR04:long_name = "number median radius mode 04" ;
float NMR05(station, time) ;
NMR05:_FillValue = NaNf ;
NMR05:mdims = 1 ;
NMR05:units = "nm" ;
NMR05:long_name = "number median radius mode 05" ;
float NMR06(station, time) ;
NMR06:_FillValue = NaNf ;
NMR06:mdims = 1 ;
NMR06:units = "nm" ;
NMR06:long_name = "number median radius mode 06" ;
float NMR07(station, time) ;
NMR07:_FillValue = NaNf ;
NMR07:mdims = 1 ;
NMR07:units = "nm" ;
NMR07:long_name = "number median radius mode 07" ;
float NMR08(station, time) ;
NMR08:_FillValue = NaNf ;
NMR08:mdims = 1 ;
NMR08:units = "nm" ;
NMR08:long_name = "number median radius mode 08" ;
float NMR09(station, time) ;
NMR09:_FillValue = NaNf ;
NMR09:mdims = 1 ;
NMR09:units = "nm" ;
NMR09:long_name = "number median radius mode 09" ;
float NMR10(station, time) ;
NMR10:_FillValue = NaNf ;
NMR10:mdims = 1 ;
NMR10:units = "nm" ;
NMR10:long_name = "number median radius mode 10" ;
float NMR12(station, time) ;
NMR12:_FillValue = NaNf ;
NMR12:mdims = 1 ;
NMR12:units = "nm" ;
NMR12:long_name = "number median radius mode 12" ;
float NMR14(station, time) ;
NMR14:_FillValue = NaNf ;
NMR14:mdims = 1 ;
NMR14:units = "nm" ;
NMR14:long_name = "number median radius mode 14" ;
float SIGMA00(station, time) ;
SIGMA00:_FillValue = NaNf ;
SIGMA00:mdims = 1 ;
SIGMA00:units = "-" ;
SIGMA00:long_name = "Aerosol mode 0 number concentration" ;
float SIGMA01(station, time) ;
SIGMA01:_FillValue = NaNf ;
SIGMA01:mdims = 1 ;
SIGMA01:units = "-" ;
SIGMA01:long_name = "Std. dev. mode 01" ;
float SIGMA02(station, time) ;
SIGMA02:_FillValue = NaNf ;
SIGMA02:mdims = 1 ;
SIGMA02:units = "-" ;
SIGMA02:long_name = "Std. dev. mode 02" ;
float SIGMA04(station, time) ;
SIGMA04:_FillValue = NaNf ;
SIGMA04:mdims = 1 ;
SIGMA04:units = "-" ;
SIGMA04:long_name = "Std. dev. mode 04" ;
float SIGMA05(station, time) ;
SIGMA05:_FillValue = NaNf ;
SIGMA05:mdims = 1 ;
SIGMA05:units = "-" ;
SIGMA05:long_name = "Std. dev. mode 05" ;
float SIGMA06(station, time) ;
SIGMA06:_FillValue = NaNf ;
SIGMA06:mdims = 1 ;
SIGMA06:units = "-" ;
SIGMA06:long_name = "Std. dev. mode 06" ;
float SIGMA07(station, time) ;
SIGMA07:_FillValue = NaNf ;
SIGMA07:mdims = 1 ;
SIGMA07:units = "-" ;
SIGMA07:long_name = "Std. dev. mode 07" ;
float SIGMA08(station, time) ;
SIGMA08:_FillValue = NaNf ;
SIGMA08:mdims = 1 ;
SIGMA08:units = "-" ;
SIGMA08:long_name = "Std. dev. mode 08" ;
float SIGMA09(station, time) ;
SIGMA09:_FillValue = NaNf ;
SIGMA09:mdims = 1 ;
SIGMA09:units = "-" ;
SIGMA09:long_name = "Std. dev. mode 09" ;
float SIGMA10(station, time) ;
SIGMA10:_FillValue = NaNf ;
SIGMA10:mdims = 1 ;
SIGMA10:units = "-" ;
SIGMA10:long_name = "Std. dev. mode 10" ;
float SIGMA12(station, time) ;
SIGMA12:_FillValue = NaNf ;
SIGMA12:mdims = 1 ;
SIGMA12:units = "-" ;
SIGMA12:long_name = "Std. dev. mode 12" ;
float SIGMA14(station, time) ;
SIGMA14:_FillValue = NaNf ;
SIGMA14:mdims = 1 ;
SIGMA14:units = "-" ;
SIGMA14:long_name = "Std. dev. mode 14" ;
float activation_supersaturation(station, time) ;
activation_supersaturation:_FillValue = NaNf ;
activation_supersaturation:mdims = 1 ;
activation_supersaturation:units = "1" ;
activation_supersaturation:long_name = "Maximum supersaturation Abdul-Razzak in vertically changing clouds" ;
float Smax_cldv_supZero(station, time) ;
Smax_cldv_supZero:_FillValue = NaNf ;
Smax_cldv_supZero:mdims = 1 ;
Smax_cldv_supZero:units = "m" ;
Smax_cldv_supZero:long_name = "0 if Smax_cldv=0, 1 if Smax_cldv>0" ;
float Smax_incld(station, time) ;
Smax_incld:_FillValue = NaNf ;
Smax_incld:mdims = 1 ;
Smax_incld:units = "1" ;
Smax_incld:long_name = "Maximum supersaturation Abdul-Razzak in cloud growing or regrowing" ;
float Smax_incld_supZero(station, time) ;
Smax_incld_supZero:_FillValue = NaNf ;
Smax_incld_supZero:mdims = 1 ;
Smax_incld_supZero:units = "m" ;
Smax_incld_supZero:long_name = "0 if Smax_incld=0, 1 if Smax_incld>0" ;
float temperature(station, time) ;
temperature:_FillValue = NaNf ;
temperature:mdims = 1 ;
temperature:units = "K" ;
temperature:long_name = "Temperature" ;
float vertical_velocity(station, time) ;
vertical_velocity:_FillValue = NaNf ;
vertical_velocity:mdims = 1 ;
vertical_velocity:units = "m/s" ;
vertical_velocity:long_name = "Diagnostic sub-grid vertical velocity" ;
float WSUBI(station, time) ;
WSUBI:_FillValue = NaNf ;
WSUBI:mdims = 1 ;
WSUBI:units = "m/s" ;
WSUBI:long_name = "Diagnostic sub-grid vertical velocity for ice" ;
float cdnc(station, time) ;
cdnc:_FillValue = NaNf ;
cdnc:mdims = 1 ;
cdnc:units = "m-3" ;
cdnc:long_name = "Average cloud water number conc" ;
float cdnc_cloud_top(station, time) ;
cdnc_cloud_top:_FillValue = NaNf ;
cdnc_cloud_top:units = "m-3" ;
cdnc_cloud_top:long_name = "Average Cloud Top droplet number" ;
double act_d_mode_1(station, time) ;
act_d_mode_1:_FillValue = NaN ;
double act_d_mode_2(station, time) ;
act_d_mode_2:_FillValue = NaN ;
double act_d_mode_4(station, time) ;
act_d_mode_4:_FillValue = NaN ;
double act_d_mode_5(station, time) ;
act_d_mode_5:_FillValue = NaN ;
double act_d_mode_6(station, time) ;
act_d_mode_6:_FillValue = NaN ;
double act_d_mode_7(station, time) ;
act_d_mode_7:_FillValue = NaN ;
double act_d_mode_8(station, time) ;
act_d_mode_8:_FillValue = NaN ;
double act_d_mode_9(station, time) ;
act_d_mode_9:_FillValue = NaN ;
double act_d_mode_10(station, time) ;
act_d_mode_10:_FillValue = NaN ;
double act_d_mode_12(station, time) ;
act_d_mode_12:_FillValue = NaN ;
double act_d_mode_14(station, time) ;
act_d_mode_14:_FillValue = NaN ;
string mode(mode) ;
double act_diam_all(mode, station, time) ;
act_diam_all:_FillValue = NaN ;
double act_diam_min(station, time) ;
act_diam_min:_FillValue = NaN ;
act_diam_min:long_name = "activation_diameter_minimum" ;
act_diam_min:description = "The minimum activation diameter of the modes." ;
act_diam_min:units = "nm" ;