writing vertical coordinate with netcdf

Hi,

 I'm in trouble to convert Regional climate data with sigma level to netcdf output. I try to define the coordinate to write the netcdf format as following example, but it's very slow. Is there anyone know how to efficiently convert binary format to netcdf, such as FORTRAN code?

   sigma = new((/29/), float)
   sigma = (/0.00,0.03,0.05,0.07,0.09,0.11,0.13,0.16,0.20 ,0.25,0.31, \
            0.37,0.43 ,0.49 ,0.55,0.61,0.67,0.73,0.78,0.826,0.862 , \
            0.890,0.912,0.932,0.952,0.970,0.984,0.994,1.000/)

   sigma_half = new((/28/), float)
   do klev = 0 , 27
   sigma_half(klev) = (sigma(klev)+sigma(klev+1))/2.
   end do
   lev = sigma_half(27:8:-1)
   lev!0 = "lev"
   lev_at_long_name = "lev"
   lev_at_units = " "
   lev_at_standard_name = "sigma level"
   lev_at_positive = "down"
   lev_at_axis = "z"

Any suggestion are much appreciated!
 
Li
----------------------------------------------
Department of Meteorology
University of Hawaii at Manoa
2525 Correa Road, HIG 375
Honolulu, HI 96822, USA

Office : 1-808-9562047
---------------------------------------------------

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Dec 09 2005 - 13:49:56 MST