On Fri, Mar 16, 2012 at 2:52 PM, Boris Vladimir Comi <gle...@comunidad.unam.mx> wrote: > Hello, I am writing to request your help in the realization of a script. I > am new to this and I'm just learning the wonderful world of python and this > has made me a little difficult. > > Briefly I commented what I intend to do: > > I detect a class of atmospheric phenomena known as Mesoscale Convective > Systems (MCS). To accomplish this, I have a database of satellites by 2004, > every half hour. These data were downloaded from the server: > unidata2.ssec.wisc.edu and data format is: 200404271515.goes12ir (Area > Format) > > The first thing to do is detect a convective system in the satellite data, > using the following characteristics: > > MCS Criteria > > Minimum Size: Continuous cold cloud shield (TIR <219 K and Must Have an > area> 34000 km ²) > > Duration: Definition of Minimum Size Must Be A Period of Exceed or ≥ 3h > > Initiation: Time When the Minimum Size is met > > Maximum Extention: Time when the continuous cloud shield (TIR <219 K) is at > its maximum size > > Termination: When the Time Minimun Size is not satisfied > > where: Temperature Infrarred is TIR > > To achieve this, first I created a script in python to identify a MCS in my > database (script in attached) > > The script is run from a linux terminal ($ python TIR.py), to run it I get > the following error: > > > File "/home/mcidasv/JYTHON/TIR.py", line 22 > count = count + 1; > ^ > SyntaxError: invalid syntax > > If anyone can help me with this script or any idea you suggest to improve > it, I would greatly appreciate. > > > > Boris Vladimir Comi Gonzalez > Universidad Nacional Autónoma de México > Grupo de Tormentas Convecivas > > _______________________________________________ > Tutor maillist - Tutor@python.org > To unsubscribe or change subscription options: > http://mail.python.org/mailman/listinfo/tutor >
Here is your code: for i in xrange(ad.getLines()): for j in xrange(ad.getElements()): if (data[0][i][j] > 199.5 and (data[0][i][j] < 200.5 count = count + 1; print "For file",name," count = ",count Notice that this line is broken: if (data[0][i][j] > 199.5 and (data[0][i][j] < 200.5 change it to this and try running again: if (data[0][i][j] > 199.5) and (data[0][i][j] < 200.5) -- Joel Goldstick _______________________________________________ Tutor maillist - Tutor@python.org To unsubscribe or change subscription options: http://mail.python.org/mailman/listinfo/tutor