[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