write_matrix.f Fortran runtime error: No such file or directory

Hi All,
I have recently upgraded to a 64-bit machine running
Fedora Core 16. I've installed the following
ncl_ncarg-5.2.1.Linux_RedHat_x86_64_nodap_gcc412.tar

All of my old ncl programs have been working fine, but today
I have a strange error.
I get the following error message
.
At line 172 of file writematrix.f
Fortran runtime error: No such file or directory

Where should write writematrix.f be? I have used locate, but can't find
the file.
The script that worked fine
on my on Linux 32-bit machine is here
Any ideas?
thanks,
Dave

; SAME as average_CTT_per_file.ncl, but read into
ave_subset_CTT_ncl_script.F90
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
files = systemfunc("ls 2011_01_3*.nc")
file_list = asciiread("file_list",-1,"string")
setfileoption("nc","SuppressClose",False)
f = addfiles(files,"r")
ListSetType (f, "join")
out_file = str_get_cols(file_list(0),0,14)
output_file = out_file + ".dat"
diri = "/home/dadams/Desktop/CPTEC_data/"
year = stringtofloat(str_get_cols(file_list,0,3))
month = stringtofloat(str_get_cols(file_list,5,6))
day = stringtofloat(str_get_cols(file_list,8,9))
hour = stringtofloat(str_get_cols(file_list,11,12))
minute = stringtofloat(str_get_cols(file_list,13,14))
new_temp = f[:]->new_temp
lat = f[0]->lat
lon = f[0]->lon
time_dimension = dimsizes(new_temp(:,0,0)); number of time steps (files)
new_temp@_FillValue = 32767
new_temp = new_temp/100.0
CTT_average = new((/1,dimsizes(new_temp(:,0,0))/),typeof(new_temp))
do n = 0, time_dimension-1
average = avg(new_temp(n,122:126,126:129))
non_missing = num(.not.ismissing(new_temp(n,122:126,126:129)))
CTT_average(0,n) = average
end do
m = dimsizes(CTT_average(0,:))
data = new ( (/m,6/), typeof(CTT_average) )
data(:,0) = year(:)
data(:,1) = month(:)
data(:,2) = day(:)
data(:,3) = hour(:)
data(:,4) = minute(:)
data(:,5) = CTT_average(0,:)
fmtx = "f5.0,2x,f3.0,2x,f3.0,2x,f3.0,2x,f3.0,2x,f9.3"
opt = True
opt@fout = diri + output_file
write_matrix (data, fmtx, opt)
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 12 23:25:06 2012