Local emission mapping to WRF-Chem

[sstikle]

Dear All,
I am trying to map local emissions to WRF-Chem. The projection of local emission is in UTM, how to convert to GCS. Here is the ncdump of the file showing some details below. I tried with cdo but hard luck. Can any one help to resolve this issue to make the nc file compatible to use in the anthro emis tool.

ncdump -h HH_final_emis.nc
netcdf HH_final_emis {
dimensions:
POINT_Y = 166 ;
transverse_mercator = UNLIMITED ; // (0 currently)
POINT_X = 203 ;
variables:
double POINT_Y(POINT_Y) ;
POINT_Y:missing_value = -9999. ;
POINT_Y:long_name = "POINT_Y" ;
POINT_Y:standard_name = "POINT_Y" ;
POINT_Y:grid_mapping = "transverse_mercator" ;
int transverse_mercator(transverse_mercator) ;
transverse_mercator:grid_mapping_name = "transverse_mercator" ;
transverse_mercator:longitude_of_central_meridian = 81. ;
transverse_mercator:latitude_of_projection_origin = 0. ;
transverse_mercator:scale_factor_at_central_meridian = 0.9996 ;
transverse_mercator:false_easting = 500000. ;
transverse_mercator:false_northing = 0. ;
double POINT_X(POINT_X) ;
POINT_X:missing_value = -9999. ;
POINT_X:long_name = "POINT_X" ;
POINT_X:standard_name = "POINT_X" ;
POINT_X:grid_mapping = "transverse_mercator" ;
double x(POINT_Y, POINT_X) ;
x:missing_value = -9999. ;
x:long_name = "x coordinate of projection" ;
x:standard_name = "projection_x_coordinate" ;
x:units = "Meter" ;
x:grid_mapping = "transverse_mercator" ;
double y(POINT_Y, POINT_X) ;
y:missing_value = -9999. ;
y:long_name = "y coordinate of projection" ;
y:standard_name = "projection_y_coordinate" ;
y:units = "Meter" ;
y:grid_mapping = "transverse_mercator" ;
double PM_2_5_Emi(POINT_Y, POINT_X) ;
PM_2_5_Emi:missing_value = -9999. ;
PM_2_5_Emi:long_name = "PM_2_5_Emi" ;
PM_2_5_Emi:standard_name = "PM_2_5_Emi" ;
PM_2_5_Emi:coordinates = "x y" ;
PM_2_5_Emi:esri_pe_string = "PROJCS[\"WGS_1984_UTM_Zone_44N\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",81.0],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]" ;
PM_2_5_Emi:grid_mapping = "transverse_mercator" ;
double OC_Emissio(POINT_Y, POINT_X) ;
OC_Emissio:missing_value = -9999. ;
OC_Emissio:long_name = "OC_Emissio" ;
OC_Emissio:standard_name = "OC_Emissio" ;
OC_Emissio:coordinates = "x y" ;
OC_Emissio:esri_pe_string = "PROJCS[\"WGS_1984_UTM_Zone_44N\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",81.0],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]" ;
OC_Emissio:grid_mapping = "transverse_mercator" ;


Best,
Suvarna

 

[sstikle]

Hi all,
Sorry for the repost.
I am still struggling with conversion of nonregular grid to regular grid with mercator projection. Grid size is 1kmx1km. I have used ARCGIS and reprojected file from WGS_UTM to WGS_GCS, but not able to change it to mercator. Anyone suggest a way to do it? It will be really great help. Thanks in advance.

Best,
Suvarna