Skip to content

Add reflectivity calculation back to WSM6 scheme #2098

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 6 commits into from
Aug 30, 2024
Merged

Add reflectivity calculation back to WSM6 scheme #2098

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
0e34287
Add reflectivity calculation back for WSM6
weiwangncar Aug 20, 2024
2e18018
To trigger a test...
weiwangncar Aug 21, 2024
964bb2d
To trigger a test...
weiwangncar Aug 22, 2024
09dc5e5
To trigger a test, again
weiwangncar Aug 24, 2024
e8edc22
Try to trigger another test...
weiwangncar Aug 29, 2024
507a8e2
Recover the original code...
weiwangncar Aug 29, 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
19 changes: 18 additions & 1 deletion phys/module_mp_wsm6.F
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
module module_mp_wsm6
use ccpp_kind_types,only: kind_phys

use mp_wsm6,only: mp_wsm6_run
use mp_wsm6,only: mp_wsm6_run,refl10cm_wsm6
use mp_wsm6_effectrad,only: mp_wsm6_effectRad_run


Expand Down Expand Up @@ -231,6 +231,23 @@ subroutine wsm6(th,q,qc,qr,qi,qs,qg,den,pii,p,delz, &
enddo
enddo

if (diagflag .and. do_radar_ref == 1) then
do i=its,ite
do k=kts,kte
t1d(k)=th(i,k,j)*pii(i,k,j)
p1d(k)=p(i,k,j)
qv1d(k)=q(i,k,j)
qr1d(k)=qr(i,k,j)
qs1d(k)=qs(i,k,j)
qg1d(k)=qg(i,k,j)
enddo
call refl10cm_wsm6 (qv1d,qr1d,qs1d,qg1d,t1d,p1d,dBZ,kts,kte)
do k = kts, kte
refl_10cm(i,k,j) = MAX(-35., dBZ(k))
enddo
enddo
endif

enddo

end subroutine wsm6
Expand Down