Updating GRIB2 labeling for MSTAV, and adding SOILL on 9 RUC levels. (#…

…740) * Updating GRIB2 labeling for MSTAV, and adding SOILL on 9 RUC levels. * Entering change log for INITPOST_NETCDF.f * Changing shortname of MSTAV to specify RUC LSM (defined at 0 m below ground)
NOAA-EMC · Jul 13, 2023 · 83e59b5 · 83e59b5
1 parent 1398ca8
commit 83e59b5
Show file tree

Hide file tree

Showing 4 changed files with 183 additions and 38 deletions.
diff --git a/parm/fv3lam_rrfs.xml b/parm/fv3lam_rrfs.xml
@@ -391,6 +391,16 @@
        <scale>3.0</scale>
        </param>
 
+       <param>
+       <shortname>SOILL_ON_DEPTH_BEL_LAND_SFC</shortname>
+       <pname>SOILL</pname>
+       <scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
+       <level>0. 1. 4. 10. 30. 60. 100. 160. 300.</level>
+       <scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
+       <level2>0. 1. 4. 10. 30. 60. 100. 160. 300.</level2>
+       <scale>3.0</scale>
+       </param>
+
        <param>
        <shortname>CNWAT_ON_SURFACE</shortname>
        <pname>CNWAT</pname>
@@ -540,13 +550,7 @@
        </param>
 
        <param>
-       <shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
-       <pname>MSTAV</pname>
-       <table_info>NCEP</table_info>
-       <scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
-       <level>0.</level>
-       <scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
-       <level2>100.</level2>
+       <shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
        <scale>3.0</scale>
        </param>
 
@@ -2306,6 +2310,16 @@
        <scale>3.0</scale>
        </param>
 
+       <param>
+       <shortname>SOILL_ON_DEPTH_BEL_LAND_SFC</shortname>
+       <pname>SOILL</pname>
+       <scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
+       <level>0. 1. 4. 10. 30. 60. 100. 160. 300.</level>
+       <scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
+       <level2>0. 1. 4. 10. 30. 60. 100. 160. 300.</level2>
+       <scale>3.0</scale>
+       </param>
+
        <param>
        <shortname>SFEXC_ON_SURFACE</shortname>
        <pname>SFEXC</pname>
@@ -2320,13 +2334,7 @@
        </param>
 
        <param>
-       <shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
-       <pname>MSTAV</pname>
-       <table_info>NCEP</table_info>
-       <scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
-       <level>0.</level>
-       <scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
-       <level2>100.</level2>
+       <shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
        <scale>3.0</scale>
        </param>
 

diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml
@@ -1508,12 +1508,10 @@
 
     <param>
        <post_avblfldidx>171</post_avblfldidx>
-       <shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
+       <shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
        <pname>MSTAV</pname>
+       <table_info>NCEP</table_info>
        <fixed_sfc1_type>depth_bel_land_sfc</fixed_sfc1_type>
-       <scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
-       <fixed_sfc2_type>depth_bel_land_sfc</fixed_sfc2_type>
-       <scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
        <scale>3.0</scale>
     </param>
 

diff --git a/parm/postxconfig-NT-fv3lam_rrfs.txt b/parm/postxconfig-NT-fv3lam_rrfs.txt
@@ -1,7 +1,7 @@
 3
 4
-237
-287
+238
+288
 PRSLEV
 32769
 ncep_nco
@@ -1905,6 +1905,43 @@ depth_bel_land_sfc
 ?
 ?
 ?
+225
+SOILL_ON_DEPTH_BEL_LAND_SFC
+?
+1
+tmpl4_0
+SOILL
+?
+?
+depth_bel_land_sfc
+1
+2
+9
+0. 1. 4. 10. 30. 60. 100. 160. 300.
+depth_bel_land_sfc
+1
+2
+9
+0. 1. 4. 10. 30. 60. 100. 160. 300.
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+1
+3.0
+0
+0
+0
+?
+?
+?
 118
 CNWAT_ON_SURFACE
 ?
@@ -2728,15 +2765,15 @@ MSTAV
 NCEP
 ?
 depth_bel_land_sfc
-1
-2
-1
-0.
-depth_bel_land_sfc
-1
-2
-1
-100.
+0
+?
+0
+?
+?
+0
+?
+0
+?
 ?
 ?
 0
@@ -12540,6 +12577,43 @@ depth_bel_land_sfc
 ?
 ?
 ?
+225
+SOILL_ON_DEPTH_BEL_LAND_SFC
+?
+1
+tmpl4_0
+SOILL
+?
+?
+depth_bel_land_sfc
+1
+2
+9
+0. 1. 4. 10. 30. 60. 100. 160. 300.
+depth_bel_land_sfc
+1
+2
+9
+0. 1. 4. 10. 30. 60. 100. 160. 300.
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+1
+3.0
+0
+0
+0
+?
+?
+?
 169
 SFEXC_ON_SURFACE
 ?
@@ -12623,15 +12697,15 @@ MSTAV
 NCEP
 ?
 depth_bel_land_sfc
-1
-2
-1
-0.
-depth_bel_land_sfc
-1
-2
-1
-100.
+0
+?
+0
+?
+?
+0
+?
+0
+?
 ?
 ?
 0

diff --git a/sorc/ncep_post.fd/INITPOST_NETCDF.f b/sorc/ncep_post.fd/INITPOST_NETCDF.f
@@ -37,6 +37,7 @@
 !> 2023-04-21 | Eric James    | Read in / calculate some fields needed for GSL p-type diagnosis for RRFS
 !> 2023-05-31 | Wen Meng      | Bug fix in qrmax initialization
 !> 2023-06-14 | Wen Meng      ! Bug fix of reading seaswtc and modification of sndepac calculation
+!> 2023-07-06 | Eric James    ! Read in SOILL on 9 levels for RRFS
 !>
 !> @author Hui-Ya Chuang @date 2016-03-04
       SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
@@ -2245,6 +2246,70 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
       enddo
      if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,4)
 
+      IF (NSOIL==9) THEN
+
+      VarName='soill5'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,sh2o(ista_2l,jsta_2l,5))
+!     mask water areas
+!$omp parallel do private(i,j)
+      do j=jsta,jend
+        do i=ista,iend
+          if (sm(i,j) /= 0.0) sh2o(i,j,5) = spval
+        enddo
+      enddo
+     if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,5)
+
+      VarName='soill6'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,sh2o(ista_2l,jsta_2l,6))
+!     mask water areas
+!$omp parallel do private(i,j)
+      do j=jsta,jend
+        do i=ista,iend
+          if (sm(i,j) /= 0.0) sh2o(i,j,6) = spval
+        enddo
+      enddo
+     if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,6)
+
+      VarName='soill7'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,sh2o(ista_2l,jsta_2l,7))
+!     mask water areas
+!$omp parallel do private(i,j)
+      do j=jsta,jend
+        do i=ista,iend
+          if (sm(i,j) /= 0.0) sh2o(i,j,7) = spval
+        enddo
+      enddo
+     if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,7)
+
+      VarName='soill8'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,sh2o(ista_2l,jsta_2l,8))
+!     mask water areas
+!$omp parallel do private(i,j)
+      do j=jsta,jend
+        do i=ista,iend
+          if (sm(i,j) /= 0.0) sh2o(i,j,8) = spval
+        enddo
+      enddo
+     if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,8)
+
+      VarName='soill9'
+      call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+      spval,VarName,sh2o(ista_2l,jsta_2l,9))
+!     mask water areas
+!$omp parallel do private(i,j)
+      do j=jsta,jend
+        do i=ista,iend
+          if (sm(i,j) /= 0.0) sh2o(i,j,9) = spval
+        enddo
+      enddo
+     if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,9)
+
+      END IF
+
 ! volumetric soil moisture using nemsio
       VarName='soilw1'
       call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &