<p><b>laura@ucar.edu</b> 2010-06-14 10:09:34 -0600 (Mon, 14 Jun 2010)</p><p>Increased maximum relative humidity to 70% to initiate convection<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-06-14 16:08:13 UTC (rev 345)
+++ branches/atmos_physics/src/core_hyd_atmos/module_test_cases.F        2010-06-14 16:09:34 UTC (rev 346)
@@ -70,8 +70,11 @@
       real (kind=RKIND), parameter :: u_perturbation = 1., pert_radius = 0.1, latitude_pert = 40., longitude_pert = 20.
       real (kind=RKIND), parameter :: theta_c = pii/4.0
       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 :: 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
@@ -372,7 +375,8 @@
          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)
@@ -380,11 +384,20 @@
                rel_hum(k,iCell) = min(rh_max,rel_hum(k,iCell))
             end do
          end do
+
+!        iCell = 5121
+!        write(6,*) '--- Relative humidity:'
+!        201 FORMAT(2i4,5(1x,e15.8))
+!        DO k = 1, nz1
+!           write(6,201) iCell,k,0.5*(pressure(k,iCell)+pressure(k+1,iCell)), &amp;
+!                        rel_hum(k,iCell)
+!        ENDDO
+!        write(6,*)
+
       else
          rel_hum(:,:) = 0.
       end if
 
-
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
       ! iteration 
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

</font>
</pre>