[ncl-talk] esmf regrid

[Buzan, Jonathan R]

Hello NCL-Talk,

I am attempting to regrid 0.5x0.5 (360x720) to 0.5’x0.5’ (17400x43200) resolution.  I am not sure what I am doing wrong.


Error:

(0)     get_src_grid_info: source lat dims = (360)
(0)     get_src_grid_info: source lon dims = (720)
(0)     get_src_grid_info: source grid type is 'rectilinear'
(0)     curvilinear_to_SCRIP: calculating grid corners...
(0)     curvilinear_to_SCRIP: no lat values are at the poles, so
(0)            calculating grid corners using
(0)            calc_SCRIP_corners_noboundaries...
(0)     calc_SCRIP_corners_noboundaries
(0)          min/max original lat: -89.75/89.75
(0)          min/max original lon: -179.75/179.75
(0)     calc_SCRIP_corners_noboundaries
(0)          min/max Extlat2d: -90.25/90.25
(0)          min/max Extlon2d: -180.25/180.25
(0)     calc_SCRIP_corners_noboundaries
(0)          min/max ExtGridCenter_lat: -90/90
(0)          min/max ExtGridCenter_lon: -180/180
(0)     =====> CPU Elapsed Time: rectilinear_to_SCRIP: 0.182972 seconds <=====
(0)     get_dst_grid_info: destination lat dims = (17400)
(0)     get_dst_grid_info: destination lon dims = (43200)
ncendef: ncid 196608: NetCDF: One or more variable sizes violate format constraints
fatal:NetCDF: Operation not allowed in define mode: error attempting to write variable (grid_dims) to file (destination_grid_file.nc)

relevant code:

	Opt                  = True
	Opt at InterpMethod     = "bilinear";"conser"
;	Opt at InterpMethod     = "neareststod";"conser"
;	Opt at InterpMethod     = "bilinear";"conser"
	Opt at WgtFileName      = wgt_FILE_NAME
	print("madeit1")
	Opt at SrcLargeFile     = "True"
	print("madeit2")
	Opt at DstLargeFile     = "True"
	print("madeit3")
	Opt at SrcGridLat     = latyi
    	Opt at SrcGridLon     = lonxi
	Opt at SrcRegional    = False
	Opt at ForceOverwrite = True
	Opt at Debug          = True
	Opt at PrintTimings   = True

;********* Desired lat lon **************

      Opt at DstGridLat = latyo
      Opt at DstGridLon = long

      cntrycodes         = ESMF_regrid(PRECTmms_targetfile,Opt)