Hi Linda:
I apologize for not checking earlier that you had replied to my email. In mpas_atmphys_driver_radiation_sw.F, and at the bottom of the subroutine radiation_sw_to_MPAS, you find the following loop:
do k = kts,kte
do i = its,ite
rthratensw(k,i) = rthratensw_p(i,k,j)
enddo
enddo
Before that loop, you could add something like:
do i = its,ite
do k = kts,kte
if(pres2_hyd_p(i,k,j) .ge. 50000. and pres2_hyd_p_p(i,k+1,j) .le. 50000.) then
rthratensw_p(i,k,j) = .... (what you wou;d like)
endif
enddo
enddo
I think that this should work. Let me know.
Laura