<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Here’s what I can get with printVarSummary(var);</div><div class=""><br class=""></div><div class="">Variable: var</div><div class="">Type: float</div><div class="">Total Size: 3588096000 bytes</div><div class=""> 897024000 values</div><div class="">Number of Dimensions: 3</div><div class="">Dimensions and sizes: [time | 7300] x [j | 384] x [i | 320]</div><div class="">Coordinates:</div><div class=""> time: [719050.5..726349.5]</div><div class=""> j: [1..384]</div><div class=""> i: [1..320]</div><div class="">Number Of Attributes: 12</div><div class=""> standard_name : sea_ice_area_fraction</div><div class=""> long_name : Sea Ice Area Fraction</div><div class=""> comment : aice_d no change, CMIP5_table_comment: fraction of grid cell covered by sea ice.</div><div class=""> units : %</div><div class=""> original_name : aice_d</div><div class=""> cell_methods : time: mean (interval: 30 days)</div><div class=""> cell_measures : area: areacello</div><div class=""> history : 2014-09-29T18:12:07Z altered by CMOR: replaced missing value flag (1e+30) with standard missing value (1e+20).</div><div class=""> missing_value : 1e+20</div><div class=""> _FillValue : 1e+20</div><div class=""> associated_files : baseURL: <a href="http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation" class="">http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation</a> gridspecFile: gridspec_seaIce_fx_CCSM4_historical_r0i0p0.nc areacello: areacello_fx_CCSM4_historical_r0i0p0.nc</div><div class=""> coordinates : lat lon</div><div class=""><br class=""></div><div class="">Thank you,</div><div class="">Tomoko</div><div class=""><br class=""></div><div class=""><br class=""></div><div><blockquote type="cite" class=""><div class="">On Nov 12, 2017, at 9:02 AM, music piano <<a href="mailto:musicpianoljy@gmail.com" class="">musicpianoljy@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hi, Tomoko<div class=""><br class=""></div><div class="">Would you try to use printVarSummary(var) to see what type of the data is ? </div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Nov 11, 2017 at 10:50 PM, Tomoko Koyama <span dir="ltr" class=""><<a href="mailto:Tomoko.Koyama@colorado.edu" target="_blank" class="">Tomoko.Koyama@colorado.edu</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br class="">
<br class="">
I am trying to regrid Sea Ice Area Fraction and Sea Ice Thickness data from some CMIP5 daily data to EASE 100km grid.<br class="">
The following script works for Sea Level Pressure, Near-Surface Temperature, Geopotential Height etc. but not for those ice parameters.<br class="">
<br class="">
The error message is “Assignment type mismatch, right hand side can't be coerced to type of left hand side” occurring at regridding;<br class="">
var_regrid = ESMF_regrid(var,Opt)<br class="">
<br class="">
Does anyone see why this might be???<br class="">
<br class="">
Thank you in advance,<br class="">
Tomoko<br class="">
<br class="">
<br class="">
+++++++++++++++++<br class="">
+ Error message +<br class="">
+++++++++++++++++<br class="">
<br class="">
Copyright (C) 1995-2015 - All Rights Reserved<br class="">
University Corporation for Atmospheric Research<br class="">
NCAR Command Language Version 6.3.0<br class="">
The use of this software is governed by a License Agreement.<br class="">
See <a href="http://www.ncl.ucar.edu/" rel="noreferrer" target="_blank" class="">http://www.ncl.ucar.edu/</a> for more details.<br class="">
(0) get_src_grid_info: source lat dims = (384)<br class="">
(0) get_src_grid_info: source lon dims = (320)<br class="">
(0) get_src_grid_info: source grid type is 'rectilinear'<br class="">
(0) curvilinear_to_SCRIP: calculating grid corners...<br class="">
(0) curvilinear_to_SCRIP: one or more lat values are at the<br class="">
(0) poles, so calculating grid corners using<br class="">
(0) calc_SCRIP_corners_boundaries.<wbr class="">..<br class="">
(0) calc_SCRIP_corners_boundaries<br class="">
(0) min/max original lat: 1/384<br class="">
(0) min/max original lon: 1/320<br class="">
fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side<br class="">
^Mfatal:["Execute.c":8575]:<wbr class="">Execute: Error occurred at or near line 770 in file $NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl<br class="">
<br class="">
^Mfatal:["Execute.c":8575]:<wbr class="">Execute: Error occurred at or near line 1744 in file $NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl<br class="">
<br class="">
^Mfatal:["Execute.c":8575]:<wbr class="">Execute: Error occurred at or near line 1822 in file $NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl<br class="">
<br class="">
^Mfatal:["Execute.c":8575]:<wbr class="">Execute: Error occurred at or near line 3851 in file $NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl<br class="">
<br class="">
^Mfatal:["Execute.c":8575]:<wbr class="">Execute: Error occurred at or near line 3983 in file $NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl<br class="">
<br class="">
^Mfatal:["Execute.c":8573]:<wbr class="">Execute: Error occurred at or near line 86<br class="">
<br class="">
++++++++++<br class="">
+ Script +<br class="">
++++++++++<br class="">
<br class="">
load "$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/gsn_code.ncl"<br class="">
load "$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/gsn_csm.ncl"<br class="">
load "$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/csm/contributed.<wbr class="">ncl"<br class="">
load "$NCARG_ROOT/lib/ncarg/<wbr class="">nclscripts/esmf/ESMF_<wbr class="">regridding.ncl"<br class="">
<br class="">
method = "bilinear"<br class="">
model = "CCSM4"<br class="">
gridtype = "latlon"<br class="">
<br class="">
rdiri = "/work/koyama_sice"<br class="">
in_diri = "/home/koyama/CMIP5/data/<wbr class="">input"<br class="">
esmf_diri = in_diri + "/" + model<br class="">
<br class="">
begin<br class="">
; Source grid<br class="">
srcGridDiri = rdiri + "/" + model<br class="">
fili = systemfunc("cd " + srcGridDiri + " ; ls sic*.nc")<br class="">
nfil = dimsizes(fili)<br class="">
<br class="">
; Destination grid<br class="">
dstGridDiri ="/home/koyama/CMIP5/data/<wbr class="">grid_coord"<br class="">
dstFileName = "NRims180x180.nc"<br class="">
dstGridFili = dstGridDiri + "/" + dstFileName<br class="">
dfile = addfile(dstGridDiri + "/" + dstFileName, "r")<br class="">
lat2d = dfile->latitude<br class="">
lon2d = dfile->longitude<br class="">
<br class="">
do ifil = 0, nfil-1<br class="">
srcGridFili = fili(ifil)<br class="">
sfile = addfile(srcGridDiri + "/" + srcGridFili, "r")<br class="">
dSizes = getfiledimsizes(sfile)<br class="">
var = sfile->sic<br class="">
<br class="">
lat1d = ndtooned(lat2d)<br class="">
lon1d = ndtooned(lon2d)<br class="">
lon1d = where(lon1d.lt.0,360.0+lon1d,<wbr class="">lon1d)<br class="">
latlon_dims = dimsizes(lat2d) ; 180x180<br class="">
<br class="">
Opt = True<br class="">
<br class="">
Opt@ForceOverwrite = True<br class="">
<br class="">
Opt@SrcFileName = esmf_diri + "/" + model +"_SCRIP.nc" ; output file names<br class="">
Opt@DstFileName = esmf_diri + "/" + "EASE_ESMF.nc"<br class="">
Opt@WgtFileName = esmf_diri + "/" + model +"_2_EASE_"+method+".nc"<br class="">
<br class="">
Opt@SrcInputFileName = srcGridFili ; optional, but good idea<br class="">
Opt@DstInputFileName = dstGridFili<br class="">
<br class="">
Opt@SkipSrcGrid = False<br class="">
Opt@SkipDstGrid = False<br class="">
Opt@SkipWgtGen = False<br class="">
Opt@DstESMF = True<br class="">
<br class="">
Indexes = ind(.not.ismissing(lon1d))<br class="">
Opt@DstGridType = "unstructured" ; Destination grid<br class="">
Opt@DstGridLat = lat1d(Indexes) ; Strip off missing data<br class="">
Opt@DstGridLon = lon1d(Indexes)<br class="">
<br class="">
Opt@RemapIndexes = True ; This is necessary to remap<br class="">
Opt@Indexes = Indexes ; regridded values back to<br class="">
Opt@IndexesDims = latlon_dims ; locations on destination grid<br class="">
<br class="">
Opt@InterpMethod = method<br class="">
<br class="">
Opt@Debug = True<br class="">
Opt@PrintTimings = True<br class="">
<br class="">
Opt@CopyVarCoords = False ; we can't copy the coords because<br class="">
; the weights on the file are only<br class="">
; for the non-missing lat/lon values.<br class="">
<br class="">
Opt@SrcGridMask = where(ismissing(var(0,:,:)),0,<wbr class="">1)<br class="">
<br class="">
var_regrid = ESMF_regrid(var,Opt) ; Regrid var<br class="">
<br class="">
copy_VarAtts(var,var_regrid)<br class="">
var_regrid@_FillValue = 1.e+20<br class="">
var_regrid@missing_value = 1.e+20<br class="">
var_regrid@_FillValue = var@_FillValue<br class="">
<br class="">
var_regrid!0 = "time"<br class="">
var_regrid!1 = "i"<br class="">
var_regrid!2 = "j"<br class="">
<br class="">
time = ispan(1,dSizes(3),1) ; CCSM4<br class="">
time!0 = "time"<br class="">
time@long_name = "time"<br class="">
time@units = "month"<br class="">
time&time = time<br class="">
<br class="">
var_regrid&time = time<br class="">
var_regrid&i = lat2d&i<br class="">
var_regrid&j = lat2d&j<br class="">
<br class="">
delete( var_regrid@lat1d )<br class="">
delete( var_regrid@lon1d )<br class="">
if ( isatt( var_regrid, "lat2d" ) ) then<br class="">
delete( var_regrid@lat2d )<br class="">
end if<br class="">
if ( isatt( var_regrid, "lon2d" ) ) then<br class="">
delete( var_regrid@lon2d )<br class="">
end if<br class="">
<br class="">
;;----------------------------<wbr class="">------------------------------<wbr class="">------------<br class="">
;; Write regridded data to file<br class="">
;;----------------------------<wbr class="">------------------------------<wbr class="">------------<br class="">
str = str_split(srcGridFili, ".")<br class="">
filo = str(0) + ".<a href="http://secure-web.cisco.com/1qp3lTKKuGOiRenNSV2mBVPtaoJvr4FX79fNvr1SHNiH45LiiRgo791rBK5WvVJ328_9O9negnRBCXH87H-TT4hFZFC1Jx2Gv7S--nVCvs6J0xREDrJzfEWPXiFW0ZsKtPsAQY3KT5kCR2tSWv9ApyU_7pLbFlMoQoVPcHXKt28S2KbzbcZLHQW2uwsMY_9X1pYYswOB1CkCT8ykYCmwxrzjA3upf-XiQssRxyO-Hcq3GROvwl6HS6cRm8ADD85jf4dqfHXti59WTN27GCubmYXCahfUApc1vwD37weD2iAORkRo21fwXdt2L5agRiKwr_942CeBi4j5MClyorMHhE3J-G1fR5XRcSN3wOaf09LYw2mDSo2Lv64Q-tgRyuNz5JQkuXP8nMW1nYnEeI56nyx4oAO94r5Pl1Wt9GdAME62yUUhdRJKL9xMorz0xxBS1/http%3A%2F%2Frgrd.nc" rel="noreferrer" target="_blank" class="">rgrd.nc</a>"<br class="">
<br class="">
if ( isfilepresent(srcGridDiri+"/"+<wbr class="">filo) .eq. True ) then<br class="">
system( "rm " + filo )<br class="">
end if<br class="">
<br class="">
setfileoption("nc","Format","<wbr class="">LargeFile")<br class="">
fo = addfile( srcGridDiri + "/" + filo,"c")<br class="">
fo->var = var_regrid<br class="">
fo->lat = lat2d<br class="">
fo->lon = lon2d<br class="">
<br class="">
;; Clean up<br class="">
delete( var )<br class="">
delete( sfile )<br class="">
delete( Opt )<br class="">
delete( filo )<br class="">
delete( fo )<br class="">
delete( var_regrid )<br class="">
<br class="">
end do<br class="">
<br class="">
delete( dfile )<br class="">
delete( lat2d )<br class="">
delete( lon2d )<br class="">
<br class="">
end<br class="">
<br class="">
<br class="">
______________________________<wbr class="">_________________<br class="">
ncl-talk mailing list<br class="">
<a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">
List instructions, subscriber options, unsubscribe:<br class="">
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/<wbr class="">mailman/listinfo/ncl-talk</a><br class="">
</blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></body></html>