Ground Up, 2D mesoscale model, Sucess!

Joined
Jun 19, 2005
Messages
884
Location
New Mexico
So, since the fall, a little project of mine is to figure out from the ground up how to build a mesoscale model. After many failed attempts, I finally can report success. The model I have built is a 2D, non-hydrostatic, fully compressible model, which includes water vapor.

The key, at least for me, is that a mesoscale model needs to built around some assumed stable condition. This inhibits a sort of gravity wave that would be produced due to the fact of finite resolution. While there are still a few kinks and unknowns, the model remains stable for several hours. Another key is that computational diffusion is required, but the amount of computational diffusion which is needed is still somewhat unknown. Too little it goes hay wire, too much it goes hay wire. So, I have built the model with the amount of computational diffusion that I can select.

Below is a run from the model, starting with an area of higher potential temperate. Within this region I have also included some water vapor. There are no initial winds. Note, the clouds are a very rough approximation based on the amount of condensation.

Clouds at 2 mins...
02mincloudswj9.jpg


Clouds at 15 mins
16mincloudseh3.jpg


Clouds at 60 mins
60mincloudsmi7.jpg


Initial Potential Temperature
thetastartlm3.jpg


2 min Potential Temperature
theta02minsdt0.jpg


15 min Potential Temperature
theta15minsiu7.jpg


60 min Potential Temperature
theta60minscd0.jpg

The values on the axes are kilometers... Ground Elevation ~1km.

Note the cloud that forms dies off, partly because there is no further source of moisture.
 
Last edited by a moderator:
Just a quick question? How did you create your model? Through a site, etc. please provide information. Thanks.
 
Robert have you considered writing about a series of blogs or articles about what you've done?

This is something I'm very curious about and really want to play with in the future when I get a bit more free time and motivation.

I know before you were using matlab, it looks like this is still the case? How is the performance?
 
I'm still using MATLAB, mainly because I am used to programming in it. Although, I have heard that MATLAB can be slower when compared to C++ (because I believe matlab is an interpretive language), I did some searching, and I'm not sure that you would get much better performance (maybe some one has a more quantitative answer here?). For the km by km grid with 25X25 points an hour simulation takes about a minute with .1s steps. Although, the time stepping can be modified. I just got it working so I'm not sure how much faster I could run it (I may want to slow it). For a reference I found Mesoscale Meteorological Modeling by Pielke to be a big help.

Oh, and I don't plan on writing about it in a blog... I'm not really a blog person. If anyone has any questions though, I'd be happy to answer them.
 
I'm still using MATLAB, mainly because I am used to programming in it. Although, I have heard that MATLAB can be slower when compared to C++ (because I believe matlab is an interpretive language), I did some searching, and I'm not sure that you would get much better performance (maybe some one has a more quantitative answer here?). For the km by km grid with 25X25 points an hour simulation takes about a minute with .1s steps. Although, the time stepping can be modified. I just got it working so I'm not sure how much faster I could run it (I may want to slow it). For a reference I found Mesoscale Meteorological Modeling by Pielke to be a big help.

Oh, and I don't plan on writing about it in a blog... I'm not really a blog person. If anyone has any questions though, I'd be happy to answer them.

Out of the products available, which one did you select?

http://www.mathworks.com/products/

And can you link model data to your model and it interpret the data available?
 
I use matlab 2007b student version I believe. I'm not using any toolboxes for the model, although the curve fitting is toolbox they have is quite nice. At school I can use 2008b.

I haven't done anything yet to quickly allow for the assimilation of other data into mine. It will still be a while if ever, before it gets there. I don't really know if I will since it is only 2d. I'm doing 2d because if it was 3d the processing time is so much longer. Note there are models online that you can download, like MM5. I have no experiance in how to use it, but I believe anyone can put it on their computer.
 
Hello Robert,

I like it very much! Before months I wrote an simple advection model. It calculated the 24h temperature forecast in few minutes. Well, the resolution was 2x2 and the code was on 10 lines. The acutual temperature on the next day was close to it:) You can look at my page for weather modeling software for Windows. Do you provide the code of your model. Is it free for use or it is under some license?

Greetings,

Angel
 
Since posting, I have been modifying the code to meet a simulation I would like to run. It has uncovered a few instabilities which are being corrected. Once I have finished I will "clean" the code up (like make it understandable to people besides me by making some comments in the code). Once this is done I would be happy to share the files, with an understanding that any resulting commercial use would require a new licensing agreement. This should take about 2 weeks. Note, I'm new to giving code out since I usually make code for my work, that no one besides me uses.
 
Sorry to those who have wanted to play with the code. I'm still trying to get it to do what I want... which in the process uncovers stability problems, like a forced boundary can amplify pressure waves at the corners. I've been working with a boundary that forces winds. This is because I'm trying to get an anvil and all I get is the shape below.

Top panel is clouds, lower is east west winds. No initial converagance in the winds, temperature gradient developed it...

screengrab.jpg


it's neat though I find the lowering in the storm occurs in the colder portion of the atmosphere.
 
Last edited by a moderator:
I'm now in the process of commenting the code and putting it into a form for some general basic use. What's been taking me so long is I have been playing with a forced boundary... this created some problems due to the computational diffusion required with the resolution I used. I'm currently working on commenting a version of the code that has no forced boundaries (except that nothing can enter or leave). However, people can edit it freely once I give out, and try to deal with the hassles of forced boundary conditions if they wish.
 
Here is a link to a video of a downburst simulation. The colors represent potential temperature. Y axis is up down, X axis is east west....

Hope you guys like it...

http://www.youtube.com/watch?v=YErU7RJx9oE

Working of improving the differencing scheme, so that some finer effects can be seen.

Also if you know of any good journal articles of downburst simulations or observations, I would be interested to see them.
 
Last edited by a moderator:
Back
Top