From 1fc671ccf61d9de08ec01761d0d0c2b43bc9a18c Mon Sep 17 00:00:00 2001 From: Jim Bresch Date: Thu, 8 Sep 2011 05:30:53 +0000 Subject: [PATCH] Back out the ruc changes in revision 603 (7/11/11). These changes will be added in a future release. git-svn-id: https://svn-wrf-wps.cgd.ucar.edu/trunk@612 86b71a92-4018-0410-97f8-d555beccfc3a --- ungrib/src/rrpr.F | 77 ------------------------------------------------------- 1 file changed, 77 deletions(-) diff --git a/ungrib/src/rrpr.F b/ungrib/src/rrpr.F index bc81865..9fbf453 100644 --- a/ungrib/src/rrpr.F +++ b/ungrib/src/rrpr.F @@ -187,7 +187,6 @@ subroutine rrpr(hstart, ntimes, interval, nlvl, maxlvl, plvl, debug_level, out_f call refw_storage(nint(level), field, ptr2d, map%nx, map%ny) nullify (ptr2d) enddo rdloop - ! ! We have reached the end of file, so time to close it. ! @@ -460,21 +459,6 @@ subroutine rrpr(hstart, ntimes, interval, nlvl, maxlvl, plvl, debug_level, out_f call put_storage(200100, 'SNOW', scr2d, map%nx, map%ny) deallocate(scr2d) endif - -! Add residual soil moisture to SOILM* if initialized from RUC model - if (index(map%source,'NOAA GSD') .ne. 0) then - call mprintf(.true.,DEBUG, & - "RRPR: Adjusting RUC soil moisture ") - if ( is_there(200100,'SOILM000') .and. & - is_there(200100,'SOILM005') .and. & - is_there(200100,'SOILM020') .and. & - is_there(200100,'SOILM040') .and. & - is_there(200100,'SOILM160') .and. & - is_there(200100,'SOILM300') ) then - call fix_ruc_soilm (map%nx, map%ny) - endif - endif - ! ! Check to see if we need to fill SOILHGT from SOILGEO. ! (From Wei Wang, 2007 June 21) @@ -901,64 +885,3 @@ subroutine fix_gfs_rh (ix, jx, plvl) deallocate (rh) deallocate (tt) end subroutine fix_gfs_rh - - -subroutine fix_ruc_soilm (ix, jx) -! This routine adds residual soil moisture if initialized fron RUC - use storage_module - implicit none - integer :: ix, jx, i, j - REAL , DIMENSION(100) :: lqmi - real, allocatable, dimension(:,:) :: soilm000, soilm005, soilm020, & - soilm040, soilm160, soilm300,soilcat - allocate(soilm000(ix,jx)) - allocate(soilm005(ix,jx)) - allocate(soilm020(ix,jx)) - allocate(soilm040(ix,jx)) - allocate(soilm160(ix,jx)) - allocate(soilm300(ix,jx)) - allocate(soilcat(ix,jx)) - call get_storage(200100, 'SOILM000', soilm000, ix, jx) - call get_storage(200100, 'SOILM005', soilm005, ix, jx) - call get_storage(200100, 'SOILM020', soilm020, ix, jx) - call get_storage(200100, 'SOILM040', soilm040, ix, jx) - call get_storage(200100, 'SOILM160', soilm160, ix, jx) - call get_storage(200100, 'SOILM300', soilm300, ix, jx) - - call get_storage(200100, 'SOILCAT', soilcat, ix, jx) - - lqmi(1:16) = & - (/0.045, 0.057, 0.065, 0.067, 0.034, 0.078, 0.10, & - 0.089, 0.095, 0.10, 0.070, 0.068, 0.078, 0.0, & - 0.004, 0.065 /) - - do j = 1, jx - do i = 1, ix - - SOILM000(i,j)=SOILM000(i,j) + lqmi(nint(soilcat(i,j))) - SOILM005(i,j)=SOILM005(i,j) + lqmi(nint(soilcat(i,j))) - SOILM020(i,j)=SOILM020(i,j) + lqmi(nint(soilcat(i,j))) - SOILM040(i,j)=SOILM040(i,j) + lqmi(nint(soilcat(i,j))) - SOILM160(i,j)=SOILM160(i,j) + lqmi(nint(soilcat(i,j))) - SOILM300(i,j)=SOILM300(i,j) + lqmi(nint(soilcat(i,j))) - enddo - enddo - call put_storage(200100, 'SOILM000', soilm000, ix, jx) - call put_storage(200100, 'SOILM005', soilm005, ix, jx) - call put_storage(200100, 'SOILM020', soilm020, ix, jx) - call put_storage(200100, 'SOILM040', soilm040, ix, jx) - call put_storage(200100, 'SOILM160', soilm160, ix, jx) - call put_storage(200100, 'SOILM300', soilm300, ix, jx) - - print *,'fix_ruc_soilm is done!' - - deallocate(soilm000) - deallocate(soilm005) - deallocate(soilm020) - deallocate(soilm040) - deallocate(soilm160) - deallocate(soilm300) - deallocate(soilcat) - -end subroutine fix_ruc_soilm - -- 2.11.4.GIT