| Info |
|---|
In this tutorial - how can you integrate your code with UAL
- how can you access data via UAL
|
Accessing data from UAL requires some modification to your code. In this part of tutorial, we will take a closer look on how to access CPO via UAL.
| Warning |
|---|
Warning Stop here for a moment. Make sure you have followed the configuration setup before proceeding any further! Configuration related tutorial is here -> Click me! <- |
3.1 Accessing data using Python
...
Exercise no. 3 (approx. 15 min)
In this exercise you will browse MDSPlus database using jTraverser application.
1. source ITMv1 script by invoking
source $ITMSCRIPTDIR/ITMv1 kepler test 4.10b > /dev/null
2. Go to example directory
cd $ITMWORK/tutorials/split-2014$TUTORIAL_DIR/cpo_basics/python/
3. Execute sample code
| Code Block |
|---|
|
shell> python ./put_cpo_array.py
shell> python ./put_cpo_slices.py |
4. Open example file
...
| Handling CPOs: put() vs. putSlice() |
|---|
| Code Block |
|---|
| language | py |
|---|
| title | $TUTORIAL_DIR/ |
|---|
|
|
...
| cpo_basics/python/put_cpo_array.py |
|
|
...
...
| import sys
from pylab import *
import ual
cpo = ual.itm( |
|
...
...
3)
cpo.create()
if not cpo.isConnected():
print 'error during itmdb entry creation'
sys.exit(1)
cpo.equilibriumArray.resize( |
|
...
...
equi = cpo.equilibriumArray |
|
...
#First fill fields which are not time-dependent.
equi.array[0].datainfo.dataprovider = 'MKO'
equi.array[0].datainfo.putdate = ' |
|
...
...
...
2016'
equi.array[0].codeparam.parameters = 'param'
# ---- a loop ----
for i in range(0, 10):
#Fill time-dependent fields
equi.array[ |
|
...
i].eqgeometry.boundary.resize(1)
equi.array[ |
|
...
i].eqgeometry.boundary[0].r = sin(arange(0,2*pi,2*pi/100)) + (1/float(100-i))
equi.array[ |
|
...
i].eqgeometry.boundary[0].z = cos(arange(0,2*pi,2*pi/100)) + (1/float(100-i))
|
|
...
#Do not forget time!!
equi.array[ |
|
...
i].time = i
# ---- a loop ----
#Save data in the database
equi.put() # <= Called outside the loop
#close the pulse file
cpo.close() |
| | Code Block |
|---|
| language | py |
|---|
| title | $TUTORIAL_DIR/cpo_basics/python/put_cpo_slices.py |
|---|
| linenumbers | true |
|---|
| import sys
from pylab import *
import ual
cpo = ual.itm(14,4)
cpo.create()
if not cpo.isConnected():
print 'error during itmdb entry creation'
sys.exit(1)
equi = cpo.equilibrium
#First fill fields which are not time-dependent.
equi.datainfo.dataprovider = 'MKO'
equi.datainfo.putdate = '20/09/2016'
equi.codeparam.parameters = 'param'
#Save time independent fields
equi.putNonTimed()
# ---- a loop ----
for i in range(0, 10):
#Fill time-dependent fields
equi.eqgeometry.boundary.resize(1)
equi.eqgeometry.boundary[0].r |
|
...
= sin(arange(0,2*pi,2*pi/100)) + (1/float(100-i))
equi. |
|
...
eqgeometry.boundary[0].z = cos(arange(0,2*pi,2*pi/100)) + (1/float(100-i))
#Do not forget time!!
equi.time = i
#Append this slice in the database
equi. |
|
...
putSlice() # <= Called inside the loop
# ---- a loop ----
#close the pulse file
cpo.close() |
|
 Image Added | VI related notice
I will use VI in every place where text files are modified. If you have any other text file editor of your choice - fell free to use it instead. |
vi $TUTORIAL_DIR/cpo_basics/python/put_cpo.py
What you can see here is a simple code that stores particular CPO into MDSPlus database using UAL.
Let's check how to read these data in Fortran.
...
2. Change directory to a demo location for this exercise
cd $ITMWORK/tutorials/split-2014/$TUTORIAL_DIR/cpo_basics/fortran
3. Open GetSlice.f90 file
vi GetSlice.f90
...
Esc Esc : q ! Enter
| Handling CPOs: get() vs. getSlice() |
|---|
| Code Block |
|---|
| title | $TUTORIAL_DIR/cpo_basics/fortran/get_cpo_array.f90 |
|---|
| linenumbers | true |
|---|
| program diagnostic
use euitm_schemas
use euitm_routines
implicit none
integer :: idx, i, arraySize
type (type_equilibrium), pointer :: eqArray(:) ! <= Array !!!
call euitm_open('euitm', 14, 4, idx)
call euitm_get(idx, 'equilibrium', eqArray)
arraySize = size(eqArray)
write (*,*) "Number of slices: ", arraySize
do i=1, arraySize
write (*,*) "Time: ", eqArray(i)%time
write (*,*) "Value of r: ", eqArray(i)%eqgeometry%boundary(1)%r(1)
write (*,*) "Value of z: ", eqArray(i)%eqgeometry%boundary(1)%z(1)
enddo
call euitm_close(idx)
end program |
| | Code Block |
|---|
| title | $TUTORIAL_DIR/cpo_basics/fortran/get_cpo_slices.f90 |
|---|
| linenumbers | true |
|---|
|
|
...
| program diagnostic
use euitm_schemas
use euitm_routines |
|
...
...
...
implicit none
integer :: idx, i, arraySize
type (type_equilibrium) :: equilibrium
real(EUITM_R8), pointer :: |
|
...
...
...
real(EUITM_R8) :: time
call euitm_open('euitm', |
|
...
...
4, idx)
call euitm_get_times(idx, 'equilibrium', |
|
...
timeVector)
arraySize = size(timeVector)
|
|
...
...
...
slices: ", arraySize
write (*,*) "Time vector: ", timeVector
do i = 1, arraySize
time = timeVector(i)
call euitm_get_slice(idx, 'equilibrium', equilibrium, time, 1)
write (*,*) "Time [", time, "]: ", equilibrium%time
write (*,*) "Value of r: ", |
|
...
equilibrium%eqgeometry%boundary(1)%r(1)
write (*,*) "Value of z: ", |
|
...
equilibrium%eqgeometry%boundary(1) |
|
...
%z(1)
enddo
call euitm_close(idx)
end program |
|
4. Compile the call euitm_disconnect()
end program
code
| Code Block |
|---|
shell> make clean
shell> make |
5. Run the code
| Code Block |
|---|
shell> ./get_cpo_array.exe
shell> ./get_cpo_slices.exe |
make clean
make
./readEquilibrium
6. You should see values that we have stored using Python based code.
...
All these topics will be covered in separate tutorial: 1.4 Using FC2K for embedding Fortran code into Kepler
