<p><b>laura@ucar.edu</b> 2011-10-17 14:41:26 -0600 (Mon, 17 Oct 2011)</p><p>corrected the dimensions of f_ice and f_rain<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F        2011-10-17 20:19:18 UTC (rev 1103)
+++ branches/atmos_physics/src/core_physics/module_driver_radiation_lw.F        2011-10-17 20:41:26 UTC (rev 1104)
@@ -29,8 +29,8 @@
subroutine allocate_radiation_lw
!=============================================================================================
- if(.not.allocated(f_ice) ) allocate(f_ice(ims:ime,jms:jme) )
- if(.not.allocated(f_rain) ) allocate(f_rain(ims:ime,jms:jme) )
+ if(.not.allocated(f_ice) ) allocate(f_ice(ims:ime,kms:kme,jms:jme) )
+ if(.not.allocated(f_rain) ) allocate(f_rain(ims:ime,kms:kme,jms:jme) )
if(.not.allocated(sfc_emiss_p) ) allocate(sfc_emiss_p(ims:ime,jms:jme) )
if(.not.allocated(snow_p) ) allocate(snow_p(ims:ime,jms:jme) )
@@ -210,11 +210,13 @@
!initialization:
do j = jts,jte
+ do k = kts,kte
do i = its,ite
- f_ice(i,j) = 0.
- f_rain(i,j) = 0.
+ f_ice(i,k,j) = 0.
+ f_rain(i,k,j) = 0.
enddo
enddo
+ enddo
do j = jts,jte
do i = its,ite
@@ -521,8 +523,6 @@
radiation_lw_select: select case (trim(radt_lw_scheme))
case ("rrtmg_lw")
- f_ice(:,:) = 0.
- f_rain(:,:) = 0.
write(0,*) '--- enter subroutine rrtmg_lwrad:'
call rrtmg_lwrad( &
rthratenlw = rthratenlw_p , lwupt = lwupt_p , lwuptc = lwuptc_p , &
Modified: branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F        2011-10-17 20:19:18 UTC (rev 1103)
+++ branches/atmos_physics/src/core_physics/module_driver_radiation_sw.F        2011-10-17 20:41:26 UTC (rev 1104)
@@ -30,8 +30,8 @@
subroutine allocate_radiation_sw
!=============================================================================================
- if(.not.allocated(f_ice) ) allocate(f_ice(ims:ime,jms:jme) )
- if(.not.allocated(f_rain) ) allocate(f_rain(ims:ime,jms:jme) )
+ if(.not.allocated(f_ice) ) allocate(f_ice(ims:ime,kms:kme,jms:jme) )
+ if(.not.allocated(f_rain) ) allocate(f_rain(ims:ime,kms:kme,jms:jme) )
if(.not.allocated(xlat_p) ) allocate(xlat_p(ims:ime,jms:jme) )
if(.not.allocated(xlon_p) ) allocate(xlon_p(ims:ime,jms:jme) )
@@ -218,8 +218,18 @@
enddo
enddo
+!initialization:
do j = jts,jte
+ do k = kts,kte
do i = its,ite
+ f_ice(i,k,j) = 0.
+ f_rain(i,k,j) = 0.
+ enddo
+ enddo
+ enddo
+
+ do j = jts,jte
+ do i = its,ite
coszr_p(i,j) = 0.
gsw_p(i,j) = 0.
swcf_p(i,j) = 0.
@@ -521,9 +531,6 @@
radiation_sw_select: select case (trim(radt_sw_scheme))
case ("rrtmg_sw")
- f_ice(:,:) = 0.
- f_rain(:,:) = 0.
-
write(0,*) '--- enter subroutine rrtmg_swrad:'
call rrtmg_swrad( &
rthratensw = rthratensw_p , swupt = swupt_p , swuptc = swuptc_p , &
</font>
</pre>