[Dart-dev] [3300] DART/trunk/ncep_obs/prep_bufr/src: Minor format change in output; outputs location in degrees now instead of radians.

[nancy at subversion.ucar.edu]

An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080409/a22fa523/attachment.html
-------------- next part --------------
Modified: DART/trunk/ncep_obs/prep_bufr/src/prepbufr.f
===================================================================
--- DART/trunk/ncep_obs/prep_bufr/src/prepbufr.f	2008-04-09 21:44:13 UTC (rev 3299)
+++ DART/trunk/ncep_obs/prep_bufr/src/prepbufr.f	2008-04-09 22:04:02 UTC (rev 3300)
@@ -48,7 +48,8 @@
 
       dimension tdata(8), udata(8), vdata(8), qdata(8), pdata(8)
       integer :: wtype, ptype, qtype, ttype
-      integer :: pc_t, pc_q, pc_u, pc_v, pc_p
+c     what are the pc_x values?
+      integer :: pc_t, pc_q, pc_u, pc_v, pc_p 
       integer :: tqm, pqm, qqm, uqm, vqm, qctype_use(max_otype)
       logical :: found, uotype, uqcflag, use_this_data_real, 
      +            use_this_data_int
@@ -374,6 +375,10 @@
             tdata(4) = ppb
             tdata(5) = tob
 
+c    in some old files this appears to be out of range
+c    and it seems to be unused in converting to an obs_seq so
+c    i feel ok setting it to something that will fit in an I2 field.
+            if (pc_t > 99) pc_t = 99
             write(lunobs, 800) tdata, ttype, tqm, subset(1:6), pc_t
 
           endif
@@ -396,6 +401,7 @@
             tdata(4) = zob
             tdata(5) = tob
 
+            if (pc_t > 99) pc_t = 99
             write(lunobs, 800) tdata, ttype, tqm, subset(1:6), pc_t
 
           endif
@@ -429,6 +435,7 @@
             qdata(5) = qob
 
             if (qoe .lt. 9.9) then  ! skip large qoe obs.
+              if (pc_q > 99) pc_q = 99
               write(lunobs, 800) qdata, qtype, qqm, subset(1:6), pc_q
             endif
 
@@ -462,6 +469,7 @@
            qdata(5) = qob
 
            if(qoe .lt. 9.9) then   ! skip large qoe obs
+             if (pc_q > 99) pc_q = 99
              write(lunobs, 800) qdata, qtype, qqm, subset(1:6), pc_q
            endif
 
@@ -485,6 +493,7 @@
             pdata(4) = zob
             pdata(5) = pob
 
+            if (pc_p > 99) pc_p = 99
             write(lunobs, 800) pdata, ptype, pqm, subset(1:6), pc_p
 
           endif
@@ -511,6 +520,8 @@
             vdata(5) = vob
             vdata(6) = uob
 
+            if (pc_u > 99) pc_u = 99
+            if (pc_v > 99) pc_v = 99
             write(lunobs, 800) udata, wtype, uqm, subset(1:6), pc_u
             write(lunobs, 800) vdata, wtype, vqm, subset(1:6), pc_v
 
@@ -539,6 +550,8 @@
             vdata(5) = vob
             vdata(6) = uob
 
+            if (pc_u > 99) pc_u = 99
+            if (pc_v > 99) pc_v = 99
             write(lunobs, 800) udata, wtype, uqm, subset(1:6), pc_u
             write(lunobs, 800) vdata, wtype, vqm, subset(1:6), pc_v
 

Modified: DART/trunk/ncep_obs/prep_bufr/src/prepbufr_03Z.f
===================================================================
--- DART/trunk/ncep_obs/prep_bufr/src/prepbufr_03Z.f	2008-04-09 21:44:13 UTC (rev 3299)
+++ DART/trunk/ncep_obs/prep_bufr/src/prepbufr_03Z.f	2008-04-09 22:04:02 UTC (rev 3300)
@@ -374,6 +374,8 @@
             tdata(4) = ppb
             tdata(5) = tob
 
+c    see comment in prepbufr.f
+            if (pc_t > 99) pc_t = 99
             write(lunobs, 800) tdata, ttype, tqm, subset(1:6), pc_t
 
           endif
@@ -396,6 +398,7 @@
             tdata(4) = zob
             tdata(5) = tob
 
+            if (pc_t > 99) pc_t = 99
             write(lunobs, 800) tdata, ttype, tqm, subset(1:6), pc_t
 
           endif
@@ -429,6 +432,7 @@
             qdata(5) = qob
 
             if (qoe .lt. 9.9) then  ! skip large qoe obs.
+              if (pc_q > 99) pc_q = 99
               write(lunobs, 800) qdata, qtype, qqm, subset(1:6), pc_q
             endif
 
@@ -462,6 +466,7 @@
            qdata(5) = qob
 
            if(qoe .lt. 9.9) then   ! skip large qoe obs
+             if (pc_q > 99) pc_q = 99
              write(lunobs, 800) qdata, qtype, qqm, subset(1:6), pc_q
            endif
 
@@ -485,6 +490,7 @@
             pdata(4) = zob
             pdata(5) = pob
 
+            if (pc_p > 99) pc_p = 99
             write(lunobs, 800) pdata, ptype, pqm, subset(1:6), pc_p
 
           endif
@@ -511,6 +517,8 @@
             vdata(5) = vob
             vdata(6) = uob
 
+            if (pc_u > 99) pc_u = 99
+            if (pc_v > 99) pc_v = 99
             write(lunobs, 800) udata, wtype, uqm, subset(1:6), pc_u
             write(lunobs, 800) vdata, wtype, vqm, subset(1:6), pc_v
 
@@ -539,6 +547,8 @@
             vdata(5) = vob
             vdata(6) = uob
 
+            if (pc_u > 99) pc_u = 99
+            if (pc_v > 99) pc_v = 99
             write(lunobs, 800) udata, wtype, uqm, subset(1:6), pc_u
             write(lunobs, 800) vdata, wtype, vqm, subset(1:6), pc_v