[Dart-dev] [3443] DART/trunk/models/cam/model_mod.f90: In
model_interpolate, (re) move the KIND_PRESSURE test on
nancy at ucar.edu
nancy at ucar.edu
Tue Jul 8 12:43:09 MDT 2008
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080708/c776ece7/attachment.html
-------------- next part --------------
Modified: DART/trunk/models/cam/model_mod.f90
===================================================================
--- DART/trunk/models/cam/model_mod.f90 2008-07-02 13:48:02 UTC (rev 3442)
+++ DART/trunk/models/cam/model_mod.f90 2008-07-08 18:43:09 UTC (rev 3443)
@@ -2773,8 +2773,10 @@
vals(2,1) = phis(lon_ind_above, lat_ind_below) / gravity_const
vals(2,2) = phis(lon_ind_above, lat_ind_above) / gravity_const
-elseif (obs_type == KIND_PRESSURE) then
- ! Calculate pressures from surface pressures and A and B coeffs.
+! Move this to the end of the block? It's no good here; short circuits GPS
+! which asks for pressures on heights
+! elseif (obs_type == KIND_PRESSURE) then
+! ! Calculate pressures from surface pressures and A and B coeffs.
elseif (vert_is_level(location)) then
! Case 1: model level specified in vertical
@@ -2823,6 +2825,8 @@
if (vstatus /= 1) call get_val(vals(2,1),x, lon_ind_above, lat_ind_below, 1, obs_type, vstatus)
if (vstatus /= 1) call get_val(vals(2,2),x, lon_ind_above, lat_ind_above, 1, obs_type, vstatus)
+! Need option for vert_is_undefined
+
end if
! lat is already converted to degrees by get_location
More information about the Dart-dev
mailing list