Skip to content

Commit

Permalink
Minor tidy-up in response to Russ Treadon's review.
Browse files Browse the repository at this point in the history
  • Loading branch information
ClaraDraper-NOAA committed Feb 7, 2024
1 parent 3001f11 commit 304e36c
Show file tree
Hide file tree
Showing 3 changed files with 7 additions and 9 deletions.
2 changes: 0 additions & 2 deletions src/enkf/statevec.f90
Original file line number Diff line number Diff line change
Expand Up @@ -159,8 +159,6 @@ subroutine init_statevec()
print *, 'ns3d: ', ns3d, ', ns2d: ', ns2d, ', nsdim: ', nsdim
endif



call getgridinfo(statefileprefixes(nstatefields/2+1), .false.)

end subroutine init_statevec
Expand Down
12 changes: 6 additions & 6 deletions src/gsi/read_prepbufr.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1981,16 +1981,16 @@ subroutine read_prepbufr(nread,ndata,nodata,infile,obstype,lunout,twindin,sis,&
pqm(k)=2 ! otherwise, type 183 will be discarded.
qm=2
tqm(k)=2
if (kx==187) obserr(3,k)=2.0
if (kx==181) obserr(3,k)=2.0
if (kx==183) obserr(3,k)=2.0
if (kx==187) obserr(3,k)=2.0_r_double
if (kx==181) obserr(3,k)=2.0_r_double
if (kx==183) obserr(3,k)=2.0_r_double
endif
if (qob .and. qm == 9 ) then
qm = 2
! qob err specified as fraction of qsat, multiplied by 10.
if (kx==187) obserr(2,k)=1.0
if (kx==181) obserr(2,k)=1.0
if (kx==183) obserr(2,k)=1.0
if (kx==187) obserr(2,k)=1.0_r_double
if (kx==181) obserr(2,k)=1.0_r_double
if (kx==183) obserr(2,k)=1.0_r_double
endif
endif
! Set usage variable
Expand Down
2 changes: 1 addition & 1 deletion src/gsi/setupq.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1016,7 +1016,7 @@ subroutine setupq(obsLL,odiagLL,lunin,mype,bwork,awork,nele,nobs,is,conv_diagsav
my_head => null()
ENDDO

endif ! i_use_2mq4b block
endif ! l_pbl_pseudo_itype
!!!!!!!!!!!!!!!!!! PBL pseudo surface obs !!!!!!!!!!!!!!!!!!!!!!!

! End of loop over observations
Expand Down

0 comments on commit 304e36c

Please sign in to comment.