From da9077d09a6e294d5b383710c6c5fdddecff46e6 Mon Sep 17 00:00:00 2001 From: "Jun.Wang" Date: Mon, 17 Feb 2020 21:36:06 +0000 Subject: [PATCH] add reproduce in g_sum call --- driver/fvGFS/atmosphere.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/driver/fvGFS/atmosphere.F90 b/driver/fvGFS/atmosphere.F90 index 39ba8c4c8..0fc5ba08e 100644 --- a/driver/fvGFS/atmosphere.F90 +++ b/driver/fvGFS/atmosphere.F90 @@ -1412,10 +1412,10 @@ subroutine atmosphere_state_update (Time, IPD_Data, IAU_Data, Atm_block, flip_vc if (IAU_Data%drymassfixer) then ! global mean total pressure and water before IAU psumb = g_sum(Atm(n)%domain,sum(Atm(n)%delp(isc:iec,jsc:jec,1:npz),dim=3),& - isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1) + isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1,reproduce=.true.) qsumb = g_sum(Atm(n)%domain,& sum(Atm(n)%delp(isc:iec,jsc:jec,1:npz)*sum(Atm(n)%q(isc:iec,jsc:jec,1:npz,1:nwat),4),dim=3),& - isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1) + isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1,reproduce=.true.) if (is_master()) then print *,'dry ps before IAU/physics',psumb+Atm(n)%ptop-qsumb endif @@ -1549,11 +1549,11 @@ subroutine atmosphere_state_update (Time, IPD_Data, IAU_Data, Atm_block, flip_vc if (IAU_Data%in_interval .and. IAU_data%drymassfixer) then ! global mean total pressure psum = g_sum(Atm(n)%domain,sum(Atm(n)%delp(isc:iec,jsc:jec,1:npz),dim=3),& - isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1) + isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1,reproduce=.true.) ! global mean total water (before adjustment) qsum = g_sum(Atm(n)%domain,& sum(Atm(n)%delp(isc:iec,jsc:jec,1:npz)*sum(Atm(n)%q(isc:iec,jsc:jec,1:npz,1:nwat),4),dim=3),& - isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1) + isc,iec,jsc,jec,Atm(n)%ng,Atm(n)%gridstruct%area_64,1,reproduce=.true.) betad = (psum - (psumb - qsumb))/qsum if (is_master()) then print *,'dry ps after IAU/physics',psum+Atm(n)%ptop-qsum