From 0f0302ac04199bdd678eaa6dec798750d3ea3590 Mon Sep 17 00:00:00 2001 From: Jim Bresch Date: Sat, 13 Jan 2018 22:12:21 -0700 Subject: [PATCH] The fix for the July 2017 gfs changes should not be applied to the gfs supplemental pressure-level files (gfsb) since they do not contain masked surface fields. Currently, ungrib crashes when fed gfsb files. A work-around would be to cat the gfsb file to the corresponding gfs file. This modified code only calls the subroutine when a gfs file contains a soil temperature field. --- ungrib/src/rrpr.F | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/ungrib/src/rrpr.F b/ungrib/src/rrpr.F index 0950667..fe4b3e5 100644 --- a/ungrib/src/rrpr.F +++ b/ungrib/src/rrpr.F @@ -758,8 +758,10 @@ subroutine rrpr(hstart, ntimes, interval, nlvl, maxlvl, plvl, & if (index(map%source,'NCEP GFS') .ne. 0 ) then call mprintf(.true.,DEBUG, & "RRPR: Adjusting GFS masked fields ") - call get_dims(200100, 'ST000010') - call fix_gfs_miss (map%nx, map%ny, 200100.) + if ( is_there(200100, 'ST000010')) then + call get_dims(200100, 'ST000010') + call fix_gfs_miss (map%nx, map%ny, 200100.) + endif endif ! Add residual soil moisture to SOILM* if initialized from the GSD RUC model or from NCEP RUC -- 2.11.4.GIT