Skip to content

Commit

Permalink
Update upp submodule for fixing in-line post issues for RRFS and HAFS (
Browse files Browse the repository at this point in the history
…#558)

* update upp revision to 8767929

* modify dx/dy calculation for RRFS with domain over north pole
  • Loading branch information
WenMeng-NOAA authored Jul 15, 2022
1 parent dc74ab1 commit b9d61f2
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 4 deletions.
11 changes: 8 additions & 3 deletions io/post_fv3.F90
Original file line number Diff line number Diff line change
Expand Up @@ -614,13 +614,18 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
call exch(gdlon)

!$omp parallel do default(none),private(i,j,ip1), &
!$omp& shared(jsta,jend_m,im,dx,gdlat,gdlon,dy,ista,iend_m)
!$omp& shared(jsta,jend_m,im,dx,gdlat,gdlon,dy,ista,iend_m,maptype,dxval,dyval,gdsdegr)
do j = jsta, jend_m
do i = ista, iend_m
ip1 = i + 1
!if (ip1 > im) ip1 = ip1 - im
dx(i,j) = erad*cos(gdlat(i,j)*dtr)*(gdlon(ip1,j)-gdlon(i,j))*dtr
dy(i,j) = erad*(gdlat(i,j+1)-gdlat(i,j))*dtr ! like A*DPH
if(maptype==207)then
dx(i,j)=erad*dxval*dtr/gdsdegr
dy(i,j)=erad*dyval*dtr/gdsdegr
else
dx(i,j) = erad*cos(gdlat(i,j)*dtr)*(gdlon(ip1,j)-gdlon(i,j))*dtr
dy(i,j) = erad*(gdlat(i,j+1)-gdlat(i,j))*dtr ! like A*DPH
endif
end do
end do
!
Expand Down

0 comments on commit b9d61f2

Please sign in to comment.