!+
!
! .. _read.jnl:
!
!
! ========
! read.jnl
! ========
!
!
! DESCRIPTION
! ===========
!
! read TROPFLUX file ${PROJECT_D}/TropFlux_19890101_20091231.nc
!
! plot the first time step of swr
!
! EXAMPLES
! ========
!
! ::
!
!  ferret> GO read
!
! or::
!
!  $ ferret -nojnl -script read.jnl
!
! SEE ALSO
! ========
!
! matlab/octave similar program : :ref:`read.m`
! 
! IDL+SAXO similar program : :ref:`read.pro`
!
! TODO
! ====
!
! usage of ${TROPFLUX_D}
!
! make it usable for real end-user
!
! learn ferret
!
! clean metafiles
!
! pb with Fprint on zeus probably because of an mtt problem
!
! pb with ferret on adonis
!
! pb with underscore in dataset label on zeus
!
! pas de floatstr function sur zeus ?!
! cf http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/handing-string-data-symbols/STRING-FUNCTIONS?searchterm=floatstr
!
! pb missing value
! EVOLUTIONS
! ==========
!
! - fplod 20110106T164602Z aedon.locean-ipsl.upmc.fr (Darwin)
!
!   * add plot
!
! - fplod 20110106T090443Z aedon.locean-ipsl.upmc.fr (Darwin)
!
!   * creation thanks to ACP
!
!-
!
CANCEL DATA/ALL
CANCEL SYMBOL/ALL
!
! pour creer des metafile convertible en ps
set mode metafile 
! pour avoir des traits plus gros
go bold 

!++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2010/msg00690.html
!let mm = 1
!say `mm,zw=2`
! ++ http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00771.html
!repeat/range=1:2/name=m (use "http://dodsserver/sst_`m`.nc"; show dat)
!repeat/range=1:12/name=mm ( show sym mm; say `mm,zw=2`) 
!repeat/range=1:12/name=mm ( show sym mm; LET CMM B`mm,zw=2`A; LIST CMM) 

! pour les log
! http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2009/msg00036.html
show symbol session*
let date = "($session_date)"
let datetime =  "($session_date) ($session_time)"
list date
list datetime

! DEFINE VARIABLE yyyymmdd_b = "19890101"
! DEFINE VARIABLE yyyymmdd_e = "20091231"
!++ DEFINE VARIABLE cb = {spawn:printf '%8.8d' `b` "}

DEFINE SYMBOL yyyymmdd_b = 19890101
DEFINE SYMBOL yyyymmdd_e = 20091231

! ++ pb floatstr si yyyymmdd_b est un entier DEFINE VARIABLE fullfilename = "/usr/zeus/incas/fplod/tropflux_d"+"TropFlux_" + floatstr(yyyymmdd_b,"(I8.8)") + "_" + floatstr(yyyymmdd_e,"(I8.8)") + ".nc"
!DEFINE VARIABLE filename = "TropFlux_" + yyyymmdd_b + "_" + yyyymmdd_e + ".nc"
DEFINE SYMBOL filename = TropFlux_($yyyymmdd_b)_($yyyymmdd_e).nc
! ++ $PROJECT_D defini ?
!DEFINE SYMBOL tropflux_d = /usr/zeus/incas/fplod/tropflux_d/
!DEFINE SYMBOL fullfilename = ($tropflux_d)($filename)
!DEFINE VARIABLE tropflux_d = {spawn:"echo $PROJECT_ID"}
DEFINE SYMBOL tropflux_d = {spawn:"echo $PROJECT_ID"}
SHOW SYM tropflux_d
! DEFINE SYMBOL fullfilename = "$tropflux_d($filename)"
DEFINE SYMBOL fullfilename = `tropflux_d`($filename)
CANCEL SYMBOL filename
! open 
! test if the file exists
!++DEFINE VARIABLE a = {"dummy", SPAWN:"ls `fullfilename`"} 
!++DEFINE VARIABLE nfiles = `a,RETURN=IEND`
!++CANCEL a
!++LIST nfiles
!++IF `nfiles EQ 2` THEN EXIT
!++CANCEL nfiles

! ++ use fullfilename
!use "/usr/zeus/incas/fplod/tropflux_d/TropFlux_19890101_20091231.nc"
!use "`fullfilename`"
use "($fullfilename)"
!
! plot the first timestep of swr
shade/l=1 swr
!
! plot world map
go land
!
! create a Postscript file
sp Fprint -l cps -p portrait -o TropFlux_19890101_20091231.ferret.ps metafile.plt
!
! create a GIF file
frame/file=TropFlux_19890101_20091231.ferret.gif
!
!++exit