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

Output WRf Solar

vanessalmeida24

New member
Good morning, I'm running WRF Solar 4.3 and I'm not getting the output for variable FARMS" "W m-2" SWDDIR2 SWDDNI2 SWDDNI2
SWDDIF2 only outputs SWDOWN2 how should I proceed? Thanks
 
&time_control
start_year = 2004, 2004, 2007,
start_month = 01, 01, 05,
start_day = 01, 01, 31,
start_hour = 00, 00, 00,
end_year = 2004, 2004, 2007,
end_month = 01, 01, 06,
end_day = 31, 31, 09,
end_hour = 23, 23, 12,
interval_seconds = 3600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000,1000,1000,
restart = .false.,
restart_interval = 9999,
ignore_iofields_warning = .true.,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput1 = 2
io_form_auxhist2 = 2
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360,360,360,
io_form_auxinput4 = 2,
debug_level = 0,
/

&domains
time_step = 30,
use_adaptive_time_step = .true.,
step_to_output_time = .true.,
target_cfl = 0.60, 0.60, 0.60,
target_hcfl = 0.84, 0.84, 0.84,
max_step_increase_pct = 5, 51, 51,
starting_time_step = 90, 40, 13,
max_time_step = 150, 45, 15,
min_time_step = 30, 10, 3,
adaptation_domain = 1,
max_dom = 2,
e_we = 101, 106, 163,
e_sn = 101, 106, 127,
s_vert = 1,
e_vert = 45, 45, 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 = 9000, 3000,
dy = 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 35, 61,
j_parent_start = 1, 25, 49,
parent_grid_ratio = 1, 5, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 0,
/

&p &physics
mp_physics = 6, 6, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 3, 3, 5,
sf_sfclay_physics = 5, 5, 2,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 5, 5, 2,
bl_mynn_edmf = 0, 0,0,
shcu_physics = 5,5,5,
bldt = 0, 0, 0,
bl_mynn_tkeadvect = .true.,
cu_physics = 1, 1, 1,
cu_rad_feedback = .false.,
cudt = 0, 0, 0,
kfeta_trigger = 2,
isfflx = 1,
iz0tlnd = 1,
icloud = 1,
icloud_bl = 0,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0,
num_land_cat = 21,
aer_opt = 1,
swint_opt = 2,
sst_update = 1,
usemonalb = .true.,
ensdim = 144,
/

&fdda
grid_fdda = 2, 2, 2,
gfdda_inname = "wrffdda_d<domain>",
gfdda_end_h = 300, 300, 300,
gfdda_interval_m = 60, 60, 60,
fgdt = 0, 0, 0,
if_no_pbl_nudging_uv = 0, 0, 0,
if_no_pbl_nudging_t = 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1,
if_zfac_uv = 1, 1, 1, 0,
k_zfac_uv = 20, 20, 20, 0,
if_zfac_t = 1, 1, 1, 0,
k_zfac_t = 20, 20, 20, 0,
if_zfac_q = 1, 1, 1, 0,
k_zfac_q = 20, 20, 20, 0,
guv = 0.0003, 0.000075, 0.000075,
gt = 0.0003, 0.000075, 0.000075,
gq = 0.0003, 0.000075, 0.000075,
xwavenum = 14,
ywavenum = 10,
if_ramping = 0,
dtramp_min = 60.0,
io_form_gfdda = 2,
/

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

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

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&diags
solar_diagnostics = 1,
/
 
radt = 3, 3, 5,
From the namelist info file
radt (max_dom) = 30, ! minutes between radiation physics calls recommend 1 min per km of dx (e.g. 10 for 10 km); use the same value for all nests.

I will say that WRF-Solar I am not super familiar with. The best person who might be able to help is the admins or the project scientist listed on this page

However, I also noticed that you are mixing physics schemes on your different domains. In my experience when running wrf-chem or wrf by itself that sometimes leads to issues when the schemes aren't matching. Not sure if this is applicable to Solar.

I would try the radt change first and see what happens.

Also in the WRF folder there is a file called README.namelist located in the /run folder. That can also help with the physics package descriptions.

Sorry I can't be of much more help.
 
Please take a look at Registry.EM_COMMON, and look at the variables you want to output, for example swddif2, you can find the lines like:
state real swddif2 ij misc 1 - rd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

Please change the oin eta:
state real swddif2 ij misc 1 - rhd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

And do the same change for other variables.

Then please recompile the code by ./clean -a , ./configure and ./compile eem_real

I suppose you should have those variables in your output. Please let me know if you still have any issue.
 
Please take a look at Registry.EM_COMMON, and look at the variables you want to output, for example swddif2, you can find the lines like:
state real swddif2 ij misc 1 - rd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

Please change the oin eta:
state real swddif2 ij misc 1 - rhd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

And do the same change for other variables.

Then please recompile the code by ./clean -a , ./configure and ./compile eem_real

I suppose you should have those variables in your output. Please let me know if you still have any issue.
This worked for me! Thanks !
 
Dear haiquingsong, did you get SWDDIF2 in the output with the given namelist or did you modify it, can you suggest me some reference for all these options. With gratitude !!
Follow Chen Ming:
Please take a look at Registry.EM_COMMON, and look at the variables you want to output, for example swddif2, you can find the lines like:
state real swddif2 ij misc 1 - rd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

Please change the oin eta:
state real swddif2 ij misc 1 - rhd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

And do the same change for other variables.

Then please recompile the code by ./clean -a , ./configure and ./compile eem_real



Then,modify your namelist.input
add this:
&diags
solar_diagnostics = 1,
 
Follow Chen Ming:
Please take a look at Registry.EM_COMMON, and look at the variables you want to output, for example swddif2, you can find the lines like:
state real swddif2 ij misc 1 - rd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

Please change the oin eta:
state real swddif2 ij misc 1 - rhd "SWDDIF2" "Shortwave surface downward diffuse irradiance from FARMS" "W m-2" ""

And do the same change for other variables.

Then please recompile the code by ./clean -a , ./configure and ./compile eem_real



Then,modify your namelist.input
add this:
&diags
solar_diagnostics = 1,
Thank you for your response, I could get swddif2 (timeseries) after adding tslist, but not in wrfout files. I followed the same procedure as stated, but I will try it once more and update. Thank you!
 
Top