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
This archive was generated by hypermail 2.2.0 : Sun Dec 11 2005 - 21:55:06 MST