<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 < 0.5) then
lon = lon + 360.0
call latlon_to_ij(proj, lat, lon, x, y)
+ else if (x >= 360.5) then
+ lon = lon - 360.0
+ call latlon_to_ij(proj, lat, lon, x, y)
end if
+if (y < 1.0) y = 1.0
+if (y > 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 < 0.5) then
lon = lon + 360.0
call latlon_to_ij(proj, lat, lon, x, y)
+ else if (x >= 2500.5) then
+ lon = lon - 360.0
+ call latlon_to_ij(proj, lat, lon, x, y)
end if
+if (y < 1.0) y = 1.0
+if (y > 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 < 0.5) then
lon = lon + 360.0
call latlon_to_ij(proj, lat, lon, x, y)
+ else if (x >= 2500.5) then
+ lon = lon - 360.0
+ call latlon_to_ij(proj, lat, lon, x, y)
end if
+if (y < 1.0) y = 1.0
+if (y > 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>