Re: [ferret_users] Embed Ferret in Fortran program : A Solution

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi -
Good solution

To keep from having to have the sleep command in there, you could put in 
a read statement that would cause the program to pause while it waits 
for your input.  Replace the sleep statement with:

  print *, ' enter 1 to continue'
  read *, n

Then you can look at the plot for just a moment, or have it in front of 
you if you see something odd.

Haoguo Hu wrote:
> Hello, Ferrer-eers,
>
> Thanks you all for the reply. I got a solution on Glen Carver's idea, 
> that is, let Fortran write the script first, and then call Ferret.
> I'm running s sea ice ocean coupled model and I want to check the 
> model result of the first a few steps. Here is the solution:
>
> !-------Begin Fortran Program-----------------------
>       Program main  !Ice Ocean model
>       include 'iceocean.h'
>       call model_initial()
>
>       do istep = 1, iend
>          call pom()  !Princeton Ocean Model
>          call Ice()  !Sea Ice Model
>          call output( result_iceocean.nc )
>              open(10,file='plot.jnl')     !Call Ferre plotting
>           write (10,*) 'use "result_iceocean.nc"'
>               write (10,*) 'fill/k=1/l=',istep, 'etw' !surface elevation
>               write (10,*) 'vector/over/l=', istep, 'ua,va'  !velocity
>           write (10,*) 'spawn sleep 2s'     !Pause for viewing
>          close(10)                      !End call Ferre plotting
>       call system("ferret -nojnl -script plot.jnl")
>        enddo  !do istep = 1, iend
>
>        end  !end Program main  !Ice Ocean model
> !-------End Fortran Program-----------------------
>
> I do not like the Ferret window keep opening and closing, and it also 
> waste 2 seconds on sleeping (without sleeping you can see nothing), 
> but it is something better than nothing.
>
> I'm not sure whether or not there is some way Fortran and Ferret can 
> know each other, interactive communicating will make the Ferret more 
> powerful.
>
>
> Cheers
>
>
> Haoguo
>
>
> Glenn Carver wrote:
> > Hi,
> >
> > You can't make the ferret window stay open when you try to use ferret 
> > like this. Each call to system() in that loop invokes a separate 
> > command for Ferret.  You end up with 5 instances of Ferret running.
> >
> > You would be better off using the loop to write the journal 
> > instructions and then at the end of the program invoke Ferret (or do 
> > it yourself once the fortran code has finished).
> >
> > e.g.
> >
> > program main
> > integer :: istep
> > open(10,file='myplot.jnl')
> > do istep = 1, 5
> >     write (10,*) 'fill/k=1,l=',istep,' sst'
> > end do
> > call system('ferret -nojnl -script myplot.jnl')
> > end program
> >
> > Hope that helps.
> >
> >   Glenn
> >
> >
> >
> > On 19 Nov 2008, at 4:41pm, Haoguo Hu wrote:
> >
> > > Dear Ferret-eers,
> > >
> > > I try to embed Ferret in Fortran program when debugging, here is an 
> > > example:
> > >
> > > !----Begin Fortran program----------------------------
> > >        program main  !Fortran call Ferret
> > >     integer istep
> > >           do istep = 1,5
> > >          !----------
> > >             call system("ferret -nojnl -script plot.jnl")
> > >           enddo
> > >        end
> > > !----End Fortran program-------------------------------
> > >
> > >
> > > plot.jnl
> > > !----Begin plot.jnl------------------------------------
> > >  use "coads_climatology"
> > >  fill/k=1/l=1 sst
> > >  spawn sleep 10s  !pause for viewing and changing the L=2,3,..
> > > !----End plot.jnl--------------------------------------
> > >
> > > I can change the "L=1,2,3,.." in plot.jnl manually.
> > > The questions are:
> > > 1. How to make "L=istep" happen automatically?
> > > 2. How to let the Ferret window open only once and the plotting keep 
> > > going with "istep"?
> > >
> > >
> > > Thanks
> > >
> > >
> > > HH