!+
!
! .. _average_demo.jnl:
!
! ================
! average_demo.jnl
! ================
!
! DESCRIPTION
! ===========
!
! plot sst anomaly from ferret demo dataset
! 
! no other interest than learning ferret and be able to work on tropflux 
! dataset (see :ref:`average_pb.jnl`)
!
! EXAMPLES
! ========
!
! ::
!
!  ferret> GO average_demo
!
! or::
!
!  $ ferret -nojnl -script average_demo.jnl
!
! SEE ALSO
! ========
!
! TODO
! ====
!
! make it work
!
! save figure
!
! save data to be compare to equivalent nco or cdo or .. implementation
!
! do the same with levitus, coads (ie a well known dataset)
! see http://www.oc.nps.edu/~bird/oc3030/ferret/ferret_notes.html
! ocean_atlas_subset.nc seems to be the closest dataset provided by ferret
! to tropflux one because of xt,x,y,z axis of sst variable
!
! EVOLUTIONS
! ==========
!
! $Id$
!
! $URL$
!
! - fplod 20111104T150732Z aedon.locean-ipsl.upmc.fr (Darwin)
!
!   * creation to try to answer to [TropFlux:7] see :ref:`average_pb.jnl`
!     and learning ferret
!
!-
CANCEL DATA/ALL
CANCEL SYMBOL/ALL
!
! pour creer des metafile convertible en ps
set mode metafile
! pour avoir des traits plus gros
go bold
!
! pour les log
! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html
let date = "($session_date)"
let datetime = "($session_date) ($session_time)"
!
SET DATA ocean_atlas_subset
!
! have a look to data structure
SHOW DATA/VARIABLES
SHOW GRID TEMP
! have a look to memory
SHOW MEMORY
!
! original statement in [tropflux:7]
! list/t=01-may-2001:30-jun-2001 TEMP[x=66:78@ave,y=6:18@ave]
SET REGION/z=1/x=66:78/y=6:18/t=16-jan-0001:17-apr-0001
list TEMP[z=1,x=@ave,y=@ave]
LET sst_ave1 = TEMP[z=1,x=@ave,y=@ave]
SHOW GRID sst_ave1
LIST sst_ave1
LET sst_ave1 = TEMP[z=1,x=@ave,y=@ave,t=@ave]
SHOW GRID sst_ave1
LIST sst_ave1
PLOT sst_ave1
LET/TITLE="TEMP anomaly x=66:78" sst_anom1 = TEMP - sst_ave1
SHOW DATA/VARIABLES
!++ PLOT sst_anom1
!
! original statement in [tropflux:7]
! list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave]
list/t=01-may-2001:30-jun-2001 sst[x=66:70@ave,y=6:18@ave]
LET sst_ave2 = sst[x=66:70@ave,y=6:18@ave,t=01-may-2001:30-jun-2001]
LET/TITLE="TEMP anomaly x=66:70" sst_anom2 = sst - sst_ave2
PLOT sst_anom2
!
! save data
! ++ LIST/file=average_demo.output/format=(20E11.3)/order=xy/L=7 sst_ave1 sst_ave2
!
! save figure++
!++exit