Google Reader is dead!

Google has finally pulled the plug on the Reader. It is time to move on another feed reader. In the meantime,  here is the last screenshot from the Trends page:

More than 45k read items since October 2008. Looking forward 450k more to read :)

Thanks to Mihai Parparita's Python scripts I archived all there is to archive from the Reader.

Happy news reading!!!

Issues with GDL Python module

1-) Importing GDL segfaults in IPython

$ ii # this is alias for ipython --pylab

Python 2.7, IPython 0.14.dev

I1 import GDL

I2 Segmentation fault (core dumped)

however, works fine in Python shell

>>> import GDL

>>> GDL.findgen(5)

>>> GDL.function('findgen',5)

array([ 0.,  1.,  2.,  3.,  4.], dtype=float32)

possibly readline module issue? Interestingly, works fine in Spyder's IPython console.

2-) Source compiled (with --enable-python_module option) GDL Python module fails to work

 >>> import GDL

Traceback (most recent call last):

  File "", line 1, in

ImportError: ./GDL.so: invalid ELF header

possibly, due to wrong lib arch use, during compilation?

3-) GDL's dindgen yields wrong result while accessed through its Python module

>>> GDL.function('dindgen',10)

array([  0.00000000e+000,   1.00000000e+000,   2.00000000e+000,

         3.00000000e+000,   4.00000000e+000,   0.00000000e+000,

         5.92878775e-323,   3.30571061e-313,   6.94380683e-310,

         4.03179200e-313])

Posted this issue on dindgen with GDL python module

Calling IDL, GDL, and FL from Python

Thanks to Anthony Smith's pIDLy module, accessing IDL, and its clones GDL and FL from within Python becomes really easy. 

Launch an IPython session, then typing these example lines to see pidly in action:

I1 import pidly

I2 idl = pidly.IDL()

I3 fl = pidly.IDL('fl', idl_prompt='FL> ')

I4 gdl = pidly.IDL('gdl', idl_prompt='GDL> ')

I5 idl.findgen(10)

O5 array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.], dtype=float32)

I6 fl.findgen(10)

O6 array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.], dtype=float32)

I7 gdl.findgen(10)

O7 array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.], dtype=float32)

I8 np.arange(10, dtype='float32')

O8 array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.], dtype=float32)

Testing these with:

I2 pidly.__version__

O2 '0.2.4+'

I3 np.__version__

O3 '1.8.0.dev-82c0bb8'

$ idl

IDL Version 8.1 (linux x86_64 m64). (c) 2011, ITT Visual Information Solutions

$ gdl    # the version that ships with Fedora 16

GDL - GNU Data Language, Version 0.9.2

$ fl      # latest development snapshot (10/Oct/2012)

Fawlty Language 0.79.18 (linux amd64 m64) Copyright (c) 2000-2012

Example access to an IDL procedure:

I5 idl.pro('cgplot', np.sin(np.linspace(0, 6*np.pi, 1000)))

I6 idl.pro('iplot', np.sin(np.linspace(0, 6*np.pi, 1000)))

For more information, just read the quick usage guide on pidly's web-site, or clone the source from https://github.com/anthonyjsmith/pIDLy

Worldcloud for ccnworks

Here is a worldcloud created for my code repository on:

http://code.google.com/p/ccnworks/

First get the commit message history into a text file (the last commit is r287):

svn log >> test.txt

Then with a bit of VIM trick and manual step, remove unwanted lines:

sort u


Select of all the text, and import into

http://www.wordle.net/create

Tadaaa, here comes the result :)

Indeed, we are data scientists here in atmospheric science!

Alternative wordcloud implementation comes from Andreas Mueller. I saw his A Wordcloud in Python titled posting via Planet SciPy. In order to run his wordcloud.py I needed to change unsigned int[:,:] integral_image -> unsigned long[:,:] integral_image in query_integral_image.pyx to get past 'ValueError: Buffer dtype mismatch, expected 'unsigned int' but got 'unsigned long''.

Then feeding the same log data I get the following image, which is nicer looking than the wordle's image :)

Python videos

# from AMS

Reading, writing and visualizing GRIB and netCDF data with Python (invited) (91st American Meteorological Society Annual Meeting)

Porting radar simulation software to Python: a case study in the benefits of Python (91st American Meteorological Society Annual Meeting)

WxMAP2—a python-based system for model diagnosis and weather forecast tool development (91st American Meteorological Society Annual Meeting)

Using CDAT for model development and analysis of large data sets (invited) (91st American Meteorological Society Annual Meeting)

Community Tools for the Analysis and Visualization of Geoscientific Data (invited) (91st American Meteorological Society Annual Meeting)

DV3D: 3D Climate data visualization using python and VTK (91st American Meteorological Society Annual Meeting)

HootPy: A Console-Based Processing and Visualization Package for Meteorological Models and Observations (91st American Meteorological Society Annual Meeting)

Python Adoption at the Naval Research Laboratory (91st American Meteorological Society Annual Meeting)

# PyCon10 Videos

Introduction to Python Profiling

PyCon 2010:Python's Dusty Corners (#182)
PyCon 2010:Keynote: Relentlessly Pursuing Opportunities With Python, or why the AIs will Spare Us Al
PyCon 2010:The Mighty Dictionary (#55)
PyCon 2010:Keynote: State of CPython
PyCon 2010:Deconstruction of an Object (#62)
The Ring of Python (#189)
Debating 'til Dawn: Topics to keep you up all night (#183)

# PyCon2009 videos

Eggs and Buildout Deployment in Python (Part 001)
Eggs and Buildout Deployment in Python (Part 002)
Eggs and Buildout Deployment in Python (Part 003)

Introduction to Python Profiling
Leveraging ReStructuredText in your Projects
The State of the Python Community: Leading the Python tribe
Keynote: Guido van Rossum
How Python is Developed

An Introduction to Object-Oriented Programming (Part 001)
An Introduction to Object-Oriented Programming (Part 002)
An Introduction to Object-Oriented Programming (Part 003) 

Python 401: Some Advanced Topics (Part 001)
Python 401: Some Advanced Topics (Part 002) 
Python 401: Some Advanced Topics (Part 003)

Python 101 (Part 001)
Python 101 (Part 002)

Python 102 (Part 001)
Python 102 (Part 002)
Python 102 (Part 003)

# various 

sagedays29: Fernando Perez -- Scientific Computing Using Python

SciPy 2010 - Travis Oliphant - Moving Forward from the Last Decade of SciPy

SciPy 2010 - David Beazley - Python Concurrency
Getting Started with Numpy Arrays in Python

PyAtl November 2009: "Python Data Structures" by Brandon Rhodes

Google's Python Class

Python for Scientific and High Performance Computing (1 of 4)

What is Python? by Wesley Chun

Advanced Python or Understanding Python

The Importance of Programming Literacy (PyCon 2007)

Google I/O 2008 - Painless Python Part 1 of 2
Google I/O 2008 - Painless Python Part 2 of 2 

Sage: Unifying Mathematical Software, video part 1, video part 2

A Demonstration of the 'IPython' Interactive Shell

Various open-source videos

What's under the hat? A sneak peek at Fedora 13 by Jesse Keating
Git er done with GIT by Jesse Keating

Keynote: MeeGo: A Free & Standard Linux OS for the Mobile Industry

Google I/O 2009 - The Myth of the Genius Programmer
Google I/O 2009 - Do You Believe in the Users?
Google I/O 2008 - Project Hosting on Google Code
Google I/O 2008 - Open Source Projects and Poisonous People
Google I/O 2008 - Open Source is Magic
What's In It for Me? Benefits from Open Sourcing Code
Learning from StackOverflow.com

Paul Graham - Great Hackers (2004)
Paul Graham - Open Source Software Development (2005)


# resources

http://docs.sfepy.org/euroscipy-tutorials/tutorials/index.html
Euroscipy 2010 tutorials
RCE 48: NumPy
The Pragmatic Bookshelf | List of Tips
Code Like a Pythonista: Idiomatic Python
NumPy MedKit by Stéfan van der Walt (.tar.gz)
PEP 3099 -- Things that will Not Change in Python 3000
Building Skills in Python

NumPy 1.5 Beginner's Guide review

This new NumPy book provides an easy start to NumPy for those who want to churn numbers in Python programming language, whether you don't have much programming background or you are switching from your favorite language. The author of the book uses quite a friendly tone throughout the book. Most of the important aspects of NumPy is well covered with well explained examples. The examples provided are step by step explained, starting from the basic array/matrix creation to more complex tasks like signal analysis and linear algebra related calculations.

To get best out of this book, it is recommended that you would try out the examples and challenge yourself with the exercises given in have a go sections of the book. If you are feeling lazy for some reason, you can get the source code from the book's page. I finished perusing the book in about a week (having a few years of NumPy experience has definitely role in this). For me reading and understanding through stock market related examples were a bit boring, but if you are in to financial business you might enjoy putting NumPy in your toolstack with the help of this book. Notably, besides the basic NumPy beginners chapters, the book extends into basic Matplotlib and SciPy lands. A lot of examples use Matplotlib to create plots and better illustrate the operations at hand (e.g. curve fitting, statistical distributions) It is a bit surprising to see MaskedArray module didn't get any mention in the book. However, with the basics you gained, it should be fairly easy to start experimenting with masked array functions of NumPy.

Overall, if you are looking for a book to get started in NumPy in about 200 pages, you might give this one a chance. It is available in both print and electronic formats. Further on, you can try their advanced matplotlib and and Sage books, if you are willing to enhance your scientific Python skills.

Final Note: I received a review copy from the publisher. Thanks to Packt publishing for their contributions to open-source literature and providing me a free copy of the book.

Virgin Adaları'ndan dönüş

Buradan oraya gidiş oradan buraya gelişten daha kolaydı nedense. Yükseklik farkından olsa gerek herhalde. Taa 2000 üstü metrelerden deniz seviyesine varmak nispeten daha kolayca oldu. O kadar bisikletimi götürmeyi istesem de adaya, bu isteğimi bir türlü gerçekleştiremedim. Halbuki onca zamandır da bisikletin içinde geldiği kutusunu saklamaktayım, olurda bir gün bir adaya giderim, olurda bisikletimi götürmek isterim diye. Ne adada bisiklete binebildim, ne de adanın çevresini koşarak ya da kanoyla kürek çekerek dolaşma hayalim gerçek oldu. Artık bu hayaller başka bir zamana kaldı. Karayipler’de 7000’i aşkın irili ufaklı ada varmış. Bir gün yine denk gelir adalar etrafında bulutları çalışmak isteyen ve bulutlardaki yağış oluşumunu anlayama çalışan birileri.

Bu aşağıdaki ilginç şekilli ada (her ne kadar konunun genelinden biraz uzak olsa da) KK adlı okur ve arkadaşa gidiyor. Gerçi her ne kadar adanın tapusu bende olmasa da, kendisinin adanın yer bilgisini bulması karşılığında bir sonraki sefere görüşmemizde bir şişe “Cruzan rum” garantilediğini bilmesini isterim. KK, ister desen tanıma kullanabilir, isterse resme ait bilgilerinden uçağın o anki yerini tespit etmeye çalışabilir, ya da olmadı eşe dosta sorabilir, tercih tamamen kendisinin.

Bugün ayın 28’i olmuş. 17 gün ne çabukta gelip geçmiş. Arda kalan günlerde gezip görüleni ve çalışılıp işleneni kısaca özetlemek için bir güzel sözsüz şarkıyla eşlik edilmiş bir resim şovu hazırladım. Şarkı yaklaşık 4 dakika 28 saniye sürüyor. Saniyeler güne tuttu, hadi öyleyse izlenen resim sayısı da adada geçen gün sayısı ile şarkının dakikasından türeyip 17x4 = 68 olsun. Tıklayıp keyifli seyirler dilerim (bu kapanışta çok fazla AA ile seyahat etmekten bir etkilenme olsa gerek :)

Virgin Adaları Günlüğü: -3. Gün

Bugün adaya gelişimin 2. haftası.  Bunca zaman geçmesine rağmen ilk kez C-130'de uçuş fırsatını yakaladım.  Çoktan beri aksatılmış günlüğün bu bölümünü yazarken 10. görevin eve dönüşüne yaklaşık 1.5 saat var. İçerisi bolca gürültülü, kulaklık takmaksızın, kulakları tıkamaksızın yolculuk etmek pek mümkün değil. Gerçi bu uçağın amacı yolculuk etmekten ziyade, bulutların içine uçup onların özelliklerini ölçmek. BUZ-T isminden de anlaşılacağı üzere asıl gaye bulutlardaki buzu ölçmek (gerek miktar, gerek şekil, gerekse de yapısal olarak) ve buzun yağış oluşumundaki işlevini daha detaylıca anlamaya çalışmak.

Uçak bulutları ziyaret ederken ayrıcalıklı ön koltuklarda oturmadıkça, dışarıyı görmek pek mümkün değil. Bugünde ileriye bakan kamera arızalı olduğu için ekranlara kamera resimleri yansımadı --uçağın altına bakan kameradan gelenler hariç.  Ziyareti söz konusu olan bulutlar epeyce gelişmiş olduklarından içerlerine doğru dikey hava hareketi epey güçleniyor. Bu sebeple emniyet (bu kelime önemli) kemerlerini açıp dolaşmak pek akıl karı bir iş değil. Bana göre “Havaların Kralı” ile bulutları ziyaret daha eğlenceliydi çünkü içerdeki herkes dışarıdaki durumu rahatça görebiliyordu, ama bu uçakta bulutları anlamamıza yardımcı olan çok daha fazla sayıda araç gereç mevcut.  Mesela bu ucuşta kaydedilen veri sayısı yaklaşık olarak 500 ve 6 saatten uzun süren yolculuğun ardından her saniyede bir kez örneklenmiş verilerden oluşan dosyanın boyutu yaklaşık 130MB civarında.

Sağ yanımda iki tane büyük hava tüpü, önümde bulut çekirdeği sayma ve ölçme cihazı, sol yanımda önündeki cihazlarla ve bilgisayarla meşgul biri, dişlerimin arasında az önce yediğim elmadan sadece geriye kalan sapı, başımda içinden bulut alakalı konuşmalar geçen kocaman bir kulaklık ve, ve aklımda bir türlü çıkmak bilmeyen bir güzel S, daha hala 1 saatlik yol var geri dönüşümüze...

Virgin Adaları Günlüğü - 2. Gün

Dün pek bir değişik etkinliği olmayan bir gündü. Gözlükler yine buğulandı, garip kuşlar yine ötüp durdu, deniz sakinliğini korudu, bilgisayar başında yapılan çalışmalar devam etti, 5. araştırma ucuşun ardından 36 kişi bir odada bir araya gelip uçuş değerlendirmesinde bulundu ama gece yağmur yoktu.  Günün önemli olaylarından biri uçuşlara katılabilmem için gerekli kısa güvenlik eğitimini tamamlamış olmamdı. Öğlen bunu yazarken gürültüsü bol, içinde ve dışında bir dünya ölçüm araç ve gereci olduğundan bahsettiğim uçağın içindeyiz. Bugün uçuş yok ama içerde olağan koşuşturmacalar sürmekte, bulutlarında istemesiyle bir-iki güne uçuşa gidiş yakındır.

Virgin Islands Diary - Day 1

Today I have once again realized that, I can't think the same way in Turkish as I do in English and vice versa. Yesterday night, and this morning I tried to translate my first diary posting to English, but to no avail. Neither I was able to list the words as nicely I did in Turkish nor my sleepy eyes were able to fight more. So I have decided to alternate these island blogs in between languages. Writing Turkish for even days and English for odd days unlike I originally planned to write dual language for all days. Unlike me, Erden Eruç can handle dual language blogging (as well as circumnavigating the earth) very well with using the beautiful intricacies of both languages. Three years of staying in an English speaking country apparently is not enough for me, still lots of grammar mistakes, still asking toothbrush instead of toothpaste and filling in arrival dates instead of departure dates in forms. Of course, fatigue from long travel and not getting proper sleep, high humidity and trying to focus on writing and improving an important research related text and the rash of trying to response for those who wonders about me using this turtle speed internet connection. Tell me please, where is my grammar fairy?

Day 1 goes fast indeed. Waking up 7:30 and trying to cope with the high humidity and heat aren't as easy as it looks. Think of it like living in your bathroom while your hot shower is on. Luckily these rooms have air conditioning to make life easier. It feels like changing climates from going inside and out and outside in. I don’t remember glasses getting fogged in summer, like you step inside a warmer place from a colder place. It was a busy day today. Finally we have created an almost page long draft that describes our efforts about numerical cloud parcel model experiments that will go in a report. It took many hours indeed to finalize those few paragraphs. Each sentence triggered a new question and each question lead to a new thought provoking discussion. In our case, even jokes turned into deep intellectual conversations and sometimes deep intellectual conversations into jokes. Anyone out can tell me what "exact non-equilibrium thermodynamics" really means? :) or how Thales' deductive reasoning could play role in explaining your thought process while delving deep into the explanation of internal cloud model interactions? Probably floating on the ocean while sipping a cold drink might cause similar triggers in one's mind. Still yet to try if this claim is true or not. Students of science occupation tend to verify things, don't they?

Even though I wasn’t able to fly today, I had chance to connect to the C-130 data system to see the measurements in real-time. This is great since without being in the aircraft you can observe how things are looking out there, except experiencing some cloud bumpiness and sharp research maneuvers. Today they performed the 4th research flight. The previous 3 research flights have already been processed and shared. Meaning I have more data to study and understand. I like variety yes, and I have atmospheric dataset collections from many different locations of the earth. When I want to go the South Pacific or Arabia or even to the central US (actually, I just go there instead of all these buzz :D), I just browse into my folders and have a fresh look at what those probes tell me, what the radar reveal, and how their imagers capture clouds like. More research is on the way, so better tigthen my mind.

Those flurries have visited us again at around the same time like they did last night. Free nightly washing service. No wonder this tropical island is so green :)

Virgin Adaları Günlüğü - 0. Gün

Bir kez daha SciPy konferansına gitmek isterken, yolum başka bir yöne düşüyor.  Gerçi bu kez gidilen yer okyanus aşırı değil, aksine okyanus içinde bir yerde: Virgin adalarında.  Bir dünya dolusu atmosferik ölçüm cihazı bulunan bir uçakla buzun yağış oluşumundaki etkisini araştırmak daha cazip geldi sanırım bu sefer :) Hem C-130 ile ilk defa uçacak olmak hem de birkaç zamandır okuduğum makalelerin yazarları ile tanışacak olmak bu işi gerçekten çok heyecan verici kılıyor.

Adaya ulaşması epey zaman aldı. İki uçak ve bir otobüste geçen, beklemelerle birlikte 15-16 saat süren bir yolculuk oldu.  Sonrasında uykulu gözlerle gecenin köründe yazılan bu yazı ve yazılması gereken bir paragrafı bitirmenin telaşı.  Keşke sadece bitmeyenler bunlar olsa: geride SciPy11 zamanında tamamlamayı planladığım basit bir bulut modeli, yolculuk sırasında bitiririm dediğim “Termodinamik Yasalarına Çok Kısa Bir Giriş” isimli sesli kitap (bilin bakalım 0. gün ismi nereden çıktı), ve diferensiyal denklemlerin hızlı bir tekrarını yapmamı sağlayan görüntülü ders notlarının son birkaç bölümü hala beni beklemekte.

Öyle ya da böyle bugünden itibaren ikinci saha araştırmama merhaba diyorum. Daha önce Arabistan’da bulutlara uçmuştuk, simdi sıra Karayipler’de gezinen bulutlarda.  Kendime ve takım arkadaşlarıma ölçüme müsait, bol bulutlu, tropik fırtınasız günler diliyorum.  Beni bekleyen ya da bir sebepten beni bekleyemeyen herkese de çok selam şimdiden.  Gecenin bu geç saatinde, dışarısı 25 °C’ye yakın sıcak iken, birden başlayan sağanak yağmur buradaki zamanın çok hareketli geçeceğini haber ediyor sanki.  Biz burada BUZ-T dolayısıyla meşgulken sizlere de buzu eksik olmayan günler diliyorum :)

Records number of citations

It has been a little more than a year since I mentioned about "Record number of equations" in a book about cloud microphysics. This time, the record holder is in a different category, that is the number of references used in a journal article. In her 2003 review article*, Prof. Finlayson-Pitts references a total of 367 articles.

If you are not much in the chemical aspects of the sea salt aerosol but want to learn more about the physical and measurement sides of the subject then you can consult to your one-of-the-high-citing books (over 1800 references) from Lewis and Schwarts (2005)**.

With number of my local article collection reaching over 400 and citation listings reaching up to upper 2 and sometimes 3 digits in most of them, scientific literacy has a great place in the life of a PhD student in atmospheric sciences.

References:
*Finlayson-Pitts, B. J., 2003: The tropospheric chemistry of sea salt: A molecular-level view of the chemistry of NaCl and NaBr. Chem. Rev., 103, 4801-4822.

**Lewis, E. R., and S. E. Schwartz, 2005: Sea Salt Aerosol Production: Mechanisms, Methods, Measurements, and Models. American Geophysical Union, Washington, DC.

Precipitation Diagram

from: http://weatherfaqs.org.uk/node/160

Aerosol effects figure in IPCC Climate Change 2007 report

Figure 2.10. Schematic diagram showing the various radiative mechanisms associated with cloud effects that have been identified as significant in relation to aerosols (modified from Haywood and Boucher, 2000). The small black dots represent aerosol particles; the larger open circles cloud droplets. Straight lines represent the incident and reflected solar radiation, and wavy lines represent terrestrial radiation. The filled white circles indicate cloud droplet number concentration (CDNC). The unperturbed cloud contains larger cloud drops as only natural aerosols are available as cloud condensation nuclei, while the perturbed cloud contains a greater number of smaller cloud drops as both natural and anthropogenic aerosols are available as cloud condensation nuclei (CCN). The vertical grey dashes represent rainfall, and LWC refers to the liquid water content.

IPCC: Climate Change 2007: The Physical Science Basis, in: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge, UK and New York, NY, USA, 2007.

Seminar on clouds

Yesterday, I gave a seminar talk to our department faculty and students, titled as "An Aerosol - Cloud Condensation Nuclei and Cloud Droplet Concentration Closure Study" that is summarizing a majority of my thesis related work on aerosol and cloud micro-physical research, starting and finishing with the presentation slides shown above and at the very end of this post.

One of the most exciting parts of the talk was after I finished demonstrating my research results and showed a few of the technicalities of the research. First of all, thanks go to Dr. Jefferson R. Snider for providing his IDL-based parcel model code and being patient with me while I was re-writing the model in Python. Both versions of the code could be accessed through ccnworks by browsing the source tab, and looking under the thesis folder or alternatively, directly contacting to Dr. Snider to get the latest updated version of the model.

When I was working on the model conversion, at one point I had stuck getting similar outputs from the two languages using the same initial conditions, and with the almost identically progressing code. (Some might object on the definitions of two conditions' being "the same" and "almost identical" in the land of high-precision floating point arithmetic.) Continuing in my numerical instability problem, I had created a simple animation to better demonstrate the issue I was experiencing. As some eyes might easily catch, the droplets in parcel model reaching to so bigger sizes comparing to the results produced by IDL, solely based on condensational growth theory. It took a while to find and correct the variable that was updating its original reference while it was supposed to be operating on a copy of the content instead. This issue and a few other discrepancies were fixed in around May 2010. As a personal important todo note at this point is to restructure the parcel model code so that it behaves as an external library, thus making all the important thermodynamical and cloud micro-physics related functions greatly re-usable. This is when I figure passing dictionaries (or changing the structure of the code) to outer scope of the original file that they are declared in.

The final technicality demonstration comes from the work of Dobashi et. al., (2008) (Their "Feedback Control of Cumuliform Cloud Formation based on Computational Fluid Dynamics" titled work.) Those are by far the most realistic cloud appearances I have ever seen in a physical cloud formation processes applied simulation study. Although cloud modellers are usually interested in studying temporal and spatial statistical properties of clouds rather than synthesizing their realistic shapes and appearances, such a similar work would have a great educational use in visualising and thus better demonstrating the physical processes governing the cloud formation and further precipitation development. I can say that the authors should have prototyped their simulations using a C-like language in order to make low level access to CPU/GPU because of the high computation demands, but there are fast computation techniques available in Python that would make this undertaking possible to be realized be it simulating a research flight or playing in between sequences of a simulation to probe important properties of the formation.

For the curious, I share high-quality PDF version of the slides. Note that I had to take out a couple pages from this original presentation in order to fit 45 minutes allocated presentation time, and it was unfortunate that one of the expected audience was not in the classroom to see a few indistinguishably minute research unrelated points. All analyses and plots in my slides are performed/created using Python and its scientific tool-stack with some additional annotation help from OpenOffice on the slides using DejaVu-Sans font; for the first time with full consistency. There is one exception to the pure Python code, that is when the Python based model calculations taking up five to ten times more than the original IDL-based code, and that is where exactly Cython magic were sought and applied to boost million-times executed functions. Here, I greet the community once more, for making the scientific computation freely, publicly, and easily available to everyone without putting any restrictions in both code and verbal discussions level as how this core work should really be in its very essence within the science endeavor.

I also salute Matthew Turk and his work for that fantastic presentation style even though I couldn't copy his great architectural design throughout my slides. Yet not him, I have other masters around me for letting me embrace Tufte's sayings :)

My work on clouds is on-going and leading me to a PhD level research after successful completion of my Masters thesis within the next month. Feel free to say hi or ask more clouds if you share any research similarities.

Now, it is time for some abstraction: Youth Uprising - Rarefaction

Please keep up imagination and creativity flames are burning all along...

Vim magic2 - search and replace in visual selection

Since I thought making generic analysis and plotting routines would have been harder for 8-9 different source airborne data, I am showing some pure uniqueness and laziness of variable name selections on the given screenshot. Like in many cases VIM is here to my rescue :) First we switch to line selection mode by Shift-V combination and highlight the area of interest and then hitting ":" and continuing typing by s/06/08/g in the command insertion area. The magic is done --all the "06" instances are replaced by "08"s readying the code for our April08 case.

Vim magic: Synchronized scrolling

As shown in this screen capture, I have two minor modifications to my original closure-fitting code. The original code runs using the lab PCASP size ranges and excluding the first PCASP channels data, "zair-fit" version uses the field calibrated PCASP size bins, and still excluding the first two channels. The final file "zair-fit2" runs the code with the size bins and including the pseudo-problematic first two channel PCASP concentration data.

VIM provides me a very practical way to synchronized scrolling of all open windows so that I can easily compare the results from those well-formatted table structure outputs. Solution is simple: execute the following line in each window as shown in this blog entry.

:set scrollbind

Damlalar her yerde

Evet, hepsini içmek serbest hatta kendi sıvınızda dahil. Olur ya kurak bir yerde mahsur kalarsınız diye Les Stroud’un Kalahari çölündeki bölümünü izleyebilirsiniz ya da burda kalıp okumaya devam ederek kendinize banyo yaparken soracak yeni sorular edinebilirsiniz.

Duş almaktan söz ediyoruz ya, ilk soru şöyle belirdi: sıcak bir duşta, pek bir havalandırması olmayan bir odanın içinde, yukarıya lamba tarafına doğru bakınca havada asılı danseden bir sürü parçacık göze çarpıyor. Gerçektende havada asılı duran parçacıklar, üfleyince bile uçuşuveriyorlar. Soruyorum, bu parçacıklar ne ölçüde aktifleşmiş? İçerideki ortam hangi oranda neme doymuştur? Sonraki soru bir şekilde banyo duvarına ya da perdesine sıçramış parçacıkların büyümesini, süzülüşünü ve diğerlerine karışıp akıp kaybolmasını izlerken aklıma geliyor. Çok sayıda küçük damlacık ve az sayıda büyük parçacık. Gamma ve Marshall-Pallmer dağılımları ne kadar doğru tanımlamakta bu parçacık dizilimini? (Diğer dağılımlara Kolmogorov–Smirnov testi yapmak isteyen?) Son sorularsa en ilginci. Bu kez damlacıkların üstümden sıçramasını izlerken düşünüyorum ve seyre dalıyorum. Çarpan yüzeyden seken damlacıklar (bu durumda benim omzum) ne şekilde parçalanıyor ve nasıl bir yörünge takip ediyor, ne zaman sonra çiseleyerek yere doğru düşüyor ve ne kadarı buharlaşıverip ya da yan duvarda kümeleşen parçacıklara yapışıp etkisiz hale geliyor?

Ya bu kapta yoğunlaşan parçacıklara ne demeli? Nedir gerçekte onları bu şekilde dengede tutan? Nasıl bir etkileşimle öyle değişik şekiller alabilmekteler ve ne kadar büyümeleri gerekir tutundukları zeminden aşağıya kaymaya başlayıncaya kadar.

Bazıları Percy W. Bridgman’ın “Eğer ölçemiyorsak büyük ihtimalle ne demek istediğimizi bilmiyoruz” sözünü fazlasıyla dikkate alaraktan ya da Meinrat O. Andreae’nin bir konuşmasında alıntıladığı gibi “Gözlemleyerek öğrenmeyecek şey yoktur” sözünün büyüsüne fazlasıyla kendilerini kaptırdıklarından, yukarıdaki soruların bazılarının yanıtlarına ulaşmışlar gibi. Tabii bu bazıları dediklerim akışkanlar mekaniği ve yağış fiziği konularında ileri seviyede bilgisi olup dahası bu minik dünyalardaki hareketlilikleri çok yüksek hızlarda izleyebilecek araç-gerece sahip kişilerdir.

Böyle bir uğraşın sonucunda ortaya çıkanlar ise sadece bir su damlacığını analiz ederek yağmur damlalarının genel dağılımı hakkında bir görüye ulaşılabileceğini gösteriyor. (“Bilinen bir konu üstüne çalışarak, yeni bir olguyu keşfetme.” Bu alıntının nereden geldiğini çözemediyseniz sizi James Burke’nin “Evrenin değistiği gün” belgeselinin ilk bölümüne alalım.) Devam edersek, düşen yağmur parçacıklarının parçalanmaya başlamadan önce teorik olarak ulaşabilecekleri çapın en fazla 6-7 mm dolayında olabileceği (bazı istisnaları ile birlikte). En ilginci ise, her ne kadar Galileo zamanından beri bilinen ve yağış ile ilgili kaynaklarda damlaların belli bir hızın üstüne gidemeyeceği bilgisinin (en azından benim böyle bilişim) yerel olarak doğruluğunu yitirişi. Parçalanış esnasında anne damladan kopan parçacıkların tekrar dengeye ulaşıncaya kadar limit hızdan aşırı hızlarda hareket ediyor oluşu. Detaylarını merak edenler şu makaleyle devam edebilirler,  ya da şu bağlantıdan görüntülenmiş bir parçacığın dansını seyredebilirler.

Buradan karşısı daha mavi

    Yok yok, buranın karşısı değil. Bu taraf geçen haftanın macerası idi. Dağın maceracısı çokmuş: kimi normal yoldan çıkıp, bisikletle patikalardan inermiş, kimileri de dağdan aşağıya atıyor kendini paraşüte güvenip :) 

    Evet işte o tepelerin ardı deniz. Şehir bir dağ ile bir tepe arasında alabildiğine uzanmış. Ne güzel bir gün öncesinden karşı taraf kendini açıkça gösteriyordu, fakat oradan dağa bakınca pek açık bir görüntü almak mümkün olmadı. Kirlilikten mi desem, buluttan mı yoksa bulut olamamış parçacıklardan mı desem dağ kendini göstermedi ertesi günü. Yandan uzanan el bana gidilecek tarafları gösteriyor: Gündoğdu, oradan aşağı Kurşunlu, sonra yan yan Mudanya.

(Son cümle tamamlanamadan bir bölünme oldu, annem tarafından ekmek ve yufka almaya gitmem gerekçesiyle, yazı bekler patatesler kızarıyor neticesinde :) Ben yolda, acaba diyorum şimdi Texas'daki SciPy10 eğlencesinde olmak nasıl olurdu diye söyleniyorum kendi kendime, ya da Oregon'daki Bulut Fiziği konferansında geziniyor olup, poster sunmak, makalelerde ve kitaplarda adlarını gördüğümüz kişilerle karşılaşıp iki laf yapmak hiç fena olmazdı.)

    Karşısı yalnızca daha mavi değil, ayrıca daha da yeşil. Yol boyu isteyene zeytinlik, sahibine incirlik, yoldan geçenine çamlıklar dizili. Kimi dut vardı, kimi erik, şeftali, armut yemek isteyenine. Ara ara çeşmelerde vardı yol boyunda, ama su içmek isteyeni dağa çıkmalı...

    Bu resmi yukarıdakinin önüne koyunca körfeze baktığımız anlaşılıyordur sanırım. Aslında tepeden denize inince kanoya tekneye binip karşı tarafa kürek çekmek vardı. Belki ileride; tepeler yeşilini, deniz de mavisini korumaya devam ederse...

    Her inişin bir çıkışı varmış. Gerçi bisikletler güzel güzel poz verse de, elde taşındılar epeyce. Özellikle bu çıkış dağa çıkışın çok benzeri, yalnız bir farkla; burası 5-10 dakikada tükenip yine aşağıya akıyor, fakat dağ yolunda buna benzer en az 100 tanesi aşılmayı bekleyen. 

    Bu sefer N900'ün şarjı tam doluydu yola çıkmadan. Bizim gezi yaklaşık 10 saat sürsede hem resimler hem de yol koordinatları bu kez bir aksiliğe uğramadılar. Bu yukarıdaki Runsaturday sitesinden; GPS olunca gezi daha başka bir hal oluyor. Gerçi yol boyu halinden memnun olmayanlar olsa da, bol iniş-çıkışlı ve bol manzaralı, yine son dönüş kısmı geceye kalan eğlenceli bir bisiklet turu oldu.