Forgot to include the list!
-------- Original Message --------
Subject: Re: [ncl-talk] write_matrix.f Fortran runtime error: No such
file or directory
Date: Mon, 16 Apr 2012 12:36:24 -0600
From: David B. Reusch <dreusch@ees.nmt.edu>
Organization: New Mexico Tech
To: David Adams <dave.k.adams@gmail.com>
Hi Dave,
I think what you want here is:
cdo -f nc import_binary test.ctl ofile.nc
The error you're getting indicates you're missing the required CDO
operator ("import_binary" in this case). The ctl file is enough to
identify the input binary file.
Cheers,
Dave
David Adams said the following on 4/16/12 12:25 PM:
> Hi Dennis et al.
> I have install cdo operators and following the syntax you suggested.
> My binary file is
> S10216956_201111201830
> and my ctl file is test.ctl
>
> DSET ^S10216956_201111201830
> OPTIONS yrev byteswapped
> UNDEF 0
> TITLE Temp. Brilho para America do Sul
> XDEF 1870 LINEAR -100.00 0.04
> YDEF 1714 LINEAR -56.00 0.04
> ZDEF 01 LEVELS 1 1
> TDEF 1 LINEAR 17:00Z21mar2010 30mn
> VARS 1
> cls 0 -1,40,2,1 'Temp.
> ENDVARS
>
>
> Then I used the following command:
> cdo -f nc S10216956_201111201830 test.ctl ofile.nc <http://ofile.nc>
>
> which gives this error
> Error (operatorInqModID) : operator missing! S10216956_201111201830 is
> a file on disk!
>
> Can you help here?
>
> thanks,
> Dave
>
>
>
>
>
> On Mon, Apr 16, 2012 at 10:23 AM, Dennis Shea <shea@ucar.edu
> <mailto:shea@ucar.edu>> wrote:
>
> [1] You can read binary files in NCL. NCL does not require nc.
>
> [2] When reading binary in NCL, you can preset how NCL handles binary.
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml
>
> setfileoption("bin","ReadByteOrder","LittleEndian")
>
> [3] The Examples URL has a section on reading/converting "GrADS files"
> (a binary file with a text file describing how the file
> was written.)
>
>
> http://www.ncl.ucar.edu/Applications/
> See Miscellaneous at bottom
>
> If you insall the Climate Data Operators, their
> command line utilitry is the way to go.
>
>
> cdo -f nc import_binary file.ctl ofile.nc <http://ofile.nc>
>
>
> +++++++++++++++++++++++++++++++++++++++++++++++++++++
>
> On 04/16/2012 11:00 AM, David Adams wrote:
>
> Hi Mary,
> yes, that worked. I downloaded
> ncl_ncarg-6.0.0.Linux_Debian_x86_64_nodap_gcc445.tar
> and I was able to use write_matrix.
>
> On a related issue, I am trying to create netcdf files (so
> that I can
> use NCL)
> from a GrADS ctl file of binary GOES 12 Channel 4 data.
> It seems to be a 64-bit problem, because the data values are
> fine when I
> process them in GrADS on my 32-bit machine. I tried their
> byteswapped
> and non-byteswapped option, but I just get two different strange
> values. I know the 64-bit value are "correct" in that the
> patterns of
> the brightness temperature are identical to the patterns of
> the correct
> brightness temperature with my 32-bit machine.
>
> For example, one particular value at a point is
> 112
> for byteswapped
> It becomes
> 65458
> for non-byteswapped. The correct value on my 32-bit machine is
> 28850 (i.e., 288.5 K)
>
> any ideas?
>
> thanks in advance,
> Dave
>
>
>
>
>
>
>
>
>
> On Mon, Apr 16, 2012 at 7:06 AM, Mary Haley <haley@ucar.edu
> <mailto:haley@ucar.edu>
> <mailto:haley@ucar.edu <mailto:haley@ucar.edu>>> wrote:
>
> Dave,
>
> I would recommend downloading an NCL binary that was
> compiled with
> the highest version of gcc possible, which I believe is
> 4.4.5 for
> NCL V6.0.0.
>
> --Mary
>
> On Apr 13, 2012, at 9:49 PM, David Adams wrote:
>
> > Hi Mary et al.
> > I have the following version of gcc
> >
> > gcc version 4.6.3 20120306 (Red Hat 4.6.3-2) (GCC)
> >
> > Which binaries should I download given this version of gcc?
> >
> > thanks,
> > Dave
> >
> >
> >
> > On Fri, Apr 13, 2012 at 7:14 AM, Mary Haley <haley@ucar.edu
> <mailto:haley@ucar.edu>
> <mailto:haley@ucar.edu <mailto:haley@ucar.edu>>> wrote:
> > Hi Dave,
> >
> > This is a very strange error indeed. "writematrix.f" is the
> underlying Fortran code for the "write_matrix"
> > procedure. Both this code and the C wrapper are "built into" the
> "ncl" executable, so there shouldn't
> > be any issue with using it.
> >
> > I wonder if maybe it's a compatibility issue with the version of
> NCL you downloaded for your system.
> >
> > I would recommend upgrading to NCL V6.0.0, and trying one of
> those binaries instead.
> >
> > What version of gcc are you running? This might affect which
> binary you download.
> >
> > gcc --version
> >
> > For NCL V6.0.0, we offer binaries that were compiled with gcc
> 4.3.2, gcc 4.4.5, and gcc 4.1.2.
> > I tend to like the Debian binaries, as they seem to work well
> across many flavors of Linux.
> > I've had trouble with using the RedHat binaries on other
> non-RedHat systems.
> >
> > --Mary
> >
> >
> > On Apr 12, 2012, at 11:24 PM, David Adams wrote:
> >
> > > 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
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- Associate Research Professor of Climatology Dept of Earth and Environmental Science MSEC 304; 801 Leroy Place New Mexico Tech Socorro, NM 87801 -- Associate Research Professor of Climatology Dept of Earth and Environmental Science MSEC 304; 801 Leroy Place New Mexico Tech Socorro, NM 87801
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Apr 16 12:38:31 2012
This archive was generated by hypermail 2.1.8 : Fri Apr 20 2012 - 16:21:18 MDT