thegooglecodearchive / ddscat

Automatically exported from code.google.com/p/ddscat
0 stars 0 forks source link

The VTR file with Poynting vector can not be visualized: Segmentation fault (core dumped) #13

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago
Hello,

I would like to plot E field as vectors, but trying to modify DDPOSTPROCESS.f90 
accordingly I encountered an error, namely Segmentation fault (core dumped). 
Testing exaples provided I found out that the convertion to the VTR format does 
not work even for Poynting vector (ddpostprocess executed with IVTR = 3) 
causing the same error. As the problem can be seen already with the basic 
example, likely it originates from some issues in writing procedure for vectors 
in DDPOSTPROCESS.f90, but unfortunately it is beyond my understanding.

What steps will reproduce the problem?

1. Start the example provided in ELLIPSOID_NEARFLD_B with NRFLD = 2. Note that 
the file w000r000k000.EB1 and 2 are 86.1 MB each instead of 133 MB as mentioned 
in the user guide for 7.3 version (p. 83)
2. Proceed with ddpostprocess with IVTR = 3 expecting to get a VTR file with 
Poynting vectors for futher visualization (the opportunity is listed in 
DDPOSTPROCESS.f90, line 85)
3. Obtain a Segmentation fault (core dumped) in the log. The resulting VTR file 
contains neither scalars nor vectors (according to Mayavi2 messages when it is 
uploaded as a source file)

What is the expected output? What do you see instead?

The VTR file with Poynting vector components is created but can not be read by 
Mayavi2 properly.The error Segmentation fault (core dumped) is listed in the 
log. The file is around 81 MB compared with VTR created with IVTR = 1 or 2

What version of the product are you using? On what operating system?

I run DDSCAT7.3 on Ubuntu 12.04 LTS 64 bit 

Original issue reported on code.google.com by aalyamk...@gmail.com on 27 Nov 2013 at 10:38