<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)), &amp;
-!                        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, &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,*) ' 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,*) '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
+
 ! 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>