About pres2hybrid of NCL

From: ±èÀοø <rladlsdnjs88_at_nyahnyahspammersnyahnyah>
Date: Wed Jun 13 2012 - 20:42:09 MDT

Hello.
My name is In-won Kim. Currently I am Ph.D. student.
Now I am trying to interpolate data on constant pressure levels to hybrid
levels using the pres2hybrid to make initial of CAM model.
But, there are some errors.

*Error message*
*=====================================================================*
*fatal:pres2hybrid: The three rightmost dimensions of 'xi' must be nlevi x
nlat x nlon*
*fatal:Execute: Error occurred at or near line 17 in file test1.ncl*

I don't understand why this error occurs, because my data has
3-dimension(26 layers) like this.

        float TMP_P0_L100_GLL0(lv_ISBL0, lat_0, lon_0) ;
                TMP_P0_L100_GLL0:initial_time = "06/06/2011 (00:00)" ;
                TMP_P0_L100_GLL0:forecast_time = 0 ;
                TMP_P0_L100_GLL0:level_type = "Isobaric surface (Pa)" ;

TMP_P0_L100_GLL0:parameter_template_discipline_category_number = 0, 0, 0, 0
;
                TMP_P0_L100_GLL0:parameter_discipline_and_category =
"Meteorological products, Temperature" ;
                TMP_P0_L100_GLL0:grid_type = "Latitude/longitude" ;
                TMP_P0_L100_GLL0:_FillValue = -999.f ;
                TMP_P0_L100_GLL0:units = "K" ;
                TMP_P0_L100_GLL0:long_name = "Temperature" ;
                TMP_P0_L100_GLL0:production_status = "Operational products"
;
                TMP_P0_L100_GLL0:center = "US National Weather Service -
NCEP (WMC)" ;

*Could you explain how to solve or give me the sample data(possible to
interpolate)?*
*
*
*
*
Thanks.

Sincerely yours.

*Detailed code*
*======================================================================*
    begin

      f = addfile("gfs.t00z.pgrb2f00.nc", "r")
      tp = f->TMP_P0_L100_GLL0
      fP = addfile("./gfs_sfc_pressure.nc","r")
      psfc = fP->PRES_P0_L1_GLL0
      p0 = 100000.

      fC = addfile("./cami_1980-01-01_0.47x0.63_L26_c071226_org.nc","r")
      hya = fC->hyam
      hyb = fC->hybm

      plev = f->lv_ISBL0
      plev = plev*100
      plev@units = "Pa"

     th = pres2hybrid(plev,psfc,p0,tp,hya,hyb,0)

     end

-- 
-- 
In-Won Kim
Intergrated Climate System Modeling Lab.
Department of Environmental Atmospheric Sciences
599-1 Daeyeon-dong, Nam-gu
Pukyong National University,
Busan, S.Korea
Tel. +82-51-629-6643
Fax +82-51-629-7991
Received on Wed Jun 13 20:42:16 2012

This archive was generated by hypermail 2.1.8 : Fri Jun 15 2012 - 14:51:31 MDT