Skip to content

Fix Lifted Index calculation in RRFS - #1630

Open
BenjaminBlake-NOAA wants to merge 7 commits into
NOAA-EMC:developfrom
BenjaminBlake-NOAA:feature/lftx
Open

BenjaminBlake-NOAA wants to merge 7 commits into
NOAA-EMC:developfrom
BenjaminBlake-NOAA:feature/lftx

Conversation

@BenjaminBlake-NOAA

@BenjaminBlake-NOAA BenjaminBlake-NOAA commented Oct 21, 2026 •

Copy link
Copy Markdown
Collaborator

This PR resolves issue #1629

A user reported that the values for Lifted Index (UPP ID 030) in the RRFS output have an apparent temperature offset, and thus are not consistent with the expected temperature-difference scale. In CLDRAD.f, it was discovered that TFRZ (273.15 K) is being added to the Lifted Index. Here are the values of Lifted Index before and after my code changes:

143:717243998:vt=2026032700:500-1000 mb:36 hour fcst:LFTX Surface Lifted Index [K]:
    ndata=14452641:undef=0:mean=281.216:min=262.3:max=314.9
143:717243998:vt=2026032700:500-1000 mb:36 hour fcst:LFTX Surface Lifted Index [K]:
    ndata=14452641:undef=0:mean=8.05802:min=-10.8:max=41.8

@WenMeng-NOAA @ShellieRowe-NOAA I have RRFS run directories here on WCOSS2 and Ursa if you'd like to take a look:

WCOSS2: /lfs/h2/emc/vpppg/noscrub/Benjamin.Blake/upp-pr1630-lftx/rrfs_2026032512
Ursa: /scratch4/NCEPDEV/stmp/Benjamin.Blake/upp-pr1630-intel/rrfs_2026032512

LFTX is in the 2DFLD and NATLEV GRIB2 files.

Tagging @AndrewBenjamin-NOAA for awareness.

@BenjaminBlake-NOAA BenjaminBlake-NOAA added bug Something isn't working Baseline Change The baselines of the UPP regression tests are changed. RRFS labels Oct 21, 2026
@WenMeng-NOAA WenMeng-NOAA linked an issue Oct 21, 2026 that may be closed by this pull request
@ShellieRowe-NOAA

ShellieRowe-NOAA commented Oct 21, 2026 via email

Copy link
Copy Markdown

@BenjaminBlake-NOAA

Copy link
Copy Markdown
Collaborator Author

@ShellieRowe-NOAA Thanks for testing and reviewing the changes! I think we can ignore the flagged specific humidity field since all of the values are undefined.

@WenMeng-NOAA

Copy link
Copy Markdown
Collaborator

@BenjaminBlake-NOAA Your changes look good to me and my test results are consistent with yours.

@WenMeng-NOAA

Copy link
Copy Markdown
Collaborator

@gspetro-NOAA This PR is ready for the UPP RTs on RDHPCS. Please refer to @BenjaminBlake-NOAA 's testing:
/scratch4/NCEPDEV/stmp/Benjamin.Blake/upp-pr1630-intel/rrfs_2026032512

The only baseline change is the rrfs test.

@BenjaminBlake-NOAA

Copy link
Copy Markdown
Collaborator Author

@WenMeng-NOAA @ShellieRowe-NOAA Here are several plots illustrating the differences in the lifted index.

Before changes (all values are greater than 260K):
image

After changes:
image

And here is the Best (4 layer) Lifted Index which compares well with the Lifted Index:
image

@gspetro-NOAA

Copy link
Copy Markdown
Collaborator

Just an update: Baselines have been updated successfully for Ursa (llvm & intel classic) and Hercules, and everything passes on rerun. For some reason, the baseline creation run on Orion keeps timing out. I'm hoping with my latest increase in wallclock time, this next run will pass, and then the test time will be normal once the baselines are updated, but I have to wait for the run to complete.

@BenjaminBlake-NOAA

Copy link
Copy Markdown
Collaborator Author

@gspetro-NOAA Sounds good, thanks!

@gspetro-NOAA gspetro-NOAA left a comment

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I updated the BL files for rrfs and all tests are now passing on all relevant systems. I unfortunately forgot to push the Orion logs from before I updated the baselines, so I pushed the ones from after. The rt_verify.pl printout from before the BL update showed the rrfs test running significantly overtime on Orion, but after the baseline update, RTs ran normally.

Since I overwrote those logs, here is the printout:

#-----------------
There were no errors detected in the work directory
#-----------------

Checking timing results for all regresssion tests...

#-----------------
gefsv13 is 67   seconds UNDER time
nmmb    is 91   seconds UNDER time
hfip    is 85   seconds UNDER time
gfs     is 741  seconds UNDER time
mpas    is 82   seconds UNDER time
rrfs    is 2027 seconds OVER  time
hafs    is 17   seconds UNDER time
sfs     is 53   seconds UNDER time
rap     is 16   seconds UNDER time
gcafs   is 23   seconds UNDER time
gefsv12 is 36   seconds UNDER time
hrrr    is 256  seconds UNDER time
aqm     is 15   seconds UNDER time
3drtma  is 66   seconds UNDER time
aigfs   is 17   seconds UNDER time
#-----------------

Checking for baseline files that have changed...

#-----------------
rrfs:
	1) /work2/noaa/epic/gpetro/orion/RTs/upp-rts/1630/ci/rundir/upp-ORION/rrfs_2026032512/NATLEV36.tm00
[sudo] password for gpetro: 
	   ***New baseline file copied successfully***
	2) /work2/noaa/epic/gpetro/orion/RTs/upp-rts/1630/ci/rundir/upp-ORION/rrfs_2026032512/2DFLD36.tm00
	   ***New baseline file copied successfully***

This branch has not been deployed

No deployments
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Baseline Change The baselines of the UPP regression tests are changed. bug Something isn't working RRFS

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Lifted Index has incorrect values in RRFS

4 participants