[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20040504: Using sfcfil to create surface file from text data
- Subject: 20040504: Using sfcfil to create surface file from text data
- Date: 04 May 2004 16:13:33 -0600
Charley,
Your PARM line in your ascii file (and your SFPRMF input list) does
not need STID.
Your sfedit is looking for 4 parameters, and having to extend into the
next line to find it.
Steve Chiswell
Unidata User Support
On Tue, 2004-05-04 at 15:18, Charley Knoderer wrote:
> I have been trying to import air quality data into a surface gempak
> file using the following script:
>
> #!/bin/csh
> #
> #Get air quality data from ftp and creates airq gem file
> #
> mkdir /tmp/airq
> cd /tmp/airq
> #
> #lftp -u shcoperator,Navi7r webserver3.sonomatech.com << _FTP_EOF_
> #cd /outgoing/forecasts
> #get PMFN.dat
> #_FTP_EOF_
> #
> # Temporary copy of AirQ.dat to test file format.
> cp /home/gempak/Charley/AirQ.dat /tmp/airq/AirQ.dat
> #
> #
> #Define Date
> set YY=`date -u +%y`
> set MM=`date -u +%m`
> set DD=`date -u +%d`
> set DATE=$YY$MM$DD
> #set GEMFILE=/home/gempak/metdat/gempak/airq/${DATE}_pmfn.gem
> set GEMFILE=/home/gempak/metdat/gempak/airq/${DATE}_airq.gem
> #set TMPGEMFILE=${DATE}_pmfn.gem
> set TMPGEMFILE=${DATE}_airq.gem
> #
> sfcfil << SFCFIL_EOF
>
> SFOUTF = $TMPGEMFILE
> SFPRMF = PMFN;PM10;OZON;STID
> STNFIL = /home/gempak/nawips5.6j/gempak/tables/stns/airqstations.tbl
> SHIPFL = NO
> TIMSTN = 24/800
> r
>
> exit
>
> SFCFIL_EOF
> gpend
> sfedit << SFEDIT_EOF
>
> #SFEFIL = /tmp/airq/PMFN.dat
> SFEFIL = /tmp/airq/AirQ.dat
> SFFILE = $TMPGEMFILE
> r
>
> exit
>
> SFEDIT_EOF
> gpend
> mv ${TMPGEMFILE} ${GEMFILE}
> #rm /tmp/airq/*
> mv /tmp/airq/* /home/gempak/Charley/
> rmdir /tmp/airq
>
>
>
> The data file looks like this:
>
> PARM = PMFN;PM10;OZON;STID
>
> STID DATTIM PMFN PM10 OZON
> A0000981 040504/0000 131 159 134
> A0000981 040504/0100 115 19 85
> A0000981 040504/0200 179 88 148
> A0000981 040504/0300 62 192 66
> A0000981 040504/0400 77 2 65
> A0000981 040504/0500 6 58 185
> A0000981 040504/0600 53 31 176
> A0000981 040504/0700 0 60 134
> A0000981 040504/0800 135 195 138
> A0000981 040504/0900 142 53 126
> A0000981 040504/1000 37 166 163
> A0000981 040504/1100 173 156 9
> A0000981 040504/1200 128 140 99
> A0000981 040504/1300 74 167 125
> A0000981 040504/1400 111 168 166
> A0000981 040504/1500 22 45 49
> A0000981 040504/1600 170 111 139
> A0000981 040504/1700 21 17 15
> A0000981 040504/1800 143 103 154
> A0000981 040504/1900 68 14 47
> A0000981 040504/2000 37 114 12
> A0000985 040504/0000 191 6 57
> A0000985 040504/0100 161 144 122
> A0000985 040504/0200 164 70 99
> A0000985 040504/0300 187 110 13
> A0000985 040504/0400 64 152 59
> A0000985 040504/0500 19 187 39
> A0000985 040504/0600 163 110 90
> A0000985 040504/0700 15 100 52
> A0000985 040504/0800 138 1 149
> A0000985 040504/0900 1 149 42
> A0000985 040504/1000 40 116 105
> A0000985 040504/1100 127 49 109
> A0000985 040504/1200 67 96 3
> A0000985 040504/1300 106 197 81
> A0000985 040504/1400 155 84 44
> A0000985 040504/1500 199 52 115
> A0000985 040504/1600 44 174 106
> A0000985 040504/1700 180 192 48
> A0000985 040504/1800 197 118 153
> A0000985 040504/1900 122 34 21
> A0000995 040504/0000 66 168 71
> A0000995 040504/0100 8 101 24
> A0000995 040504/0200 123 187 53
> A0000995 040504/0300 2 183 185
> A0000995 040504/0400 33 51 78
> A0000995 040504/0500 87 190 194
> A0000995 040504/0600 42 73 148
> A0000995 040504/0700 4 64 4
> A0000995 040504/0800 175 90 193
> A0000995 040504/0900 189 45 68
> A0000995 040504/1000 58 11 111
> A0000995 040504/1100 124 134 35
> A0000995 040504/1200 11 133 196
> A0000995 040504/1300 178 140 38
> A0000995 040504/1400 158 198 145
> A0000995 040504/1500 146 73 91
> A0000995 040504/1600 53 117 3
> A0000995 040504/1700 143 97 26
> A0000995 040504/1800 186 96 10
> A0000995 040504/1900 36 47 115
> A0000995 040504/2000 95 131 146
> A0000999 040504/0000 107 21 12
> A0000999 040504/0100 126 55 74
> A0000999 040504/0200 195 56 163
> A0000999 040504/0300 58 55 109
> A0000999 040504/0400 145 133 129
> A0000999 040504/0500 61 68 118
> A0000999 040504/0600 82 60 26
> A0000999 040504/0700 6 197 79
> A0000999 040504/0800 2 149 121
> A0000999 040504/0900 163 97 119
> A0000999 040504/1000 52 178 137
> A0000999 040504/1100 143 117 193
> A0000999 040504/1200 62 129 0
> A0000999 040504/1300 133 70 132
> A0000999 040504/1400 46 93 127
> A0000999 040504/1500 166 81 17
> A0000999 040504/1600 76 136 7
> A0000999 040504/1700 193 103 86
> A0000999 040504/1800 17 33 194
> A0000999 040504/1900 174 150 102
> A0000999 040504/2000 172 160 220
>
> The problem I am having is that the surface gem file being created
> does not contain all of the data. It seems to be skipping every other
> hour. For example, two stations will have all data for the 1400 hour
> and then the other two stations will have data for 1500 hours, then it
> switches back the two stations that had data at 1400 hours for the
> 1600 hour? Any ideas as to why this is happening?
>
> -Charley
>
> Charley Knoderer
> Meteorologist 1
> Sonoma Technology, Inc.
> 1360 Redwood Way, Suite C
> Petaluma, CA 94954
> (707) 665-9900
> FAX: (707) 665-9800
> www.sonomatech.com
>