diff --git a/sorc/ncep_post.fd/SURFCE.f b/sorc/ncep_post.fd/SURFCE.f index e3ca85ab8..492f418e3 100644 --- a/sorc/ncep_post.fd/SURFCE.f +++ b/sorc/ncep_post.fd/SURFCE.f @@ -46,6 +46,7 @@ !> 2023-10-04 | W Meng | Fix mismatched IDs from 526-530 !> 2023-10-05 | E James | Correcting bug fix in GSL precip type for RRFS (was using 1000x 1h pcp) !> 2024-01-23 | E James | Using consistent snow ratio SR from history files throughout GSL precip type diagnosis. +!> 2024-01-30 | A Jensen | Comment out graupel precipitation warning. !> !> @note !> USAGE: CALL SURFCE @@ -5153,11 +5154,17 @@ SUBROUTINE SURFCE totprcp = (RAINC_BUCKET(I,J) + RAINNC_BUCKET(I,J))*1.e-3 ENDIF snowratio = 0.0 - if(graup_bucket(i,j)*1.e-3 > totprcp.and.graup_bucket(i,j)/=spval)then - print *,'WARNING - Graupel is higher that total precip at point',i,j - print *,'totprcp,graup_bucket(i,j),snow_bucket(i,j),rainnc_bucket',& - totprcp,graup_bucket(i,j),snow_bucket(i,j),rainnc_bucket(i,j) - endif +!-- This following warning message prints too often and is being commented out by +!-- Anders Jensen on 30 Jan 2024. I think that this warning message prints only when +!-- graupel alone is reaching the surface. Total precipitation is interpolated +!-- and precipitation from individual hydrometeor categories is not. Thus, when +!-- total precipitation equals graupel precipitation and total precipitation is +!-- interpolated and graupel precipitation is not, the two values may not be equal. +! if(graup_bucket(i,j)*1.e-3 > totprcp.and.graup_bucket(i,j)/=spval)then +! print *,'WARNING - Graupel is higher than total precip at point',i,j +! print *,'totprcp,graup_bucket(i,j)*1.e-3,snow_bucket(i,j),rainnc_bucket',& +! totprcp,graup_bucket(i,j)*1.e-3,snow_bucket(i,j),rainnc_bucket(i,j) +! endif ! --------------------------------------------------------------- ! Minimum 1h precipitation to even consider p-type specification