Check my LES run - 2D output

Check my LES run - 2D output#

The below notebook allows to compare a selected set of large-eddy simulation.
All requested 2D variables are here read and plotted. “load_sims_2d” reports variable that are missing. Maps will not show, if a variable is missing or if it contains missing values.
In case you find inconsistencies in the benchmark set or wish to expand the code, please contact Ann Fridlind (ann.fridlind@nasa.gov), Timothy Juliano (tjuliano@ucar.edu), and Florian Tornow (ft2544@columbia.edu).

%run functions_plotting.py

## select LES run(s) of interest
sim_keyword = 'FixN_2'  

## select specfic time
time_vec = [18.0]

vars_mean_list = pd.read_excel('https://docs.google.com/spreadsheets/d/1Vl8jYGviet7EtXZuQiitrx4NSkV1x27aJAhxxjBb9zI/export?gid=1756539842&format=xlsx',
                              sheet_name='2D')


pd.set_option('display.max_rows', None)
vars_mean_list

## short names
var_vec_2d = np.asarray(vars_mean_list.variable_id[vars_mean_list.dimensions == 'time, x, y'])

## long names
var_vec_2d_long = np.asarray(vars_mean_list.standard_name[vars_mean_list.dimensions == 'time, x, y'])

## units
var_vec_2d_units = np.asarray(vars_mean_list.units[vars_mean_list.dimensions == 'time, x, y'])

## load all simulations located in subfolders of the given directory
df_col_2d = load_sims_2d('/data/project/comble-mip/output_les/',var_vec_2d,t_shift=-2,times=time_vec,keyword=sim_keyword)

Loading variables: f(time,x,y)
/data/project/comble-mip/output_les/wrf/WRF_Lx25_dx100_FixN_2D.nc
NaN values in pr
NaN values in lwp
NaN values in iwp
NaN values in opt
NaN values in alb
NaN values in olr11
/data/project/comble-mip/output_les/scale-sdm/scale-sn14_Lx25_dx100_FixN_2d.nc
...adjusting x and y values
/data/project/comble-mip/output_les/dharma/DHARMA_Lx125_dx100_FixN_2D.nc
No field for hfss
No field for hfls
No field for ustar
/data/project/comble-mip/output_les/dharma/DHARMA_Lx25_dx100_FixN_2D.nc
No field for ustar
/data/project/comble-mip/output_les/uclales-salsa/UCLALES-SALSA_Lx25_dx100_FixN_2D.nc
...adjusting x and y values

%run functions_plotting.py
plot_2d_map(df_col_2d,var_vec_2d,times=time_vec)

../../_images/f6e0d37b1c011042452ae5a5c582e3a0eb3d65f59a02097c1e55bcd2e0f69adc.png

../../_images/27d1e55fa2955f963eb42647a4b1973fb348c6ebe99c3b336a87dbf48e1cd01e.png

../../_images/f47dcaea90e769e187ba23d171d4192f340a341e070728817de3f7aeba8cd872.png

../../_images/8f0d04e05996b77aaf9165ea2be083ef51a98f5a9531e0d710e7c5ca295f0c73.png

../../_images/babee63da42f1fda0f3945924ab8df08953f1cad2e1648afe79cbcbd3683ccd3.png

../../_images/c5f8d672a6079cf3271d88839ceca8d2e69b4c491aa3dde1c483d369d9ba5e35.png

../../_images/635c5ca759bcab3405324610447d10f598690b61e1f1ccb7831a4cc7dd0cc9d3.png

../../_images/1a421bc44a00363f8ac6bd372458ab654c36a5e88ae58b4324cc376a703ac009.png

../../_images/8e007a3435a5bee02b69da31a5cb3500a50d85cc9128c0384ca5dce6ed30ee80.png

../../_images/1e5c6347cbc7fe93c8eaeb2c80b8801aaf15d2b7f65082521febdbf4f5aa7ec3.png

../../_images/96c8e844ee05b3a2a5b5e7cc2c3b3737d263f3b40a620daa888dcaf76a68efd6.png

<Figure size 1000x600 with 0 Axes>