Enable output of VEGMIN/VEGMAX (#1021)

* Updates to enable output of VEGMIN and VEGMAX from RRFS. * Adding change log entry.
NOAA-EMC · Aug 26, 2024 · 897fcb0 · 897fcb0
1 parent fc85241
commit 897fcb0
Show file tree

Hide file tree

Showing 3 changed files with 107 additions and 2 deletions.
diff --git a/parm/postxconfig-NT-rrfs.txt b/parm/postxconfig-NT-rrfs.txt
@@ -1,7 +1,7 @@
 4
 5
 3
-229
+231
 278
 PRSLEV
 32769
@@ -20237,6 +20237,90 @@ surface
 ?
 ?
 ?
+726
+VEG_MIN_ON_SURFACE
+?
+1
+tmpl4_0
+VEGMIN
+NCEP
+?
+surface
+0
+?
+0
+?
+?
+0
+?
+0
+?
+?
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+0
+0.0
+0
+0.0
+1
+3.0
+0
+0
+0
+?
+?
+?
+729
+VEG_MAX_ON_SURFACE
+?
+1
+tmpl4_0
+VEGMAX
+NCEP
+?
+surface
+0
+?
+0
+?
+?
+0
+?
+0
+?
+?
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+0
+0.0
+0
+0.0
+1
+3.0
+0
+0
+0
+?
+?
+?
 262
 INST_CSDSF_ON_SURFACE
 ?

diff --git a/parm/rrfs_postcntrl.xml b/parm/rrfs_postcntrl.xml
@@ -3198,6 +3198,16 @@
        <scale>-3.0</scale>
        </param>
 
+       <param>
+       <shortname>VEG_MIN_ON_SURFACE</shortname>
+       <scale>3.0</scale>
+       </param>
+
+       <param>
+       <shortname>VEG_MAX_ON_SURFACE</shortname>
+       <scale>3.0</scale>
+       </param>
+
        <param>
        <shortname>INST_CSDSF_ON_SURFACE</shortname>
        <pname>CSDSF</pname>

diff --git a/sorc/ncep_post.fd/INITPOST_NETCDF.f b/sorc/ncep_post.fd/INITPOST_NETCDF.f
@@ -54,6 +54,7 @@
 !> 2024-04-23 | Eric James    | Updating smoke emissions to be 3D variable (ebu_smoke)
 !> 2024-05-01 | Eric James    | set "prec_acc_dt1" as 15 min for RRFS
 !> 2024-05-09 | Eric James    | Enable reading of clear-sky downwelling shortwave irradiance
+!> 2024-08-20 | Eric James    | Enable reading of seasonal max/min green vegetation fraction
 !>
 !> @author Hui-Ya Chuang @date 2016-03-04
 !----------------------------------------------------------------------
@@ -92,7 +93,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
               uz0, vz0, ptop, htop, pbot, hbot, ptopl, pbotl, ttopl, ptopm, pbotm, ttopm,       &
               ptoph, pboth, pblcfr, ttoph, runoff, tecan, tetran, tedir, twa, maxtshltr,        &
               mintshltr, maxrhshltr, fdnsst, acgraup, graup_bucket, acfrain, frzrn_bucket,      &
-              snow_acm, snow_bkt, snownc, graupelnc, qrmax, swddif, swddni, xlaixy,             &
+              snow_acm, snow_bkt, snownc, graupelnc, qrmax, swddif, swddni, xlaixy, shdmin, shdmax,&
               minrhshltr, dzice, smcwlt, suntime, fieldcapa, htopd, hbotd, htops, hbots,        &
               cuppt, dusmass, ducmass, dusmass25, ducmass25, aswintoa,rel_vort_maxhy1,          &
               maxqshltr, minqshltr, acond, sr, u10h, v10h,refd_max, w_up_max, w_dn_max,         &
@@ -2796,6 +2797,16 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
       call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
       spval,VarName,xlaixy)
 
+! leaf area index
+      VarName='shdmin'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,shdmin)
+
+! leaf area index
+      VarName='shdmax'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,shdmax)
+
 ! time averaged incoming sfc uv-b using getgb
       VarName='duvb_ave'
       call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &