Re: [ferret_users] time-series and mark

["Karl M. Smith" <karl.smith@xxxxxxxx>]

Hi Dessy,

If you list the first couple of values of your tvar (list tvar[l=1:5]) 
you will see that it is only reading every other line of your input file.
To skip just the first line on your input file, use /skip=1.  The '/' in 
your format causes it to skip a line each time it reads the next line of 
data.

Karl


dessy berlianty wrote:
> Dear Ferret Users,
>
>  
>
> I have question about plotting my time series data (110 days in hourly).
>
> Hope someone can give some help for my basic question.
>
> For information, following is description of my data (format:  
> f7.2,f12.1,4f12.5), i.e.:
>
>  
>
> HOURS:       SECONDS:    v1          v2          v3           v4       
>
>    1.00      3600.0    -0.87476    -0.64783    -0.86129    -0.90325
>
>    2.00      7200.0    -0.75127    -0.77024    -0.70260    -0.70213
>
>    3.00     10800.0    -0.85790    -0.70214    -0.84490    -0.86482
>
>    4.00     14400.0    -0.48004    -0.48037    -0.49582    -0.51512
>
> …
>
> 2639.00   9500400.0     0.73257     0.73720     0.73312     0.73499
>
> 2640.00   9504000.0    -0.04312     0.24171    -0.03135    -0.06120
>
>  
>
> Others information is about my script:
>
>  
>
> !----START SCRIPT---
>
> set data/save
>
> !define axis, start 11 Feb 2010 01:00 - 1 Jun 2010 00:00
>
> define axis/t="11-FEB-2010:01:00":"1-JUN-2010:00:00":1/unit=hour taxis
>
>  
>
> define grid/t=taxis tgrid
>
> !to skip header, use '/'
>
> file/format=(/f7.2,f12.1,4f12.5)/var="tvar,hvar,v1,v2,v3,v4"/grid=tgrid/order=t 
> abcd.ver
>
>  
>
> set data abcd.ver
>
> show data
>
>  
>
> !only interest to plot v3
>
> plot v3
>
>  
>
> !untuk convert dari ASCII ke netcdf
>
> save/clobber/file=elev.cdf "tvar,hvar,v1,v2,v3,v4"
>
> cancel data abcd.ver
>
>  
>
> can data/all
>
>  
>
> use elev.cdf
>
> set data elev.cdf
>
> show data
>
> set win/new
>
> plot v4
>
> !----END SCRIPT, dessy---
>
>  
>
>  
>
> Question:
>
>    1. After I run the script in Ferret window, for show data of
>       abcd.ver, I think the results is ok, like this:
>
> Name
>
> 	
>
> Title
>
> 	
>
> I
>
> 	
>
> J
>
> 	
>
> K
>
> 	
>
> L
>
> tvar
>
> 	
>
> tvar
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:2640
>
> hvar
>
> 	
>
> hvar
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:2640
>
> v1
>
> 	
>
> v1
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:2640
>
> v2
>
> 	
>
> v2
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:2640
>
> v3
>
> 	
>
> v3
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:2640
>
>  
>
> But the problem is when I run script to show data of elev.cdf:
>
> 1> .\elev.cdf (default)
>
> Name
>
> 	
>
> Title
>
> 	
>
> I
>
> 	
>
> J
>
> 	
>
> K
>
> 	
>
> L
>
> tvar
>
> 	
>
> tvar
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:1320
>
> hvar
>
> 	
>
> hvar
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:1320
>
> v1
>
> 	
>
> v1
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:1320
>
> v2
>
> 	
>
> v2
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:1320
>
> v3
>
> 	
>
> v3
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> ...
>
> 	
>
> 1:1320
>
>  
>
> It shows that L=1320 (or 55 days) and causing the x-axis is only 
> between 11 Feb 2010 - 6 April 2010 (it's only 55 days).
>
> I have no idea to change the x-axis become 11 Feb 2010 - 31 May 2010 
> (110 days).
>
>  
>
>    2. Does anyone have idea or suggestion how to expand graph from
>       long data to short data?
>    3. How to save my cdf data (after convert from ASCII) to my
>       directory folder?
>    4. How to make mark position in my domain area? I already try run
>       mp_demo.jnl, but I need more example about it.
>
>  
>
> That's all my questions, thanks in advance for all any suggestion and 
> help.
>
>  
>
> Sincerely yours,
>
>  
>
> dessy
>
>  
>
>  

--
Karl M. Smith, Ph.D.
Research Scientist and Software Developer
UW/JISAO and NOAA/PMEL/TMAP
karl.smith@xxxxxxxx
(206) 526-4806

"The contents of this message are mine personally and do
not necessarily reflect any position of the Government
or the National Oceanic and Atmospheric Administration."