hiccups
program and work-dipom is basically no different than running little problems to solve. Plays little more time. This month was marked by two larger ones. The first problem was to make it very difficult to locate. Fortunately
now the whole code in a version control system called Subversion. Many people use open source Subversion Kham
so that multiple people can work and do not come in your way. Very practical.
The second problem which I have now is not technical but mathematical programming. This point is probably in every major work ahead of time? Was all for nothing I have done? No! At least that is my opinion.
heard It is so simple. I must find out where a particle was previously a time-step. I know where it is now, and I know the 3D velocity field. So it's all there, right? But just the beginning of the story. , it happens no chemical transformation and no diffusion. But that is precisely the problem.
You think you could just carry the particles that fall on the ice. Just at the temperature formula which omitted? But no, there is no diffusion, there are sharp transitions from one property to another. For example, if two glaciers flow together, that's an ice and thus the particle is 2000 years old and the other 40000th At one point I then got a transfer of 2,000 to 4,000 Exactly at the interface I have on a jump. This leads to some problems. At the temperature there is something not because it is no matter how old the ice, the temperature there has not cracked. To detect the gene-history of the particles, such jumps should be possible. The origin and the time where it has landed on the ice must be maintained. It is the simplest formula, like this: The change is 0! The problem was solved until 2002 was the first time for ice sheet. The idea was to make it so so in 2002. See the picture below. It basically
two ways how to implement a model. It shares its entire area so as to Greenland in a grid, say 20x20 km and 100 in the high layers. In the Euler method one has a fixed grid and looking at what happened. In the Lagrange model you look at a particle and it follows all the time. In this case, as it were, then bends the bars. Both options have advantages and disadvantages. My
ice sheet has a fixed grid. So now it comes: there is the semi-Lagrangian method. This is a mixture of both. One uses a fixed grid and tracked in each time step, only the particle that hits a lattice point as the endpoint. Because change is not with me so the property should it enough if I know where the particles come and how the property was there the previous time step.
The difficulty is now to know where it was PRIOR particles. In a 3D grid with a different velocity at every point within a grid cell! And an ever-changing speed. Yes when that does not work then the whole method does not, or will the inaccurate, so inaccurate and everything becomes less accurate with each time step. If we accurately compute time just 400,000 years then add up the deviations.
is unfortunately nothing in the article about accurate extreeeem only a simple formula: the place where it was before the same is the place where it's gone now minus the back down ....
So I searched more articles the semi-lagrange use and after a day I also found some where it is described in more detail, with all the numerical odds and ends you have to consider too. So I'm confident it will work and I will meet again on Monday with my supervisor.
We think that is more complex and the more complicated it is. Make a simple task at first is evolving into an interesting and complex problem. We lost very quickly to details. To every detail, there are dozens of articles and sometimes you just think everything is too complicated and it is no longer the easiest things safe. So it was good that I've written here once so I know what itself the original objective was.
but then starts to really, because then you can make virtual ice cores and compare them with real! About isotopes and for which one uses it I'll write something later.
Well, now, a slightly longer post but then I want if I am then in a week in Graz not always explain what I make. I'm looking forward to daham, so until then.
PS In August I'll go to Iceland! The hike, etc., not from the AWI.
0 comments:
Post a Comment