Skip to content

Commit

Permalink
Updating a continue statement in fv_nudge (NOAA-EMC#153)
Browse files Browse the repository at this point in the history
* update in fv_nudge to fix 666 label logic
  • Loading branch information
laurenchilutti committed Oct 7, 2021
1 parent 14a6bf5 commit c2decd7
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions tools/fv_nudge.F90
Original file line number Diff line number Diff line change
Expand Up @@ -782,7 +782,7 @@ subroutine ps_nudging(dt, factor, factor_nwp, npz, ak, bk, ps_obs, mask, tm, ps,
if ( kmax < km ) call mpp_error(FATAL,'==> KMAX must be larger than km')

do j=js,je
do 666 i=is,ie
do i=is,ie
do k=1, km+1
pk0(k) = (ak0(k) + bk0(k)*ps_obs(i,j))**kappa
enddo
Expand All @@ -800,8 +800,8 @@ subroutine ps_nudging(dt, factor, factor_nwp, npz, ak, bk, ps_obs, mask, tm, ps,
pt0 = tm(i,j)/(pk0(km+1)-pk0(km))*(kappa*(pn0(km+1)-pn0(km)))
pst = pk0(km+1) + (gz0(i,j)-phis(i,j))/(cp_air*pt0)
endif
666 continue ! i-loop
ps_dt(i,j) = pst**(1./kappa) - ps(i,j)
666 ps_dt(i,j) = pst**(1./kappa) - ps(i,j)
enddo ! i-loop
enddo ! j-loop

if( nf_ps>0 ) call del2_scalar(ps_dt, del2_cd, 1, nf_ps, bd, npx, npy, gridstruct, domain)
Expand Down

0 comments on commit c2decd7

Please sign in to comment.