<p><b>laura@ucar.edu</b> 2010-07-23 15:36:08 -0600 (Fri, 23 Jul 2010)</p><p>Added print statements to calculation of initial relative humidity and water vapor<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-23 21:34:17 UTC (rev 423)
+++ branches/atmos_physics/src/core_hyd_atmos/module_test_cases.F        2010-07-23 21:36:08 UTC (rev 424)
@@ -70,10 +70,10 @@
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
+! 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
@@ -370,6 +370,7 @@
!
! initialization for moisture
!
+ 201 format(i4,8(1x,e15.8))
if (config_mp_physics /= 0) then
do iCell=1,grid % nCells
@@ -382,20 +383,14 @@
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
-! 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
- else
rel_hum(:,:) = 0.
+
end if
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -436,7 +431,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
@@ -449,14 +444,24 @@
end do
end do
- write(6,*) 'ptop = ',ptop,' zt = ',geopotential(nz,1)/gravity
+ 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,*) ' 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,*) '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
+
! 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>