height-radius plot

[JEFFREY S GALL]

I was hoping someone could help me with this problem that I have been stuck on for a while.

I want to make a height-radius contour plot of temperature at 1 given time, i.e. have

radius from 0 to 150km on the x-axis and height z=39m to 18000m on the y-axis.



I open the file using:



       OPEN(541,FILE='0.3sf_grads.dat',FORM='UNFORMATTED',

     &   ACCESS='DIRECT',RECL=KX*IX*4)



where IX=32 and KX=30, that is 32 points in the x-dir and 30 points in the y-dir.



I then write the radius, height and temperature to the file:


       IREC=1

       WRITE(541,REC=IREC) ((R(K,I),K=1,KX),I=1,IR)

       IREC=IREC+1

       WRITE(541,REC=IREC) ((Z(K,I),K=1,KX),I=1,IR)

       IREC=IREC+1

       WRITE(541,REC=IREC) ((T(K,I),K=1,KX),I=1,IR)

       IREC=IREC+1



It should be noted that R(K,1)=0.0 and R(K,IR)=155.0km.  Also, Z(1,I)=39.0m and

Z(KX,I)=18000m.



Now, I want to create the corresponding .ctl file (0.3sf_grads.ctl) such that temperature

is contoured on the r-z place.  I have checked my tempearture, radius, and height values

and I know they are correct.  I would like the distance (radius) in km to be plotted on

the x-axis and height in m on the y-axis.



Any help in fixing this problem would be most appreciated, as I havent been able to find

a solution to the problem.  Thanks.



Jeff






-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070207/e63ac1c2/attachment.html