Re: Navigation

• To: Charlie Zender <zender@xxxxxxxxxxxxx>
• Subject: Re: Navigation
• From: Ana Gusmao <angusmao@xxxxxxxxxxxxxxxxxxx>
• Date: Fri, 3 Sep 1999 08:53:30 -0300 (EST)

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
> 

• Follow-Ups:
  • Re: Navigation
    • From: Charlie Zender

• References:
  • Re: Navigation
    • From: Charlie Zender