1)Add 0-1km helicity in GFS produt list; 2) bug fix in cloud fraction (…

…#578) and VIL.
NOAA-EMC · Oct 25, 2022 · 331bc78 · 331bc78
1 parent ac6bfbd
commit 331bc78
Show file tree

Hide file tree

Showing 11 changed files with 51 additions and 43 deletions.
diff --git a/parm/global_1x1_paramlist_g2 b/parm/global_1x1_paramlist_g2
@@ -167,6 +167,7 @@ HGT:PV=2e-06 (Km^2/kg/s) surface
 HGT:surface
 HGT:tropopause
 HINDEX:surface
+HLCY:1000-0 m above ground
 HLCY:3000-0 m above ground
 HPBL:surface
 ICAHT:max wind

diff --git a/parm/global_1x1_paramlist_g2.f000 b/parm/global_1x1_paramlist_g2.f000
@@ -159,6 +159,7 @@ HGT:PV=2e-06 (Km^2/kg/s) surface
 HGT:surface
 HGT:tropopause
 HINDEX:surface
+HLCY:1000-0 m above ground
 HLCY:3000-0 m above ground
 HPBL:surface
 ICAHT:max wind

diff --git a/parm/postcntrl_gfs.xml b/parm/postcntrl_gfs.xml
@@ -314,8 +314,8 @@
 
     <param>
        <shortname>HLCY_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
-       <level>3000.</level>
-       <level2>0.</level2>
+       <level>3000. 1000.</level>
+       <level2>0. 0.</level2>
        <scale>4.0</scale>
     </param>
 
@@ -445,7 +445,7 @@
     </param>
 
     <param>
-       <shortname>VIS_ON_SURFACE</shortname>
+       <shortname>GSD_VIS_ON_SURFACE</shortname>
        <scale>6.0</scale>
     </param>
 

diff --git a/parm/postcntrl_gfs_f00.xml b/parm/postcntrl_gfs_f00.xml
@@ -314,8 +314,8 @@
 
     <param>
        <shortname>HLCY_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
-       <level>3000.</level>
-       <level2>0.</level2>
+       <level>3000. 1000.</level>
+       <level2>0. 0.</level2>
        <scale>4.0</scale>
     </param>
 
@@ -360,7 +360,7 @@
     </param>
 
     <param>
-       <shortname>VIS_ON_SURFACE</shortname>
+       <shortname>GSD_VIS_ON_SURFACE</shortname>
        <scale>6.0</scale>
     </param>
 

diff --git a/parm/postcntrl_gfs_f00_two.xml b/parm/postcntrl_gfs_f00_two.xml
@@ -314,8 +314,8 @@
 
     <param>
        <shortname>HLCY_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
-       <level>3000.</level>
-       <level2>0.</level2>
+       <level>3000. 1000.</level>
+       <level2>0. 0.</level2>
        <scale>4.0</scale>
     </param>
 
@@ -360,7 +360,7 @@
     </param>
 
     <param>
-       <shortname>VIS_ON_SURFACE</shortname>
+       <shortname>GSD_VIS_ON_SURFACE</shortname>
        <scale>6.0</scale>
     </param>
 

diff --git a/parm/postcntrl_gfs_two.xml b/parm/postcntrl_gfs_two.xml
@@ -314,8 +314,8 @@
 
     <param>
        <shortname>HLCY_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
-       <level>3000.</level>
-       <level2>0.</level2>
+       <level>3000. 1000.</level>
+       <level2>0. 0.</level2>
        <scale>4.0</scale>
     </param>
 
@@ -445,7 +445,7 @@
     </param>
 
     <param>
-       <shortname>VIS_ON_SURFACE</shortname>
+       <shortname>GSD_VIS_ON_SURFACE</shortname>
        <scale>6.0</scale>
     </param>
 

diff --git a/parm/postxconfig-NT-GFS-F00-TWO.txt b/parm/postxconfig-NT-GFS-F00-TWO.txt
@@ -1767,13 +1767,13 @@ HLCY
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-3000.
+2
+3000. 1000.
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-0.
+2
+0. 0.
 ?
 ?
 0
@@ -2015,9 +2015,9 @@ surface
 ?
 ?
 ?
-180
-VIS_ON_SURFACE
-?
+410
+GSD_VIS_ON_SURFACE
+GSD_visibility on surface
 1
 tmpl4_0
 VIS

diff --git a/parm/postxconfig-NT-GFS-F00.txt b/parm/postxconfig-NT-GFS-F00.txt
@@ -1766,13 +1766,13 @@ HLCY
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-3000.
+2
+3000. 1000.
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-0.
+2
+0. 0.
 ?
 ?
 0
@@ -2014,9 +2014,9 @@ surface
 ?
 ?
 ?
-180
-VIS_ON_SURFACE
-?
+410
+GSD_VIS_ON_SURFACE
+GSD_visibility on surface
 1
 tmpl4_0
 VIS

diff --git a/parm/postxconfig-NT-GFS-TWO.txt b/parm/postxconfig-NT-GFS-TWO.txt
@@ -1767,13 +1767,13 @@ HLCY
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-3000.
+2
+3000. 1000.
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-0.
+2
+0. 0.
 ?
 ?
 0
@@ -2607,9 +2607,9 @@ entire_atmos_single_lyr
 ?
 ?
 ?
-180
-VIS_ON_SURFACE
-?
+410
+GSD_VIS_ON_SURFACE
+GSD_visibility on surface
 1
 tmpl4_0
 VIS

diff --git a/parm/postxconfig-NT-GFS.txt b/parm/postxconfig-NT-GFS.txt
@@ -1766,13 +1766,13 @@ HLCY
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-3000.
+2
+3000. 1000.
 spec_hgt_lvl_above_grnd
 0
 ?
-1
-0.
+2
+0. 0.
 ?
 ?
 0
@@ -2606,9 +2606,9 @@ entire_atmos_single_lyr
 ?
 ?
 ?
-180
-VIS_ON_SURFACE
-?
+410
+GSD_VIS_ON_SURFACE
+GSD_visibility on surface
 1
 tmpl4_0
 VIS

diff --git a/sorc/ncep_post.fd/MDLFLD.f b/sorc/ncep_post.fd/MDLFLD.f
@@ -44,7 +44,8 @@
 !!   20-11-10  J MENG - USE UPP_PHYSICS MODULE
 !!   21-04-01  J MENG - COMPUTATION ON DEFINED POINTS ONLY
 !!   21-07-07  J MENG - 2D DECOMPOSITION
-!!   22-09-22  L Zhang- ADD NO3 and NH4 output for UFS-Aerosols model
+!!   22-09-22  L Zhang - ADD NO3 and NH4 output for UFS-Aerosols model
+!!   22-10-20  W Meng - Bug fix for cloud fraction and vertically integrated liquid
 !!
 !! USAGE:    CALL MDLFLD
 !!   INPUT ARGUMENT LIST:
@@ -1217,8 +1218,11 @@ SUBROUTINE MDLFLD
 !$omp parallel do private(i,j)
                DO J=JSTA,JEND
                  DO I=ista,iend
-                   IF(abs(CFR(I,J,LL)-SPVAL) > SMALL)                   &
-     &                 GRID1(I,J) = CFR(I,J,LL)*H100
+                   IF(abs(CFR(I,J,LL)-SPVAL) > SMALL) THEN
+                       GRID1(I,J) = CFR(I,J,LL)*H100
+                   ELSE
+                       GRID1(I,J) = SPVAL
+                   ENDIF
                  ENDDO
                ENDDO
                CALL BOUND(GRID1,D00,H100)
@@ -3051,6 +3055,8 @@ SUBROUTINE MDLFLD
               if(zint(i,j,l) < spval .and.zint(i,j,l+1)<spval.and.DBZ(I,J,L)<spval) then
                 GRID1(I,J)=GRID1(I,J)+0.00344* &
                 (10.**(DBZ(I,J,L)/10.))**0.57143*(ZINT(I,J,L)-ZINT(I,J,L+1))/1000.
+              else
+                GRID1(I,J)=spval
               endif
             ENDDO
           ENDDO