diff --git a/parm/fv3lam.xml b/parm/fv3lam.xml
index f0287bfa5..38667958e 100755
--- a/parm/fv3lam.xml
+++ b/parm/fv3lam.xml
@@ -1762,13 +1762,6 @@
4.0
-
- LTNG_ON_SURFACE
- LTNG
- NCEP
- 1.0
-
-
VIL_ON_ENTIRE_ATMOS
VIL
@@ -3404,6 +3397,13 @@
3.0
+
+ GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+ LTNG
+ NCEP
+ 4.0
+
+
SDEN_ON_SURFACE
SDEN
diff --git a/parm/fv3lam_post_avblflds.xml b/parm/fv3lam_post_avblflds.xml
index ff331dbc0..552d9bb63 100644
--- a/parm/fv3lam_post_avblflds.xml
+++ b/parm/fv3lam_post_avblflds.xml
@@ -5189,7 +5189,10 @@
702
GSD_MAX_LTG_THREAT1_ON_ENTIRE_ATMOS
+ tmpl4_8
LTGTHREAT1
+ NCEP
+ MAX
entire_atmos
3.0
@@ -5197,7 +5200,10 @@
703
GSD_MAX_LTG_THREAT2_ON_ENTIRE_ATMOS
+ tmpl4_8
LTGTHREAT2
+ NCEP
+ MAX
entire_atmos
3.0
@@ -5205,7 +5211,10 @@
704
GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+ tmpl4_8
LTNG
+ NCEP
+ MAX
entire_atmos
3.0
diff --git a/parm/fv3lam_rrfs.xml b/parm/fv3lam_rrfs.xml
index dc09f31f9..908143b7c 100755
--- a/parm/fv3lam_rrfs.xml
+++ b/parm/fv3lam_rrfs.xml
@@ -1774,13 +1774,6 @@
4.0
-
- LTNG_ON_SURFACE
- LTNG
- NCEP
- 1.0
-
-
VIL_ON_ENTIRE_ATMOS
VIL
@@ -3612,6 +3605,12 @@
4.0
+
+ GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+ LTNG
+ 3.0
+
+
SDEN_ON_SURFACE
SDEN
diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml
index 591f1be5b..3149c1b5d 100755
--- a/parm/post_avblflds.xml
+++ b/parm/post_avblflds.xml
@@ -6552,7 +6552,10 @@
702
GSD_MAX_LTG_THREAT1_ON_ENTIRE_ATMOS
+ tmpl4_8
LTGTHREAT1
+ NCEP
+ MAX
entire_atmos
3.0
@@ -6560,7 +6563,10 @@
703
GSD_MAX_LTG_THREAT2_ON_ENTIRE_ATMOS
+ tmpl4_8
LTGTHREAT2
+ NCEP
+ MAX
entire_atmos
3.0
@@ -6568,7 +6574,10 @@
704
GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+ tmpl4_8
LTNG
+ NCEP
+ MAX
entire_atmos
3.0
diff --git a/parm/postxconfig-NT-3drtma.txt b/parm/postxconfig-NT-3drtma.txt
index d480e0ba7..a35a3ba50 100644
--- a/parm/postxconfig-NT-3drtma.txt
+++ b/parm/postxconfig-NT-3drtma.txt
@@ -3759,10 +3759,10 @@ equil_lvl
GSD_MAX_LTG_THREAT1_ON_ENTIRE_ATMOS
?
1
-tmpl4_0
+tmpl4_8
LTGTHREAT1
-?
-?
+NCEP
+MAX
entire_atmos
0
?
@@ -3796,10 +3796,10 @@ entire_atmos
GSD_MAX_LTG_THREAT2_ON_ENTIRE_ATMOS
?
1
-tmpl4_0
+tmpl4_8
LTGTHREAT2
-?
-?
+NCEP
+MAX
entire_atmos
0
?
@@ -3833,10 +3833,10 @@ entire_atmos
GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
?
1
-tmpl4_0
+tmpl4_8
LTNG
-?
-?
+NCEP
+MAX
entire_atmos
0
?
diff --git a/parm/postxconfig-NT-fv3lam.txt b/parm/postxconfig-NT-fv3lam.txt
index 5881aa4a1..80d5f117d 100644
--- a/parm/postxconfig-NT-fv3lam.txt
+++ b/parm/postxconfig-NT-fv3lam.txt
@@ -1,6 +1,6 @@
2
-220
-264
+221
+263
PRSLEV
32769
ncep_nco
@@ -9341,43 +9341,6 @@ surface
?
?
?
-444
-LTNG_ON_SURFACE
-lightning
-1
-tmpl4_0
-LTNG
-NCEP
-?
-surface
-0
-?
-0
-?
-?
-0
-?
-0
-?
-?
-?
-0
-0.0
-0
-0.0
-?
-0
-0.0
-0
-0.0
-1
-1.0
-0
-0
-0
-?
-?
-?
581
VIL_ON_ENTIRE_ATMOS
entire atmosphere Vertically Integrated Liquid (kg/m-2)
@@ -17904,6 +17867,43 @@ surface
?
?
?
+704
+GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+?
+1
+tmpl4_8
+LTNG
+NCEP
+MAX
+entire_atmos
+0
+?
+0
+?
+?
+0
+?
+0
+?
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+1
+4.0
+0
+0
+0
+?
+?
+?
1006
SDEN_ON_SURFACE
?
diff --git a/parm/postxconfig-NT-fv3lam_rrfs.txt b/parm/postxconfig-NT-fv3lam_rrfs.txt
index 098022236..a953f30e7 100644
--- a/parm/postxconfig-NT-fv3lam_rrfs.txt
+++ b/parm/postxconfig-NT-fv3lam_rrfs.txt
@@ -1,7 +1,7 @@
3
4
-241
-289
+242
+288
PRSLEV
32769
ncep_nco
@@ -9638,43 +9638,6 @@ surface
?
?
?
-444
-LTNG_ON_SURFACE
-lightning
-1
-tmpl4_0
-LTNG
-NCEP
-?
-surface
-0
-?
-0
-?
-?
-0
-?
-0
-?
-?
-?
-0
-0.0
-0
-0.0
-?
-0
-0.0
-0
-0.0
-1
-1.0
-0
-0
-0
-?
-?
-?
581
VIL_ON_ENTIRE_ATMOS
entire atmosphere Vertically Integrated Liquid (kg/m-2)
@@ -19607,6 +19570,43 @@ top_of_atmos
?
?
?
+704
+GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
+?
+1
+tmpl4_8
+LTNG
+NCEP
+MAX
+entire_atmos
+0
+?
+0
+?
+?
+0
+?
+0
+?
+?
+?
+0
+0.0
+0
+0.0
+?
+0
+0.0
+0
+0.0
+1
+3.0
+0
+0
+0
+?
+?
+?
1006
SDEN_ON_SURFACE
?
diff --git a/sorc/ncep_post.fd/INITPOST_NETCDF.f b/sorc/ncep_post.fd/INITPOST_NETCDF.f
index e6220ff87..3d3419fb8 100644
--- a/sorc/ncep_post.fd/INITPOST_NETCDF.f
+++ b/sorc/ncep_post.fd/INITPOST_NETCDF.f
@@ -28,6 +28,7 @@
!> 2022-12-07 | Wen Meng | Read AOD from AQM model
!> 2022-12-23 | Eric Aligo | Read six winter weather diagnostics from model
!> 2023-01-30 | Sam Trahan | Read cldfra or cldfra_bl, whichever is available
+!> 2023-03-02 | Sam Trahan | Read lightning threat index fields
!>
!> @author Hui-Ya Chuang @date 2016-03-04
SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
@@ -66,7 +67,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
alwoutc,alwtoac,aswoutc,aswtoac,alwinc,aswinc,avgpotevp,snoavg, &
ti,aod550,du_aod550,ss_aod550,su_aod550,oc_aod550,bc_aod550,prate_max,maod,dustpm10, &
dustcb,bccb,occb,sulfcb,sscb,dustallcb,ssallcb,dustpm,sspm,pp25cb,pp10cb,no3cb,nh4cb,&
- pwat, ebb, hwp, aqm_aod550
+ pwat, ebb, hwp, aqm_aod550, ltg1_max,ltg2_max,ltg3_max
use soil, only: sldpth, sllevel, sh2o, smc, stc
use masks, only: lmv, lmh, htm, vtm, gdlat, gdlon, dx, dy, hbm2, sm, sice
use physcons_post, only: grav => con_g, fv => con_fvirt, rgas => con_rd, &
@@ -1002,6 +1003,24 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
if(debugprint)print*,'sample ',VarName,' =',hwp(isa,jsa)
endif
+! lightning threat index 1
+ VarName='ltg1_max'
+ call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+ spval,VarName,ltg1_max(ista_2l,jsta_2l))
+ if(debugprint)print*,'sample ',VarName,' =',ltg1_max(isa,jsa)
+
+! lightning threat index 2
+ VarName='ltg2_max'
+ call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+ spval,VarName,ltg2_max(ista_2l,jsta_2l))
+ if(debugprint)print*,'sample ',VarName,' =',ltg2_max(isa,jsa)
+
+! lightning threat index 3
+ VarName='ltg3_max'
+ call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
+ spval,VarName,ltg3_max(ista_2l,jsta_2l))
+ if(debugprint)print*,'sample ',VarName,' =',ltg3_max(isa,jsa)
+
! surface pressure
VarName='pressfc'
call read_netcdf_2d_para(ncid3d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
diff --git a/sorc/ncep_post.fd/MDL2AGL.f b/sorc/ncep_post.fd/MDL2AGL.f
index 55d97b07c..82a0fc551 100644
--- a/sorc/ncep_post.fd/MDL2AGL.f
+++ b/sorc/ncep_post.fd/MDL2AGL.f
@@ -17,6 +17,7 @@
!! 21-04-01 J MENG - computation on defined points only
!! 21-07-26 W Meng - Restrict computation from undefined grids
!! 21-10-14 J MENG - 2D DECOMPOSITION
+!! 2023-03-02 S TRAHAN - copy lightning threat index 3 element-by-element
!!
!! USAGE: CALL MDL2P
!! INPUT ARGUMENT LIST:
@@ -771,8 +772,15 @@ SUBROUTINE MDL2AGL
fld_info(cfld)%tinvstat = 1
endif
fld_info(cfld)%ntrange = 1
- datapd(1:iend-ista+1,1:jend-jsta+1,cfld)=GRID1(ista:iend,jsta:jend)
+ do j=1,jend-jsta+1
+ jj = jsta+j-1
+ do i=1,iend-ista+1
+ ii = ista+i-1
+ datapd(i,j,cfld) = GRID1(ii,jj)
+ enddo
+ enddo
endif
+ !datapd(1:iend-ista+1,1:jend-jsta+1,cfld)=GRID1(ista:iend,jsta:jend)
END IF
!--- GSD Updraft Helicity