Hey Stephen: I don't know if you got this email
Daniel Fritz
prodicalboxer at YAHOO.COM
Sat Feb 9 17:30:35 EST 2008
Stephen,
I made those necessary corrections but the only
thing after I get my fw file, I typed up the control
file from my terminal on Florence2000, I opened the
ctl file and then ran track4 and it is saying that
htere is an error on line 104 and 170 and talks about
a non numeric arg for......here it is
Draw wind streamlines and plot Florence
track...
Non-numeric args to numeric operation
Error occurred on line 170
In file /home2/fritz/script/FLORENCEWHIRLPOOL.gs
Error occurred on line 104
In file /home2/fritz/script/FLORENCEWHIRLPOOL.gs
This does not make sense whats so ever because this
same script worked for floy99 and Gert99......
does the line (in the track4 script)have to change
where it says timeinfo 9, does it change to timeinfo
11, I changed it and it still game me the same
thing??? Is something wrong with me creating the
binary file????? I thought I was doing it right...It
seemed like everything was getting ready to run
fine....the script asked for the storm and
everything.......then when I type in Florence, it
gives me an error......
here is the script to be converted to fw
DSET c:\pcgrads\lib\dev\dat\wind.sep99.fw
TITLE SEP1999 WIND (SIG995) NCEP-NCAR REANALYSIS
UNDEF -9.96921e+36
XDEF 31 Linear 255 2.5
YDEF 21 Linear 0 2.5
ZDEF 1 LEVELS 1000
TDEF 60 LINEAR 00z01Sep1999 12hr
VARS 2
ugrd10m 0 -999 sig995 u-wind (m/s)
vgrd10m 0 -999 sig995 v-wind (m/s)
ENDVARS
geosci5.ncat.edu{fritz}1005: cat FWRITE2
# This sample writes out unwnd and vwnd variables to
part of the N Atlantic
# reads all timestep
datefile1='/home2/fritz/uwnd.sig995.2000.nc'
datefile2='/home2/fritz/vwnd.sig995.2000.nc'
outfile='/home2/fritz/wind.Florence2000.fw'
# open source files, set map environment
'sdfopen 'datefile1
'sdfopen 'datefile2
'set lat 0 50';'set lon 255 330'
time1='00z01Sep2000';time2='12z30Sep2000'
'set time 'time1' 'time2
'q dim'
timeinfo=sublin(result,5)
t1=subwrd(timeinfo,11);t2=subwrd(timeinfo,13)
say 'trange = 't1' 't2
'set gxout fwrite'
'set fwrite 'outfile
t=t1
while(t<=t2)
'set t 't
'd uwnd'
'd vwnd.2'
t=t+1
endwhile
'disable fwrite'
I know this is right by your standards.......I know
that this is a smal problem, because if this works for
two hurricanes, I know that this will work for the
other ones......but note: this is the first hurricane
that I have done fwrite statement and ctl .files for
binary files so I could be doing something wrong....
DANNY :):):):):)
--- smcmillan at planalytics.com
--- Stephen R McMillan <smcmillan at PLANALYTICS.COM>
wrote:
> Danny,
> See my detailed email from earlier this evening. It
> explains exactly how
> to get the dimension environment that you need for
> your control file. It
> was the one with two files attached, including the
> sample control file
> with step-by-step comments. Good luck!
> Stephen
>
>
>
> Daniel Fritz <prodicalboxer at YAHOO.COM>
> Sent by: GRADSUSR at LIST.CINECA.IT
> 02/07/2008 11:22 PM
> Please respond to
> GRADSUSR at LIST.CINECA.IT
>
>
> To
> GRADSUSR at LIST.CINECA.IT
> cc
>
> Subject
> Re: converting binary file fw to control file ctl
>
>
>
>
>
> Hey Eric,
> could you go alittle more into detail because i am
> not quite understanding...do you mean my lat and lon
> dimension environment.....how do i know what var is
> if
> I am extracting data then writting the script for
> the
> fwrite????
> DANNY
> --- Eric Altshuler <ela at COLA.IGES.ORG> wrote:
>
> > Danny,
> >
> > You need to know how the fwrite file was written,
> > i.e. what was the dimension environment when the
> > 'display VAR' command was given with gxout set to
> > 'fwrite', which variables were displayed, etc.
> There
> > is no way to create a ctl file for a binary
> dataset
> > without knowing how the dataset was written, since
> > it has no metadata. If you created the file in a
> > grads session or script, you need to figure out
> what
> > is actually in that file.
> >
> > Eric
> >
> > ----- Original Message -----
> > From: "Daniel Fritz" <prodicalboxer at YAHOO.COM>
> > To: GRADSUSR at LIST.CINECA.IT
> > Sent: Thursday, February 7, 2008 10:41:31 PM
> > (GMT-0500) America/New_York
> > Subject: converting binary file fw to control file
> > ctl
> >
> > Hello Everyone,
> > from using gxout fwrite in grads
> > there come a binary file example: filename.fx ,
> how
> > do you convert the fw binary file to a control
> > file????? I don't quite understand from reading
> the
> > Grads manual......
> > can anyone help with this???
> > It's 10:30pm here, I will be up for
> > awhile..please feel free to email...
> > DANNY
> >
> > --- Adam Sobel <ahs129 at COLUMBIA.EDU> wrote:
> >
> > > Thanks Arlindo, took a bit of redoing to get the
> > ctl
> > > file to work with
> > > "open" but it
> > > seems to have fixed the problem now!
> > >
> > > Adam
> > >
> > > Arlindo da Silva wrote:
> > > > On Feb 7, 2008 12:15 AM, Adam Sobel
> > > <ahs129 at columbia.edu
> > > > <mailto:ahs129 at columbia.edu>> wrote:
> > > >
> > > > I have a set of netcdf files, each is 1
> day.
> > > I have a .ctl file
> > > > so I can
> > > > view a bunch of them together in sequence
> > > using xdfopen. The ctl file
> > > > is below.
> > > >
> > > >
> > > >
> > > > If I set t to the whole time range (in
> this
> > > case set t 1 12) I can
> > > > view
> > > > the variables u or v just fine, e.g. make
> > > animations for the whole
> > > > period, looks good. But then if I do
> > > >
> > > > define rvort=hcurl(u,v)
> > > >
> > > > I get a segmentation fault
> > > >
> > > >
> > > > This is a well known problem in the v1.8/v1.9
> > > series (but not v1.7b9
> > > > and new v2). There is a memory leak in the
> > > sdfopen/xdfopen part of the
> > > > code that causes the amount of memory to
> > increase
> > > with the size of the
> > > > timeseries until you ran out of memory. (On a
> > > linux system you can
> > > > watch this happen by opening another window
> and
> > > executing "top"). I
> > > > usually experience this problem with much
> longer
> > > data sets, but all
> > > > depends on your system load and horizontal
> > > resolution.
> > > >
> > > > I have a patched a version of GrADS for use at
> > > NASA/GMAO which is
> > > > built with the sdfopen/xdfopen code from
> v1.7b9
> > > sdf/xdfopen code,
> > > > with some of the bugfixes up to v1.9b4. Since
> it
> > > *may* not address all
> > > > NetCDF files supported by v1.9b4, COLA
> > recommended
> > > that this patch not
> > > > be included in v1.9.0-rc1. Their recommended
> > > action is to prepare a
> > > > ctl file with "DTYPE netcdf" and use the
> "open"
> > > command. Another
> > > > option is the just released v2, but this
> version
> > > is still in active
> > > > development and is not as stable as v1.9 yet.
> > BTW,
> > > this problem is
> > > > common to netcdf and hdf files alike.
> > > >
> > > > Arlindo
> > > >
> > > >
> > > > --
> > > > Arlindo da Silva
> > > > dasilva at alum.mit.edu
> > <mailto:dasilva at alum.mit.edu>
> > >
> > >
> > > --
> > >
> >
>
----------------------------------------------------------------
> > > Adam H. Sobel
> > > Department of Applied Physics and Applied
> > > Mathematics
> > > Department of Earth and Environmental Sciences
> > > Columbia University
> > > 500 W. 120th Street, Room 217
> > > New York, NY 10027 USA
> > > tel: 212-854-6587
> > > fax: 212-854-8257
> > > e-mail: ahs129 at columbia.edu
> > > web: http://www.columbia.edu/~ahs129/home.html
> > >
> > > UNTIL AUGUST 2008, ON SABBATICAL at
> > > Bureau of Meteorology Research Centre
> > > GPO Box 1289
> > > Melbourne VIC 3001
> > > Australia
> > > tel: +61-(0)3-9669-4645
> > >
> >
>
----------------------------------------------------------------
> > >
> >
> >
> >
> >
> >
>
____________________________________________________________________________________
> > Looking for last minute shopping deals?
>
=== message truncated ===
____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
More information about the gradsusr
mailing list