Skip to content

Use Qv or specific humidity in interpolation if they are available #1959

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
weiwangncar merged 5 commits into from
Jan 24, 2024
Merged

Use Qv or specific humidity in interpolation if they are available #1959

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
bf96db6
use Qv or specific humidity in vertical interpolation if they are ava…
weiwangncar Dec 27, 2023
7fbd23b
add a namelist use_sh_mr to allow user to choose whether to use SH or…
weiwangncar Jan 13, 2024
1b6990a
rename namelist use_sh_mr to use_sh_qv
weiwangncar Jan 13, 2024
1d46c01
correct a logic
weiwangncar Jan 13, 2024
0e6d007
change the logical variable to use_sh_qv
weiwangncar Jan 13, 2024
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
1 change: 1 addition & 0 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2282,6 +2282,7 @@ rconfig integer interp_method_type namelist,domains 1 2
rconfig logical aggregate_lu namelist,domains 1 .false. irh "aggregate_lu" "T/F aggregate the grass, shrubs, trees in LU"
rconfig logical rh2qv_wrt_liquid namelist,domains 1 .true. irh "rh2qv_wrt_liquid" "T = rh=>Qv assumes RH wrt liquid water, F = allows ice"
rconfig integer rh2qv_method namelist,domains 1 1 irh "rh2qv_method" "1=old MM5 method, 2=new WMO method"
rconfig logical use_sh_qv namelist,domains 1 .false. irh "use_sh_qv" "T/F whether to use SH or mixing ratio in input"
rconfig real qv_max_p_safe namelist,domains 1 10000 irh "qv_max_p_safe" "Threshhold pressure, Qv > flag set to value" "Pa"
rconfig real qv_max_flag namelist,domains 1 1.E-5 irh "qv_max_flag" "Qv flag for max" "kg kg{-1}"
rconfig real qv_max_value namelist,domains 1 3.E-6 irh "qv_max_value" "Qv value for max" "kg kg{-1}"
Expand Down
32 changes: 32 additions & 0 deletions dyn_em/module_initialize_real.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1181,6 +1181,7 @@ SUBROUTINE init_domain_rk ( grid &
END IF

! Some data sets do not provide a 3d geopotential height field.
! This calculation is more accurate if the data is bottom-up.

IF ( grid%ght_gc(i_valid,grid%num_metgrid_levels/2,j_valid) .LT. 1 ) THEN
DO j = jts, MIN(jte,jde-1)
Expand Down Expand Up @@ -1239,6 +1240,15 @@ SUBROUTINE init_domain_rk ( grid &
END DO
END IF

IF ( flag_sh .EQ. 1 ) THEN
DO j = jts, min(jde-1,jte)
DO i = its, min(ide-1,ite)
IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE
grid%q2(i,j)=grid%qv_gc(i,1,j)
END DO
END DO
END IF

! The requested ptop for real data cases.

p_top_requested = grid%p_top_requested
Expand Down Expand Up @@ -1708,6 +1718,23 @@ SUBROUTINE init_domain_rk ( grid &
ims , ime , jms , jme , kms , kme , &
its , ite , jts , jte , kts , kte )

! when specific humidity is available, qv_gc is computed from sh_gc
IF (config_flags%use_sh_qv .and. (flag_sh .eq. 1 .or. flag_qv .eq. 1)) THEN
CALL vert_interp ( grid%qv_gc , grid%pd_gc , moist(:,:,:,P_QV) , grid%pb , &
grid%hgtmaxw , grid%hgttrop , grid%pmaxw , grid%ptrop , &
grid%pmaxwnn , grid%ptropnn , &
0 , 0 , &
config_flags%maxw_horiz_pres_diff , config_flags%trop_horiz_pres_diff , &
config_flags%maxw_above_this_level , &
num_metgrid_levels , 'Q' , &
interp_type , lagrange_order , extrap_type , &
lowest_lev_from_sfc , use_levels_below_ground , use_surface , &
zap_close_levels , force_sfc_in_vinterp , grid%id , &
ids , ide , jds , jde , kds , kde , &
ims , ime , jms , jme , kms , kme , &
its , ite , jts , jte , kts , kte )
END IF

! If this is theta being interpolated, AND we have extra levels for temperature,
! convert those extra levels (trop and max wind) to potential temp.

Expand Down Expand Up @@ -1778,6 +1805,8 @@ SUBROUTINE init_domain_rk ( grid &
its , ite , jts , jte , kts , kte )
END IF

! do not compute qv from RH if flag_sh or flag_qv = 1, or use_sh_qv = F
IF ( flag_sh .ne. 1 .or. flag_qv .ne. 1 .or. .not. config_flags%use_sh_qv ) THEN
IF ( config_flags%rh2qv_method .eq. 1 ) THEN
CALL rh_to_mxrat1(grid%u_1, grid%v_1, grid%p , moist(:,:,:,P_QV) , &
config_flags%rh2qv_wrt_liquid , &
Expand All @@ -1799,6 +1828,7 @@ SUBROUTINE init_domain_rk ( grid &
ims , ime , jms , jme , kms , kme , &
its , ite , jts , jte , kts , kte-1 )
END IF
END IF

IF ( .NOT. config_flags%interp_theta ) THEN
CALL t_to_theta ( grid%t_2 , grid%p , p00 , &
Expand Down Expand Up @@ -4042,6 +4072,7 @@ SUBROUTINE init_domain_rk ( grid &
ims , ime , jms , jme , kms , kme , &
its , ite , jts , jte , kts , kte )

IF ( flag_sh .ne. 1 .or. flag_qv .ne. 1 .or. .not. config_flags%use_sh_qv ) THEN
IF ( config_flags%rh2qv_method .eq. 1 ) THEN
CALL rh_to_mxrat1(grid%u_1, grid%v_1, grid%p_hyd , moist(:,:,:,P_QV) , &
config_flags%rh2qv_wrt_liquid , &
Expand All @@ -4063,6 +4094,7 @@ SUBROUTINE init_domain_rk ( grid &
ims , ime , jms , jme , kms , kme , &
its , ite , jts , jte , kts , kte-1 )
END IF
END IF

! Compute pressure similarly to how computed within model, with final Qv.

Expand Down
3 changes: 2 additions & 1 deletion run/README.namelist
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -251,7 +251,8 @@ Namelist variables specifically for the WPS input for real:
rh2qv_method = 1, ! which method to use to computer mixing ratio from RH:
default is option 1, the old MM5 method; option 2 uses a WMO
recommended method (WMO-No. 49, corrigendum, August 2000) -
there is a difference between the two methods though small
use_sh_qv = .false., ! whether to use specific humidity or mixing ratio data from input
recommended if input data has high vertical resolution
interp_theta = .false. ! If set to .false., it will vertically interpolate temperature
instead of potential temperature, which may reduce bias when
compared with input data
Expand Down