@@ -73,19 +73,19 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, lextop,&
7373 use surface_perturbation, only: cdfnor,ppfbet
7474
7575 ! For Thompson MP
76- use module_mp_thompson, only: calc_effectRad_thmpsn = > calc_effectRad, &
77- Nt_c_l_thmpsn = > Nt_c_l, &
78- Nt_c_o_thmpsn = > Nt_c_o, &
79- re_qc_min_thmpsn = > re_qc_min, &
80- re_qc_max_thmpsn = > re_qc_max, &
81- re_qi_min_thmpsn = > re_qi_min, &
82- re_qi_max_thmpsn = > re_qi_max, &
83- re_qs_min_thmpsn = > re_qs_min, &
84- re_qs_max_thmpsn = > re_qs_max
76+ use module_mp_thompson, only: calc_effectRad_thompson = > calc_effectRad, &
77+ Nt_c_l_thompson = > Nt_c_l, &
78+ Nt_c_o_thompson = > Nt_c_o, &
79+ re_qc_min_thompson = > re_qc_min, &
80+ re_qc_max_thompson = > re_qc_max, &
81+ re_qi_min_thompson = > re_qi_min, &
82+ re_qi_max_thompson = > re_qi_max, &
83+ re_qs_min_thompson = > re_qs_min, &
84+ re_qs_max_thompson = > re_qs_max
8585 use module_mp_thompson_make_number_concentrations, only: &
86- make_IceNumber_thmpsn = > make_IceNumber, &
87- make_DropletNumber_thmpsn = > make_DropletNumber, &
88- make_RainNumber_thmpsn = > make_RainNumber
86+ make_IceNumber_thompson = > make_IceNumber, &
87+ make_DropletNumber_thompson = > make_DropletNumber, &
88+ make_RainNumber_thompson = > make_RainNumber
8989
9090 use module_mp_tempo_params, only: &
9191 ty_tempo_cfg, &
@@ -792,13 +792,13 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, lextop,&
792792 qs_mp (i,k) = tracer1(i,k,ntsw)/ (1 .- qvs)
793793 if (nint (slmsk(i)) == 1 ) then
794794 if (imp_physics == imp_physics_thompson) then
795- nc_mp (i,k) = Nt_c_l_thmpsn * orho(i,k)
795+ nc_mp (i,k) = Nt_c_l_thompson * orho(i,k)
796796 else
797797 nc_mp (i,k) = Nt_c_l_tempo* orho(i,k)
798798 endif
799799 else
800800 if (imp_physics == imp_physics_thompson) then
801- nc_mp (i,k) = Nt_c_o_thmpsn * orho(i,k)
801+ nc_mp (i,k) = Nt_c_o_thompson * orho(i,k)
802802 else
803803 nc_mp (i,k) = Nt_c_o_tempo* orho(i,k)
804804 endif
@@ -923,14 +923,14 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, lextop,&
923923 do i= 1 ,im
924924 if ((ltaerosol .or. mraerosol) .and. qc_mp(i,k)>1.e-12 .and. nc_mp(i,k)<100 .) then
925925 if (imp_physics == imp_physics_thompson) then
926- nc_mp(i,k) = make_DropletNumber_thmpsn (qc_mp(i,k)* rho(i,k), nwfa(i,k)* rho(i,k)) * orho(i,k)
926+ nc_mp(i,k) = make_DropletNumber_thompson (qc_mp(i,k)* rho(i,k), nwfa(i,k)* rho(i,k)) * orho(i,k)
927927 else
928928 nc_mp(i,k) = make_DropletNumber_tempo(qc_mp(i,k)* rho(i,k), nwfa(i,k)* rho(i,k)) * orho(i,k)
929929 endif
930930 endif
931931 if (qi_mp(i,k)>1.e-12 .and. ni_mp(i,k)<100 .) then
932932 if (imp_physics == imp_physics_thompson) then
933- ni_mp(i,k) = make_IceNumber_thmpsn (qi_mp(i,k)* rho(i,k), tlyr(i,k)) * orho(i,k)
933+ ni_mp(i,k) = make_IceNumber_thompson (qi_mp(i,k)* rho(i,k), tlyr(i,k)) * orho(i,k)
934934 else
935935 ni_mp(i,k) = make_IceNumber_tempo(qi_mp(i,k)* rho(i,k), tlyr(i,k)) * orho(i,k)
936936 endif
@@ -946,18 +946,18 @@ subroutine GFS_rrtmg_pre_run (im, levs, lm, lmk, lmp, n_var_lndp, lextop,&
946946 ! it will raise the low limit from 5 to 10, but the high limit will remain 125.
947947
948948 if (imp_physics == imp_physics_thompson) then
949- call calc_effectRad_thmpsn (tlyr(i,:), plyr(i,:)* 100 ., qv_mp(i,:), qc_mp(i,:), &
949+ call calc_effectRad_thompson (tlyr(i,:), plyr(i,:)* 100 ., qv_mp(i,:), qc_mp(i,:), &
950950 nc_mp(i,:), qi_mp(i,:), ni_mp(i,:), qs_mp(i,:), &
951951 effrl(i,:), effri(i,:), effrs(i,:), islmsk, 1 , lm )
952952 ! Scale Thompson's effective radii from meter to micron
953953 do k= 1 ,lm
954- effrl(i,k) = MAX (re_qc_min_thmpsn , MIN (effrl(i,k), re_qc_max_thmpsn ))* 1.e6
955- effri(i,k) = MAX (re_qi_min_thmpsn , MIN (effri(i,k), re_qi_max_thmpsn ))* 1.e6
956- effrs(i,k) = MAX (re_qs_min_thmpsn , MIN (effrs(i,k), re_qs_max_thmpsn ))* 1.e6
954+ effrl(i,k) = MAX (re_qc_min_thompson , MIN (effrl(i,k), re_qc_max_thompson ))* 1.e6
955+ effri(i,k) = MAX (re_qi_min_thompson , MIN (effri(i,k), re_qi_max_thompson ))* 1.e6
956+ effrs(i,k) = MAX (re_qs_min_thompson , MIN (effrs(i,k), re_qs_max_thompson ))* 1.e6
957957 end do
958- effrl(i,lmk) = re_qc_min_thmpsn * 1.e6
959- effri(i,lmk) = re_qi_min_thmpsn * 1.e6
960- effrs(i,lmk) = re_qs_min_thmpsn * 1.e6
958+ effrl(i,lmk) = re_qc_min_thompson * 1.e6
959+ effri(i,lmk) = re_qi_min_thompson * 1.e6
960+ effrs(i,lmk) = re_qs_min_thompson * 1.e6
961961 else
962962 call calc_effectRad_tempo(t1d= tlyr(i,:), p1d= plyr(i,:)* 100 ., qv1d= qv_mp(i,:), qc1d= qc_mp(i,:), &
963963 nc1d= nc_mp(i,:), qi1d= qi_mp(i,:), ni1d= ni_mp(i,:), qs1d= qs_mp(i,:), &
0 commit comments