Affected Releases:
CESM1.2.z
CESM1.1.z
CESM1.0.z
(Bugzilla 2054)
Po-lun Ma reported a bug in COSP simulator.
I think there might be a bug in COSP. I am looking at CESM1.2.1: cesm1_2_1/models/atm/cam/src/physics/cam/cospsimulator_intr.F90
Line 1969-1984 reads in liquid and ice mixing ratios for COSP calculation, but line 1974 seems to use convective liquid mixing ratio as ice.
do k=1,pver
do i=1,ncol
if (cld(i,k) .gt. 0._r8) then
!! note: convective mixing ratio is the sum of shallow and deep convective clouds in CAM5
mr_ccliq(i,k) = sh_cldliq(i,k) + dp_cldliq(i,k)
mr_ccice(i,k) = sh_cldliq(i,k) + dp_cldliq(i,k)
mr_lsliq(i,k)=state%q(i,k,ixcldliq) ! mr_lsliq, mixing_ratio_large_scale_cloud_liquid, state only includes stratiform (kg/kg)
mr_lsice(i,k)=state%q(i,k,ixcldice) ! mr_lsice - mixing_ratio_large_scale_cloud_ice, state only includes stratiform (kg/kg)
else
mr_ccliq(i,k) = 0._r8
mr_ccice(i,k) = 0._r8
mr_lsliq(i,k) = 0._r8
mr_lsice(i,k) = 0._r8
end if
end do
end do
CESM1.2.z
CESM1.1.z
CESM1.0.z
(Bugzilla 2054)
Po-lun Ma reported a bug in COSP simulator.
I think there might be a bug in COSP. I am looking at CESM1.2.1: cesm1_2_1/models/atm/cam/src/physics/cam/cospsimulator_intr.F90
Line 1969-1984 reads in liquid and ice mixing ratios for COSP calculation, but line 1974 seems to use convective liquid mixing ratio as ice.
do k=1,pver
do i=1,ncol
if (cld(i,k) .gt. 0._r8) then
!! note: convective mixing ratio is the sum of shallow and deep convective clouds in CAM5
mr_ccliq(i,k) = sh_cldliq(i,k) + dp_cldliq(i,k)
mr_ccice(i,k) = sh_cldliq(i,k) + dp_cldliq(i,k)
mr_lsliq(i,k)=state%q(i,k,ixcldliq) ! mr_lsliq, mixing_ratio_large_scale_cloud_liquid, state only includes stratiform (kg/kg)
mr_lsice(i,k)=state%q(i,k,ixcldice) ! mr_lsice - mixing_ratio_large_scale_cloud_ice, state only includes stratiform (kg/kg)
else
mr_ccliq(i,k) = 0._r8
mr_ccice(i,k) = 0._r8
mr_lsliq(i,k) = 0._r8
mr_lsice(i,k) = 0._r8
end if
end do
end do