[Dart-dev] [4773] DART/trunk/observations/NCEP/prep_bufr/src: Add missing override for land moisture errors. Obs whic

[nancy at ucar.edu]

Revision: 4773
Author:   nancy
Date:     2011-03-02 16:53:26 -0700 (Wed, 02 Mar 2011)
Log Message:
-----------
Add missing override for land moisture errors.  Obs which NCEP
does not assimilate come in with errors 1e10 so they will have
no effect; for surface temperature & winds we have been resetting
the obs error based on a namelist item so that the obs can be
used; the moisture namelist override was accidently omitted.
This has been missing since nov 2007.

Modified Paths:
--------------
    DART/trunk/observations/NCEP/prep_bufr/src/prepbufr.f
    DART/trunk/observations/NCEP/prep_bufr/src/prepbufr_03Z.f

-------------- next part --------------
Modified: DART/trunk/observations/NCEP/prep_bufr/src/prepbufr.f
===================================================================
--- DART/trunk/observations/NCEP/prep_bufr/src/prepbufr.f	2011-03-02 20:51:48 UTC (rev 4772)
+++ DART/trunk/observations/NCEP/prep_bufr/src/prepbufr.f	2011-03-02 23:53:26 UTC (rev 4773)
@@ -563,10 +563,10 @@
      &        use_this_data_int(pqm,qctype_use,inum_qctype) .and.
      &        qob .lt. 1.0e9                                ) then
 
-           qoe_rh = qoe
-
            es   = es0 * (tob / t0) ** fact1 * exp(fact2*(fact3-1./tob))
            qsat = eps * es / (pob - omeps * es)
+           if (qoe .gt. 1.e9) qoe = land_moist_error
+           qoe_rh = qoe
            qoe  = max(0.1, qoe * qsat * 1000.0) ! to g/kg, set min value
 
            if( .not. use_this_data_int(tqm,qctype_use,inum_qctype)) then

Modified: DART/trunk/observations/NCEP/prep_bufr/src/prepbufr_03Z.f
===================================================================
--- DART/trunk/observations/NCEP/prep_bufr/src/prepbufr_03Z.f	2011-03-02 20:51:48 UTC (rev 4772)
+++ DART/trunk/observations/NCEP/prep_bufr/src/prepbufr_03Z.f	2011-03-02 23:53:26 UTC (rev 4773)
@@ -582,10 +582,10 @@
      &        use_this_data_int(pqm,qctype_use,inum_qctype) .and.
      &        qob .lt. 1.0e9                                ) then
 
-           qoe_rh = qoe
-
            es   = es0 * (tob / t0) ** fact1 * exp(fact2*(fact3-1./tob))
            qsat = eps * es / (pob - omeps * es)
+           if (qoe .gt. 1.e9) qoe = land_moist_error
+           qoe_rh = qoe
            qoe  = max(0.1, qoe * qsat * 1000.0) ! to g/kg, set min value
 
            if( .not. use_this_data_int(tqm,qctype_use,inum_qctype)) then