From ff6f39aa7053b5a37af03ef17a453161671933b8 Mon Sep 17 00:00:00 2001
From: KaiWang-NOAA <109170734+KaiWang-NOAA@users.noreply.github.com>
Date: Wed, 11 Jun 2025 14:57:47 -0400
Subject: [PATCH] A temporary solution to fix erroneous values of LAI from
 noah-mp

Fix erroneous values of LAI coming from Noah-MP LSM (over snow/ice land use types) that breaks the AQM for some grid cells in the NA domain when testing the GFSv17 CCPP package with Noah-MP LSM for AQM
---
 src/shr/aqm_methods.F90 | 15 ++++++++++++++-
 1 file changed, 14 insertions(+), 1 deletion(-)

diff --git a/src/shr/aqm_methods.F90 b/src/shr/aqm_methods.F90
index 9406d6fd..4d888bab 100644
--- a/src/shr/aqm_methods.F90
+++ b/src/shr/aqm_methods.F90
@@ -1221,7 +1221,20 @@ LOGICAL FUNCTION  XTRACT3 ( FNAME, VNAME,                           &
       case ("HFX")
         p2d => stateIn % hfx
       case ("LAI")
-        p2d => stateIn % xlai
+!        p2d => stateIn % xlai
+!     a temporary solution provided by Patrick Campbell from ARL for GFSv17 CCPP and
+!     AQMv8 coupling
+        k = 0
+        do r = row0, row1
+          do c = col0, col1  
+            k = k + 1
+               if ( stateIn % xlai(c,r) > 10.0 ) then  !zero out erroneously high values
+                 buffer(k) = 0.0 * stateIn % xlai(c,r) !from Noah-MP option, which should
+               else                                    !max out at ~ 4.5 from tables and
+                 buffer(k) = stateIn % xlai(c,r)       !~ 10 from satellite input data
+               end if
+          end do 
+        end do   
       case ("LH")
         p2d => stateIn % lh
       case ("PRSFC")