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 !
 
Top