[Dart-dev] [4519] DART/trunk/models/tiegcm/model_mod.f90: Now obs_val is set to DART missing value (missing_r8) when there is a problem in computing the interpolation (istatus /= 0). In the l

[nancy at ucar.edu]

Revision: 4519
Author:   tmatsuo
Date:     2010-10-12 16:44:19 -0600 (Tue, 12 Oct 2010)
Log Message:
-----------
Now obs_val is set to DART missing value (missing_r8) when there is a problem in computing the interpolation (istatus /= 0).  In the last version obs_val was set to 0 when istatus /= 0. 

Modified Paths:
--------------
    DART/trunk/models/tiegcm/model_mod.f90

-------------- next part --------------
Modified: DART/trunk/models/tiegcm/model_mod.f90
===================================================================
--- DART/trunk/models/tiegcm/model_mod.f90	2010-10-12 21:01:21 UTC (rev 4518)
+++ DART/trunk/models/tiegcm/model_mod.f90	2010-10-12 22:44:19 UTC (rev 4519)
@@ -17,7 +17,7 @@
 !------------------------------------------------------------------
 
 ! DART Modules
-use        types_mod, only : r8, digits12
+use        types_mod, only : r8, digits12, missing_r8
 use time_manager_mod, only : time_type, set_calendar_type, set_time_missing, &
                              set_time, get_time, print_time, &
                              set_date, get_date, print_date, & 
@@ -383,13 +383,10 @@
    end do
    obs_val = lat_fract * a(2) + (1.0 - lat_fract) * a(1)
 else
-   obs_val = 0.
+   obs_val = missing_r8
 endif
 
-!write(11,*,access='APPEND') lon, lat, lon_fract, lat_fract, a(1), a(2), obs_val
-!print*, 'model_interpolate', lon, lat, height,obs_val
 
-
 end subroutine model_interpolate
 
 
@@ -460,7 +457,6 @@
 
 val = exp(frac_lev * log(val_bottom)  +  (1.0 - frac_lev) * log(val_top))
 
-!write(11,*,access='APPEND') lat_index, lon_index, lev_top, lev_bottom, frac_lev, val_top, val_bottom, val
 
 end subroutine get_val