Re: About pres2hybrid of NCL

Hi In-won Kim,

The fourth argument of pres2hybrid, "xi", must have rightmost dimensions nlevi x nlat x nlon, where nlevi is the length of "p" (the first argument), and nlat and nlon must be the rightmost two dimensions of "ps" (the second argument).

You should do "printVarSummary" on all of these variables to make sure they meet these size requirements.

--Mary

On Jun 13, 2012, at 8:42 PM, 김�� wrote:

> 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
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 14 11:14:51 2012