<p><b>laura@ucar.edu</b> 2011-02-09 09:50:49 -0700 (Wed, 09 Feb 2011)</p><p>changed array dimensions for long- and short-wave radiation fluxes when an additional layer is added above the model top<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_lw.F
===================================================================
--- branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_lw.F        2011-02-09 16:45:51 UTC (rev 728)
+++ branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_lw.F        2011-02-09 16:50:49 UTC (rev 729)
@@ -11413,7 +11413,8 @@
! Layer longwave fluxes (including extra layer above model top)
! Vertical ordering is from bottom to top (W m-2)
- REAL, DIMENSION( ims:ime, kms:kme+2, jms:jme ), &
+! REAL, DIMENSION( ims:ime, kms:kme+2, jms:jme ), &
+ REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), &
OPTIONAL, INTENT(OUT) :: &
LWUPFLX,LWUPFLXC,LWDNFLX,LWDNFLXC
Modified: branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_sw.F
===================================================================
--- branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_sw.F        2011-02-09 16:45:51 UTC (rev 728)
+++ branches/atmos_physics/src/core_physics/physics_wrf/module_ra_rrtmg_sw.F        2011-02-09 16:50:49 UTC (rev 729)
@@ -9663,7 +9663,8 @@
! Layer shortwave fluxes (including extra layer above model top)
! Vertical ordering is from bottom to top (W m-2)
- REAL, DIMENSION( ims:ime, kms:kme+2, jms:jme ), &
+! REAL, DIMENSION( ims:ime, kms:kme+2, jms:jme ), &
+ REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), &
OPTIONAL, INTENT(OUT) :: &
SWUPFLX,SWUPFLXC,SWDNFLX,SWDNFLXC
@@ -9822,6 +9823,9 @@
!------------------------------------------------------------------
+ write(0,*)
+ write(0,*) '--- enter subroutine rrtmg_swrad:'
+
!-----CALCULATE SHORT WAVE RADIATION
!
! All fields are ordered vertically from bottom to top
@@ -9849,7 +9853,6 @@
! clat(i) = xxlat
coszrs = sin(xxlat) * sin(declin) + cos(xxlat) * cos(declin) * cos(hrang)
coszr(i,j) = coszrs
-
! Set flag to prevent shortwave calculation when sun below horizon
if (coszrs.le.0.0) dorrsw = .false.
! Perform shortwave calculation if sun above horizon
</font>
</pre>