@@ -1694,7 +1694,7 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
1694
1694
do k= 1 ,levs
1695
1695
tvflg(k)= one ! initialize as sensible
1696
1696
do j= 1 ,20
1697
- if (.not. isnan(vtcd) .and. tpc(k,j)==vtcd) tvflg(k)= zero ! reset flag if virtual
1697
+ if (tpc(k,j)==vtcd) tvflg(k)= zero ! reset flag if virtual
1698
1698
if (tpc(k,j)>= bmiss) exit ! end of stack
1699
1699
end do
1700
1700
end do
@@ -1713,7 +1713,7 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
1713
1713
exit
1714
1714
end if
1715
1715
end if
1716
- if (.not. isnan(vtcd) .and. tpc(k,j)==vtcd) then
1716
+ if (tpc(k,j)==vtcd) then
1717
1717
obsdat(3 ,k)= tobaux(1 ,k,j+1 )
1718
1718
qcmark(3 ,k)= min (tobaux(2 ,k,j+1 ),qcmark_huge)
1719
1719
tqm(k)= nint (qcmark(3 ,k))
@@ -1737,7 +1737,7 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
1737
1737
tobs4q(1 ,k)= obsdat(3 ,k) ! temp obs read in prepbufr
1738
1738
tqm4q(1 ,k)= tqm(k)
1739
1739
do j= 1 ,20
1740
- if (.not. isnan(vtcd) .and. tpc(k,j)==vtcd) tvflg4q(1 ,k)= zero ! reset flag if virtual
1740
+ if (tpc(k,j)==vtcd) tvflg4q(1 ,k)= zero ! reset flag if virtual
1741
1741
if (tpc(k,j)>= bmiss) exit ! end of stack
1742
1742
end do
1743
1743
end do
@@ -3413,7 +3413,7 @@ subroutine sonde_ext(obsdat,tpc,qcmark,obserr,drfdat,levsio,kx,vtcd)
3413
3413
do k= 1 ,levs
3414
3414
tvflg(k)= one ! initialize as sensible
3415
3415
do j= 1 ,20
3416
- if (.not. isnan(vtcd) .and. tpc(k,j)==vtcd) tvflg(k)= zero ! reset flag if virtual
3416
+ if (tpc(k,j)==vtcd) tvflg(k)= zero ! reset flag if virtual
3417
3417
if (tpc(k,j)>= bmiss) exit ! end of stack
3418
3418
end do
3419
3419
end do
0 commit comments