1. new publications
·
paper submitted to PoP by Michele
Romanelli et al on the microstability analysis of FTU pellet injected
plasmas to be read here.
· Talk on a general overview of all the stabilizing impact of few key parameters (Zeff, s, q, Ti/Te, a, density peaking, fraction of trapped) given by myself at the European-TTF Topic-meeting on electron transport in Milano Feb. 18-20 2004.
2. quite a few upgrades in the interfacing tools
· plotkize.m is now ploting a figure to test the validity of the balooning representation in 2 ways:
- by comparing the mode widths with the distance between resonating surfaces, knowing that the mode width should be more or equal to the distance between resonant surfaces. This might be a problem at low shear, so with this figure you can quantify what it means a magnetic shear "too" low...
- Also it is testing that the distance between resonant surfaces is smaller compare to the gradient lengths. This can only be a problem when you do artificial scans and increase the normalized gradients too much, or in real, real steep ITBs.
· now you will be able to save the ExB shearing rate (using the new mdskize.m and jettokize.m) and then plot it with the ITG growth rates using plotkize_omeb.m.
- If you read from the mdsplus database, then mdskize.m will calculate the ExB shearing rate using the neoclassical formula for the poloidal velocity, and reading "Vrot" to evaluate the toroidal velocity and "Bpol" to have the poloidal B. Then the formula used to calculate the ExB shearing rate is the local one, i.e. on the low field side. This is unconsistent with the fact that the "Bpol" we read is a flux average value. If anyone has an idea to do it better, please let me know and please be careful I might have done something wrong, so check it carefully! I need your feedback...
- If you read data from JETTO, jettokize.m is reading "omeb" which is defined as the ExB shearing rate. I don't know which formula it is there for V_poloidal and for the ExB shearing rate.
- For both cases, to plot the resulting ExB shearing rate together with the ITG-TEM growth rates, just use plotkize_omeb.m, it will then read the new file called omExB_11111_060.mat which is a matlab file where r/a and the shearing rate are saved.
3. one upgrade in the code itself
· As many of you have already noticed, sometimes kinezero crashes, and it is always for the same reason: an integration from the NAG for the passing particles is required to have a precision of 2% but sometimes it does not converge. So the precision has been changed to 6%, I have checked that the results were unchanged in a couple of cases, you should do a check on one of your cases as well. And it will avoid to crash so often. To apply this, just copy the new fonctdeomc.f90 in your kinezero directory and type 'make clean' then 'make'. If kinezero still crashes, then just change slightly the radii and wave numbers you are choosing to avoid calculating again the same crashing case, for example 15 radial points instead of 16.
4. reminder
· don't forget that when the collisions are much higher than the typical trapped electron drift frequency you should use kizewote (with plotwote.m and plotwote_omeb.m), which is a version of Kinezero without Trapped Electrons.