<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Hi</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">I need to understand a fundamental feature of the ESMF "conserve" regridding: does it conserve fluxes or absolute values?</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">I've got a <b>population count</b> and <b>population density</b> data on a native 1deg x 1deg grid. Ultimately, I need a population count on a 1.9deg x 2.5deg grid. I tried to achieve this in two ways:</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">1. Regridding the population count directly: This gives me a grid sum of 1.67 billion instead of 7.96 billion (world population in 2020) as obtained from the native grid sum. Strangely enough, if I multiply 1.67 billion with 1.9*2.5, the result gets pretty close to 7.96 billion but I'm unable to interpret why that might be the case.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">2. Regridding population density data and multiplying by a matrix of gridpoint areas (I have obtained this separately) to get total population. In this case, the grid sum turns out to be 13.5 billion which is significantly higher than the actual world population of 7.96 billion again.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">So, my first question is, which is the right way to regrid population when using conservative remapping? And further, are there any obvious issues in my approach?</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Here's my code to generate weights file:</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><font color="#0000ff"><i>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"<br><br>begin<br>    sfile = addfile("/work/users/tan/gpw-v4-population-density-rev11_totpop_1_deg_nc/<a href="http://gpw_v4_population_density_rev11_1_deg.nc">gpw_v4_population_density_rev11_1_deg.nc</a>","r")<br>    dfile = addfile("/data/sync/modelinput/CAMS-Emissions-tbu/<a href="http://HTAP_transient_1.9x2.5_emis_NO_20201217.nc">HTAP_transient_1.9x2.5_emis_NO_20201217.nc</a>","r")<br><br>    varname = "Population Density, v4.11 (2000, 2005, 2010, 2015, 2020): 1 degree" ; A TERRIBLE VARIABLE NAME!<br>    pop     = sfile->$varname$<br>    pop2020 = pop(4,:,:) ; EXTRACTING 2D POPULATION DATA FOR THE YEAR 2020<br><br>    nox_cams = dfile->glseas(0,:,:)<br><br>    poplat = pop2020&latitude<br>    poplon = pop2020&longitude<br><br>    camslat = nox_cams&lat<br>    camslon = nox_cams&lon<br><br><br>   Opt                = True<br>   Opt@ForceOverwrite = True<br>   Opt@InterpMethod   = "conserve"<br>   Opt@SrcGridLon        = poplon<br>   Opt@SrcGridLat        = poplat<br>   Opt@DstGridLon        = camslon<br>   Opt@DstGridLat        = camslat<br><br>   pop_regrid = ESMF_regrid(pop2020,Opt)<br>   printVarSummary(pop_regrid)<br><br>end<br></i></font></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Here's the code to generate the regridded data using the newly created weights file:</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><font color="#0000ff"><i>begin<br><br>sfile = addfile("/work/users/tan/gpw-v4-population-density-rev11_totpop_1_deg_nc/<a href="http://gpw_v4_population_density_rev11_1_deg.nc">gpw_v4_population_density_rev11_1_deg.nc</a>","r")<br>dfile = addfile("/work/users/tan/gpw-v4-population-density-rev11_totpop_1_deg_nc/<a href="http://popdensity2020_regridded.nc">popdensity2020_regridded.nc</a>","c")<br>varname = "Population Density, v4.11 (2000, 2005, 2010, 2015, 2020): 1 degree" <br>pop     = sfile->$varname$<br>pop2020 = pop(4,:,:)<br><br>WF = "/work/users/tan/nclscripts/pop-regridfiles/<a href="http://weights_file.nc">weights_file.nc</a>"<br><br>pop2020_regrid = ESMF_regrid_with_weights(pop2020,WF,False) ;REGRIDDING<br>dfile->popdensity2020 = pop2020_regrid  ;STORING<br>print("Regridded!")</i></font></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><font color="#0000ff"><i><br>end<br></i></font></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Thanks a lot.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">best regards,</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Tabish</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><font face="trebuchet ms, sans-serif" color="#666666">-------------------------------<br><span style="background-color:rgb(255,255,255)">Dr Tabish Ansari</span></font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Research Associate </font></div><div><font face="trebuchet ms, sans-serif" color="#666666">Air Quality Modelling Group</font></div><div><span style="background-color:rgb(255,255,255)"><span style="font-weight:normal"><font face="trebuchet ms, sans-serif" color="#666666">IASS-Potsdam </font></span></span></div><div><font face="trebuchet ms, sans-serif" color="#666666">Germany</font></div></div></div></div></div>