<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 (&quot;rrtmg_lw&quot;)
-       f_ice(:,:)  = 0.
-       f_rain(:,:) = 0.
        write(0,*) '--- enter subroutine rrtmg_lwrad:'
        call rrtmg_lwrad( &amp;
                 rthratenlw = rthratenlw_p , lwupt     = lwupt_p    , lwuptc    = lwuptc_p  , &amp;

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 (&quot;rrtmg_sw&quot;)
-       f_ice(:,:)  = 0.
-       f_rain(:,:) = 0.
-
        write(0,*) '--- enter subroutine rrtmg_swrad:'
        call rrtmg_swrad( &amp;
                 rthratensw = rthratensw_p , swupt     = swupt_p    , swuptc    = swuptc_p  , &amp;

</font>
</pre>