Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Changes to output variables at 160 m and 320 m above ground #693

Merged
merged 6 commits into from
Apr 19, 2023
Merged

Changes to output variables at 160 m and 320 m above ground #693

Show file tree
Hide file tree
Changes from 5 commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
90ce1d7
Changes to the extinction from smoke and dust; now reading in "ext550…
EricJames-NOAA Apr 13, 2023
77c5b30
Enabling output of 160m and 320m winds (and other variables) for rene…
EricJames-NOAA Apr 13, 2023
b6699f1
Merge remote-tracking branch 'origin/develop' into ext550_changes
EricJames-NOAA Apr 17, 2023
9ccd977
Updating change logs.
EricJames-NOAA Apr 17, 2023
502ac43
Cleaning up references to AODTOT, which is no longer needed.
EricJames-NOAA Apr 18, 2023
74f3be1
Updating change log to indicate aodtot is removed.
EricJames-NOAA Apr 18, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 5 additions & 5 deletions parm/fv3lam_rrfs.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1158,7 +1158,7 @@
<param>
<shortname>TMP_ON_SPEC_HGT_LVL_ABOVE_GRND_FDHGT</shortname>
<pname>TMP</pname>
<level>30. 50. 80. 100.</level>
<level>30. 50. 80. 100. 160. 320.</level>
<scale>3.0</scale>
</param>

Expand All @@ -1172,7 +1172,7 @@
<param>
<shortname>UGRD_ON_SPEC_HGT_LVL_ABOVE_GRND_FDHGT</shortname>
<pname>UGRD</pname>
<level>30. 50. 80. 100.</level>
<level>30. 50. 80. 100. 160. 320.</level>
<scale>4.0</scale>
</param>

Expand All @@ -1186,7 +1186,7 @@
<param>
<shortname>VGRD_ON_SPEC_HGT_LVL_ABOVE_GRND_FDHGT</shortname>
<pname>VGRD</pname>
<level>30. 50. 80. 100.</level>
<level>30. 50. 80. 100. 160. 320.</level>
<scale>4.0</scale>
</param>

Expand All @@ -1200,14 +1200,14 @@
<param>
<shortname>SPFH_ON_SPEC_HGT_LVL_ABOVE_GRND_FDHGT</shortname>
<pname>SPFH</pname>
<level>30. 50. 80. 100.</level>
<level>30. 50. 80. 100. 160. 320.</level>
<scale>5.0</scale>
</param>

<param>
<shortname>PRES_ON_SPEC_HGT_LVL_ABOVE_GRND_FDHGT</shortname>
<pname>PRES</pname>
<level>30. 50. 80. 100.</level>
<level>30. 50. 80. 100. 160. 320.</level>
<scale>3.0</scale>
</param>

Expand Down
20 changes: 10 additions & 10 deletions parm/postxconfig-NT-fv3lam_rrfs.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6171,8 +6171,8 @@ TMP
spec_hgt_lvl_above_grnd
0
?
4
30. 50. 80. 100.
6
30. 50. 80. 100. 160. 320.
?
0
?
Expand Down Expand Up @@ -6245,8 +6245,8 @@ UGRD
spec_hgt_lvl_above_grnd
0
?
4
30. 50. 80. 100.
6
30. 50. 80. 100. 160. 320.
?
0
?
Expand Down Expand Up @@ -6319,8 +6319,8 @@ VGRD
spec_hgt_lvl_above_grnd
0
?
4
30. 50. 80. 100.
6
30. 50. 80. 100. 160. 320.
?
0
?
Expand Down Expand Up @@ -6393,8 +6393,8 @@ SPFH
spec_hgt_lvl_above_grnd
0
?
4
30. 50. 80. 100.
6
30. 50. 80. 100. 160. 320.
?
0
?
Expand Down Expand Up @@ -6430,8 +6430,8 @@ PRES
spec_hgt_lvl_above_grnd
0
?
4
30. 50. 80. 100.
6
30. 50. 80. 100. 160. 320.
?
0
?
Expand Down
2 changes: 0 additions & 2 deletions sorc/ncep_post.fd/ALLOCATE_ALL.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -615,7 +615,6 @@ SUBROUTINE ALLOCATE_ALL()
allocate(mean_frp(ista_2l:iend_2u,jsta_2l:jend_2u))
allocate(ebb(ista_2l:iend_2u,jsta_2l:jend_2u))
allocate(hwp(ista_2l:iend_2u,jsta_2l:jend_2u))
allocate(aodtot(ista_2l:iend_2u,jsta_2l:jend_2u))
!Initialization
!$omp parallel do private(i,j)
do j=jsta_2l,jend_2u
Expand Down Expand Up @@ -643,7 +642,6 @@ SUBROUTINE ALLOCATE_ALL()
mean_frp(i,j)=spval
ebb(i,j)=spval
hwp(i,j)=spval
aodtot(i,j)=spval
enddo
enddo
allocate(smoke(ista_2l:iend_2u,jsta_2l:jend_2u,lm,nbin_sm))
Expand Down
12 changes: 3 additions & 9 deletions sorc/ncep_post.fd/CLDRAD.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -73,6 +73,7 @@
!> 2023-02-10 | Eric James | Removing neighbourhood check from GSL exp2 ceiling diagnostic
!> 2023-02-23 | Eric James | Adding coarse PM from RRFS, and using AOD from FV3 for RRFS
!> 2023-04-04 | Li(Kate Zhang) | Add namelist optoin for CCPP-Chem (UFS-Chem) model.
!> 2023-04-17 | Eric James | Getting rid of special treatment for RRFS AOD (use RAP/HRRR approach)
!>
!> @author Russ Treadon W/NP2 @date 1993-08-30
SUBROUTINE CLDRAD
Expand All @@ -97,7 +98,7 @@ SUBROUTINE CLDRAD
ALWINC, ALWTOAC, SWDDNI, SWDDIF, SWDNBC, SWDDNIC, &
SWDDIFC, SWUPBC, LWDNBC, LWUPBC, SWUPT, &
TAOD5502D, AERSSA2D, AERASY2D, MEAN_FRP, EBB, HWP, &
AODTOT, LWP, IWP, AVGCPRATE, &
LWP, IWP, AVGCPRATE, &
DUSTCB,SSCB,BCCB,OCCB,SULFCB,DUSTPM,SSPM,aod550, &
du_aod550,ss_aod550,su_aod550,oc_aod550,bc_aod550, &
PWAT,DUSTPM10,MAOD,NO3CB,NH4CB,aqm_aod550
Expand Down Expand Up @@ -446,16 +447,9 @@ SUBROUTINE CLDRAD
! TOTAL COLUMN AOD (TAOD553D FROM HRRR-SMOKE)
!
IF (IGET(735) > 0) THEN
IF (MODELNAME == 'RAPR') THEN
IF (MODELNAME == 'RAPR' .OR. MODELNAME == 'FV3R') THEN
CALL CALPW(GRID1(ista:iend,jsta:jend),19)
CALL BOUND(GRID1,D00,H99999)
ELSE IF (MODELNAME == 'FV3R') THEN
GRID1=SPVAL
DO J=JSTA,JEND
DO I=ISTA,IEND
if (AODTOT(I,J) < SPVAL) GRID1(I,J) = AODTOT(I,J)
WenMeng-NOAA marked this conversation as resolved.
Show resolved Hide resolved
ENDDO
ENDDO
ENDIF
if(grib == "grib2" )then
cfld = cfld + 1
Expand Down
5 changes: 3 additions & 2 deletions sorc/ncep_post.fd/CTLBLK.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ module CTLBLK_mod
! 2023-03-21 Jesse Meng - Add slrutah_on option to use U Utah SLR
! 2023-04-04 Li(Kate Zhang) Add namelist optoin for CCPP-Chem (UFS-Chem)
! and 2D diag. output (d2d_chem) for GEFS-Aerosols and CCPP-Chem model.
! 2023-04-17 Eric James - Adding 160 and 320 m above ground to HTFD for RRFS output.
!-----------------------------------------------------------------------
!
implicit none
Expand All @@ -29,7 +30,7 @@ module CTLBLK_mod
end type
integer, parameter :: komax=70
integer, parameter :: LSMDEF=46 ! default number of p levels
integer,PARAMETER :: NFD=18,NBND=6
integer,PARAMETER :: NFD=20,NBND=6
REAL, PARAMETER :: QMIN = 1.E-15
!
integer :: novegtype ! max number of veg type
Expand Down Expand Up @@ -115,7 +116,7 @@ module CTLBLK_mod
integer, parameter :: nbin_sm = 1 ! smoke
!
! SET FD LEVEL HEIGHTS IN GEOPOTENTAL METERS.
DATA HTFD / 20.E0,30.E0,40.E0,50.E0,80.E0,100.E0,305.E0,457.E0,610.E0, &
DATA HTFD / 20.E0,30.E0,40.E0,50.E0,80.E0,100.E0,160.E0,305.E0,320.E0,457.E0,610.E0, &
914.E0,1524.E0,1829.E0,2134.E0,2743.E0,3658.E0,4572.E0, &
6000.E0,7010.E0/
!
Expand Down
1 change: 0 additions & 1 deletion sorc/ncep_post.fd/DEALLOCATE.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -209,7 +209,6 @@ SUBROUTINE DE_ALLOCATE
deallocate(mean_frp)
deallocate(ebb)
deallocate(hwp)
deallocate(aodtot)
deallocate(smoke)
deallocate(fv3dust)
deallocate(coarsepm)
Expand Down
30 changes: 11 additions & 19 deletions sorc/ncep_post.fd/INITPOST_NETCDF.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -31,8 +31,9 @@
!> 2023-02-23 | Eric James | Read coarse PM and aodtot from RRFS
!> 2023-03-02 | Sam Trahan | Read lightning threat index fields
!> 2023-03-22 | WM Lewis | Read RRFS effective radii (EFFRL, EFFRI, EFFRS)
!> !> 2023-04-04 |Li(Kate Zhang) |Add namelist optoin for CCPP-Chem(UFS-Chem)
!> 2023-04-04 |Li(Kate Zhang) |Add namelist optoin for CCPP-Chem(UFS-Chem)
! and 2D diag. output (d2d_chem) for GEFS-Aerosols and CCPP-Chem model.
!> 2023-04-17 | Eric James | Read in unified ext550 extinction for RRFS
WenMeng-NOAA marked this conversation as resolved.
Show resolved Hide resolved
!>
!> @author Hui-Ya Chuang @date 2016-03-04
SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
Expand Down Expand Up @@ -73,7 +74,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, aodtot, aqm_aod550, ltg1_max,ltg2_max,ltg3_max
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, &
Expand Down Expand Up @@ -200,11 +201,10 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
real, allocatable :: div3d(:,:,:)
real(kind=4),allocatable :: vcrd(:,:)
real :: dum_const
real, allocatable :: extsmoke(:,:,:), extdust(:,:,:)
real, allocatable :: ext550(:,:,:)

if (modelname == 'FV3R') then
allocate(extsmoke(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
allocate(extdust(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
allocate(ext550(ista_2l:iend_2u,jsta_2l:jend_2u,lm))
endif

!***********************************************************************
Expand Down Expand Up @@ -506,13 +506,13 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
end if
if(me==0)print*,'nhcas= ',nhcas
if (nhcas == 0 ) then !non-hydrostatic case
nrec=19
nrec=18
allocate (recname(nrec))
recname=[character(len=20) :: 'ugrd','vgrd','spfh','tmp','o3mr', &
'presnh','dzdt', 'clwmr','dpres', &
'delz','icmr','rwmr', &
'snmr','grle','smoke','dust', &
'coarsepm','smoke_ext','dust_ext']
'coarsepm','ext550']
else
nrec=8
allocate (recname(nrec))
Expand Down Expand Up @@ -874,9 +874,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
call read_netcdf_3d_para(ncid3d,im,jm,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,recname(17),coarsepm(ista_2l,jsta_2l,1,1),lm)
call read_netcdf_3d_para(ncid2d,im,jm,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,recname(18),extsmoke(ista_2l,jsta_2l,1),lm)
call read_netcdf_3d_para(ncid2d,im,jm,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,recname(19),extdust(ista_2l,jsta_2l,1),lm)
spval,recname(18),ext550(ista_2l,jsta_2l,1),lm)
endif

! calculate CWM from FV3 output
Expand Down Expand Up @@ -1024,11 +1022,6 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,hwp(ista_2l,jsta_2l))
if(debugprint)print*,'sample ',VarName,' =',hwp(isa,jsa)
! total aod
VarName='aodtot'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,aodtot(ista_2l,jsta_2l))
if(debugprint)print*,'sample ',VarName,' =',aodtot(isa,jsa)
endif

! lightning threat index 1
Expand Down Expand Up @@ -2472,8 +2465,8 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
do l = 1, lm
do j = jsta_2l, jend_2u
do i = ista_2l, iend_2u
if(extsmoke(i,j,l)<spval.and.extdust(i,j,l)<spval)then
taod5503d ( i, j, l) = extsmoke ( i, j, l ) + extdust ( i, j, l )
if(ext550(i,j,l)<spval)then
taod5503d ( i, j, l) = ext550 ( i, j, l )
dz = ZINT( i, j, l ) - ZINT( i, j, l+1 )
aextc55 ( i, j, l ) = taod5503d ( i, j, l ) / dz
endif
Expand All @@ -2486,8 +2479,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
end do
end do
end do
deallocate(extsmoke)
deallocate(extdust)
deallocate(ext550)
end if

!$omp parallel do private(i,j)
Expand Down
2 changes: 1 addition & 1 deletion sorc/ncep_post.fd/VRBLS2D_mod.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,7 @@ module vrbls2d
,SWNORMmean(:,:),SNFDEN(:,:),SNDEPAC(:,:),SWDDNI(:,:),SWDDIF(:,:) &
,SWDNBC(:,:),SWDDNIC(:,:),SWDDIFC(:,:), SWUPBC(:,:), SWUPT(:,:) &
,TAOD5502D(:,:),AERASY2D(:,:),AERSSA2D(:,:),MEAN_FRP(:,:),EBB(:,:) &
,HWP(:,:),AODTOT(:,:),LWP(:,:),IWP(:,:) &
,HWP(:,:),LWP(:,:),IWP(:,:) &
! add new fields for GFS
,SFCUX(:,:),SFCVX(:,:),SFCUXI(:,:), SFCVXI(:,:),AVGALBEDO(:,:),AVGCPRATE(:,:) &
,AVGPREC(:,:),PTOP(:,:),PBOT(:,:),AVGCFRACH(:,:) &
Expand Down