Skip to content

Updates to module_SGSCloud_RadPre.F90 and .meta (from Tanya) #42

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
DomHeinzeller merged 2 commits into from
Jun 29, 2020
Merged

Updates to module_SGSCloud_RadPre.F90 and .meta (from Tanya) #42

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
46d0200
Merge pull request #3 from NOAA-GSD/gsd/develop
joeolson42 Jun 26, 2020
3e214f7
Updates to module_SGSCloud_RadPre.F90 and .meta (from Tanya)
joeolson42 Jun 26, 2020
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
72 changes: 34 additions & 38 deletions physics/module_SGSCloud_RadPre.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@ subroutine sgscloud_radpre_run( &
flag_init,flag_restart, &
do_mynnedmf, &
qc, qi, qv, T3D, P3D, &
qr, qs, &
qr, qs, qg, &
qci_conv, &
imfdeepcnv, imfdeepcnv_gf, &
qc_save, qi_save, &
Expand Down Expand Up @@ -68,7 +68,7 @@ subroutine sgscloud_radpre_run( &
& nlay, imp_physics, imp_physics_gfdl
logical, intent(in) :: flag_init, flag_restart, do_mynnedmf
real(kind=kind_phys), dimension(im,levs), intent(inout) :: qc, qi
real(kind=kind_phys), dimension(im,levs), intent(inout) :: qr, qs
real(kind=kind_phys), dimension(im,levs), intent(inout) :: qr, qs, qg
! qci_conv only allocated if GF is used
real(kind=kind_phys), dimension(:,:), intent(inout) :: qci_conv
real(kind=kind_phys), dimension(im,levs), intent(in) :: T3D,delp, &
Expand Down Expand Up @@ -117,22 +117,20 @@ subroutine sgscloud_radpre_run( &
if ( qi(i,k) > 1E-7 .OR. qc(i,k) > 1E-7 ) then
es = min( p3d(i,k), fpvs( t3d(i,k) ) ) ! fpvs and prsl in pa
qsat = max( QMIN, eps * es / (p3d(i,k) + epsm1*es) )
rhgrid = max( 0., min( 0.95, qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) ! g/kg
rhgrid = max( 0., min( 1., qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) + qr(i,k) + qs(i,k) + qg(i,k) ! g/kg
clwt = 1.0e-6 * (p3d(i,k)*0.00001)

if (h2oliq > clwt) then
onemrh= max( 1.e-10, 1.0-rhgrid )
tem1 = min(max((onemrh*qsat)**0.49,0.0001),1.0) !jhan
tem1 = 100.0 / tem1
value = max( min( tem1*(h2oliq), 50.0 ), 0.0 )
value = max( min( tem1*(h2oliq-clwt), 50.0 ), 0.0 )
tem2 = sqrt( sqrt(rhgrid) )

clouds1(i,k) = max( tem2*(1.0-exp(-value)), 0.0 )
endif
!clouds1(i,k)=(1.-exp(-coef_alph*h2oliq/ &
! & ((1.-rhgrid)*qsat*1000.0)**coef_gamm))*(rhgrid**coef_p)
!clouds1(i,k)=max(0.0,MIN(1.,clouds1(i,k)))

endif
enddo
enddo
Expand Down Expand Up @@ -213,27 +211,20 @@ subroutine sgscloud_radpre_run( &

es = min( p3d(i,k), fpvs( t3d(i,k) ) ) ! fpvs and prsl in pa
qsat = max( QMIN, eps * es / (p3d(i,k) + epsm1*es) )
rhgrid = max( 0., min( 0.95, qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) ! g/kg
rhgrid = max( 0., min( 1., qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) + qr(i,k) + qs(i,k) + qg(i,k) ! g/kg
clwt = 1.0e-6 * (p3d(i,k)*0.00001)

if (h2oliq > clwt) then
onemrh= max( 1.e-10, 1.0-rhgrid )
tem1 = min(max((onemrh*qsat)**0.49,0.0001),1.0) !jhan
tem1 = 100.0 / tem1
value = max( min( tem1*(h2oliq), 50.0 ), 0.0 )
value = max( min( tem1*(h2oliq-clwt), 50.0 ), 0.0 )
tem2 = sqrt( sqrt(rhgrid) )

clouds1(i,k) = max( tem2*(1.0-exp(-value)), 0.0 )
endif

!es = min( p3d(i,k), fpvs( t3d(i,k) ) ) ! fpvs and prsl in pa
!qsat = max( QMIN, eps * es / (p3d(i,k) + epsm1*es) )
!rhgrid = max( 0., min( 0.95, qv(i,k)/qsat ) )
!h2oliq=1000.0*( qc(i,k) + qi(i,k) ) ! g/kg
!clouds1(i,k)=(1.-exp(-coef_alph*h2oliq/ &
! & ((1.-rhgrid)*qsat*1000.0)**coef_gamm))*(rhgrid**coef_p)
!clouds1(i,k)=max(0.0,MIN(1.,clouds1(i,k)))
endif
enddo
enddo
Expand Down Expand Up @@ -265,27 +256,32 @@ subroutine sgscloud_radpre_run( &
if(qi(i,k)>1.E-8)clouds5(i,k)=max(173.45 + 2.14*Tc, 20.)
endif

! Xu-Randall (1996) cloud fraction
es = min( p3d(i,k), fpvs( t3d(i,k) ) ) ! fpvs and prsl in pa
qsat = max( QMIN, eps * es / (p3d(i,k) + epsm1*es) )
rhgrid = max( 0., min( 0.95, qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) ! g/kg
clwt = 1.0e-6 * (p3d(i,k)*0.00001)

if (h2oliq > clwt) then
onemrh= max( 1.e-10, 1.0-rhgrid )
tem1 = min(max((onemrh*qsat)**0.49,0.0001),1.0) !jhan
tem1 = 100.0 / tem1
value = max( min( tem1*(h2oliq), 50.0 ), 0.0 )
tem2 = sqrt( sqrt(rhgrid) )

clouds1(i,k) = max( tem2*(1.0-exp(-value)), 0.0 )
if ( do_mynnedmf .or. (imp_physics == imp_physics_gfdl) ) then
!print *,'MYNN PBL or GFDL MP cldcov used'
else
clouds1(i,k) = 0.0
endif
!print*,"XuRandla- cf:",clouds1(i,k)," rh:",rhgrid," qt:",h2oliq
!print*,"XuRandlb- clwt:",clwt," qsat:",qsat," p:",p3d(i,k)
endif
!print *,'GF with Xu-Randall cloud fraction'
! Xu-Randall (1996) cloud fraction
es = min( p3d(i,k), fpvs( t3d(i,k) ) ) ! fpvs and prsl in pa
qsat = max( QMIN, eps * es / (p3d(i,k) + epsm1*es) )
rhgrid = max( 0., min( 1.00, qv(i,k)/qsat ) )
h2oliq = qc(i,k) + qi(i,k) + qr(i,k) + qs(i,k) + qg(i,k) ! g/kg
clwt = 1.0e-6 * (p3d(i,k)*0.00001)

if (h2oliq > clwt) then
onemrh= max( 1.e-10, 1.0-rhgrid )
tem1 = min(max((onemrh*qsat)**0.49,0.0001),1.0) !jhan
tem1 = 100.0 / tem1
value = max( min( tem1*(h2oliq-clwt), 50.0 ), 0.0 )
tem2 = sqrt( sqrt(rhgrid) )

clouds1(i,k) = max( tem2*(1.0-exp(-value)), 0.0 )
else
clouds1(i,k) = 0.0
endif
!print*,"XuRandla- cf:",clouds1(i,k)," rh:",rhgrid," qt:",h2oliq
!print*,"XuRandlb- clwt:",clwt," qsat:",qsat," p:",p3d(i,k)
endif ! not MYNN PBL or GFDL MP
endif ! qci_conv
enddo
enddo
endif ! imfdeepcnv_gf
Expand Down
9 changes: 9 additions & 0 deletions physics/module_SGSCloud_RadPre.meta
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,6 +106,15 @@
kind = kind_phys
intent = inout
optional = F
[qg]
standard_name = graupel_mixing_ratio
long_name = graupel mixing ratio wrt dry+vapor (no condensates)
units = kg kg-1
dimensions = (horizontal_dimension,vertical_dimension)
type = real
kind = kind_phys
intent = inout
optional = F
[qci_conv]
standard_name = convective_cloud_condesate_after_rainout
long_name = convective cloud condesate after rainout
Expand Down