Project

General

Profile

why the value is much smaller after remapbil and how to correct it ?

Added by allen zhang almost 6 years ago

I have an regional emission inventory of 0.002deg spatial resolution, whose header file is shown as follows:

netcdf INDY.Total.2012.hourly {
dimensions:
Lon = 560 ;
Lat = 521 ;
Time = 8760 ;
variables:
double Lon(Lon) ;
Lon:units = "degree" ;
Lon:long_name = "Lon" ;
double Lat(Lat) ;
Lat:units = "degree" ;
Lat:long_name = "Lat" ;
int Time(Time) ;
Time:units = "hour" ;
Time:long_name = "Time" ;
double Carbon\ Emission(Time, Lat, Lon) ;
Carbon\ Emission:units = "Kilogram" ;
Carbon\ Emission:_FillValue = NaN ;
}

and I want to match it to a WRF domain, the grid definition is attached as grid03.
For original emission inventory file, I have modified something so that it can be processed by remapbil

netcdf \2012_hourly_emission_new {
dimensions:
Lon = 560 ;
Lat = 521 ;
time = UNLIMITED ; // (8760 currently)
variables:
double Lon(Lon) ;
Lon:standard_name = "longitude" ;
Lon:long_name = "Lon" ;
Lon:units = "degree_east" ;
double Lat(Lat) ;
Lat:standard_name = "latitude" ;
Lat:long_name = "Lat" ;
Lat:units = "degree_north" ;
double time(time) ;
time:standard_name = "time" ;
time:units = "hours since 2012-01-01 00:00:00" ;
time:calendar = "proleptic_gregorian" ;
double Carbon\ Emission(time, Lat, Lon) ;
Carbon\ Emission:units = "Kilogram" ;
Carbon\ Emission:_FillValue = NaN ;
Carbon\ Emission:missing_value = NaN ;

// global attributes:
:CDI = "Climate Data Interface version 1.6.2 (http://code.zmaw.de/projects/cdi)" ;
:Conventions = "CF-1.4" ;
:history = "Sat Jun 9 23:54:34 2018: ncatted -O -a axis,Lon,d,, -a axis,Lat,d,, 2012_hourly_emission_new.nc\n",
"Sat Jun 09 23:12:08 2018: cdo settaxis,2012-01-01,00:00:00,1hour -settunits,hours 2012_hourly_emission.nc 2012_hourly_emission_new.nc\n",
"Sat Jun 9 22:00:51 2018: ncatted -a units,Lon,m,c,degree_east -a units,Lat,m,c,degree_north 2012_hourly_emission.nc\n",
"Sat Jun 9 15:44:52 2018: ncatted -a units,Lon,m,c,degree_east -a units,Lat,m,c,degree_notth 2012_hourly_emission.nc" ;
:CDO = "Climate Data Operators version 1.6.2 (http://code.zmaw.de/projects/cdo)" ;

After remapbil, we get the emission file corresponding to the WRF domain,

netcdf emission_d03_2012_hourly_remapbil {
dimensions:
x = 87 ;
y = 87 ;
time = UNLIMITED ; // (8760 currently)
variables:
double XLONG ;
XLONG:standard_name = "longitude" ;
XLONG:long_name = "longitude" ;
XLONG:units = "degree_east" ;
XLONG:_CoordinateAxisType = "Lon" ;
double XLAT ;
XLAT:standard_name = "latitude" ;
XLAT:long_name = "latitude" ;
XLAT:units = "degree_north" ;
XLAT:_CoordinateAxisType = "Lat" ;
double time(time) ;
time:standard_name = "time" ;
time:units = "hours since 2012-01-01 00:00:00" ;
time:calendar = "proleptic_gregorian" ;
double Carbon\ Emission(time, y, x) ;
Carbon\ Emission:units = "Kilogram" ;
Carbon\ Emission:coordinates = "XLONG XLAT" ;
Carbon\ Emission:_FillValue = NaN ;
Carbon\ Emission:missing_value = NaN ;

// global attributes:
:CDI = "Climate Data Interface version 1.6.2 (http://code.zmaw.de/projects/cdi)" ;
:Conventions = "CF-1.4" ;
:history = "Sat Jun 09 23:56:50 2018: cdo remapbil,grid03 2012_hourly_emission_new.nc emission_d03_2012_hourly_remapbil.nc\n",
"Sat Jun 9 23:54:34 2018: ncatted -O -a axis,Lon,d,, -a axis,Lat,d,, 2012_hourly_emission_new.nc\n",
"Sat Jun 09 23:12:08 2018: cdo settaxis,2012-01-01,00:00:00,1hour -settunits,hours 2012_hourly_emission.nc 2012_hourly_emission_new.nc\n",
"Sat Jun 9 22:00:51 2018: ncatted -a units,Lon,m,c,degree_east -a units,Lat,m,c,degree_north 2012_hourly_emission.nc\n",
"Sat Jun 9 15:44:52 2018: ncatted -a units,Lon,m,c,degree_east -a units,Lat,m,c,degree_notth 2012_hourly_emission.nc" ;
:CDO = "Climate Data Operators version 1.6.2 (http://code.zmaw.de/projects/cdo)" ;

But for data of emission,it seems to be strange since the original emission are between 0 and 169857 Kg but the after are between 0 and 2271 Kg
I don't know how to correct it or may be this process is wrong , can anybody give me some advice?
I also attached the sinfov file of original and after some process for references since the emission file is too large to upload.


Replies (1)

RE: why the value is much smaller after remapbil and how to correct it ? - Added by allen zhang almost 6 years ago

ADD some information
the griddes of original emission file is

#
  1. gridID 2 #
    gridtype = generic
    gridsize = 291760
    xname = Lon
    xlongname = Lon
    xunits = degree
    yname = Lat
    ylongname = Lat
    yunits = degree
    xsize = 560
    ysize = 521
    xfirst = -86.694
    xinc = 0.002
    yfirst = 40.3788
    yinc = -0.002
    ~

since generic gridtype cannot be remapilled, and I have modified it and get

#
  1. gridID 2 #
    gridtype = lonlat
    gridsize = 291760
    xname = Lon
    xlongname = Lon
    xunits = degree_east
    yname = Lat
    ylongname = Lat
    yunits = degree_north
    xsize = 560
    ysize = 521
    xfirst = -86.694
    xinc = 0.002
    yfirst = 40.3788
    yinc = -0.002
    ~
    (1-1/1)