From 1a711f256e3cb195bc5683575084c71bf1dc2348 Mon Sep 17 00:00:00 2001 From: Tracy Date: Wed, 5 Feb 2020 16:01:59 -0700 Subject: [PATCH 1/2] Add if-statement for imp_physics when calling GSD VIS (dependency is that imp_physics=28) --- sorc/ncep_post.fd/MDLFLD.f | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/sorc/ncep_post.fd/MDLFLD.f b/sorc/ncep_post.fd/MDLFLD.f index 30e879f16..f7863566c 100644 --- a/sorc/ncep_post.fd/MDLFLD.f +++ b/sorc/ncep_post.fd/MDLFLD.f @@ -3331,6 +3331,7 @@ SUBROUTINE MDLFLD ! --- GSD VISIBILITY ! IF (IGET(410).GT.0) THEN + IF (IMP_PHYSICS == 28) THEN CALL CALVIS_GSD(CZEN,VIS) DO J=JSTA,JEND DO I=1,IM @@ -3346,7 +3347,10 @@ SUBROUTINE MDLFLD fld_info(cfld)%lvl=LVLSXML(1,IGET(410)) datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend) endif + ELSE + print*, 'GSD Visibility available for imp_physics=28, current imp_physics=',imp_physics ENDIF + ENDIF ! ! --- RADAR REFLECT - 1km ! From 3b79e89b6fb4bf07f10180a4ce25907e7fdccbaa Mon Sep 17 00:00:00 2001 From: Tracy Date: Thu, 6 Feb 2020 11:19:01 -0700 Subject: [PATCH 2/2] Fix the way that the parameter 'rg' is read in, should be read in from params.F --- sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f b/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f index 8623d144b..84d4b38b5 100644 --- a/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f +++ b/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f @@ -2122,7 +2122,7 @@ REAL FUNCTION EFFR(pmid,t,q,qqw,qqi,qqr,f_rimef, nlice, nrain, & ! JASON OTKIN AND WILLIAM LEWIS ! 09 DECEMBER 2014 - use params_mod, only: pi, rd, d608 + use params_mod, only: pi, rd, d608, rg implicit none @@ -2174,7 +2174,7 @@ REAL FUNCTION EFFR(pmid,t,q,qqw,qqi,qqr,f_rimef, nlice, nrain, & real, parameter :: min_c=2.e-6, min_r=20.e-6, min_i=4.e-6,min_s=20.e-6, min_g=20.e-6 real, parameter :: max_c=1.e-2, max_r=1.e-2, max_i=1.e-3,max_s=2.e-2, max_g=5.e-0 - real :: rg, am_g, bm_g, mu_g + real :: am_g, bm_g, mu_g real :: cgg(3), cge(3), oge1, obmg, ogg1, ogg2 double precision :: no_exp, no_min, lm_exp, lamg, lamc, lamr, lami, lams