<div dir="ltr"><div dir="ltr"><div>Bonjour Adriana,</div><div><br></div><div> The <b>area_conserve_remap </b>and <b>area_conserve_remap_Wrap</b> have been <span style="color:rgb(255,0,0)"><b>deprecated</b></span>. <br></div><div>The<a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml"><b> ESMF regrid</b></a> library is the recommended replacement. It uses a superior methodology.<br></div><div><br></div><div>The <b>source</b> and <b>destination</b> grids are both <b>rectilinear</b>. <br></div><div><br></div><div>I modified one of the rectilinear-to-rectilinear grid examples.</div><div><br></div><div>Attached is a sample script that uses ESMF conservative remapping. It creates:</div><div><br></div><div>[1] An ESMF weight file that can be used on subsequent regrid variables.</div><div>[2] A netCDF file containing the regridded emissions variable.</div><div>[3] A panel plot showing the source and regridded variables at an arbitrary level and time. <br></div><div> NOTE: Volcano emission(s) are point source(s) so the plots won't show much.</div><div><br></div><div>You can copy and modify the attached file to interpolate (regrid) to different resolutions.</div><div><br></div><div>====</div><div><br></div><div><b>%></b> ncl adriana.esmf_regrid_1.ncl<br></div><div><br></div><div>Cheers</div><div>D<br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 28, 2019 at 8:05 AM <<a href="mailto:adriana.bossolasco@cnrs-orleans.fr">adriana.bossolasco@cnrs-orleans.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all!<br>
<br>
I write because I don't understand what happen now with my regridding. I <br>
am trying to regrid an emission file from 0.9x1.9 lat/lon to a grid <br>
1.9x2.5 lat/lon.<br>
I did with successful before, but now I have a file with more verticals <br>
levels and doesn't work and I get the following Error message:<br>
<br>
Error in "cremapbin":<br>
Could not map global lat array into grid array<br>
Aborted (core dumped)<br>
<br>
I will copy part of my code:<br>
<br>
<br>
;*************************************************************<br>
; read<br>
;*************************************************************<br>
<br>
srcFileName = <br>
"<a href="http://emissions-cmip6_so4_a1_contvolcano_vertical_2000-2015.nc" rel="noreferrer" target="_blank">emissions-cmip6_so4_a1_contvolcano_vertical_2000-2015.nc</a>"<br>
f = addfile(srcFileName, "r")<br>
<br>
x = f->emiss_volcanoes<br>
<br>
printVarSummary(x)<br>
printMinMax(x, True)<br>
;*************************************************************<br>
; get dimension sizes and other information<br>
;*************************************************************<br>
dimx = dimsizes( x )<br>
ntim = dimx( 0 )<br>
nlev = dimx( 1 )<br>
nlat = dimx( 2 ) ; [0.9 x 1.9 ]<br>
mlon = dimx( 3 ) ;<br>
<br>
latS = x&lat(0) ; south extent of input grid<br>
latN = x&lat(nlat-1) ; north extent<br>
<br>
dimx = dimsizes( x )<br>
ntim = dimx( 0 )<br>
nlev = dimx( 1 )<br>
nlat = dimx( 2 ) ;<br>
mlon = dimx( 3 ) ;<br>
<br>
latS = x&lat(0) ; south extent of input grid<br>
latN = x&lat(nlat-1) ; north extent<br>
<br>
<br>
nlati = 192 ; global extent of input grid <br>
[180/0.5]<br>
<br>
;*************************************************************<br>
; perform conservative remapping to two different grid resolutions<br>
;*************************************************************<br>
<br>
opt = True<br>
<br>
<br>
; ---<br>
NLAT25x19= 96 ; RES = "2.5x1.9"<br>
MLON25x19= 144<br>
LAT25x19 = latGau (NLAT25x19, "lat", "latitude" , "degrees_north")<br>
LON25x19 = lonGlobeFo(MLON25x19, "lon", "longitude", "degrees_east" )<br>
<br>
opt@NLATo = NLAT25x19<br>
LAT_REG25x19= LAT25x19({latS:latN})<br>
X25x19 = area_conserve_remap_Wrap (x&lon, x&lat, x ,LON25x19, <br>
LAT_REG25x19, opt)<br>
<br>
<br>
<br>
printVarSummary(X25x19)<br>
printMinMax(X25x19, True)<br>
<br>
<br>
And I get this Summary of my variables<br>
<br>
Variable: x<br>
Type: float<br>
Total Size: 1274019840 bytes<br>
318504960 values<br>
Number of Dimensions: 4<br>
Dimensions and sizes: [time | 192] x [altitude | 30] x [lat | 192] x <br>
[lon | 288]<br>
Coordinates:<br>
time: [91325..97138]<br>
altitude: [0.125..7.375]<br>
lat: [-90..90]<br>
lon: [ 0..358.75]<br>
Number Of Attributes: 3<br>
units : molecules/cm3/s<br>
long_name : Continuously outgassing volcano SO4 = (0.025*SO2)*0.5 <br>
emissions<br>
molecular_weight : 115<br>
(0)<br>
(0) Continuously outgassing volcano SO4 = (0.025*SO2)*0.5 emissions <br>
(molecules/cm3/s) : min=0 max=2.19543e+06<br>
<br>
<br>
Variable: latS<br>
Type: float<br>
Total Size: 4 bytes<br>
1 values<br>
Number of Dimensions: 1<br>
Dimensions and sizes: [1]<br>
Coordinates:<br>
Number Of Attributes: 2<br>
units : degrees_north<br>
long_name : Latitude<br>
(0) -90<br>
<br>
<br>
Variable: latN<br>
Type: float<br>
Total Size: 4 bytes<br>
1 values<br>
Number of Dimensions: 1<br>
Dimensions and sizes: [1]<br>
Coordinates:<br>
Number Of Attributes: 2<br>
units : degrees_north<br>
long_name : Latitude<br>
(0) 90<br>
Error in "cremapbin":<br>
Could not map global lat array into grid array<br>
Aborted (core dumped)<br>
<br>
<br>
<br>
I really appreciate any help. I tried a lot of things but nothing work.<br>
<br>
For the other hand I run the same code for other file and works. The <br>
other file have 8 altitude levels and this 30. So I was thinking maybe <br>
is to big the file, but I try to select just two o three levels and I <br>
get the same error.<br>
<br>
Thanks<br>
<br>
Adriana<br>
-- <br>
Adriana Bossolasco<br>
LPC2E UMR 7328 (CNRS, Université d'Orléans)<br>
3A avenue de la Recherche Scientifique<br>
45071 Orléans Cedex 2 - FRANCE<br>
Tel +33 (0)2 38 25 52 85<br>
e-mail: <a href="mailto:adriana.bossolasco@cnrs-orleans.fr" target="_blank">adriana.bossolasco@cnrs-orleans.fr</a><br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>