Skip to content

Commit

Permalink
Change names of phy_f3d variables in physics/GFS_rrtmg_pre.*
Browse files Browse the repository at this point in the history
  • Loading branch information
@climbfuji
climbfuji committed Oct 7, 2020
1 parent cf34191 commit a33954b
Show file tree
Hide file tree
Showing 2 changed files with 38 additions and 40 deletions.
64 changes: 31 additions & 33 deletions physics/GFS_rrtmg_pre.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -26,9 +26,9 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
ltaerosol, lgfdlmprad, uni_cld, effr_in, do_mynnedmf, lmfshal, &
lmfdeep2, fhswr, fhlwr, solhr, sup, eps, epsm1, fvirt, &
rog, rocp, con_rd, xlat_d, xlat, xlon, coslat, sinlat, tsfc, slmsk, &
prsi, prsl, prslk, tgrs, sfc_wts, phy_f3d_mg_cld, phy_f3d_reffr, &
phy_f3d_cnvw, phy_f3d_cnvc, qgrs, aer_nm, & !inputs from here and above
coszen, coszdg, phy_f3d_leffr, phy_f3d_ieffr, phy_f3d_seffr, &
prsi, prsl, prslk, tgrs, sfc_wts, mg_cld, effrr_in, &
cnvw_in, cnvc_in, qgrs, aer_nm, & !inputs from here and above
coszen, coszdg, effrl_inout, effri_inout, effrs_inout, &
clouds1, clouds2, clouds3, clouds4, clouds5, & !in/out from here and above
kd, kt, kb, mtopa, mbota, raddt, tsfg, tsfa, de_lgth, alb1d, delp, dz, & !output from here and below
plvl, plyr, tlvl, tlyr, qlyr, olyr, gasvmr_co2, gasvmr_n2o, gasvmr_ch4,&
Expand Down Expand Up @@ -101,18 +101,16 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &

real(kind=kind_phys), dimension(:,:), intent(in) :: prsi, prsl, prslk, &
tgrs, sfc_wts, &
phy_f3d_mg_cld, &
phy_f3d_reffr, &
phy_f3d_cnvw, &
phy_f3d_cnvc
mg_cld, effrr_in, &
cnvw_in, cnvc_in

real(kind=kind_phys), dimension(:,:,:), intent(in) :: qgrs, aer_nm

real(kind=kind_phys), dimension(:), intent(inout) :: coszen, coszdg

real(kind=kind_phys), dimension(:,:), intent(inout) :: phy_f3d_leffr, &
phy_f3d_ieffr, &
phy_f3d_seffr
real(kind=kind_phys), dimension(:,:), intent(inout) :: effrl_inout, &
effri_inout, &
effrs_inout
real(kind=kind_phys), dimension(im,lm+LTP), intent(inout) :: clouds1, &
clouds2, clouds3, &
clouds4, clouds5
Expand Down Expand Up @@ -683,18 +681,18 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
do k=1,lm
k1 = k + kd
do i=1,im
cldcov(i,k1) = phy_f3d_mg_cld(i,k)
effrl(i,k1) = phy_f3d_leffr(i,k)
effri(i,k1) = phy_f3d_ieffr(i,k)
effrr(i,k1) = phy_f3d_reffr(i,k)
effrs(i,k1) = phy_f3d_seffr(i,k)
cldcov(i,k1) = mg_cld(i,k)
effrl(i,k1) = effrl_inout(i,k)
effri(i,k1) = effri_inout(i,k)
effrr(i,k1) = effrr_in(i,k)
effrs(i,k1) = effrs_inout(i,k)
enddo
enddo
else
do k=1,lm
k1 = k + kd
do i=1,im
cldcov(i,k1) = phy_f3d_mg_cld(i,k)
cldcov(i,k1) = mg_cld(i,k)
enddo
enddo
endif
Expand All @@ -720,10 +718,10 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
do k=1,lm
k1 = k + kd
do i=1,im
effrl(i,k1) = phy_f3d_leffr(i,k)
effri(i,k1) = phy_f3d_ieffr(i,k)
effrr(i,k1) = phy_f3d_reffr(i,k)
effrs(i,k1) = phy_f3d_seffr(i,k)
effrl(i,k1) = effrl_inout(i,k)
effri(i,k1) = effri_inout(i,k)
effrr(i,k1) = effrr_in(i,k)
effrs(i,k1) = effrs_inout(i,k)
! if(me==0) then
! if(effrl(i,k1)> 5.0) then
! write(6,*) 'rad driver:cloud radii:',kdt, i,k1, &
Expand Down Expand Up @@ -783,9 +781,9 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
do k=1,lm
k1 = k + kd
do i=1,im
phy_f3d_leffr(i,k) = effrl(i,k1)
phy_f3d_ieffr(i,k) = effri(i,k1)
phy_f3d_seffr(i,k) = effrs(i,k1)
effrl_inout(i,k) = effrl(i,k1)
effri_inout(i,k) = effri(i,k1)
effrs_inout(i,k) = effrs(i,k1)
enddo
enddo
else ! all other cases
Expand All @@ -806,16 +804,16 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
!GJF: this is not consistent with GFS_typedefs,
! but it looks like the Zhao-Carr-PDF scheme is not in the CCPP
deltaq(i,k1) = 0.0!Tbd%phy_f3d(i,k,5) !GJF: this variable is not in phy_f3d anymore
cnvw (i,k1) = phy_f3d_cnvw(i,k)
cnvc (i,k1) = phy_f3d_cnvc(i,k)
cnvw (i,k1) = cnvw_in(i,k)
cnvc (i,k1) = cnvc_in(i,k)
enddo
enddo
elseif ((npdf3d == 0) .and. (ncnvcld3d == 1)) then ! same as imp_physics=99
do k=1,lm
k1 = k + kd
do i=1,im
deltaq(i,k1) = 0.0
cnvw (i,k1) = phy_f3d_cnvw(i,k)
cnvw (i,k1) = cnvw_in(i,k)
cnvc (i,k1) = 0.0
enddo
enddo
Expand Down Expand Up @@ -904,18 +902,18 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &

elseif(imp_physics == imp_physics_wsm6 .or. imp_physics == imp_physics_fer_hires) then
if (kdt == 1) then
phy_f3d_leffr(:,:) = 10.
phy_f3d_ieffr(:,:) = 50.
phy_f3d_seffr(:,:) = 250.
effrl_inout(:,:) = 10.
effri_inout(:,:) = 50.
effrs_inout(:,:) = 250.
endif

call progcld5 (plyr,plvl,tlyr,qlyr,qstl,rhly,tracer1, & ! --- inputs
xlat,xlon,slmsk,dz,delp, &
ntrac-1, ntcw-1,ntiw-1,ntrw-1, &
ntsw-1,ntgl-1, &
im, lmk, lmp, uni_cld, lmfshal, lmfdeep2, &
cldcov(:,1:LMK),phy_f3d_leffr(:,:), &
phy_f3d_ieffr(:,:), phy_f3d_seffr(:,:), &
cldcov(:,1:LMK),effrl_inout(:,:), &
effri_inout(:,:), effrs_inout(:,:), &
dzb, xlat_d, julian, yearlen, &
clouds,cldsa,mtopa,mbota, de_lgth, alpha) ! --- outputs

Expand Down Expand Up @@ -946,8 +944,8 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, &
ntrac-1, ntcw-1,ntiw-1,ntrw-1, &
ntsw-1,ntgl-1, &
im, lmk, lmp, uni_cld, lmfshal, lmfdeep2, &
cldcov(:,1:LMK), phy_f3d_leffr(:,:), &
phy_f3d_ieffr(:,:), phy_f3d_seffr(:,:), &
cldcov(:,1:LMK), effrl_inout(:,:), &
effri_inout(:,:), effrs_inout(:,:), &
dzb, xlat_d, julian, yearlen, &
clouds, cldsa, mtopa ,mbota, de_lgth, alpha) ! --- outputs
endif ! MYNN PBL or GF
Expand Down
14 changes: 7 additions & 7 deletions physics/GFS_rrtmg_pre.meta
View file
Original file line number Diff line number Diff line change
Expand Up @@ -617,7 +617,7 @@
kind = kind_phys
intent = in
optional = F
[phy_f3d_mg_cld]
[mg_cld]
standard_name = cloud_fraction_for_MG
long_name = cloud fraction used by Morrison-Gettelman MP
units = frac
Expand All @@ -626,7 +626,7 @@
kind = kind_phys
intent = in
optional = F
[phy_f3d_reffr]
[effrr_in]
standard_name = effective_radius_of_stratiform_cloud_rain_particle_in_um
long_name = effective radius of cloud rain particle in micrometers
units = um
Expand All @@ -635,7 +635,7 @@
kind = kind_phys
intent = in
optional = F
[phy_f3d_cnvw]
[cnvw_in]
standard_name = convective_cloud_water_mixing_ratio_in_phy_f3d
long_name = convective cloud water mixing ratio in the phy_f3d array
units = kg kg-1
Expand All @@ -644,7 +644,7 @@
kind = kind_phys
intent = in
optional = F
[phy_f3d_cnvc]
[cnvc_in]
standard_name = convective_cloud_cover_in_phy_f3d
long_name = convective cloud cover in the phy_f3d array
units = frac
Expand Down Expand Up @@ -689,7 +689,7 @@
kind = kind_phys
intent = inout
optional = F
[phy_f3d_leffr]
[effrl_inout]
standard_name = effective_radius_of_stratiform_cloud_liquid_water_particle_in_um
long_name = eff. radius of cloud liquid water particle in micrometer
units = um
Expand All @@ -698,7 +698,7 @@
kind = kind_phys
intent = inout
optional = F
[phy_f3d_ieffr]
[effri_inout]
standard_name = effective_radius_of_stratiform_cloud_ice_particle_in_um
long_name = eff. radius of cloud ice water particle in micrometer
units = um
Expand All @@ -707,7 +707,7 @@
kind = kind_phys
intent = inout
optional = F
[phy_f3d_seffr]
[effrs_inout]
standard_name = effective_radius_of_stratiform_cloud_snow_particle_in_um
long_name = effective radius of cloud snow particle in micrometers
units = um
Expand Down

0 comments on commit a33954b

Please sign in to comment.