I was hoping someone could help me with this problem that I have been stuck on
for a while.&nbsp; <br />
I want to make a height-radius contour plot of temperature at 1 given time, i.e.
have <br />
radius from 0 to 150km on the x-axis and height z=39m to 18000m on the
y-axis.<br />
<br />
I open the file using:<br />
<br />
&nbsp; &nbsp; &nbsp;
&nbsp;OPEN(541,FILE=&#39;0.3sf_grads.dat&#39;,FORM=&#39;UNFORMATTED&#39;,<br />
&nbsp; &nbsp; &nbsp;&amp;&nbsp; &nbsp;ACCESS=&#39;DIRECT&#39;,RECL=KX*IX*4)<br
/>
<br />
where IX=32 and KX=30, that is 32 points in the x-dir and 30 points in the
y-dir.<br />
<br />
I then write the radius, height and temperature to the file:<br />
<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; IREC=1 <br />
&nbsp; &nbsp; &nbsp; &nbsp;WRITE(541,REC=IREC) ((R(K,I),K=1,KX),I=1,IR)<br />
&nbsp; &nbsp; &nbsp; &nbsp;IREC=IREC+1<br />
&nbsp; &nbsp; &nbsp; &nbsp;WRITE(541,REC=IREC) ((Z(K,I),K=1,KX),I=1,IR)<br />
&nbsp; &nbsp; &nbsp; &nbsp;IREC=IREC+1<br />
&nbsp; &nbsp; &nbsp; &nbsp;WRITE(541,REC=IREC) ((T(K,I),K=1,KX),I=1,IR)<br />
&nbsp; &nbsp; &nbsp; &nbsp;IREC=IREC+1<br />
<br />
It should be noted that R(K,1)=0.0 and R(K,IR)=155.0km.&nbsp; Also, Z(1,I)=39.0m
and <br />
Z(KX,I)=18000m.<br />
<br />
Now, I want to create the corresponding .ctl file (0.3sf_grads.ctl) such that
temperature <br />
is contoured on the r-z place.&nbsp; I have checked my tempearture, radius, and
height values <br />
and I know they are correct.&nbsp; I would like the distance (radius) in km to
be plotted on <br />
the x-axis and height in m on the y-axis.<br />
<br />
Any help in fixing this problem would be most appreciated, as I havent been able
to find <br />
a solution to the problem.&nbsp; Thanks. <br />
<br />
Jeff<p>
<br />
<br />
<br />
</p>