vis5d-users Mailing List for Vis5d+

Sorry - I should have used TCL syntax... - no parens, if you use parens you
will get the error. It should be:

vis5d_save_scene $ctx "lamps1.wrl" VIS5D_VRML

I thought the mixkit stuff was pretty verbose... I haven't used it in sometime - let me get back to you on that one.

jeff

Chris Moore wrote:
> 
> Hi Jeff/All,
> 
> I seem to get an error with vis5d_save_scene:
> 
>   Tcl_EvalFile result: invalid command name
> "vis5d_save_scene(0,lamps1.wrl,VIS5D_VRML)"
> 
> Here's some of my code:
> 
>   set format VIS5D_VRML
>   for {set time 1} {$time <= 5} {set time [expr $time + 1]} {
>         set fname lamps$time.wrl
>         vis5d_save_scene($ctx,$fname,$format)
>   }
> 
> Other bits of my Tcl script seem to work fine... turning on/off isosurfaces,
> etc.  Using vis5d+-1.2.1 on a Dec Alpha.
> 
> Also: how do I tell if mixkit is working?  Any way of getting a polygon count
> for an isosurface and compare to a version compiled without mixkit?
> 
> Hope things are well in Boulder.
> 
> Cheers,
> Chris
> 
> ______________________________________________________________
> Christopher Moore                             tel 206.526.6779
> Research Scientist, Oceanography              fax 206.526.6744
> University of Washington/JISAO/NOAA-PMEL
> --------------------------------------------------------------
> 
> "Jeff W. Boote" wrote:
> >
> > vis5d_save_scene($ctx,$fname,$format)
> >
> > $fname = filename to save
> > $format = VIS5D_VRML
> >
> > (It is written such that other scene formats can be added by adding new
> > flags for the $format variable.)
> >
> > jeff

-- 
Jeff W. Boote  <bo...@uc...>      The more that you read, the more
Software Engineer                    things you will know.  The more
VETS/SCD/NCAR                        that you learn, the more places
Boulder, Colorado, USA               you'll go.        Dr. Suess

Hi Jeff/All,

I seem to get an error with vis5d_save_scene:

  Tcl_EvalFile result: invalid command name    
"vis5d_save_scene(0,lamps1.wrl,VIS5D_VRML)"

Here's some of my code:

  set format VIS5D_VRML
  for {set time 1} {$time <= 5} {set time [expr $time + 1]} {
        set fname lamps$time.wrl
        vis5d_save_scene($ctx,$fname,$format)
  }

Other bits of my Tcl script seem to work fine... turning on/off isosurfaces,
etc.  Using vis5d+-1.2.1 on a Dec Alpha.

Also: how do I tell if mixkit is working?  Any way of getting a polygon count
for an isosurface and compare to a version compiled without mixkit?

Hope things are well in Boulder.

Cheers,
Chris

______________________________________________________________
Christopher Moore                             tel 206.526.6779
Research Scientist, Oceanography              fax 206.526.6744
University of Washington/JISAO/NOAA-PMEL
--------------------------------------------------------------


"Jeff W. Boote" wrote:
> 
> vis5d_save_scene($ctx,$fname,$format)
> 
> $fname = filename to save
> $format = VIS5D_VRML
> 
> (It is written such that other scene formats can be added by adding new
> flags for the $format variable.)
> 
> jeff

I'm glad I could help!

Dan

Daniel M. McCormick, TSgt, USAF
Lead Programmer, Fine Scale Models Team
Meteorological Models Branch
HQ Air Force Weather Agency
Offutt Air Force Base, Nebraska
(402) 294-5558
DSN: 271-5558


-----Original Message-----
From: Jim Edwards [mailto:jim...@al...]
Sent: Thursday, November 29, 2001 5:32 AM
To: Giovanni Leoncini
Cc: McCormick Daniel M TSgt HQ AFWA/DNXM; vis...@ss...
Subject: Re: halleluja!!



I try.  Thanks Dan for stepping forward with the right answer.

Jim

Giovanni Leoncini wrote:

> Hi Daniel! ;)
>
> you did it! you fixed it! ;)))
>
> thanks a lot, also to Jim, of course!!
>
> Giovanni
>
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1

I try.  Thanks Dan for stepping forward with the right answer.

Jim

Giovanni Leoncini wrote:

> Hi Daniel! ;)
>
> you did it! you fixed it! ;)))
>
> thanks a lot, also to Jim, of course!!
>
> Giovanni
>
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1

Dear Vis5d users,
	I have written a C program that converts the global ~1 Km
elevation dataset known as "GLOBE"
	(http://www.ngdc.noaa.gov/seg/topo/globe.shtml)
to the vis5d topography format on an arbitrary area of the world. Since I
have noticed this is a common problem, I decided to release the program to
the public. If you are interested, you can download it from the grib2v5d
download area:

	http://sourceforge.net/project/showfiles.php?group_id=14299
	
Please read the instruction in the README and in the Makefile files before
compiling. Notice that, in order to appreciate the increase in the
topography resolution, you may have to run vis5d with the -hirestopo
option. If you think the program is useful it could be included in
vis5d distribution.
	best wishes, Davide


-- 
__________________________________________________________
Davide Cesari	Servizio Meteorologico Regionale      __
 tel.      (39) 051/284541                            ||\
 fax       (39) 051/284664                            |||\
 e-mail    DC...@sm...                    |||/
 Address:  ARPA-SMR, Viale Silvani 6, 40122 Bologna, Italy		
__________________________________________________________

vis5d_save_scene($ctx,$fname,$format)

$fname = filename to save
$format = VIS5D_VRML

(It is written such that other scene formats can be added by adding new
flags for the $format variable.)

jeff



Chris Moore wrote:
> 
> Now that the NCAR enhancements are a part of the regular distribution,
> can anyone tell me how I can use "Save Scene as VRML" in a Tcl script?
> 
> Bill Hibbard's Vis5d Tcl script tutorial (obviously) doesn't cover the
> new option.  Is this possible?
> 
> Cheers,
> Chris
> 
> ______________________________________________________________
> Christopher Moore                             tel 206.526.6779
> Research Scientist, Oceanography              fax 206.526.6744
> University of Washington/JISAO/NOAA-PMEL
> --------------------------------------------------------------

-- 
Jeff W. Boote  <bo...@uc...>      The more that you read, the more
Software Engineer                    things you will know.  The more
VETS/SCD/NCAR                        that you learn, the more places
Boulder, Colorado, USA               you'll go.        Dr. Suess

Now that the NCAR enhancements are a part of the regular distribution,
can anyone tell me how I can use "Save Scene as VRML" in a Tcl script?

Bill Hibbard's Vis5d Tcl script tutorial (obviously) doesn't cover the
new option.  Is this possible?

Cheers,
Chris


______________________________________________________________
Christopher Moore                             tel 206.526.6779
Research Scientist, Oceanography              fax 206.526.6744
University of Washington/JISAO/NOAA-PMEL
--------------------------------------------------------------

Hi Daniel! ;)


you did it! you fixed it! ;)))

thanks a lot, also to Jim, of course!!


Giovanni




__________________________________________________
Do You Yahoo!?
Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
http://geocities.yahoo.com/ps/info1

Hi Giovanni - it looks like the problem resides in the value for
proj_args(3). 
This represents the row of the nearest pole (in this case, north), if the
pole 
were located on a larger, superimposed grid of the same scale.  Looking at
the 
corner points of the grid:

                             North Pole
                          90.000N, 95.000W
                            (53,178.745)





     54.536N, 152.856W  ----------------------  57.290N, 49.385W
         (1,65)                                     (93,65)
           |                                           |
           |                                           |
           |                                           |
     12.190N, 133.459W  ----------------------  14.335N, 65.091W
         (1,1)                                       (93,1)

Since the pole is at 90N, even with the curvature of the earth in the LC
grid, 
we can see that it would not, in fact, reside on this particular grid.
Vis5d 
starts counting the rows at the top, not the bottom.  As a result, for the
North Pole, the following equation comes in to play:

pole_row = (rows - 1) - pole_row

So you should end up with the following:

nr = 65
nc = 93

projection = 2

proj_args(1) = 25.0
proj_args(2) = 25.0
proj_args(3) = -114.745
proj_args(4) = 53.0
proj_args(5) = 95.0     // this one may be -95.0
proj_args(6) = 81.2705

I hope this helps.

Regards,
Dan McCormick

Giovanni Leoncini wrote:

> Hi All! ;)
> 
> 
> I'm trying to convert some eta netcdf files into
> vis5d. The AWIPS 211 eta files are on lambert
> conformal
> (http://www.ncep.noaa.gov/NCO/PMB/products/grids.html#GRID211)
> and in my code I'm passing, to v5dCreate, these
> parameter for the lambert projection:
> 
>       proj_args(1) = 25
>       proj_args(2) = 25
>       proj_args(3) = 53.00
>       proj_args(4) = 178.745
>       proj_args(5) = -95
>       proj_args(6) = 81.2705
> 
> the met fields are consistent with topography, but
> everything is rotated 90 deg clockwise, and I see
> hawaii north of california. Playing with these values
> (i.e. switching 3 and 4 or using 95 or 265 instead of
> -95) I'm only able to see different parts of the
> globe!
> 
> If I switch grid and colums I get the fields right,
> but the rockies still go from west to east!! 
> 
> any clue?
> 
> thanks a lot!
> 
> Giovanni 
> 
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1
> 
> 

What topo file are you using?  Is it possible that the data is ordered
wrong there?  If you can create a minimal data set which exhibits the
problem I will try to look at it here.    Don=B4t mail it to the list,
mail it directly to me or put it on an
anonymous ftp site and send me the address.

Jim


Giovanni Leoncini wrote:

> Jim,
>
> I'm not quite sure I understood how you reorder the
> array here,
>
> > #define DATA(R,C)  data[R*cols+C]
> > #define TEMP(R,C)  temp[C*rows+rows-R-1]
> to me this sounds like a rotation, that brings a (rows
> X cols) array into another (rows X cols) array, am I
> right?
>
> my netcdf array is flipped along longitude, I ordered
> it as vis5d requires, but the topography is still
> turned. I mean the sw corner in terms of lat lon is
> the upper left corner of vis5d, the nw is the upper
> right, and so on, no matter how I orient my data.
>
> I also tried to reproduce the rotation you suggested
> (and a few variations), but I only changed the data,
> not the topography in the background.
>
> I ran out of ideas. is anybody else using lambert
> conformal?
>
> thanks a lot!!
>
> Giovanni
>
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1

Jim,


I'm not quite sure I understood how you reorder the
array here, 

> #define DATA(R,C)  data[R*cols+C]
> #define TEMP(R,C)  temp[C*rows+rows-R-1]
to me this sounds like a rotation, that brings a (rows
X cols) array into another (rows X cols) array, am I
right? 

my netcdf array is flipped along longitude, I ordered
it as vis5d requires, but the topography is still
turned. I mean the sw corner in terms of lat lon is
the upper left corner of vis5d, the nw is the upper
right, and so on, no matter how I orient my data.

I also tried to reproduce the rotation you suggested
(and a few variations), but I only changed the data,
not the topography in the background.

I ran out of ideas. is anybody else using lambert
conformal? 

thanks a lot!!

Giovanni



__________________________________________________
Do You Yahoo!?
Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
http://geocities.yahoo.com/ps/info1

hello every body
can somone help me please to define a function in a
C-program that allaow the cursor to return to a new line in
a file?
what I've done is that :

#define NL(n) { int l; for(l=3D0; l<(n); l++) fputc(cfp,"\n");
}

where cfp design the name of the file and NL(n) will write
"n" new lines in cfp.
The problem is that it print ^A and not a new line !!!
I've changed fputc by fputs and than by fprintf but the
problem still exist

the compiler that I've used is cc and the machine is an IBM
RS6000

Thanks
______________________________________________________
Bo=EEte aux lettres - Caramail - http://www.caramail.com

Giovanni,

I think the problem is that data is stored in the netcdf file in a
different order than it is in a vis5d file.
You need to reorder the data in the arrays after you read it in.

For example when reading a grads file the following rotation is applied
to the data (in horizontal slices) after it is read:  (from
read_grads_i.c)

#define DATA(R,C)  data[R*cols+C]
#define TEMP(R,C)  temp[C*rows+rows-R-1]

   for (i=0;i<rows;i++) {
      for (j=0;j<cols;j++) {

         if (DATA(i,j)==missing_value) {
            TEMP(i,j) = MISSING;
         }
         else {
            TEMP(i,j) = DATA(i,j);
         }
      }
   }

   memcpy( data, temp, rows*cols*sizeof(float) );




Jim



Giovanni Leoncini wrote:

> Hi All! ;)
>
> I'm trying to convert some eta netcdf files into
> vis5d. The AWIPS 211 eta files are on lambert
> conformal
> (http://www.ncep.noaa.gov/NCO/PMB/products/grids.html#GRID211)
> and in my code I'm passing, to v5dCreate, these
> parameter for the lambert projection:
>
>       proj_args(1) = 25
>       proj_args(2) = 25
>       proj_args(3) = 53.00
>       proj_args(4) = 178.745
>       proj_args(5) = -95
>       proj_args(6) = 81.2705
>
> the met fields are consistent with topography, but
> everything is rotated 90 deg clockwise, and I see
> hawaii north of california. Playing with these values
> (i.e. switching 3 and 4 or using 95 or 265 instead of
> -95) I'm only able to see different parts of the
> globe!
>
> If I switch grid and colums I get the fields right,
> but the rockies still go from west to east!!
>
> any clue?
>
> thanks a lot!
>
> Giovanni
>
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1

Hi All! ;)


I'm trying to convert some eta netcdf files into
vis5d. The AWIPS 211 eta files are on lambert
conformal
(http://www.ncep.noaa.gov/NCO/PMB/products/grids.html#GRID211)
and in my code I'm passing, to v5dCreate, these
parameter for the lambert projection:

      proj_args(1) = 25
      proj_args(2) = 25
      proj_args(3) = 53.00
      proj_args(4) = 178.745
      proj_args(5) = -95
      proj_args(6) = 81.2705

the met fields are consistent with topography, but
everything is rotated 90 deg clockwise, and I see
hawaii north of california. Playing with these values
(i.e. switching 3 and 4 or using 95 or 265 instead of
-95) I'm only able to see different parts of the
globe!

If I switch grid and colums I get the fields right,
but the rockies still go from west to east!! 

any clue?

thanks a lot!

Giovanni 

__________________________________________________
Do You Yahoo!?
Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
http://geocities.yahoo.com/ps/info1

Jim,

Do you know where is example data file that will be available to be
converted to v5d format? These data files could be downloaded in web.

Thank you,
Jianyong

Several example conversion programs can be found in the util directory.
The documentation also includes a chapter
titled =B4 Putting Your Data Into Vis5D=B4

Jim



Jianyong Bai wrote:

> Do you know how to convert my own data into v5d format? Is there anybod=
y
> who has some examples to explain how?

Do you know how to convert my own data into v5d format? Is there anybody
who has some examples to explain how?

Dear Davide,

Using character*10 seems to be imperative, 
your suggestion solved my problem.  
Thank you very much.

Kind regards,
Marcus


On Mon, 19 Nov 2001, Davide Cesari wrote:

> Dear Marcus,
> 	I regularly compile grib2v5d (you can get the sources at
> http://grib2v5d.sf.net) with pgf90 Portland compiler and it correctly
> generates variable names on Linux (actually I don't use a DATA statement
> but I read variable names from a NAMELIST file); maybe your problem is
> that the strings composing the vector passed to v5dcreate are not 10
> characters long but less (I guess 7)? I think this is a requirement for
> v5dcreate, the example in the manual reads:
> 
>          character*10 varname(MAXVARS)
> 
> 	If the problem is not there, maybe some lines of code would be
> useful to understand better what's going on.
> 	best wishes,
> 		Davide
> -- 
> __________________________________________________________
> Davide Cesari	Servizio Meteorologico Regionale      __
>  tel.      (39) 051/284541                            ||\
>  fax       (39) 051/284664                            |||\
>  e-mail    DC...@sm...                    |||/
>  Address:  ARPA-SMR, Viale Silvani 6, 40122 Bologna, Italy		
> __________________________________________________________

Dear Marcus,
	I regularly compile grib2v5d (you can get the sources at
http://grib2v5d.sf.net) with pgf90 Portland compiler and it correctly
generates variable names on Linux (actually I don't use a DATA statement
but I read variable names from a NAMELIST file); maybe your problem is
that the strings composing the vector passed to v5dcreate are not 10
characters long but less (I guess 7)? I think this is a requirement for
v5dcreate, the example in the manual reads:

         character*10 varname(MAXVARS)

	If the problem is not there, maybe some lines of code would be
useful to understand better what's going on.
	best wishes,
		Davide


On Mon, 19 Nov 2001, Marcus Koehler wrote:

> Dear group members,
> 
> I attempt to create v5d files using Portland Group Fortran on a Linux
> machine.  I use v5dcreatesimple for a regular spaced grid.
> 
> However, the variable names seem to be not properly identified.  They
> are given as a character field in a data statement.  But instead of
> recognising the variables, P, U, V, T, Theta, PV, and SpecH,
> it appears that the strings are being "glued" together.
> 
> See here the output from v5dinfo:
> 
> File format: v5d  version: 4.3
> Compression:  1 byte per gridpoint.
> header size=11804
> sizeof(v5dstruct)=169068
> 
> NumVars = 7
> Var  Name       Units      Rows  Cols  Levels LowLev  MinVal       
>   1  P      U   Pa          36    72    30      0     1000             
>   2      V      m/s         36    72    30      0     -40.8064   
>   3   T      T  m/s         36    72    30      0     -57.8602   
>   4  eta  PV    K           36    72    30      0     185.118   
>   5    SpecH    K           36    72    30      0     239.123   
>   6  a          PVU         36    72    30      0     -207.815   
>   7   m/s       g/g         36    72    30      0     0.000458169
> 
> etc...
> 

-- 
__________________________________________________________
Davide Cesari	Servizio Meteorologico Regionale      __
 tel.      (39) 051/284541                            ||\
 fax       (39) 051/284664                            |||\
 e-mail    DC...@sm...                    |||/
 Address:  ARPA-SMR, Viale Silvani 6, 40122 Bologna, Italy		
__________________________________________________________

Dear group members,

I attempt to create v5d files using Portland Group Fortran on a Linux
machine.  I use v5dcreatesimple for a regular spaced grid.

However, the variable names seem to be not properly identified.  They
are given as a character field in a data statement.  But instead of
recognising the variables, P, U, V, T, Theta, PV, and SpecH,
it appears that the strings are being "glued" together.

See here the output from v5dinfo:

File format: v5d  version: 4.3
Compression:  1 byte per gridpoint.
header size=11804
sizeof(v5dstruct)=169068

NumVars = 7
Var  Name       Units      Rows  Cols  Levels LowLev  MinVal       
  1  P      U   Pa          36    72    30      0     1000             
  2      V      m/s         36    72    30      0     -40.8064   
  3   T      T  m/s         36    72    30      0     -57.8602   
  4  eta  PV    K           36    72    30      0     185.118   
  5    SpecH    K           36    72    30      0     239.123   
  6  a          PVU         36    72    30      0     -207.815   
  7   m/s       g/g         36    72    30      0     0.000458169

etc...


I had a similar problem with v5dsetunits and adding a '\0' to each
unit string took care of it.  However this doesn't work with the
variable names.

I would be very grateful for suggestions on how to solve this problem. 

Best regards
Marcus

Thanks Steven and others:

OSF has problems using limit with tcsh!?  I ran csh and used:

limit datasize unlimited

and now we have a working version of Vis5d using DEC's OpenGL libraries!

Thanks again for all the work getting the NCAR version re-united with
the Vis5d+ strain.

Cheers,
Chris

______________________________________________________________
Christopher Moore                             tel 206.526.6779
Research Scientist, Oceanography              fax 206.526.6744
University of Washington/JISAO/NOAA-PMEL
--------------------------------------------------------------

Try using -mbs 0, which should remove any memory limit.  (We should
probably change this to be the default...a limit of 32MB is ridiculous
these days.)

Steven

On Thu, 15 Nov 2001, Chris Moore wrote:
> I get this error message even when using -mbs 34 (or -mbs 75).  I should
> have plenty of memory available.  I can run any number of slices,
> contours, volume renderings, trajectories, wind vectors, etc, while
> animating in stereo mode.  But not a single small isosurface.

I compiled the latest version on a DEC Alpha with OpenGL.  The good news
is that it runs in stereo using OpenGL (thanks NCAR), and that the
installation was painless (thanks Sourceforge)! 

But... I get memory problems.  Heres the output when I try to render a
simple isosurface:

>/home/elegua/cmoore/ 150 % vis5d canyon.v5d 
>Setting visual to 24 bit truecolor with default color map
>Setting visual to 24 bit truecolor with default color map
>Stereo Mode Enabled
>adding dtx to ctx 0 0 0x40bb1020
>Cache size: 90 grids 24 4
> Hint... To run Vis5D more efficiently try setting canyon.v5d to '-mbs >34'
> You do not have enough memory to create isosurfaces.
> You do not have enough memory to create isosurfaces.
> You do not have enough memory to create isosurfaces.
> You do not have enough memory to create isosurfaces.
> ... (deleted repeating lines)

I get this error message even when using -mbs 34 (or -mbs 75).  I should
have plenty of memory available.  I can run any number of slices,
contours, volume renderings, trajectories, wind vectors, etc, while
animating in stereo mode.  But not a single small isosurface.

I tried the LAMPS.v5d data set: same error.

Any ideas?

Yours,
Chris

______________________________________________________________
Christopher Moore                             tel 206.526.6779
Research Scientist, Oceanography              fax 206.526.6744
University of Washington/JISAO/NOAA-PMEL
--------------------------------------------------------------

Hi,
Here's a source for map data, including coastline and politicical boundaries.
http://rimmer.ngdc.noaa.gov/coast/getcoast.html

In the portion of the form titled "Coastline data base"  there are choices
including political boundaries.  It returns ASCII data files of longitude,
latitude.

Abdoulaye SARR wrote:

> Hi all,
>
> Could someone please tell me where can I find a map with political boundaries
> to be used with vis5d.
>
> Thanks

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA  (206)526-6246

For political boundaries, you may find the following map of value:
http://viewer.digitalearth.gov/viewer.cgi?context=3D100557681921049271.xm=
l

This map includes two layers with political boundaries: states and nation=
al.

The access method at the link above uses the OGC Web Map service=20
(WMS), but you may be able to use the source data in vis5D.  Other=20
visualization programs have integrated the WMS access protocol to=20
include maps in the context of the more powerful visualization client.

George P.

At 10:50 AM +0000 11/12/01, Jim Edwards wrote:
>I have been searching for a good electronic source of world political bo=
undary
>information for a long time, I=B4ve yet to come up with anything=20
>satisfactory.  If
>you do come up with the information there is a sample program util/makem=
ap.c
>which explains how to get it into a format suitable to vis5d.  If people=
 are
>using and want to contribute maps other than the two (OUTLSUPW  OUTLUSAM=
 )
>supplied with vis5d I would consider adding those to the distribution.
>
>Jim
>
>
>
>Abdoulaye SARR wrote:
>
>>  Hi all,
>>
>>  Could someone please tell me where can I find a map with political=20
>>boundaries
>>  to be used with vis5d.
>>
>>  Thanks

--=20
________________________
George Percivall

(301) 286-4073  Office
(410) 971-0382  Mobile
per...@gs...