GGPLOT2     HEALTH     R     EMISSIONS     STATISTICS     GIS   WRF Open Air Air Pollution  

martes, 25 de septiembre de 2012

Calpuff View: Fix(Read62) Delete spurious level

Algunas de las cosas que he aprendido las he ido publicando en mi modesto blog, pero ahora también publicaré mis dudas por que estoy intentando aprender Calpuff View, por si alguien me puede ayudar =)


Estoy en la parte de MET Procesor, con la meteorología superficial y de altura. No estoy bajando datos del re-analysis, sino de radio sondeos. Al procesar [Process Read62] me aparece una advertencia con el siguiente mensaje:
"One or more of the UP.DAT files created by READ62 filed de Wuality Assurance test. CALMENT will not be able to process these files. Please use Fix(READ62) to correct the issues"


Luego, al aplicar diferentes opciones de Fix(Read62) no llego a nada

Finalmente al leer el archivo Read62 List file dice:



                          READ62 OUTPUT SUMMARY
                   VERSION:  5.54         LEVEL:  070627      




 STARTING DATE:                ENDING DATE:

                YEAR = 2008                  YEAR = 2008
          JULIAN DAY =   32            JULIAN DAY =  335
                HOUR =    0 (GMT)            HOUR =   23 (GMT)


 PRESSURE LEVELS EXTRACTED:

                     SURFACE TO   700. MB


 INPUT FILE FORMAT (1=TD6201,2=NCDC CD-ROM):   2
 OUTPUT FILE FORMAT (1=/ DELIMITED,2=COMMA DELIMITED):   1


 ALT. SOUNDING FILE FOR SUBSTITUTIONS IS NOT USED 


 DATA LEVEL ELIMINATED IF HEIGHT MISSING ?         F

 DATA LEVEL ELIMINATED IF TEMPERATURE MISSING ?    F

 DATA LEVEL ELIMINATED IF WIND DIRECTION MISSING ? F

 DATA LEVEL ELIMINATED IF WIND SPEED MISSING ?     F


 MISSING PROFILE DATA EXTRAPOLATED TO TOP ?         T

 Last valid data must be above pressure (mb):      850.0

 MISSING PROFILE DATA EXTRAPOLATED TO SURFACE ?     T

 First valid data must be below height (m AGL):      200.0


 FILENAMES: 
     Control file:          read62.inp                                                                                                                          
     Input upper air file:  2STOD_~1.FSL                                                                                                                        
     Output upper air file: UP_99999.DAT                                                                                                                        
     Output list file:      UP_99999.LST                                                                                                                        

 -------------------------------------------------
 NOTICE: Starting year in control file sets the
         expected century for the simulation.  All
         YY years are converted to YYYY years in
         the range:  1958 2057
 -------------------------------------------------


 FSL Station ID used:         85586
             WBAN ID:       99999
              WMO ID:       85586


 Temperature values used in range checks: 
     TMIN = 175.0
     TMAX = 322.0

 Pressure values used in range checks: 
     PMIN =    0.0
     PMAX = 1040.0


 THE FOLLOWING SOUNDINGS HAVE BEEN PROCESSED:

       YEAR   MONTH   DAY   JULIAN DAY   HOUR (GMT)   NO. LEVELS EXTRACTED

       2008     2       1        32          0                 69
  ->->->Data at top of sounding is missing
  ->->->First level in sounding is above surface
  ->->->Elevation missing
  ->->->Invalid or missing pressure (< PMIN or > PMAX)
       2008     2       1        32         12                 79
  ->->->Data at top of sounding is missing
  ->->->First level in sounding is above surface
  ->->->Elevation missing
  ->->->Invalid or missing pressure (< PMIN or > PMAX)

Y el archivo Lakes List File sugiere:

 LAKES TOOL FOR REPAIRING UPPER AIR DATA
 R62CLEAN    VERSION:  5.5    03/31/2004

 ------
 USAGE:
 ------

    R62CLEAN5 [MODE] [UPDAT] [SUBDAT]
    MODE   (INTEGER) : Control mode
    UPDAT  (CHAR *)  : Upper air file to be repaired 
    SUBDAT (CHAR *)  : Substitution UP.DAT (only for MODE=3)

 -------------
 OUTPUT FILES:
 -------------

   * UPRn.DAT 
   * UPRn.LST 
   * n = MODE number

 ---------
 SETTINGS:
 ---------

    CONTROL MODE:     2
    -------------------------------------------
    DATA           IN/OUT     UNIT    FILE NAME
    -------------------------------------------
    UP.DAT           IN        11     UP_99999.DAT
    UPRn.DAT         OUT       21     UP_99999r2.dat
    UPRn.LST         OUT        7     UP_99999r2.lst
    -------------------------------------------

 ==============================================================================
 UPPER AIR DATA TO BE PROCESSED:


 Starting date:                Ending date:

                Year = 2008                  Year = 2008
          Julian day =   32            Julian day =  335
                Hour =    0 (GMT)            Hour =   23 (GMT)

 Topmost pressure level:  700. MB
 Input/Output file format (1=/ delimited,2=comma delimited):   1
 Upper air station ID       : 85586

 ------------------------------------------------------------------------------
 REPAIR PHASE (I) - REPAIR (P,Z) LEVELS
 ------------------------------------------------------------------------------

   * Delete spurious level which
       o contains unrealistic values of P and Z, or
       o is missing both T and wind vector (missing WD, or WS, or both)
   * Missing Z is derived from P according to hydrostatic equation.
   * Equal (anomalous) P or Z at successive levels are corrected according to 
     hydrostatic equation.



 THE FOLLOWING SOUNDINGS FROM <<UPDAT>> HAVE BEEN PROCESSED:

       YEAR   MONTH   DAY   JULIAN DAY   HOUR (GMT)   NO. LEVELS EXTRACTED

       2008     2      32         1          0                 69

      

Entonces debería borrar el nivel espurio (falso). ¿Colo lo podría hacer? .... Se aceptan sugerencias....