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

Add tref #419

Merged
merged 6 commits into from
Dec 2, 2021
Merged

Add tref #419

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
c1cf9df
Add foundation temperature per request from coupling group
Nov 15, 2021
0601de0
add founfation T to inline post control file
Nov 16, 2021
6ac0f9f
add founfation T to inline post control file flux file part instead
Nov 16, 2021
3cc969a
Change sea ice threshold for identifying sea ice from .1 to 0
Nov 23, 2021
7c0984a
Resolved merge conflict by removing RQSTFLD.F
Nov 30, 2021
dbc6f48
made updates to only use new ice fraction threshold to identify sea
Dec 2, 2021
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
8 changes: 8 additions & 0 deletions parm/post_avblflds.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4719,6 +4719,14 @@
<scale>3.0</scale>
</param>

<param>
<post_avblfldidx>549</post_avblfldidx>
<shortname>FDNSSTMP_ON_SURFACE</shortname>
<pname>FDNSSTMP</pname>
<fixed_sfc1_type>surface</fixed_sfc1_type>
<scale>4.0</scale>
</param>

<param>
<post_avblfldidx>551</post_avblfldidx>
<shortname>CSNOW_ON_SURFACE</shortname>
Expand Down
5 changes: 5 additions & 0 deletions parm/postcntrl_gfs_f00_two.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -963,6 +963,11 @@
<scale>4.0</scale>
</param>

<param>
<shortname>FDNSSTMP_ON_SURFACE</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>TSOIL_ON_DEPTH_BEL_LAND_SFC</shortname>
<scale_fact_fixed_sfc1>2 2 2 2</scale_fact_fixed_sfc1>
Expand Down
5 changes: 5 additions & 0 deletions parm/postcntrl_gfs_flux.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,6 +54,11 @@
<scale>4.0</scale>
</param>

<param>
<shortname>FDNSSTMP_ON_SURFACE</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>TSOIL_ON_DEPTH_BEL_LAND_SFC</shortname>
<scale_fact_fixed_sfc1>2 2 2 2</scale_fact_fixed_sfc1>
Expand Down
5 changes: 5 additions & 0 deletions parm/postcntrl_gfs_flux_f00.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,6 +54,11 @@
<scale>4.0</scale>
</param>

<param>
<shortname>FDNSSTMP_ON_SURFACE</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>TSOIL_ON_DEPTH_BEL_LAND_SFC</shortname>
<scale_fact_fixed_sfc1>2 2 2 2</scale_fact_fixed_sfc1>
Expand Down
5 changes: 5 additions & 0 deletions parm/postcntrl_gfs_two.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1215,6 +1215,11 @@
<scale>4.0</scale>
</param>

<param>
<shortname>FDNSSTMP_ON_SURFACE</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>TSOIL_ON_DEPTH_BEL_LAND_SFC</shortname>
<scale_fact_fixed_sfc1>2 2 2 2</scale_fact_fixed_sfc1>
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-GFS-F00-TWO.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
2
46
47
150
GFSPRS
0
Expand Down Expand Up @@ -5842,6 +5842,43 @@ surface
?
?
?
549
FDNSSTMP_ON_SURFACE
?
1
tmpl4_0
FDNSSTMP
?
?
surface
0
?
0
?
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
116
TSOIL_ON_DEPTH_BEL_LAND_SFC
?
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-GFS-FLUX-F00.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
1
46
47
GFSPRS
0
ncep_nco
Expand Down Expand Up @@ -275,6 +275,43 @@ surface
?
?
?
549
FDNSSTMP_ON_SURFACE
?
1
tmpl4_0
FDNSSTMP
?
?
surface
0
?
0
?
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
116
TSOIL_ON_DEPTH_BEL_LAND_SFC
?
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-GFS-FLUX.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
1
111
112
GFSPRS
0
ncep_nco
Expand Down Expand Up @@ -275,6 +275,43 @@ surface
?
?
?
549
FDNSSTMP_ON_SURFACE
?
1
tmpl4_0
FDNSSTMP
?
?
surface
0
?
0
?
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
116
TSOIL_ON_DEPTH_BEL_LAND_SFC
?
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-GFS-TWO.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
2
111
112
197
GFSPRS
0
Expand Down Expand Up @@ -7581,6 +7581,43 @@ surface
?
?
?
549
FDNSSTMP_ON_SURFACE
?
1
tmpl4_0
FDNSSTMP
?
?
surface
0
?
0
?
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
116
TSOIL_ON_DEPTH_BEL_LAND_SFC
?
Expand Down
2 changes: 2 additions & 0 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 @@ -969,6 +969,7 @@ SUBROUTINE ALLOCATE_ALL()
allocate(tetran(im,jsta_2l:jend_2u))
allocate(tedir(im,jsta_2l:jend_2u))
allocate(twa(im,jsta_2l:jend_2u))
allocate(fdnsst(im,jsta_2l:jend_2u))
!Initialization
!$omp parallel do private(i,j)
do j=jsta_2l,jend_2u
Expand Down Expand Up @@ -1018,6 +1019,7 @@ SUBROUTINE ALLOCATE_ALL()
tetran(i,j)=spval
tedir(i,j)=spval
twa(i,j)=spval
fdnsst(i,j)=spval
enddo
enddo
!
Expand Down
1 change: 1 addition & 0 deletions sorc/ncep_post.fd/DEALLOCATE.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -385,6 +385,7 @@ SUBROUTINE DE_ALLOCATE
deallocate(tetran)
deallocate(tedir)
deallocate(twa)
deallocate(fdnsst)
! GSD
deallocate(rainc_bucket)
deallocate(rainc_bucket1)
Expand Down
20 changes: 18 additions & 2 deletions sorc/ncep_post.fd/FIXED.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,8 @@ SUBROUTINE FIXED

!
use vrbls3d, only: pint
use vrbls2d, only: albedo, avgalbedo, albase, mxsnal, sst, ths, epsr, ti
use vrbls2d, only: albedo, avgalbedo, albase, mxsnal, sst, ths, epsr, ti&
, fdnsst
use masks, only: gdlat, gdlon, sm, sice, lmh, lmv
use params_mod, only: small, p1000, capa
use lookup_mod, only: ITB,JTB,ITBQ,JTBQ
Expand Down Expand Up @@ -109,7 +110,7 @@ SUBROUTINE FIXED
DO I = 1,IM
GRID1(I,J) = SPVAL
IF(SM(I,J) /= SPVAL) GRID1(I,J) = 1. - SM(I,J)
IF(SICE(I,J) /= SPVAL .AND. SICE(I,J) > 0.1) GRID1(I,J) = 0.
IF(SICE(I,J) /= SPVAL .AND. SICE(I,J) > 0.0) GRID1(I,J) = 0.
WenMeng-NOAA marked this conversation as resolved.
Show resolved Hide resolved
! if(j==jm/2)print*,'i,mask= ',i,grid1(i,j)
ENDDO
ENDDO
Expand Down Expand Up @@ -331,6 +332,21 @@ SUBROUTINE FIXED
datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
endif
ENDIF
!
! FOUNDATION TEMPERAURE.
IF (IGET(549)>0) THEN
!$omp parallel do private(i,j)
DO J=JSTA,JEND
DO I=1,IM
GRID1(I,J) = FDNSST(I,J)
ENDDO
ENDDO
if(grib=='grib2') then
cfld=cfld+1
fld_info(cfld)%ifld=IAVBLFLD(IGET(549))
datapd(1:im,1:jend-jsta+1,cfld)=GRID1(1:im,jsta:jend)
endif
ENDIF

! EMISSIVIT.
IF (IGET(248)>0) THEN
Expand Down
9 changes: 7 additions & 2 deletions sorc/ncep_post.fd/INITPOST_GFS_NETCDF_PARA.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -57,7 +57,7 @@ SUBROUTINE INITPOST_GFS_NETCDF_PARA(ncid3d)
smstot, ivgtyp, isltyp, sfcevp, sfcexc, acsnow, acsnom, sst, thz0, qz0, &
uz0, vz0, ptop, htop, pbot, hbot, ptopl, pbotl, ttopl, ptopm, pbotm, ttopm, &
ptoph, pboth, pblcfr, ttoph, runoff, tecan, tetran, tedir, twa, maxtshltr, &
mintshltr, maxrhshltr, &
mintshltr, maxrhshltr, fdnsst, &
minrhshltr, dzice, smcwlt, suntime, fieldcapa, htopd, hbotd, htops, hbots, &
cuppt, dusmass, ducmass, dusmass25, ducmass25, aswintoa,rel_vort_maxhy1, &
maxqshltr, minqshltr, acond, sr, u10h, v10h,refd_max, w_up_max, w_dn_max, &
Expand Down Expand Up @@ -1021,7 +1021,12 @@ SUBROUTINE INITPOST_GFS_NETCDF_PARA(ncid3d)
enddo
if(debugprint)print*,'sample ',VarName,' = ',ths(isa,jsa)


! foundation temperature
VarName='tref'
call read_netcdf_2d_para(ncid2d,im,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,fdnsst)
if(debugprint)print*,'sample ',VarName,' = ',fdnsst(isa,jsa)

! GFS does not have time step and physics time step, make up ones since they
! are not really used anyway
! NPHS=1.
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 @@ -82,7 +82,7 @@ module vrbls2d
,avgesnow(:,:),avgpotevp(:,:),avgprec_cont(:,:),avgcprate_cont(:,:)&
,ti(:,:),aod550(:,:),du_aod550(:,:),ss_aod550(:,:),su_aod550(:,:) &
,bc_aod550(:,:),oc_aod550(:,:),landfrac(:,:),paha(:,:),pahi(:,:) &
,tecan(:,:),tetran(:,:),tedir(:,:),twa(:,:)
,tecan(:,:),tetran(:,:),tedir(:,:),twa(:,:),fdnsst(:,:)
integer, allocatable :: IVGTYP(:,:),ISLTYP(:,:),ISLOPE(:,:) &
,IEQL(:,:)

Expand Down