Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Reduce computational patch size for intermediate domain for vertical nesting #37

Merged
merged 2 commits into from
Nov 14, 2016
Merged

Reduce computational patch size for intermediate domain for vertical nesting #37

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
9a355e0
educe computational patch size for intermediate domain for vertical n…
davegill Nov 8, 2016
e3faea2
Fix my oops on indexing, squash this
davegill Nov 8, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
22 changes: 11 additions & 11 deletions external/RSL_LITE/module_dm.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4217,17 +4217,17 @@ SUBROUTINE force_domain_em_part2 ( grid, ngrid, pgrid, config_flags &

! Uncouple the variables moist and t_2 that are used to calculate ph_2

DO j = jps,jpe
DO i = ips,ipe
DO k=kps,kpe-1
DO j = MAX(jds,jps),MIN(jde-1,jpe)
DO i = MAX(ids,ips),MIN(ide-1,ipe)
DO k=kds,kde-1
grid%t_2(i,k,j) = grid%t_2(i,k,j)/(grid%mub(i,j) + grid%mu_2(i,j))
moist(i,k,j,P_QV) = moist(i,k,j,P_QV)/(grid%mub(i,j) + grid%mu_2(i,j))
END DO
END DO
END DO

DO j = jps, jpe
DO i = ips,ipe
DO j = MAX(jds,jps),MIN(jde-1,jpe)
DO i = MAX(ids,ips),MIN(ide-1,ipe)

DO k = 1, kpe-1
grid%pb(i,k,j) = ngrid%znu(k)*grid%mub(i,j)+ngrid%p_top
Expand Down Expand Up @@ -4271,8 +4271,8 @@ SUBROUTINE force_domain_em_part2 ( grid, ngrid, pgrid, config_flags &
ALLOCATE( p (ips:ipe, kps:kpe, jps:jpe) )
ALLOCATE( al(ips:ipe, kps:kpe, jps:jpe) )

DO j = jps, jpe
DO i = ips, ipe
DO j = MAX(jds,jps),MIN(jde-1,jpe)
DO i = MAX(ids,ips),MIN(ide-1,ipe)

! Integrate the hydrostatic equation (from the RHS of the bigstep vertical momentum
! equation) down from the top to get the pressure perturbation. First get the pressure
Expand Down Expand Up @@ -4339,15 +4339,15 @@ SUBROUTINE force_domain_em_part2 ( grid, ngrid, pgrid, config_flags &
DEALLOCATE(al)

! Couple the variables moist and t_2, and the newly calculated ph_2
DO j = jps, jpe
DO i = ips,ipe
DO j = MAX(jds,jps),MIN(jde-1,jpe)
DO i = MAX(ids,ips),MIN(ide-1,ipe)
DO k=kps,kpe
grid%ph_2(i,k,j) = grid%ph_2(i,k,j)*(grid%mub(i,j) + grid%mu_2(i,j))
END DO
END DO
END DO
DO j = jps, jpe
DO i = ips,ipe
DO j = MAX(jds,jps),MIN(jde-1,jpe)
DO i = MAX(ids,ips),MIN(ide-1,ipe)
DO k=kps,kpe-1
grid%t_2(i,k,j) = grid%t_2(i,k,j)*(grid%mub(i,j) + grid%mu_2(i,j))
moist(i,k,j,P_QV) = moist(i,k,j,P_QV)*(grid%mub(i,j) + grid%mu_2(i,j))
Expand Down