[ferret_users] time-series and mark

[dessy berlianty <dessyberlianty@xxxxxxxxx>]

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