On Thu, 2 Sep 1999, Charlie Zender wrote:
Hi Charlie,
I am one of the users of ccm3.6 and I am trying to
construct a set of boundary conditions for SST
that is the average of to years set of SST.
I've being trying to do this using the CDL language.
But I was not well suceeded in doing it.
The " rdsst " program to read the sst_jan1978_sep1993.nc
input data of ccm3 is down below.
I am having lots of problem witht the "headers "
So the program doesn't work.
Is there any simple way
I can read the sst_jan1978_sep1993.nc
and then generate an average of two year set from this
sst_jan1978_sep1993.nc. ?
Regards and Thanks
Ana
Ana Maria Gusmao
University of Sao Paulo
Sao Paulo, Brasil
Tel (55) (021)(11)81814713
=================================================
Program rdsst
c
parameter ( mth=12)
real hedr(339),dati(128),flati,floni
real sst(128,64)
integer hedi(93)
character*8 hedc(213)
character*22 TVBDS
real flatd(64)
c
c
c Set sea-ice surface temperature. See CCM3.2 source code
TVBDS = 'sst_jan1978_sep1993.nc' ! 12 months of SST climatology
c
c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
open (unit=51,file=TVBDS,status='old',form='unformatted')
c OROMASK defines Land=1, Sea-Ice=2 and Ocean=0
c +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
c
do 50 mo=1,mth
C - Now read headers and data from NCAR SST's -
read(51) hedi(1),(hedi(i),i=2,hedi(1))
read(51) (hedc(i),i=1,hedi(31))
read(51) (hedr(i),i=1,hedi(32))
write(6,*)' Headers from tvbds data set read'
do 40 j=1,64
flatd(j) = hedr(hedi(34)+j-1)
40 continue
do 43 j=1,64
read(51) flati,floni,(dati(i),i=1,128)
c Keep the original data set for NCAR sst's
do i=1,128
sst(i,j) = dati(i)
enddo
c -----------------------------------------------------------
43 continue
50 continue
write (6,*)' Termine'
stop
end
> Hi Daran,
>
> NCO is a set of command line utilities that do useful
> things to netCDF files. ncks is one such operator:
>
> ncks -H -v temperature -d lon,100.0 -d lat,40.0 -d lev,700.0 data.nc
>
> , e.g., prints value of temperature nearest 700mb nearest Denver.
>
> See http://www.cgd.ucar.edu/cms/nco
>
> Charlie
> --
> Through September 10:
> Charlie Zender Voice: (303) 497-1445, FAX: 497-1400
> NCAR ACD & CGD E-mail: zender@xxxxxxxxxxxxx
> P.O. Box 3000 URL: http://www.cgd.ucar.edu/cms/zender
> Boulder CO 80307-3000 PGP: finger -l zender@xxxxxxxxxxxxxxxxxxxxx
>
> After September 10:
> Charlie Zender zender@xxxxxxx (949) 824-2987/FAX-3256
> www.ess.uci.edu/zender
> Earth System Science 216 PSRF, University of California, Irvine, CA
> 92697-3100
>
