<p><b>duda</b> 2011-08-12 13:42:44 -0600 (Fri, 12 Aug 2011)</p><p>BRANCH COMMIT<br>
<br>
Add fixes to lat and lon range for interpolation of SNOALB, GREENFRAC, <br>
and ALBEDO12M in real-data initialization.<br>
<br>
<br>
M    src/core_init_nhyd_atmos/module_test_cases.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2011-08-12 19:38:58 UTC (rev 937)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2011-08-12 19:42:44 UTC (rev 938)
@@ -2531,7 +2531,12 @@
             if (x &lt; 0.5) then
                lon = lon + 360.0
                call latlon_to_ij(proj, lat, lon, x, y)
+            else if (x &gt;= 360.5) then
+               lon = lon - 360.0
+               call latlon_to_ij(proj, lat, lon, x, y)
             end if
+if (y &lt; 1.0) y = 1.0
+if (y &gt; 179.0) y = 179.0
 !            grid % snoalb % array(iCell) = interp_sequence(x, y, 1, maxsnowalb, 1, 360, 1, 180, 1, 1, -1.e30, interp_list, 1)
             grid % snoalb % array(iCell) = interp_sequence(x, y, 1, maxsnowalb, 1, 360, 1, 180, 1, 1, 0., interp_list, 1)
          else
@@ -2597,7 +2602,12 @@
             if (x &lt; 0.5) then
                lon = lon + 360.0
                call latlon_to_ij(proj, lat, lon, x, y)
+            else if (x &gt;= 2500.5) then
+               lon = lon - 360.0
+               call latlon_to_ij(proj, lat, lon, x, y)
             end if
+if (y &lt; 1.0) y = 1.0
+if (y &gt; 1249.0) y = 1249.0
             do k=1,12
                grid % greenfrac % array(k,iCell) = interp_sequence(x, y, k, vegfra, 1, 2500, 1, 1250, 1, 12, -1.e30, interp_list, 1)
             end do
@@ -2669,7 +2679,12 @@
             if (x &lt; 0.5) then
                lon = lon + 360.0
                call latlon_to_ij(proj, lat, lon, x, y)
+            else if (x &gt;= 2500.5) then
+               lon = lon - 360.0
+               call latlon_to_ij(proj, lat, lon, x, y)
             end if
+if (y &lt; 1.0) y = 1.0
+if (y &gt; 1249.0) y = 1249.0
             do k=1,12
                grid % albedo12m % array(k,iCell) = interp_sequence(x, y, k, vegfra, 1, 2500, 1, 1250, 1, 12, 0.0, interp_list, 1)
             end do

</font>
</pre>