[Dart-dev] [4405] DART/trunk/observations/MADIS/convert_madis_acars.f90: change test for rejected altitude obs to be if either equal

[nancy at ucar.edu]

Revision: 4405
Author:   nancy
Date:     2010-07-01 07:50:56 -0600 (Thu, 01 Jul 2010)
Log Message:
-----------
change test for rejected altitude obs to be if either equal
to missing value OR bad qc; the old test required both.
thanks to ryan torn for the fix.

Modified Paths:
--------------
    DART/trunk/observations/MADIS/convert_madis_acars.f90

-------------- next part --------------
Modified: DART/trunk/observations/MADIS/convert_madis_acars.f90
===================================================================
--- DART/trunk/observations/MADIS/convert_madis_acars.f90	2010-07-01 04:16:05 UTC (rev 4404)
+++ DART/trunk/observations/MADIS/convert_madis_acars.f90	2010-07-01 13:50:56 UTC (rev 4405)
@@ -37,7 +37,7 @@
 use      location_mod, only : VERTISPRESSURE
 use  obs_sequence_mod, only : obs_sequence_type, obs_type, read_obs_seq, &
                               static_init_obs_sequence, init_obs, write_obs_seq, & 
-                              init_obs_sequence, get_num_obs, & 
+                              init_obs_sequence, get_num_obs, &
                               set_copy_meta_data, set_qc_meta_data
 use        meteor_mod, only : pres_alt_to_pres, sat_vapor_pressure, &
                               specific_humidity, wind_dirspd_to_uv
@@ -195,13 +195,13 @@
 
   ! Check for duplicate observations
   do i = 1, nused
-    if ( lon(n) == lonu(i) .and. &
-         lat(n) == latu(i) .and. &
+    if ( lon(n)  == lonu(i) .and. &
+         lat(n)  == latu(i) .and. &
          tobs(n) == tobu(i) .and. & 
          palt(n) == palu(i) ) cycle obsloop
   end do
 
-  if ( palt(n) == palt_miss .and. qc_palt(n) == 0 ) cycle obsloop  
+  if ( palt(n) == palt_miss .or. qc_palt(n) == 0 ) cycle obsloop  
   pres = pres_alt_to_pres(palt(n))
 
   ! extract actual time of observation in file into oday, osec.