2D plot of filled contour with Plplot libraries using plshades

In the previous post, we show how to plot 2D data with plcont. It plots the line contour. Now, how about if we want something more colorful? we can use a filled contour plot with the command plshade. The syntax is similar. But first, let’s look at the complete program.

program simple
use plplot
implicit none
real(plflt),dimension(6) :: x,y
real(plflt),dimension(6,6)::z,x2d,y2d,z2
real(plflt)::xmin,xmax,ymin,ymax,zmin,zmax
real(plflt),dimension(10)::clevel
real::step
integer::nlevel
integer::i,j
character*1 [...]

2D Contour plot using plplot function “plcont”

In the previous post, we describe how to do a simple line plot. In this post, I will describe how to plot contour using Plplot libraries, particularly using the “plcont” function. Let’s look at the code first.

program simple
use plplot
implicit none
real(plflt),dimension(6) :: x,y
real(plflt),dimension(6,6)::z,x2d,y2d
real(plflt)::xmin,xmax,ymin,ymax,zmin,zmax
real(plflt),dimension(10)::clevel
real::step
integer::nlevel
integer::i,j
x=(/1,2,3,4,5,6/)
y=x
do i=1,6
do j=1,6
z(i,j)=x(i)*2*y(j)
x2d(i,j)=x(i)
[...]

Wiki for Corporate

my presentation during presentation skills workshop with the topic: Wiki for Corporate

Using Plplot from fortran in Kubuntu 8.04

I installed libplplot from Kubuntu Adept package manager, and chose “libplplot-fortran9″ package.
To test the code, I created the following fortran program:

program testplplot2d
use plplot
implicit none
real(plflt),dimension(6) :: x,y
real(plflt)::xmin,xmax,ymin,ymax
x=(/1,2,3,4,5,6/)
y=x**2
write(*,*) y
call plinit()
xmin=1.0
xmax=6.0
ymin=1.0
ymax=40.0
call plcol0(1)
[...]

Cleaning up “/” root directory

When I tried to install something, it turns out that my root directory usage is 100%, to check type “df -h”. I couldn’t figure out what was happening until I realize that I have updated so many Linux kernel. So I checked what is the latest Linux kernel version I have, and uninstall the rest. [...]

Installing Mayavi2 in Kubuntu/Ubuntu 8.04

Open the package manager, such as Adept or Synaptic. Add the following lines to the repository:

deb http://ppa.launchpad.net/gael-varoquaux/ppa/ubuntu hardy main
deb-src http://ppa.launchpad.net/gael-varoquaux/ppa/ubuntu hardy main

And reload the packages, search for “mayavi2″, and install. I encountered some problem (I forgot to take note the error message), it shows some import cannot be found. After I googled, it turns out [...]

Updating OpenOffice 3.1.0 to 3.1.1 in Kubuntu 8.04

I found the update from OpenOffice 3.1.0 to 3.1.1 to be pretty simple. When I run an openoffice program, there is this Update notification on the top right of the window. So I clicked and it download the update. After finished downloading, I clicked “Install”, but as expected it didn’t work.
So I went to the [...]

What does !DEC$ mean in fortran

I encounter this symbol !DEC$ in fortran code, and tried to figure out what it means. It’s rather difficult to search using google, but after a while, I got a good explanation. It is basically similar to preprocessor in C/C++, you may read further from here:
http://h21007.www2.hp.com/portal/download/files/unprot/Fortran/docs/lrm/lrm0588.htm

SOPRA Database format

This is the format of SOPRA database provided by FreeSnell manual.
The first number indicates the spectral unit (1 = eV, 2 = um, 3= 1/cm, 4 = nm). The second and the third numbers indicate the spectral range and the fourth one the number of intervals. The N and K values are then stored in [...]