Skip to content

Commit

Permalink
Remove warnings from grib1_utils.F90
Browse files Browse the repository at this point in the history
  • Loading branch information
GeorgeGayno-NOAA committed Dec 16, 2022
1 parent 6e8c131 commit 67181ae
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions sorc/chgres_cube.fd/grib2_util.F90
Original file line number Diff line number Diff line change
Expand Up @@ -56,7 +56,7 @@ subroutine rh2spfh(rh_sphum,p,t)
!print *, 'es = ', es
e = rh * es / 100.0
!print *, 'e = ', e
rh_sphum = 0.622 * e / p
rh_sphum = real((0.622 * e / p),kind=esmf_kind_r4)
!print *, 'q = ', sphum

!if (P .eq. 100000.0) THEN
Expand Down Expand Up @@ -110,7 +110,7 @@ subroutine rh2spfh_gfs(rh_sphum,p,t)
do i=1,i_input
ES = MIN(FPVSNEW(T(I,J)),P)
QC(i,j) = CON_EPS*ES/(P+CON_EPSM1*ES)
rh_sphum(i,j) = rh(i,j)*QC(i,j)/100.0
rh_sphum(i,j) = real((rh(i,j)*QC(i,j)/100.0),kind=esmf_kind_r4)
end do
end do

Expand Down Expand Up @@ -169,7 +169,7 @@ elemental function fpvsnew(t)
c1xpvs=1.-xmin*c2xpvs
! xj=min(max(c1xpvs+c2xpvs*t,1.0),real(nxpvs,krealfp))
xj=min(max(c1xpvs+c2xpvs*t,1.0),float(nxpvs))
jx=min(xj,float(nxpvs)-1.0)
jx=nint(min(xj,float(nxpvs)-1.0))
x=xmin+(jx-1)*xinc

tr=con_ttp/x
Expand Down

0 comments on commit 67181ae

Please sign in to comment.