!+
!
! .. _read.jnl:
!
!
! ========
! read.jnl
! ========
!
!
! DESCRIPTION
! ===========
!
! read TROPFLUX file ${TROPFLUX_D}/TropFlux_19890101_20091231.nc
!
! plot the first time step of swr
!
! EXAMPLES
! ========
!
! ferret> CANCEL DATA/ALL
! ferret> GO read
!
! SEE ALSO
! ========
!
! matlab/octave similar program
!
! :ref:`read.m`
! 
!
! TODO
! ====
!
! usage of $TROFLUX_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
!
! EVOLUTIONS
! ==========
!
! - fplod 20110106T164602Z aedon.locean-ipsl.upmc.fr (Darwin)
!
!   * add plot
!
! - fplod 20110106T090443Z aedon.locean-ipsl.upmc.fr (Darwin)
!
!   * creation thanks to ACP
!
!-
!
! pour creer des metafile convertible en ps
set mode metafile 
! pour avoir des traits plus gros
go bold 

LET yyyymmdd_b = "19890101"
LET yyyymmdd_e = "20091231"

! ++ pb floatstr si yyyymmdd_b est un entier LET fullfilename = "/usr/zeus/incas/fplod/tropflux_d"+"TropFlux_" + floatstr(yyyymmdd_b,"(I8.8)") + "_" + floatstr(yyyymmdd_e,"(I8.8)") + ".nc"
LET filename = "TropFlux_" + yyyymmdd_b + "_" + yyyymmdd_e + ".nc"
LET fullfilename = "/usr/zeus/incas/fplod/tropflux_d/" + filename
LIST fullfilename
! open 
! ++ use fullfilename
use "/usr/zeus/incas/fplod/tropflux_d/TropFlux_19890101_20091231.nc"
!
! 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