Skip to content

WRFDA for racm_soa_vbs_da (chem_opt=108) in the WRF-Chem model #1575

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
liujake merged 6 commits into from
Nov 19, 2021
Merged

WRFDA for racm_soa_vbs_da (chem_opt=108) in the WRF-Chem model #1575

Show file tree
Hide file tree
Changes from 5 commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
c1ce4d2
Aerosol species for chem_cv_options == 108 are added for racm_soa_vbs…
syha Jun 21, 2021
5c0cef1
Major bug fixes for mallocation issues with be%cv and be%cv_mz_chemic
syha Aug 5, 2021
68131d2
Clean up
syha Nov 2, 2021
f8a8d18
The writing format for qcstat_conv_01 is now updated in da_print_qcst…
syha Aug 12, 2021
fe88dfa
Add var/da/da_statistics/da_print_qcstat.inc and create var/run/namel…
syha Nov 3, 2021
9f44dbb
Duplicated namelist parameters for scaling12 and var12 are removed fr…
syha Nov 16, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions Registry/registry.var
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -300,6 +300,7 @@ rconfig real var_scaling8 namelist,wrfvar7 max_outer_iteratio
rconfig real var_scaling9 namelist,wrfvar7 max_outer_iterations 1.0 - "var_scaling9" "" ""
rconfig real var_scaling10 namelist,wrfvar7 max_outer_iterations 1.0 - "var_scaling10" "" ""
rconfig real var_scaling11 namelist,wrfvar7 max_outer_iterations 1.0 - "var_scaling11" "" ""
rconfig real var_scaling12 namelist,wrfvar7 45*max_outer_iterations 1.0 - "var_scaling12" "" ""
Comment thread
syha marked this conversation as resolved.
Outdated
rconfig real len_scaling1 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling1" "" ""
rconfig real len_scaling2 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling2" "" ""
rconfig real len_scaling3 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling3" "" ""
Expand All @@ -311,6 +312,7 @@ rconfig real len_scaling8 namelist,wrfvar7 max_outer_iteratio
rconfig real len_scaling9 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling9" "" ""
rconfig real len_scaling10 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling10" "" ""
rconfig real len_scaling11 namelist,wrfvar7 max_outer_iterations 1.0 - "len_scaling11" "" ""
rconfig real len_scaling12 namelist,wrfvar7 45*max_outer_iterations 1.0 - "len_scaling12" "" ""
rconfig real je_factor namelist,wrfvar7 1 1.0 - "je_factor" "" ""
rconfig real power_truncation namelist,wrfvar7 1 1.0 - "power_truncation" "" ""
rconfig logical def_sub_domain namelist,wrfvar8 1 .false. - "def_sub_domain" "" ""
Expand Down Expand Up @@ -390,6 +392,7 @@ rconfig real max_vert_var8 namelist,wrfvar13 1 99.0 - "ma
rconfig real max_vert_var9 namelist,wrfvar13 1 99.0 - "max_vert_var9" "" ""
rconfig real max_vert_var10 namelist,wrfvar13 1 99.0 - "max_vert_var10" "" ""
rconfig real max_vert_var11 namelist,wrfvar13 1 99.0 - "max_vert_var11" "" ""
rconfig real max_vert_var12 namelist,wrfvar13 45*1 99.0 - "max_vert_var12" "" ""
rconfig real max_vert_var_alpha namelist,wrfvar13 1 99.0 - "max_vert_var_alpha" "" ""
rconfig real psi_chi_factor namelist,wrfvar13 1 1.0 - "psi_chi_factor" "" ""
rconfig real psi_t_factor namelist,wrfvar13 1 1.0 - "psi_t_factor" "" ""
Expand Down
131 changes: 126 additions & 5 deletions Registry/registry.var_chem
View file
Open in desktop

Large diffs are not rendered by default.

6 changes: 1 addition & 5 deletions var/da/da_chem_sfc/da_ao_stats_chem_sfc.inc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,10 +35,6 @@ subroutine da_ao_stats_chem_sfc (stats_unit, iv, re, ob)

do ichem=PARAM_FIRST_SCALAR, num_chemic_surf
nchem=0
! stats%maximum%chem(ichem) = maxmin_type (missing_r, 0, 0)
! stats%minimum%chem(ichem) = maxmin_type (missing_r, 0, 0)
! stats%average = residual_chem_sfc_type(0.0, 0.0, 0.0, 0.0, 0.0)
! stats%rms_err = stats%average

do n=1, iv%info(chemic_surf)%nlocal
if (iv%info(chemic_surf)%proc_domain(1,n)) then
Expand Down Expand Up @@ -67,7 +63,7 @@ subroutine da_ao_stats_chem_sfc (stats_unit, iv, re, ob)
write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of AO for chem'
if (chemicda_opt == 1 .or. chemicda_opt == 2) then
call da_print_stats_chem_sfc(stats_unit, nchem, stats)
else if (chemicda_opt == 3 .or. chemicda_opt == 4 .or. chemicda_opt == 5) then
else if (chemicda_opt >= 3) then
call da_print_stats_chem_sfc(stats_unit, nchem, stats, ichem)
end if
end if
Expand Down
5 changes: 2 additions & 3 deletions var/da/da_chem_sfc/da_calculate_grady_chem_sfc.inc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,9 +18,8 @@ subroutine da_calculate_grady_chem_sfc(iv, re, jo_grad_y)
do n=1, iv%info(chemic_surf)%nlocal
if (iv%chemic_surf(n)%chem(ichem)%qc < obs_qc_pointer) re%chemic_surf(n)%chem(ichem) = 0.0

if (iv%chemic_surf(n)%chem(ichem)%qc < obs_qc_pointer) re%chemic_surf(n)%chem(ichem) = 0.0

jo_grad_y%chemic_surf(n)%chem(ichem) = -re%chemic_surf(n)%chem(ichem) / (iv%chemic_surf(n)%chem(ichem)%error * iv%chemic_surf(n)%chem(ichem)%error)
jo_grad_y%chemic_surf(n)%chem(ichem) = -re%chemic_surf(n)%chem(ichem) / &
(iv%chemic_surf(n)%chem(ichem)%error * iv%chemic_surf(n)%chem(ichem)%error)
Comment thread
syha marked this conversation as resolved.
end do
end do

Expand Down
11 changes: 10 additions & 1 deletion var/da/da_chem_sfc/da_chem_sfc.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,16 @@ module da_chem_sfc
p_chem_ic_so4_a01, p_chem_ic_so4_a02, p_chem_ic_so4_a03, p_chem_ic_so4_a04, p_chem_ic_no3_a01, p_chem_ic_no3_a02, p_chem_ic_no3_a03, p_chem_ic_no3_a04,&
p_chem_ic_nh4_a01, p_chem_ic_nh4_a02, p_chem_ic_nh4_a03, p_chem_ic_nh4_a04, p_chem_ic_cl_a01, p_chem_ic_cl_a02, p_chem_ic_cl_a03, p_chem_ic_cl_a04, &
p_chem_ic_na_a01, p_chem_ic_na_a02, p_chem_ic_na_a03, p_chem_ic_na_a04, p_chem_ic_oin_a01, p_chem_ic_oin_a02, p_chem_ic_oin_a03, p_chem_ic_oin_a04, &
p_chem_ic_so2, p_chem_ic_no2, p_chem_ic_o3, p_chem_ic_co
p_chem_ic_so2, p_chem_ic_no2, p_chem_ic_o3, p_chem_ic_co, &
p_chem_ic_so4aj, p_chem_ic_so4ai, p_chem_ic_nh4aj, p_chem_ic_nh4ai, & ! aerosols in racm_soa_vbs_da
p_chem_ic_no3aj, p_chem_ic_no3ai, p_chem_ic_naaj, p_chem_ic_naai, &
p_chem_ic_asoa1j, p_chem_ic_asoa1i, p_chem_ic_asoa2j, p_chem_ic_asoa2i, &
p_chem_ic_asoa3j, p_chem_ic_asoa3i, p_chem_ic_asoa4j, p_chem_ic_asoa4i, &
p_chem_ic_bsoa1j, p_chem_ic_bsoa1i, p_chem_ic_bsoa2j, p_chem_ic_bsoa2i, &
p_chem_ic_bsoa3j, p_chem_ic_bsoa3i, p_chem_ic_bsoa4j, p_chem_ic_bsoa4i, &
p_chem_ic_orgpaj, p_chem_ic_orgpai, p_chem_ic_ecj, p_chem_ic_eci, &
p_chem_ic_p25j, p_chem_ic_p25i, p_chem_ic_antha, p_chem_ic_seas, &
p_chem_ic_claj, p_chem_ic_clai, p_chem_ic_soila
#endif

use da_define_structures, only : iv_type, y_type, jo_type, maxmin_type, &
Expand Down
Loading