Hi Jamie,
Another option, possibly easier, would be to create a .csh or .sh from
within FORTRAN and then do a single system command to run the .csh script.
This may reduce the I/O associated with several system commands. Below is
an example .csh script that would run nex2img (Note: the directory locations
may not be the same for your csh and Gemenviron).
scriptname: /data/nex2img.csh
------------
#!/bin/csh
source /data/gempak/Gemenviron
nex2img << NEX2IMG_EOF
GRDAREA=28.5;-105.28;49;-88.78
PROJ=MER
RADDUR=20
GFUNC=ncr
KXKY=1612;2280
CPYFIL=
RADTIM=20011121/2346
RADFRQ=
STNFIL=nexrad.tbl
RADMODE=PC
LUTFIL=upc_rad24.tbl
run
exit
NEX2IMG_EOF
gpend
-------------
Regards,
Evan Lowery
elowery@xxxxxxxxxxxx
From: gembud-bounces@xxxxxxxxxxxxxxxx
[mailto:gembud-bounces@xxxxxxxxxxxxxxxx] On Behalf Of Jamie Lahowetz
Sent: Thursday, September 03, 2009 9:53 PM
To: Stonie R. Cooper
Cc: gembud@xxxxxxxxxxxxxxxx
Subject: Re: [gembud] Fortran and gempak
Here is my test program:
PROGRAM gempak_test
IMPLICIT NONE
!##################################################
CHARACTER(200) :: command
!##################################################
command = "nex2img"
CALL system(command)
It seems that it gets stuck here.
command = "GRDAREA = 28.5;-105.28;49;-88.78"
CALL system(command)
command = "PROJ = MER"
CALL system(command)
command = "RADDUR = 20"
CALL system(command)
command = "GFUNC = ncr"
CALL system(command)
command = "KXKY = 1612;2280"
CALL system(command)
command = "CPYFIL ="
CALL system(command)
command = "RADTIM = 20011121/2346"
CALL system(command)
command = "RADFRQ ="
CALL system(command)
command = "STNFIL = nexrad.tbl"
CALL system(command)
command = "RADMODE = PC"
CALL system(command)
command = "RADFIL = /home/thor/20011121/2346.gif"
CALL system(command)
command = "LUTFIL = upc_rad24.tbl"
CALL system(command)
command = "list"
CALL system(command)
command = "run"
CALL system(command)
command = "exit"
CALL system(command)
END PROGRAM gempak_test
I'm sorry to say that I cant divulge about what the main program does to
much but it requires gif images created from raw radar data received from
NCDC. Hope this helps.
On Thu, Sep 3, 2009 at 8:32 PM, Stonie R. Cooper <stonewall@xxxxxxxxxxxxx>
wrote:
Jamie - you may have to be more specific. What does this "fortran"
program do to the gif images created from nex2img? What is your desired
result?
Your CALL statement is probably getting stuck because the neximg has
been set for continuous loop. If RADFRQ is set to anything other than 0
or nothing, it will attempt to repeat after the numerical value of the
assignment.
Stonie Cooper,
Planetary Data, Inc.
Prague, NE
402-663-6599
Jamie Lahowetz wrote:
> I am curious. I run nex2img to create gif images of radar scans. These
> are then ingested into a Fortran program. Does anyone know a way to do
> this directly from fortran? If I use CALL system statements the script
> gets stuck as soon as you start nex2img.
>
> --
> Jamie Ryan Lahowetz
> University of Nebraska - Lincoln
> Graduate Student - Geosciences
> 402.304.0766
> jlahowetz@xxxxxxxxx <mailto:jlahowetz@xxxxxxxxx>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gembud mailing list
> gembud@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe, visit:
http://www.unidata.ucar.edu/mailing_lists/
--
Jamie Ryan Lahowetz
University of Nebraska - Lincoln
Graduate Student - Geosciences
402.304.0766
jlahowetz@xxxxxxxxx
