Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

solar_diagnostics output variables not showing in WRFout

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

julliana freire

New member
:cry: :cry: :cry: :cry: :cry: :cry: :cry: :cry:
I've tried running WRF-SOLAR version 4.3 on a short simulation test to get the variables GHI, DNI, and DIF. But, none of them are shown in the wrfout file.

Following another topic in the forum (https://forum.mmm.ucar.edu/phpBB3/viewtopic.php?f=104&t=10331)
I used the reference configuration of WRF-Solar(https://ral.ucar.edu/solutions/products/wrf-solar) but again without any of the variables being output.

Below is my initial namelist.input:

&time_control
run_days = 0,
run_hours = 24,
run_minutes = 0,
run_seconds = 0,
start_year = 2017,
start_month = 10,
start_day = 01,
start_hour = 00,
end_year = 2017,
end_month = 10,
end_day = 02,
end_hour = 00,
interval_seconds = 3600
input_from_file = .true.,
history_interval = 15,
frames_per_outfile = 1000,
restart = .false.,
restart_interval = 9999,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput1 = 2
io_form_auxhist2 = 2
debug_level = 0
/

&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 170,
e_sn = 100,
s_vert = 1,
e_vert = 45,
eta_levels = 1.00000, 0.99629, 0.99257, 0.98879, 0.98486, 0.98071, 0.97622, 0.97130, 0.96585, 0.95977, 0.95299, 0.94540, 0.93692, 0.92744, 0.91686, 0.90507, 0.89195, 0.87737, 0.86120, 0.84331, 0.82356, 0.80181, 0.77793, 0.75181, 0.72335, 0.69246, 0.65911, 0.62329, 0.58506, 0.54455, 0.50195, 0.45755, 0.41175, 0.36503, 0.31802, 0.27144, 0.22617, 0.18317, 0.14344, 0.10788, 0.07710, 0.05132, 0.03028, 0.01343, 0.00000,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 5000,
dy = 5000,
grid_id = 1, 2,
parent_id = 1, 1,
i_parent_start = 1, 53,
j_parent_start = 1, 25,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
feedback = 0,
/

&physics
mp_physics = 8,
ra_lw_physics = 4,
ra_sw_physics = 4,
radt = 5,
sf_sfclay_physics = 1,
sf_surface_physics = 2,
bl_pbl_physics = 5,
bldt = 0,
bl_mynn_tkeadvect = .true.,
bl_mynn_edmf = 0,
cu_physics = 0,
cu_rad_feedback = .false.,
cudt = 0,
shcu_physics = 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
icloud_bl = 0,
num_soil_layers = 4,
sf_urban_physics = 0,
num_land_cat = 21,
aer_opt = 1,
swint_opt = 2,
usemonalb = .true.,
/

&fdda
/

&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
damp_opt = 3,
zdamp = 5000.,
dampcoef = 0.2,
non_hydrostatic = .true.,
/

&bdy_control
spec_bdy_width = 5,
specified = .true.
spec_zone = 1,
relax_zone = 4,
nested = .false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&diags
solar_diagnostics = 1,
/

Follow attached the .logs (rsl.error.0000 and rsl.out.0000)

Thanks in advance for any help you can provide.
Julliana Freire.
 

Attachments

  • rsl.error.txt
    251.6 KB · Views: 42
  • rsl.out.txt
    251.5 KB · Views: 49
GHI should be in the output. It is a variable called SWDOWN.
The direct normal irradiance (DNI) and the diffuse irradiance are exposed via the Registry. You need to edit the file Registry.EM_COMMON and add the variables SWDDNI and SWDDIF to the WRF output by adding a h in the IO part of the table for these variables.
 
Thanks! It worked well.

Also, there is another option without recompiling.
1. Edit the namelist.input file, the &time_control namelist record
iofields_filename = "myoutfields.txt”
ignore_iofields_warning = .true.
2. prepare a text file ("my_output.txt") to select io fields:
+:h:0:SWDOWN,SWDDNI,SWDDIR,SWDDIF

Both options worked for my case.

Thanks!
 
Top