[problems] satellite data -regridding ESMF

Hi users,
I'm using satellite data, but when I try to run the script to change the =
grid (esmf), returns the errors:
(0) get_src_grid_info: source lat dims = (3,4096)(0) get_src_grid_info:=
 source lon dims = (3,8192)(0) get_src_grid_info: SrcGridType and/or Sr=
cGridLat/SrcGridLon were not set.(0) Cannot determine the source =
grid type.
Variable: varType: shortTotal Size: 201326592 bytes 100663296 va=
luesNumber of Dimensions: 3Dimensions and sizes: [ncl_join | 3] x [lat | 40=
96] x [lon | 8192]Coordinates: Number Of Attributes: 28 bsst.dsp_PixelType=
 : 1 bsst.dsp_PixelSize : 2 bsst.dsp_Flag : 0 bsst.dsp_nBits : 16 bsst.=
dsp_LineSize : 0 bsst.dsp_cal_name : Temperature bsst.units : Temp bsst.=
dsp_cal_eqnNumber : 2 bsst.dsp_cal_CoeffsLength : 8 bsst.dsp_cal_coeffs :=
        ( 0.075, -3 ) bsst.scale_factor : 0.075 bsst.add_off : -3 dsp_PixelTy=
pe : 1 dsp_PixelSize : 2 dsp_Flag : 0 dsp_nBits : 16 dsp_LineSize : 0 =
dsp_cal_name : Temperature units : Temp dsp_cal_eqnNumber : 2 dsp_cal_Co=
effsLength : 8 dsp_cal_coeffs : ( 0.075, -3 ) scale_factor : 0.075 add=
_off : -3 dim_0.name : lat dim_0.long_name : latitude dim_1.name : lon =
dim_1.long_name : longitude

beginfils = systemfunc("ls *.nc") f = addfiles(fils, "r") =
                     ListSetType(f,"join") var=f[:]->bsstprint=
VarSummary(var)src_lat = f[:]->lat src_lon = f[:]->lon newfi=
le = "teste.nc"system("rm -f "+newfile)ncdf = addfile(newfile ,"c") =
; open output netCDF file setvalues NhlGetWorkspaceObjectId "wsMaximum=
Size" : 300000000end setvalues Opt = TrueOpt@SrcGridL=
at = src_lat ; source gridOpt@SrcGridLon = sr=
c_lon;Opt@SrcMask2D = maskfile->MASKSURF;newlat = fspan( -8=
9.7, 89.7,224)newlon = fspan(0.0,359.7,256)tGridType ==
 "rectilinear"Opt@DstGridLat = newlatOpt@DstGridLon =
= newlonOpt@DstFileName = "Rectilinear.nc" ; destination =
filesOpt@InterpMethod = "bilinear" Opt@ForceOverwrite = True=
Opt@PrintTimings = TrueOpt@Debug = Truevar_regrid = =
ESMF_regrid(var,Opt) ; Do the regridding var_regrid!0 ="time"va=
r_regrid!1 ="lat"var_regrid!2 = "lon"var_regrid&time = var&timeva=
r_regrid&lat = newlatvar_regrid&lat@units="degrees_north"var_regrid&lon=
 = newlonvar_regrid&lon@units="degrees_east" ncdf->time_bnds = f[:]->=
time_bndsncdf->bsst = var_regridend

---Vanúcia SchumacherMestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEACel: (31) 9978 2522 DEA: (31) =
3899 1890 =

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


Received on Mon Jul 07 09:49:42 2014