On Fri, Mar 16, 2012 at 2:58 PM, Joel Goldstick <joel.goldst...@gmail.com> wrote: > 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 - tu...@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
Sorry, also put : at end of if statements -- Joel Goldstick _______________________________________________ Image-SIG maillist - Image-SIG@python.org http://mail.python.org/mailman/listinfo/image-sig
