Re: cmip5 搇on2d,lat2d攄ata

Hi,

Did you download a newer version of the "ESMF_regridding.ncl" script from:

http://www.ncl.ucar.edu/Applications/Files/ESMF_regridding.ncl

Once you do this, you either need to copy it to $NCARG_ROOT/lib/ncarg/nclscripts/esmf/,
or, if you don't have permission to do this, then put it in the same directory as
your script, and edit your script to load this script instead:

; load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
load "./ESMF_regridding.ncl"

--Mary

On Oct 10, 2012, at 1:24 AM, 孙yan wrote:

> thank you very much,Mary Haley . your suggest is right.
> I see example ESMF_regrid_6.ncl. in ESMF_regrid_6.ncl,the time dimension is
> only one,it works well.when I use the whole time dimensions,new problem
> appears:
>
> my ncl script is:
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
>
> begin
> ;---Interpolation methods
> methods = (/"patch"/)
>
> ;---Input file
> srcFileName = "tos_Omon_MPI-ESM-P_historical_r1i1p1_185001-200512.nc"
> wgtFile = "CMIP5_2_World_" + methods + ".nc"
>
> ;---Get data and lat/lon grid from CMIP5 Grid
> sfile = addfile(srcFileName,"r")
> thetao = sfile->tos(0,:,:)
> thetao@lat2d = sfile->lat
> thetao@lon2d = sfile->lon
> printVarSummary(thetao)
>
> Opt = True
> Opt@SrcFileName = "CMIP5_SCRIP.nc" ; source file name
> Opt@DstFileName = "World1deg_SCRIP.nc" ; destination file name
> Opt@ForceOverwrite = True
>
> Opt@SrcGridCornerLat = sfile->lat_vertices ; corners are necessary
> Opt@SrcGridCornerLon = sfile->lon_vertices ; for "conserve" method
> Opt@SrcMask2D = where(.not.ismissing(thetao),1,0)
> Opt@DstGridType = "1x1" ; Destination grid
> Opt@DstTitle = "World Grid 1-degree Resolution"
> Opt@DstLLCorner = (/-89.75d, 0.00d /)
> Opt@DstURCorner = (/ 89.75d, 359.75d /)
>
> print("Generating interpolation weights from CMIP5 to")
> print("World 1 degree grid using the " + methods + " method.")
>
> Opt@WgtFileName = wgtFile
> Opt@InterpMethod = methods
>
> ;----------------------------------------------------------------------
> ; Interpolate data from CMIP5 to World 1-degree grid.
> ;----------------------------------------------------------------------
> thetao_regrid = ESMF_regrid(thetao,Opt)
> printVarSummary(thetao_regrid)
>
> end
>
>
> fatal:Eq: Dimension size, for dimension number 0, of operands does not
> match, can't continue
>
> fatal:["Execute.c":7556]:Execute: Error occurred at or near line 1071 in
> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":7556]:Execute: Error occurred at or near line 2810 in
> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":7556]:Execute: Error occurred at or near line 2921 in
> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> I don't know what the problem is. thank you.
> -------------------------------------------------------------------
> -----原始邮件-----
> From: Mary Haley
> Sent: Monday, October 01, 2012 1:49 AM
> To: 孙yan
> Cc: ncl-talk@ucar.edu
> Subject: Re: [ncl-talk] cmip5 “lon2d,lat2d”data
>
>
> On Sep 27, 2012, at 6:28 PM, 孙yan wrote:
>
>> hello,
>> (1) I use cmip5 data, many datas are “lon2d,lat2d”data.when I use
>> “rcm2rgrid” function to regrid,some region is FillValue.why?
>> (2) how to deal with “lon2d,lat2d”data?
>> thank you
>>
>
> There's not enough information to help you. It doesn't surprise me that you
> would have missing values, especially if the rectilinear grid goes outside
> the range of your curvilinear grid.
>
> Did you print out any values, or try plotting the regridded data to see if
> it looks correct?
>
> You can also visit our ESMF regridding page. See example ESMF_regrid_6.ncl:
>
> http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex6
>
> --Mary
>
>>
>> ------------------------------------------------------------------
>>
>> National Key Laboratory of Atmospheric Sciences and
>> Geophysical Fluid Dynamics (LASG),
>> Institute of Atmospheric Physics (IAP),
>> Chinese Academy of Sciences (CAS),
>>
>> -------------------------------------------------------------------
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 10 10:21:47 2012