20040928: Problems with sfmod : station not added

Christian,

Use blank characters for your STID (icao) instead of 9999
as you have done for stations that do not have an icao identifier.
Otherwise, you will be confusing the station id as the same 9999
STID for multiple stations, and therefore likely not finding the
one that matches your garea=@STNM. Also, your STNM will be the 5 digit
WMO for a synoptic bulletin.

Steve Chiswell
Unidata User Support


On Tue, 2004-09-28 at 07:48, Christian Pagé wrote:
> Hi everyone,
> 
> I have problems with SFMOD. Here is the situation:
> 1- I create a file with SFCFIL with stations I want using a station table.
> 2- I merge SYNOP data using SFMOD into the new GEMPAK file.
> 3- While I can list a particular station with SFLIST in the original
> SYNOP file, it is not merged in the new file with SFMOD, even though
> there is plenty of room for additional stations, and the station is in
> the table list. The only thing is that there is no ICAO for that
> station, only the WMO ID.
> 
> Here is the output of the commands and some embedded comments:
> 130 io [/io/httpd/meteoalerte_fr/sc/aa] % sfcfil
>  SFOUTF    Output surface file               aa.gem
>  SFPRMF    Surface parameter packing file    meteoalerte.pack
>  STNFIL    Station information file         
> /io/httpd/meteoalerte_fr/sc/france_metar.txt
>  SHIPFL    Ship data file flag               n
>  TIMSTN    Times/additional stations         70/9999
>  SFFSRC    Surface file source                
>  Parameters requested: SFOUTF,SFPRMF,STNFIL,SHIPFL,TIMSTN,SFFSRC.
>  GEMPAK-SFCFIL>r
> 
>  SFCFIL PARAMETERS: 
> 
>  New surface file:      aa.gem
>  Parameter file:        meteoalerte.pack
>  Station file:
>  Number of stations in STNFIL:     222
>  Number of additional stations:   9999
>  Total number of stations:       10221
>  Total number of times:             70
> 
>  This file will be a standard surface file.
> 
> Enter <cr> to accept parameters or type EXIT:
>  The file has been created.
>  Parameters requested: SFOUTF,SFPRMF,STNFIL,SHIPFL,TIMSTN,SFFSRC.
>  GEMPAK-SFCFIL>e
> 
> % wc -l /io/httpd/meteoalerte_fr/sc/france_metar.txt
>     222 /io/httpd/meteoalerte_fr/sc/france_metar.txt
> 
> % sfmod
>  SFFILE    Surface data file                
> /io/ldm/data/gempak/surface/syn/20040928_syn.gem
>  SFOUTF    Output surface file               aa.gem
>  DATTIM    Date/time                         all
>  DATOUT    Output date/time                   
>  AREA      Data area                         @74870
>  SFPARM    Surface parameter list            tmpc
>  Parameters requested: SFFILE,SFOUTF,DATTIM,DATOUT,AREA,SFPARM.
>  GEMPAK-SFMOD>r
> 
>  SFMOD PARAMETERS: 
> 
>  Input file:        /io/ldm/data/gempak/surface/syn/20040928_syn.gem
>  Output file:       aa.gem
>  Area:              @74870
>  Output parameters: tmpc
>  Input date/time:        Output date/time:  
>  040928/0000             040928/0000         
>  040928/0100             040928/0100         
>  040928/0200             040928/0200         
>  040928/0300             040928/0300         
>  040928/0400             040928/0400         
>  040928/0500             040928/0500         
>  040928/0600             040928/0600         
>  040928/0700             040928/0700         
>  040928/0800             040928/0800         
>  040928/0900             040928/0900         
>  040928/1000             040928/1000         
>  040928/1100             040928/1100         
>  040928/1200             040928/1200         
>  040928/1300             040928/1300         
> Enter <cr> to accept parameters or type EXIT:
>  Parameters requested: SFFILE,SFOUTF,DATTIM,DATOUT,AREA,SFPARM.
>  GEMPAK-SFMOD>e
> 
> % sflist
>  SFFILE    Surface data file                
> /io/ldm/data/gempak/surface/syn/20040928_syn.gem
>  AREA      Data area                         @74870
>  DATTIM    Date/time                         all
>  SFPARM    Surface parameter list            tmpc
>  OUTPUT    Output device/filename            t
>  IDNTYP    STNM or STID                      STID
>  Parameters requested: SFFILE,AREA,DATTIM,SFPARM,OUTPUT,IDNTYP.
>  GEMPAK-SFLIST>sffile=aa.gem
>  GEMPAK-SFLIST>area=dset
>  GEMPAK-SFLIST>r
>  PARM = TMPC
> 
>     STN    YYMMDD/HHMM      TMPC
> 
> 
>  [SFLIST -3]  No stations reporting data.
>  Parameters requested: SFFILE,AREA,DATTIM,SFPARM,OUTPUT,IDNTYP.
>  GEMPAK-SFLIST>e
> 
> Station is in the station list:
> 9999     074870 GRENOBLE/LE VERSOUD              -- FR  4522    585   220
> 
> Even if the area is FR or includes the latitude and longitude, this
> station is still not included.
> And I can list it in the original file:
> % sflist
>  SFFILE    Surface data file                
> /io/ldm/data/gempak/surface/syn/20040928_syn.gem
>  AREA      Data area                         @74870
>  DATTIM    Date/time                         all
>  SFPARM    Surface parameter list            tmpc
>  OUTPUT    Output device/filename            t
>  IDNTYP    STNM or STID                      STID
>  Parameters requested: SFFILE,AREA,DATTIM,SFPARM,OUTPUT,IDNTYP.
>  GEMPAK-SFLIST>r
>  PARM = TMPC
> 
>     STN    YYMMDD/HHMM      TMPC
>    9999    040928/0100      9.20
>    9999    040928/0200      8.80
>    9999    040928/0400      7.80
>    9999    040928/0500      8.10
>    9999    040928/0700      9.20
>    9999    040928/0800     12.70
>    9999    040928/1000     16.80
>    9999    040928/1100     18.70
>    9999    040928/1300     20.20
> 
> 
>  Parameters requested: SFFILE,AREA,DATTIM,SFPARM,OUTPUT,IDNTYP.
>  GEMPAK-SFLIST>e
> 
> Why?