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