Hi,
I have downloaded CMIP5 ocean temp. data
(thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-284912.nc)
but find that the horizontal grid is not regular grid used for GRADS
(see below)
Therefor I think I need to regrid the ocean temp. data from these
irregular grid points to the regular 1X1 grid.
I am asking if NCL can do the interpolation? If so, Could you point
out to me the key sentences used in NCL scripts?
Thanks,
Tao
--------------------
netcdf thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-284912 {
dimensions:
i = 256 ;
j = 220 ;
vertices = 4 ;
lev = 19 ;
time = UNLIMITED ; // (600 currently)
variables:
int i(i) ;
i:units = "1" ;
i:long_name = "cell index along first dimension" ;
int j(j) ;
j:units = "1" ;
j:long_name = "cell index along second dimension" ;
float lat(j, i) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude coordinate" ;
lat:units = "degrees_north" ;
lat:bounds = "lat_vertices" ;
float lat_vertices(j, i, vertices) ;
lat_vertices:units = "degrees_north" ;
double lev(lev) ;
lev:bounds = "lev_bnds" ;
lev:units = "m" ;
lev:axis = "Z" ;
lev:positive = "down" ;
lev:long_name = "ocean depth coordinate" ;
lev:standard_name = "depth" ;
float lon(j, i) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude coordinate" ;
lon:units = "degrees_east" ;
lon:bounds = "lon_vertices" ;
float lon_vertices(j, i, vertices) ;
lon_vertices:units = "degrees_east" ;
float thetao(time, lev, j, i) ;
thetao:standard_name = "sea_water_potential_temperature" ;
thetao:long_name = "Sea Water Potential Temperature" ;
thetao:units = "K" ;
thetao:cell_methods = "time: mean" ;
thetao:cell_measures = "area: areacello volume: volcello" ;
thetao:missing_value = 1.e+20f ;
thetao:_FillValue = 1.e+20f ;
thetao:associated_files = "baseURL:
http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile:
gridspec_ocean_fx_MPI-ESM-LR_piControl_r0i0p0.nc areacello:
areacello_fx_MPI-ESM-LR_piControl_r0i0p0.nc volcello:
volcello_fx_MPI-ESM-LR_piControl_r0i0p0.nc" ;
thetao:coordinates = "lat lon" ;
double time(time) ;
time:bounds = "time_bnds" ;
time:units = "days since 1850-1-1 00:00:00" ;
time:calendar = "proleptic_gregorian" ;
time:axis = "T" ;
time:long_name = "time" ;
time:standard_name = "time" ;
// global attributes:
:institution = "Max Planck Institute for Meteorology" ;
:institute_id = "MPI-M" ;
:experiment_id = "piControl" ;
:source = "MPI-ESM-LR 2011; URL:
http://svn.zmaw.de/svn/cosmos/branches/releases/mpi-esm-cmip5/src/mod;
atmosphere: ECHAM6 (REV: 4556), T63L47; land: JSBACH (REV: 4556); ocean:
MPIOM (REV: 4556), GR15L40; sea ice: 4556; marine bgc: HAMOCC (REV:
4556);" ;
:model_id = "MPI-ESM-LR" ;
:forcing = "N/A" ;
:parent_experiment_id = "N/A" ;
:parent_experiment_rip = "N/A" ;
:branch_time = 0. ;
:contact = "cmip5-mpi-esm@dkrz.de" ;
:history = "Sun Nov 27 12:44:50 2011: ncrcat
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-280912.new.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_281001-281912.new.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_282001-282912.new.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_283001-283912.new.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_284001-284912.new.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/proc/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-284912.nc\n",
"Sun Nov 27 12:43:41 2011: ncks -c -F -d
lev,1,19 -v lat,lon,lat_vertices,lon_vertices,thetao
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/raw/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-280912.nc
/Volumes/tzhang5/tzhang5/CMIP5/ocntmp/work/thetao_Omon_MPI-ESM-LR_piControl_r1i1p1_280001-280912.new.nc\n",
"Model raw output postprocessing with modelling
environment (IMDI) at DKRZ: URL:
http://svn-mad.zmaw.de/svn/mad/Model/IMDI/trunk, REV: 3315
2011-06-28T12:24:00Z CMOR rewrote data to comply with CF standards and
CMIP5 requirements." ;
:comment = "started in year 3424 of tra0066 spinup run." ;
:references = "ECHAM6: n/a; JSBACH: Raddatz et al.,
2007. Will the tropical land biosphere dominate the climate-carbon cycle
feedback during the twenty first century? Climate Dynamics, 29, 565-574,
doi 10.1007/s00382-007-0247-8; MPIOM: Marsland et al., 2003. The
Max-Planck-Institute global ocean/sea ice model with orthogonal
curvilinear coordinates. Ocean Modelling, 5, 91-127; HAMOCC:
http://www.mpimet.mpg.de/fileadmin/models/MPIOM/HAMOCC5.1_TECHNICAL_REPORT.pdf;"
;
:initialization_method = 1 ;
:physics_version = 1 ;
:tracking_id = "dfa793a7-2d2d-402b-8fe2-f6599393a311" ;
:product = "output" ;
:experiment = "pre-industrial control" ;
:frequency = "mon" ;
:creation_date = "2011-06-28T12:24:22Z" ;
:Conventions = "CF-1.4" ;
:project_id = "CMIP5" ;
:table_id = "Table Omon (27 April 2011)
340eddd4fd838d90fa9ffe1345ecbd73" ;
:title = "MPI-ESM-LR model output prepared for CMIP5
pre-industrial control" ;
:parent_experiment = "N/A" ;
:modeling_realm = "ocean" ;
:realization = 1 ;
:cmor_version = "2.5.9" ;
:NCO = "4.0.0" ;
:nco_openmp_thread_number = 1 ;
data:
i = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
20,
21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37,
38,
39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55,
56,
57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73,
74,
75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91,
92,
93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107,
108,
109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122,
123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136,
137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150,
151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164,
165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178,
179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192,
193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206,
207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220,
221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234,
235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248,
249, 250, 251, 252, 253, 254, 255, 256 ;
j = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
20,
21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37,
38,
39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55,
56,
57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73,
74,
75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91,
92,
93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107,
108,
109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122,
123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136,
137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150,
151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164,
165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178,
179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192,
193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206,
207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220 ;
lev = 6, 17, 27, 37, 47, 57, 68.5, 82.5, 100, 122.5, 150, 182.5, 220,
262.5,
310, 362.5, 420, 485, 560 ;
time = 346995.5, 347025.5, 347055.5, 347086, 347116.5, 347147, 347177.5,
347208.5, 347239, 347269.5, 347300, 347330.5, 347361.5, 347391,
347420.5,
347451, 347481.5, 347512, 347542.5, 347573.5, 347604, 347634.5,
347665,
347695.5, 347726.5, 347756, 347785.5, 347816, 347846.5, 347877,
347907.5,
347938.5, 347969, 347999.5, 348030, 348060.5, 348091.5, 348121,
348150.5,
348181, 348211.5, 348242, 348272.5, 348303.5, 348334, 348364.5,
348395,
348425.5, 348456.5, 348486.5, 348516.5, 348547, 348577.5, 348608,
348638.5, 348669.5, 348700, 348730.5, 348761, 348791.5, 348822.5,
348852,
348881.5, 348912, 348942.5, 348973, 349003.5, 349034.5, 349065,
349095.5,
349126, 349156.5, 349187.5, 349217, 349246.5, 349277, 349307.5,
349338,
349368.5, 349399.5, 349430, 349460.5, 349491, 349521.5, 349552.5,
349582,
349611.5, 349642, 349672.5, 349703, 349733.5, 349764.5, 349795,
349825.5,
349856, 349886.5, 349917.5, 349947.5, 349977.5, 350008, 350038.5,
350069,
350099.5, 350130.5, 350161, 350191.5, 350222, 350252.5, 350283.5,
350313,
350342.5, 350373, 350403.5, 350434, 350464.5, 350495.5, 350526,
350556.5,
350587, 350617.5, 350648.5, 350678, 350707.5, 350738, 350768.5,
350799,
350829.5, 350860.5, 350891, 350921.5, 350952, 350982.5, 351013.5,
351043,
351072.5, 351103, 351133.5, 351164, 351194.5, 351225.5, 351256,
351286.5,
351317, 351347.5, 351378.5, 351408.5, 351438.5, 351469, 351499.5,
351530,
351560.5, 351591.5, 351622, 351652.5, 351683, 351713.5, 351744.5,
351774,
351803.5, 351834, 351864.5, 351895, 351925.5, 351956.5, 351987,
352017.5,
352048, 352078.5, 352109.5, 352139, 352168.5, 352199, 352229.5,
352260,
352290.5, 352321.5, 352352, 352382.5, 352413, 352443.5, 352474.5,
352504,
352533.5, 352564, 352594.5, 352625, 352655.5, 352686.5, 352717,
352747.5,
352778, 352808.5, 352839.5, 352869.5, 352899.5, 352930, 352960.5,
352991,
353021.5, 353052.5, 353083, 353113.5, 353144, 353174.5, 353205.5,
353235,
353264.5, 353295, 353325.5, 353356, 353386.5, 353417.5, 353448,
353478.5,
353509, 353539.5, 353570.5, 353600, 353629.5, 353660, 353690.5,
353721,
353751.5, 353782.5, 353813, 353843.5, 353874, 353904.5, 353935.5,
353965,
353994.5, 354025, 354055.5, 354086, 354116.5, 354147.5, 354178,
354208.5,
354239, 354269.5, 354300.5, 354330.5, 354360.5, 354391, 354421.5,
354452,
354482.5, 354513.5, 354544, 354574.5, 354605, 354635.5, 354666.5,
354696,
354725.5, 354756, 354786.5, 354817, 354847.5, 354878.5, 354909,
354939.5,
354970, 355000.5, 355031.5, 355061, 355090.5, 355121, 355151.5,
355182,
355212.5, 355243.5, 355274, 355304.5, 355335, 355365.5, 355396.5,
355426,
355455.5, 355486, 355516.5, 355547, 355577.5, 355608.5, 355639,
355669.5,
355700, 355730.5, 355761.5, 355791.5, 355821.5, 355852, 355882.5,
355913,
355943.5, 355974.5, 356005, 356035.5, 356066, 356096.5, 356127.5,
356157,
356186.5, 356217, 356247.5, 356278, 356308.5, 356339.5, 356370,
356400.5,
356431, 356461.5, 356492.5, 356522, 356551.5, 356582, 356612.5,
356643,
356673.5, 356704.5, 356735, 356765.5, 356796, 356826.5, 356857.5,
356887,
356916.5, 356947, 356977.5, 357008, 357038.5, 357069.5, 357100,
357130.5,
357161, 357191.5, 357222.5, 357252.5, 357282.5, 357313, 357343.5,
357374,
357404.5, 357435.5, 357466, 357496.5, 357527, 357557.5, 357588.5,
357618,
357647.5, 357678, 357708.5, 357739, 357769.5, 357800.5, 357831,
357861.5,
357892, 357922.5, 357953.5, 357983, 358012.5, 358043, 358073.5,
358104,
358134.5, 358165.5, 358196, 358226.5, 358257, 358287.5, 358318.5,
358348,
358377.5, 358408, 358438.5, 358469, 358499.5, 358530.5, 358561,
358591.5,
358622, 358652.5, 358683.5, 358713.5, 358743.5, 358774, 358804.5,
358835,
358865.5, 358896.5, 358927, 358957.5, 358988, 359018.5, 359049.5,
359079,
359108.5, 359139, 359169.5, 359200, 359230.5, 359261.5, 359292,
359322.5,
359353, 359383.5, 359414.5, 359444, 359473.5, 359504, 359534.5,
359565,
359595.5, 359626.5, 359657, 359687.5, 359718, 359748.5, 359779.5,
359809,
359838.5, 359869, 359899.5, 359930, 359960.5, 359991.5, 360022,
360052.5,
360083, 360113.5, 360144.5, 360174.5, 360204.5, 360235, 360265.5,
360296,
360326.5, 360357.5, 360388, 360418.5, 360449, 360479.5, 360510.5,
360540,
360569.5, 360600, 360630.5, 360661, 360691.5, 360722.5, 360753,
360783.5,
360814, 360844.5, 360875.5, 360905, 360934.5, 360965, 360995.5,
361026,
361056.5, 361087.5, 361118, 361148.5, 361179, 361209.5, 361240.5,
361270,
361299.5, 361330, 361360.5, 361391, 361421.5, 361452.5, 361483,
361513.5,
361544, 361574.5, 361605.5, 361635.5, 361665.5, 361696, 361726.5,
361757,
361787.5, 361818.5, 361849, 361879.5, 361910, 361940.5, 361971.5,
362001,
362030.5, 362061, 362091.5, 362122, 362152.5, 362183.5, 362214,
362244.5,
362275, 362305.5, 362336.5, 362366, 362395.5, 362426, 362456.5,
362487,
362517.5, 362548.5, 362579, 362609.5, 362640, 362670.5, 362701.5,
362731,
362760.5, 362791, 362821.5, 362852, 362882.5, 362913.5, 362944,
362974.5,
363005, 363035.5, 363066.5, 363096.5, 363126.5, 363157, 363187.5,
363218,
363248.5, 363279.5, 363310, 363340.5, 363371, 363401.5, 363432.5,
363462,
363491.5, 363522, 363552.5, 363583, 363613.5, 363644.5, 363675,
363705.5,
363736, 363766.5, 363797.5, 363827, 363856.5, 363887, 363917.5,
363948,
363978.5, 364009.5, 364040, 364070.5, 364101, 364131.5, 364162.5,
364192,
364221.5, 364252, 364282.5, 364313, 364343.5, 364374.5, 364405,
364435.5,
364466, 364496.5, 364527.5, 364557.5, 364587.5, 364618, 364648.5,
364679,
364709.5, 364740.5, 364771, 364801.5, 364832, 364862.5, 364893.5,
364923,
364952.5, 364983, 365013.5, 365044, 365074.5, 365105.5, 365136,
365166.5,
365197, 365227.5 ;
}
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Nov 27 14:25:56 2011
This archive was generated by hypermail 2.1.8 : Wed Nov 30 2011 - 19:52:47 MST