<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 < 50000.) then
+! if (ptmp < 50000.) then
+ if (ptmp <= 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)), &
+! rel_hum(k,iCell)
+! ENDDO
+! write(6,*)
+
else
rel_hum(:,:) = 0.
end if
-
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! iteration
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
</font>
</pre>