<p><b>laura@ucar.edu</b> 2010-07-26 10:11:27 -0600 (Mon, 26 Jul 2010)</p><p>cleaned up and updated with trunk<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_hyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_hyd_atmos/module_test_cases.F        2010-07-26 16:06:18 UTC (rev 441)
+++ branches/atmos_physics/src/core_hyd_atmos/module_test_cases.F        2010-07-26 16:11:27 UTC (rev 442)
@@ -72,9 +72,6 @@
       real (kind=RKIND), parameter :: lambda_c = 3.0*pii/2.0
       real (kind=RKIND), parameter :: rh_max = 0.4       ! Maximum relative humidity
       real (kind=RKIND), parameter :: k_x = 9.           ! Normal mode wave number
-!LDF begin (05-24-2010):
-!     real (kind=RKIND), parameter :: rh_max = 0.7       ! Maximum relative humidity
-!LDF end.
 
       real (kind=RKIND), dimension(:), pointer :: rdnu, rdnw, fnm, fnp, dbn, dnu, dnw
       real (kind=RKIND), dimension(:), pointer :: surface_pressure
@@ -370,29 +367,24 @@
 !
 !     initialization for moisture 
 !
-      201 format(i4,8(1x,e15.8))
       if (config_mp_physics /= 0) then
 
          do iCell=1,grid % nCells
             do k=1,nz1
                ptmp = 0.5*(pressure(k,iCell) + pressure(k+1,iCell))
-!              if (ptmp &lt; 50000.) then
-               if (ptmp &lt;= 50000.) then
+               if (ptmp &lt; 50000.) then
                   rel_hum(k,iCell) = 0.0
                else
                   rel_hum(k,iCell) = (1.-((p0-ptmp)/50000.)**1.25)
                end if
                rel_hum(k,iCell) = min(rh_max,rel_hum(k,iCell))
-!              if(iCell == 5121) write(6,201) k,ptmp,rel_hum(k,iCell)
             end do
          end do
-
       else
-
          rel_hum(:,:) = 0.
-
       end if
 
+
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
       ! iteration 
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
@@ -431,7 +423,7 @@
                    es  = 1000.*0.6112*exp(21.8745584*(temperature(k,iCell)-273.16)/(temperature(k,iCell)-7.66))
                end if
                qvs = (287.04/461.6)*es/(ptmp-es)
-!              qvs =  380.*exp(17.27*(temperature(k,iCell)-273.)/(temperature(k,iCell)-36.))/ptmp
+!               qvs =  380.*exp(17.27*(temperature(k,iCell)-273.)/(temperature(k,iCell)-36.))/ptmp
    
                scalars(index_qv,k,iCell) = rel_hum(k,iCell)*qvs
             end do
@@ -444,24 +436,14 @@
          end do
       end do
 
-      write(6,*)
-      write(6,*) '--- end sounding:'
-      iCell = 5121
-      do k = nz1,1,-1
-         write(6,201) k,0.5*(geopotential(k,iCell)+geopotential(k+1,iCell))/gravity, &amp;
-                      0.5*(pressure(k,iCell)+pressure(k+1,iCell)),                   &amp;
-                      theta(k,iCell),temperature(k,iCell),                           &amp;
-                      rel_hum(k,iCell),scalars(index_qv,k,iCell)
-      enddo
+      write(6,*) 'ptop = ',ptop,'  zt = ',geopotential(nz,1)/gravity
 
-!     write(6,*) 'ptop = ',ptop,'  zt = ',geopotential(nz,1)/gravity
+      write(6,*) ' full sounding with moisture'
+      do k=1,nz1
+         write(6,*) k, geopotential(k,1)/gravity, 0.01*pressure(k,1), theta(k,1), &amp;
+                    theta(k,1)*(pressure(k,1)/p0)**(rgas/cp)
+      end do
 
-!     write(6,*) ' full sounding with moisture'
-!     do k=1,nz1
-!        write(6,*) k, geopotential(k,1)/gravity, 0.01*pressure(k,1), theta(k,1), &amp;
-!                   theta(k,1)*(pressure(k,1)/p0)**(rgas/cp)
-!     end do
-
 ! When initializing a scalar, be sure not to put unreasonably large values
 ! into indices in the moist class
 !      scalars(2,:,:) = 1.  ! transport test

</font>
</pre>