<p><b>ringler@lanl.gov</b> 2011-12-23 15:32:14 -0700 (Fri, 23 Dec 2011)</p><p><br>
minor updates to pv solver<br>
</p><hr noshade><pre><font color="gray">Modified: branches/pv_based_swm/mpas/src/core_pvsw/module_poisson_solver.F
===================================================================
--- branches/pv_based_swm/mpas/src/core_pvsw/module_poisson_solver.F        2011-12-22 23:17:11 UTC (rev 1277)
+++ branches/pv_based_swm/mpas/src/core_pvsw/module_poisson_solver.F        2011-12-23 22:32:14 UTC (rev 1278)
@@ -8,7 +8,7 @@
implicit none
real (KIND=RKIND), parameter :: threshold = 1.0e-9
- integer, parameter :: nIterPerGS = 2, nCheckPoisson = 25, nQuitPoisson = 5000
+ integer, parameter :: nIterPerGS = 2, nCheckPoisson = 100, nQuitPoisson = 25000
private
public :: init_poisson, poisson, derive_velocity
Modified: branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F
===================================================================
--- branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F        2011-12-22 23:17:11 UTC (rev 1277)
+++ branches/pv_based_swm/mpas/src/core_pvsw/module_time_integration.F        2011-12-23 22:32:14 UTC (rev 1278)
@@ -981,7 +981,7 @@
pv_forcing(:,:) = 0.0
rtime = s % xtime % scalar
- amp = s % pvForcingAmp % scalar
+ amp = s % pvForcingAmp % scalar !! units of 1/s^2
write(6,*) 'compute forcing : ',rtime, amp
if(config_forcing) then
@@ -997,7 +997,7 @@
do iCell = 1,grid%nCells
do k=1,grid % nVertLevels
- pv_forcing(k,iCell) = amp * (cos(lwave*yCell(iCell)+phiy) - cos(kwave*xCell(iCell)+phix)) / h(k,iCell)
+ pv_forcing(k,iCell) = amp * (cos(lwave*yCell(iCell)+phiy) - cos(kwave*xCell(iCell)+phix))
enddo
enddo
</font>
</pre>