<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>