<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 < 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)
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, &
- 0.5*(pressure(k,iCell)+pressure(k+1,iCell)), &
- theta(k,iCell),temperature(k,iCell), &
- 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), &
+ 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), &
-! 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>