Skip to content

Commit 83e59b5

Browse files
EricJames-NOAAEricJames-NOAA
authored
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)
1 parent 1398ca8 commit 83e59b5

File tree

4 files changed

+183
-38
lines changed

4 files changed

+183
-38
lines changed

‎parm/fv3lam_rrfs.xml

Lines changed: 22 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -391,6 +391,16 @@
391391
<scale>3.0</scale>
392392
</param>
393393

394+
<param>
395+
<shortname>SOILL_ON_DEPTH_BEL_LAND_SFC</shortname>
396+
<pname>SOILL</pname>
397+
<scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
398+
<level>0. 1. 4. 10. 30. 60. 100. 160. 300.</level>
399+
<scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
400+
<level2>0. 1. 4. 10. 30. 60. 100. 160. 300.</level2>
401+
<scale>3.0</scale>
402+
</param>
403+
394404
<param>
395405
<shortname>CNWAT_ON_SURFACE</shortname>
396406
<pname>CNWAT</pname>
@@ -540,13 +550,7 @@
540550
</param>
541551

542552
<param>
543-
<shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
544-
<pname>MSTAV</pname>
545-
<table_info>NCEP</table_info>
546-
<scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
547-
<level>0.</level>
548-
<scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
549-
<level2>100.</level2>
553+
<shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
550554
<scale>3.0</scale>
551555
</param>
552556

@@ -2306,6 +2310,16 @@
23062310
<scale>3.0</scale>
23072311
</param>
23082312

2313+
<param>
2314+
<shortname>SOILL_ON_DEPTH_BEL_LAND_SFC</shortname>
2315+
<pname>SOILL</pname>
2316+
<scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
2317+
<level>0. 1. 4. 10. 30. 60. 100. 160. 300.</level>
2318+
<scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
2319+
<level2>0. 1. 4. 10. 30. 60. 100. 160. 300.</level2>
2320+
<scale>3.0</scale>
2321+
</param>
2322+
23092323
<param>
23102324
<shortname>SFEXC_ON_SURFACE</shortname>
23112325
<pname>SFEXC</pname>
@@ -2320,13 +2334,7 @@
23202334
</param>
23212335

23222336
<param>
2323-
<shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
2324-
<pname>MSTAV</pname>
2325-
<table_info>NCEP</table_info>
2326-
<scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
2327-
<level>0.</level>
2328-
<scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
2329-
<level2>100.</level2>
2337+
<shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
23302338
<scale>3.0</scale>
23312339
</param>
23322340

‎parm/post_avblflds.xml

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1508,12 +1508,10 @@
15081508

15091509
<param>
15101510
<post_avblfldidx>171</post_avblfldidx>
1511-
<shortname>MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
1511+
<shortname>RUC_MSTAV_ON_DEPTH_BEL_LAND_SFC</shortname>
15121512
<pname>MSTAV</pname>
1513+
<table_info>NCEP</table_info>
15131514
<fixed_sfc1_type>depth_bel_land_sfc</fixed_sfc1_type>
1514-
<scale_fact_fixed_sfc1>2</scale_fact_fixed_sfc1>
1515-
<fixed_sfc2_type>depth_bel_land_sfc</fixed_sfc2_type>
1516-
<scale_fact_fixed_sfc2>2</scale_fact_fixed_sfc2>
15171515
<scale>3.0</scale>
15181516
</param>
15191517

‎parm/postxconfig-NT-fv3lam_rrfs.txt

Lines changed: 94 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
3
22
4
3-
237
4-
287
3+
238
4+
288
55
PRSLEV
66
32769
77
ncep_nco
@@ -1905,6 +1905,43 @@ depth_bel_land_sfc
19051905
?
19061906
?
19071907
?
1908+
225
1909+
SOILL_ON_DEPTH_BEL_LAND_SFC
1910+
?
1911+
1
1912+
tmpl4_0
1913+
SOILL
1914+
?
1915+
?
1916+
depth_bel_land_sfc
1917+
1
1918+
2
1919+
9
1920+
0. 1. 4. 10. 30. 60. 100. 160. 300.
1921+
depth_bel_land_sfc
1922+
1
1923+
2
1924+
9
1925+
0. 1. 4. 10. 30. 60. 100. 160. 300.
1926+
?
1927+
?
1928+
0
1929+
0.0
1930+
0
1931+
0.0
1932+
?
1933+
0
1934+
0.0
1935+
0
1936+
0.0
1937+
1
1938+
3.0
1939+
0
1940+
0
1941+
0
1942+
?
1943+
?
1944+
?
19081945
118
19091946
CNWAT_ON_SURFACE
19101947
?
@@ -2728,15 +2765,15 @@ MSTAV
27282765
NCEP
27292766
?
27302767
depth_bel_land_sfc
2731-
1
2732-
2
2733-
1
2734-
0.
2735-
depth_bel_land_sfc
2736-
1
2737-
2
2738-
1
2739-
100.
2768+
0
2769+
?
2770+
0
2771+
?
2772+
?
2773+
0
2774+
?
2775+
0
2776+
?
27402777
?
27412778
?
27422779
0
@@ -12540,6 +12577,43 @@ depth_bel_land_sfc
1254012577
?
1254112578
?
1254212579
?
12580+
225
12581+
SOILL_ON_DEPTH_BEL_LAND_SFC
12582+
?
12583+
1
12584+
tmpl4_0
12585+
SOILL
12586+
?
12587+
?
12588+
depth_bel_land_sfc
12589+
1
12590+
2
12591+
9
12592+
0. 1. 4. 10. 30. 60. 100. 160. 300.
12593+
depth_bel_land_sfc
12594+
1
12595+
2
12596+
9
12597+
0. 1. 4. 10. 30. 60. 100. 160. 300.
12598+
?
12599+
?
12600+
0
12601+
0.0
12602+
0
12603+
0.0
12604+
?
12605+
0
12606+
0.0
12607+
0
12608+
0.0
12609+
1
12610+
3.0
12611+
0
12612+
0
12613+
0
12614+
?
12615+
?
12616+
?
1254312617
169
1254412618
SFEXC_ON_SURFACE
1254512619
?
@@ -12623,15 +12697,15 @@ MSTAV
1262312697
NCEP
1262412698
?
1262512699
depth_bel_land_sfc
12626-
1
12627-
2
12628-
1
12629-
0.
12630-
depth_bel_land_sfc
12631-
1
12632-
2
12633-
1
12634-
100.
12700+
0
12701+
?
12702+
0
12703+
?
12704+
?
12705+
0
12706+
?
12707+
0
12708+
?
1263512709
?
1263612710
?
1263712711
0

‎sorc/ncep_post.fd/INITPOST_NETCDF.f

Lines changed: 65 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -37,6 +37,7 @@
3737
!> 2023-04-21 | Eric James | Read in / calculate some fields needed for GSL p-type diagnosis for RRFS
3838
!> 2023-05-31 | Wen Meng | Bug fix in qrmax initialization
3939
!> 2023-06-14 | Wen Meng ! Bug fix of reading seaswtc and modification of sndepac calculation
40+
!> 2023-07-06 | Eric James ! Read in SOILL on 9 levels for RRFS
4041
!>
4142
!> @author Hui-Ya Chuang @date 2016-03-04
4243
SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
@@ -2245,6 +2246,70 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
22452246
enddo
22462247
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,4)
22472248
2249+
IF (NSOIL==9) THEN
2250+
2251+
VarName='soill5'
2252+
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
2253+
spval,VarName,sh2o(ista_2l,jsta_2l,5))
2254+
! mask water areas
2255+
!$omp parallel do private(i,j)
2256+
do j=jsta,jend
2257+
do i=ista,iend
2258+
if (sm(i,j) /= 0.0) sh2o(i,j,5) = spval
2259+
enddo
2260+
enddo
2261+
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,5)
2262+
2263+
VarName='soill6'
2264+
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
2265+
spval,VarName,sh2o(ista_2l,jsta_2l,6))
2266+
! mask water areas
2267+
!$omp parallel do private(i,j)
2268+
do j=jsta,jend
2269+
do i=ista,iend
2270+
if (sm(i,j) /= 0.0) sh2o(i,j,6) = spval
2271+
enddo
2272+
enddo
2273+
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,6)
2274+
2275+
VarName='soill7'
2276+
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
2277+
spval,VarName,sh2o(ista_2l,jsta_2l,7))
2278+
! mask water areas
2279+
!$omp parallel do private(i,j)
2280+
do j=jsta,jend
2281+
do i=ista,iend
2282+
if (sm(i,j) /= 0.0) sh2o(i,j,7) = spval
2283+
enddo
2284+
enddo
2285+
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,7)
2286+
2287+
VarName='soill8'
2288+
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
2289+
spval,VarName,sh2o(ista_2l,jsta_2l,8))
2290+
! mask water areas
2291+
!$omp parallel do private(i,j)
2292+
do j=jsta,jend
2293+
do i=ista,iend
2294+
if (sm(i,j) /= 0.0) sh2o(i,j,8) = spval
2295+
enddo
2296+
enddo
2297+
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,8)
2298+
2299+
VarName='soill9'
2300+
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
2301+
spval,VarName,sh2o(ista_2l,jsta_2l,9))
2302+
! mask water areas
2303+
!$omp parallel do private(i,j)
2304+
do j=jsta,jend
2305+
do i=ista,iend
2306+
if (sm(i,j) /= 0.0) sh2o(i,j,9) = spval
2307+
enddo
2308+
enddo
2309+
if(debugprint)print*,'sample l',VarName,' = ',1,sh2o(isa,jsa,9)
2310+
2311+
END IF
2312+
22482313
! volumetric soil moisture using nemsio
22492314
VarName='soilw1'
22502315
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &

0 commit comments

Comments
 (0)