<p><b>laura@ucar.edu</b> 2012-06-15 13:46:05 -0600 (Fri, 15 Jun 2012)</p><p>now prints the max and min values of physics tendencies for rk_step=3 only<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_todynamics.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_todynamics.F        2012-06-15 19:00:14 UTC (rev 1990)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_todynamics.F        2012-06-15 19:46:05 UTC (rev 1991)
@@ -12,15 +12,16 @@
  contains
  
 !=============================================================================================
-subroutine physics_addtend(mesh, state, diag, tend, tend_physics, mass, mass_edge)
+subroutine physics_addtend(mesh, state, diag, tend, tend_physics, mass, mass_edge, rk_step)
 !=============================================================================================
 
 !input variables:
 !----------------
-type(mesh_type),intent(in):: mesh
+ type(mesh_type),intent(in):: mesh
  type(state_type),intent(in):: state
  type(diag_type),intent(in):: diag
  type(tend_physics_type),intent(inout):: tend_physics
+ integer, intent(in):: rk_step
  real(kind=RKIND),dimension(:,:),intent(in):: mass
  real(kind=RKIND),dimension(:,:),intent(in):: mass_edge
 
@@ -54,9 +55,6 @@
 !ldf end.
 
 !=============================================================================================
-!write(0,*)
-!write(0,*) '--- enter subroutine physics_add_tend:'
-
  block =&gt; mesh % block
  
  nCells      = mesh % nCells
@@ -120,13 +118,6 @@
     enddo
     enddo
  endif
-write(0,*) 'max rthblten = ',maxval(rthblten(:,1:nCellsSolve))
-write(0,*) 'min rthblten = ',minval(rthblten(:,1:nCellsSolve))
-!write(0,*) 'max rqvblten = ',maxval(rqvblten(:,1:nCellsSolve))
-!write(0,*) 'min rqvblten = ',minval(rqvblten(:,1:nCellsSolve))
-!write(0,*) 'max tend     = ',maxval(tend_scalars(tend%index_qv,:,1:nCellsSolve))
-!write(0,*) 'min tend     = ',minval(tend_scalars(tend%index_qv,:,1:nCellsSolve))
-!write(0,*)
 
 !add coupled tendencies due to convection:
  if(config_conv_deep_scheme .ne. 'off') then
@@ -154,13 +145,6 @@
     enddo
     enddo
  endif
-write(0,*) 'max rthcuten = ',maxval(rthcuten(:,1:nCellsSolve))
-write(0,*) 'min rthcuten = ',minval(rthcuten(:,1:nCellsSolve))
-!write(0,*) 'max rqvcuten = ',maxval(rqvcuten(:,1:nCellsSolve))
-!write(0,*) 'min rqvcuten = ',minval(rqvcuten(:,1:nCellsSolve))
-!write(0,*) 'max tend     = ',maxval(tend_scalars(tend%index_qv,:,1:nCellsSolve))
-!write(0,*) 'min tend     = ',minval(tend_scalars(tend%index_qv,:,1:nCellsSolve))
-!write(0,*)
 
 !add coupled tendencies due to longwave radiation:
  if(config_radt_lw_scheme .ne. 'off') then
@@ -170,8 +154,6 @@
     enddo
     enddo
  endif
-write(0,*) 'max rthratenlw = ',maxval(rthratenlw(:,1:nCellsSolve))
-write(0,*) 'min rthratenlw = ',minval(rthratenlw(:,1:nCellsSolve))
     
 !add coupled tendencies due to shortwave radiation:
  if(config_radt_sw_scheme .ne. 'off') then
@@ -181,8 +163,6 @@
     enddo
     enddo
  endif
-write(0,*) 'max rthratensw = ',maxval(rthratensw(:,1:nCellsSolve))
-write(0,*) 'min rthratensw = ',minval(rthratensw(:,1:nCellsSolve))
 
 !if non-hydrostatic core, convert the tendency for the potential temperature to a
 !tendency for the modified potential temperature:
@@ -205,6 +185,21 @@
  deallocate(theta)
  deallocate(tend_th)
 
+ if(rk_step .eq. 3) then
+    write(0,*)
+    write(0,*) '--- enter subroutine physics_addtend:'
+    write(0,*) 'max rthblten   = ',maxval(rthblten(:,1:nCellsSolve))
+    write(0,*) 'min rthblten   = ',minval(rthblten(:,1:nCellsSolve))
+    write(0,*) 'max rthcuten   = ',maxval(rthcuten(:,1:nCellsSolve))
+    write(0,*) 'min rthcuten   = ',minval(rthcuten(:,1:nCellsSolve))
+    write(0,*) 'max rthratenlw = ',maxval(rthratenlw(:,1:nCellsSolve))
+    write(0,*) 'min rthratenlw = ',minval(rthratenlw(:,1:nCellsSolve))
+    write(0,*) 'max rthratensw = ',maxval(rthratensw(:,1:nCellsSolve))
+    write(0,*) 'min rthratensw = ',minval(rthratensw(:,1:nCellsSolve))
+    write(0,*) '--- end subroutine physics_addtend'
+    write(0,*)
+ endif
+
 !formats:
  201 format(2i6,10(1x,e15.8))
  202 format(3i6,10(1x,e15.8))

</font>
</pre>