Hybrid level to pressure level

I want to convert the CCSM hybrid coordinate to pressure coordinate.So I
used the vinth2p function.I used the following scripts. My level is hybrid
level at midpoints (1000*(A+B))
So I am confused with the sentence in the description of function vinth2p
where it is saying that "If *datai* is on midlevels (interfaces), then *hyam
*/*hybm* (*hyai*/*hybi*) coefficients should be input." What is mean by *the
division hyam*/*hybm. If that is the case how I can change my scripts? or
the below script is O.K? Any help would be appreciated.*

begin

   f=addfile("file.nc","r")
   hyam=f->hyam
   hybm=f->hybm
   psurf=f->PS
   P0=1000
   lev_press=(/10,20,30,50,70,100,150,200,250,300 \
                    ,400,500,600,700,850,925,1000/)
   uwnd=f->uwnd ;time,lev,lat,lon
   u_new=vinth2p(uwnd,hyam,hybm,lev_press,psurf,1,P0,1,True)
end

Thanks in advance
sabeer

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 14 17:53:35 2011