[Tutor] TypeError: only length-1 arrays can be converted to Python scalars

2012-01-26 Thread stm atoc
Hi, I have a question regarding the TypeError. I have this script as follows in Python. I have run a program in Fortran: # coding: utf-8 from pylab import * import numpy import matplotlib.pyplot as pyplot import matplotlib.mlab as mlab t=3600 N = 100 dz = 1.0 with open("act_out.list", "r") as

[Tutor] changing coordinate

2012-01-11 Thread stm atoc
Hi: I am trying to define a new coordinate for My plot. I have a plot x axis is Concentration, Y axis, depth or Z. Z= -1000 um. and I would like to see the the trend and plot of concentration up to -300 um and see how the plot changed according to the depth. So, my main question is about how to

Re: [Tutor] unsupported operand

2011-12-18 Thread stm atoc
till if I do plot based on pyplot.pcolor, I have a nice and seems correct graph. So, actually, I am not sure what is the problem for showing the regular/linear plot. I hope I made it clear now.. Would it be possible to give me some hint please? Thank you in advance, Sue On Sun, Dec 18, 2011 a

[Tutor] unsupported operand

2011-12-17 Thread stm atoc
Hi there, I would like to define concentration array based on time and then plot it. Here is the profile: from pylab import * import numpy import array import math tmax=1*3600 t = 1. outfile=open('ourtest_out.list','w') N = 100 Conc1 = arange([0, tmax]) outfile.close() lw = 2.0 #linewidth d

Re: [Tutor] New plot over the old graph

2011-12-01 Thread stm atoc
I appreciated for the accurate response. I used step by step and it is running now. Thank you very much for your advice and guidance, Sue On Thu, Dec 1, 2011 at 10:01 PM, Andreas Perstinger wrote: > On 2011-12-01 19:20, stm atoc wrote: >> >> Thanks for all information/webs

Re: [Tutor] New plot over the old graph

2011-12-01 Thread stm atoc
least with my searching. Would it be possible helping me about? I am new at python and a beginner and I do have learn while I work with it... Thanks in advance, Sue On Thu, Dec 1, 2011 at 6:25 PM, Andreas Perstinger wrote: > On 2011-12-01 14:30, stm atoc wrote: >> >> With

[Tutor] New plot over the old graph

2011-12-01 Thread stm atoc
Hi there, I would like to make a new plot with new data over the old one (with old/first data) in the same condition (shape, dimensions) for comparison and analysis data. With your help, I have a good script from the previous discussion: ** from pylab import * import numpy import ma

Re: [Tutor] plotting in python

2011-12-01 Thread stm atoc
you so much, Sue On Thu, Dec 1, 2011 at 12:39 PM, Andreas Perstinger wrote: > [Still top-posting :-( ] > > > On 2011-12-01 11:13, stm atoc wrote: >> >> Well, I did also change the line in the python script to this: >> >> plot(Conc[0],z,'r-',label=

Re: [Tutor] plotting in python

2011-12-01 Thread stm atoc
The output of the print len(Conc[0]), len(z) is 100 3600. Now I changed Conc[0] to Conc[1], and the output is: 100 100 But still I need to see from Concentration from 0. On Thu, Dec 1, 2011 at 5:26 AM, Asokan Pichai wrote: > On Thu, Dec 1, 2011 at 2:38 AM, stm atoc wrote: >> Hi there

Re: [Tutor] Do loop in Python

2011-11-30 Thread stm atoc
it seems it works. it should have a trend to reducing values... - Show quoted text - On Tue, Nov 29, 2011 at 2:00 PM, Steven D'Aprano wrote: > stm atoc wrote: >> >> Thank you so much for your reply. It was very helpful information and >> I used it in order to improve

[Tutor] plotting in python

2011-11-30 Thread stm atoc
Hi there, I have a question regarding plotting with Python. I have the following python script: # coding: utf-8 from pylab import * import numpy filename='ourtest_out.list' fh=open(filename) line=fh.readline() fh.close z=array([ float(val) for val in line.split()[1:] ]) a = numpy.loadtxt(fi

Re: [Tutor] Do loop in Python

2011-11-29 Thread stm atoc
Thank you so much. This script and all information was totally helpful and actually helped me for the next step of my work as well. Have a great time. Sue On Fri, Nov 25, 2011 at 10:44 PM, Andreas Perstinger wrote: > On 2011-11-25 14:46, stm atoc wrote: >> >> Here is the new

Re: [Tutor] Do loop in Python

2011-11-25 Thread stm atoc
end(z) nuh.append(0.001 * exp(-0.005*(z+200.0))*dz) #turbulence diffusivity m**2/s nu.append(num + nuh[z]) but still no sign of 'nu changes' over 'z'! So, would it be possible to check that again? Thanks, Sue On Fri, Nov 25, 2011 at 12:36 PM, Steven D'Apran

Re: [Tutor] Do loop in Python

2011-11-25 Thread stm atoc
27;nu' for each step of z (as height). On the other hand, for each step, 'h' is supposed to be subtracted from 'z' (like: 200-10, 190-10...) as well, at least 10 times (which was trying to be defined as N)! I hope this is somehow clear Thanks in advance, Sue - Show q

[Tutor] Do loop in Python

2011-11-25 Thread stm atoc
Hi there, I am a new python user. I have a question regarding do loop. This is a simple program that I have written: - N=10 h=10.0 # [micrometer] z=-200.0 # [micrometer] num = 0.05 #m**2/s dz = 1.0 nuh=[] tmax=3600 dt=20. nu=[]height = arange(z*dz,0,dz) outfile=open('nu.da