Project

General

Profile

how to switch the irrgular model output to the regular data ?

Added by yong sun almost 12 years ago

Hi,

I am recently handling the ocean variable produced by the IPSL model (submitted to the CMIP5). I found it is not the irregular output.

How to convert it to the regular one? I have a try to use the operator: cdo remapbil,r360x180 -selname,thetao in.nc .out.nc

as a relsult, It does not work. the error information is :

cdo remapbil: Started child process "selname,thetao thetao_Omon_IPSL-CM5A-LR_rcp45_r1i1p1_200601-205512.nc (pipe1.1)".
Warning (scanVarAttributes) : NetCDF: Variable not found - areacello
Warning (define_all_grids) : Unsupported array structure, skipped variable thetao!

cdo(2) selname: Open failed on >thetao_Omon_IPSL-CM5A-LR_rcp45_r1i1p1_200601-205512.nc<
Unsupported file structure

To my surprise, my friends employ the same method and same input file succeed.

My CDO version is 1.5.3, Would you please give me some suggestion to solve this problem ?

Thanks in advance.

Have a nice day

YONG

the information for input data file are as follows:

dimensions:
time = UNLIMITED ; // (600 currently)
lev = 31 ;
j = 149 ;
i = 182 ;
bnds = 2 ;
vertices = 4 ;
variables:
double time(time) ;
time:bounds = "time_bnds" ;
time:units = "days since 2006-01-01 00:00:00" ;
time:calendar = "noleap" ;
time:axis = "T" ;
time:long_name = "time" ;
time:standard_name = "time" ;
double time_bnds(time, bnds) ;
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" ;
double lev_bnds(lev, bnds) ;
int j(j) ;
j:units = "1" ;
j:long_name = "cell index along second dimension" ;
int i(i) ;
i:units = "1" ;
i:long_name = "cell index along first dimension" ;
float lat(j, i) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude coordinate" ;
lat:units = "degrees_north" ;
lat:bounds = "lat_vertices" ;
float lon(j, i) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude coordinate" ;
lon:units = "degrees_east" ;
lon:bounds = "lon_vertices" ;
float lat_vertices(j, i, vertices) ;
lat_vertices:units = "degrees_north" ;
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:original_name = "votemper" ;
thetao:original_units = "degC" ;
thetao:history = "2011-07-03T21:01:03Z altered by CMOR: Converted units from \'degC\' to \'K\'. 2011-07-03T21:01:03Z altered by CMOR: replaced missing value flag (9.96921e+36) with standard missing value (1e+20)." ;
thetao:cell_methods = "time: mean (interval: 30 minutes)" ;
thetao:cell_measures = "area: areacello thetao:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_ocean_fx_IPSL-CM5A-LR_rcp45_r0i0p0.nc areacello: areacello_fx_IPSL-CM5A-LR_rcp45_r0i0p0.nc volcello: volcello_fx_IPSL-CM5A-LR_rcp45_r0i0p0.nc" ;
thetao:coordinates = "lat lon" ;

// global attributes:
:institution = "IPSL (Institut Pierre Simon Laplace, Paris, France)" ;
:institute_id = "IPSL" ;
:experiment_id = "rcp45" ;
:source = "IPSL-CM5A-LR (2010) : atmos : LMDZ4 (LMDZ4_v5, 96x95x39); ocean : ORCA2 (NEMOV2_3, 2x2L31); seaIce : LIM2 (NEMOV2_3); ocnBgchem : PISCES (NEMOV2_3); land : ORCHIDEE (orchidee_1_9_4_AR5)" ;
:model_id = "IPSL-CM5A-LR" ;
:forcing = "Nat,Ant,GHG,SA,Oz,LU,SS,Ds,BC,MD,OC,AA" ;
:parent_experiment_id = "historical" ;
:parent_experiment_rip = "r1i1p1" ;
:branch_time = 2005. ;
:contact = "ipsl-cmip5 at ipsl.jussieu.fr Data manager : Sebastien Denvil" ;
:comment = "This 21th century simulation has been forced by prescribed concentration following the rcp 4.5 scenario." ;
:references = "Model documentation and further reference available here : http://icmc.ipsl.fr" ;
:initialization_method = 1 ;
:physics_version = 1 ;
:tracking_id = "4789822a-019c-44b6-bf82-5d172a7a1df4" ;
:product = "output" ;
:experiment = "RCP4.5" ;
:frequency = "mon" ;
:creation_date = "2011-07-03T21:01:18Z" ;
:history = "2011-07-03T21:01:18Z CMOR rewrote data to comply with CF standards and CMIP5 requirements." ;
:Conventions = "CF-1.4" ;
:project_id = "CMIP5" ;
:table_id = "Table Omon (31 January 2011) d2d6beec2b8fea5bbed33920a6e08bbe" ;
:title = "IPSL-CM5A-LR model output prepared for CMIP5 RCP4.5" ;
:parent_experiment = "historical" ;
:modeling_realm = "ocean" ;
:realization = 1 ;
:cmor_version = "2.7.1" ;
data:

time = 15.5, 45, 74.5, 105, 135.5, 166, 196.5, 227.5, 258, 288.5, 319, 
349.5, 380.5, 410, 439.5, 470, 500.5, 531, 561.5, 592.5, 623, 653.5, 684,
714.5, 745.5, 775, 804.5, 835, 865.5, 896, 926.5, 957.5, 988, 1018.5,
1049, 1079.5, 1110.5, 1140, 1169.5, 1200, 1230.5, 1261, 1291.5, 1322.5,
1353, 1383.5, 1414, 1444.5, 1475.5, 1505, 1534.5, 1565, 1595.5, 1626,
1656.5, 1687.5, 1718, 1748.5, 1779, 1809.5, 1840.5, 1870, 1899.5, 1930,
1960.5, 1991, 2021.5, 2052.5, 2083, 2113.5, 2144, 2174.5, 2205.5, 2235,
2264.5, 2295, 2325.5, 2356, 2386.5, 2417.5, 2448, 2478.5, 2509, 2539.5,
2570.5, 2600, 2629.5, 2660, 2690.5, 2721, 2751.5, 2782.5, 2813, 2843.5,
2874, 2904.5, 2935.5, 2965, 2994.5, 3025, 3055.5, 3086, 3116.5, 3147.5,
3178, 3208.5, 3239, 3269.5, 3300.5, 3330, 3359.5, 3390, 3420.5, 3451,
3481.5, 3512.5, 3543, 3573.5, 3604, 3634.5, 3665.5, 3695, 3724.5, 3755,
3785.5, 3816, 3846.5, 3877.5, 3908, 3938.5, 3969, 3999.5, 4030.5, 4060,
4089.5, 4120, 4150.5, 4181, 4211.5, 4242.5, 4273, 4303.5, 4334, 4364.5,
4395.5, 4425, 4454.5, 4485, 4515.5, 4546, 4576.5, 4607.5, 4638, 4668.5,
4699, 4729.5, 4760.5, 4790, 4819.5, 4850, 4880.5, 4911, 4941.5, 4972.5,
5003, 5033.5, 5064, 5094.5, 5125.5, 5155, 5184.5, 5215, 5245.5, 5276,
5306.5, 5337.5, 5368, 5398.5, 5429, 5459.5, 5490.5, 5520, 5549.5, 5580,
5610.5, 5641, 5671.5, 5702.5, 5733, 5763.5, 5794, 5824.5, 5855.5, 5885,
5914.5, 5945, 5975.5, 6006, 6036.5, 6067.5, 6098, 6128.5, 6159, 6189.5,
6220.5, 6250, 6279.5, 6310, 6340.5, 6371, 6401.5, 6432.5, 6463, 6493.5,
6524, 6554.5, 6585.5, 6615, 6644.5, 6675, 6705.5, 6736, 6766.5, 6797.5,
6828, 6858.5, 6889, 6919.5, 6950.5, 6980, 7009.5, 7040, 7070.5, 7101,
7131.5, 7162.5, 7193, 7223.5, 7254, 7284.5, 7315.5, 7345, 7374.5, 7405,
7435.5, 7466, 7496.5, 7527.5, 7558, 7588.5, 7619, 7649.5, 7680.5, 7710,
7739.5, 7770, 7800.5, 7831, 7861.5, 7892.5, 7923, 7953.5, 7984, 8014.5,
8045.5, 8075, 8104.5, 8135, 8165.5, 8196, 8226.5, 8257.5, 8288, 8318.5,
8349, 8379.5, 8410.5, 8440, 8469.5, 8500, 8530.5, 8561, 8591.5, 8622.5,
8653, 8683.5, 8714, 8744.5, 8775.5, 8805, 8834.5, 8865, 8895.5, 8926,
8956.5, 8987.5, 9018, 9048.5, 9079, 9109.5, 9140.5, 9170, 9199.5, 9230,
9260.5, 9291, 9321.5, 9352.5, 9383, 9413.5, 9444, 9474.5, 9505.5, 9535,
9564.5, 9595, 9625.5, 9656, 9686.5, 9717.5, 9748, 9778.5, 9809, 9839.5,
9870.5, 9900, 9929.5, 9960, 9990.5, 10021, 10051.5, 10082.5, 10113,
10143.5, 10174, 10204.5, 10235.5, 10265, 10294.5, 10325, 10355.5, 10386,
10416.5, 10447.5, 10478, 10508.5, 10539, 10569.5, 10600.5, 10630,
10659.5, 10690, 10720.5, 10751, 10781.5, 10812.5, 10843, 10873.5, 10904,
10934.5, 10965.5, 10995, 11024.5, 11055, 11085.5, 11116, 11146.5,
11177.5, 11208, 11238.5, 11269, 11299.5, 11330.5, 11360, 11389.5, 11420,
11450.5, 11481, 11511.5, 11542.5, 11573, 11603.5, 11634, 11664.5,
11695.5, 11725, 11754.5, 11785, 11815.5, 11846, 11876.5, 11907.5, 11938,
11968.5, 11999, 12029.5, 12060.5, 12090, 12119.5, 12150, 12180.5, 12211,
12241.5, 12272.5, 12303, 12333.5, 12364, 12394.5, 12425.5, 12455,
12484.5, 12515, 12545.5, 12576, 12606.5, 12637.5, 12668, 12698.5, 12729,
12759.5, 12790.5, 12820, 12849.5, 12880, 12910.5, 12941, 12971.5,
13002.5, 13033, 13063.5, 13094, 13124.5, 13155.5, 13185, 13214.5, 13245,
13275.5, 13306, 13336.5, 13367.5, 13398, 13428.5, 13459, 13489.5,
13520.5, 13550, 13579.5, 13610, 13640.5, 13671, 13701.5, 13732.5, 13763,
13793.5, 13824, 13854.5, 13885.5, 13915, 13944.5, 13975, 14005.5, 14036,
14066.5, 14097.5, 14128, 14158.5, 14189, 14219.5, 14250.5, 14280,
14309.5, 14340, 14370.5, 14401, 14431.5, 14462.5, 14493, 14523.5, 14554,
14584.5, 14615.5, 14645, 14674.5, 14705, 14735.5, 14766, 14796.5,
14827.5, 14858, 14888.5, 14919, 14949.5, 14980.5, 15010, 15039.5, 15070,
15100.5, 15131, 15161.5, 15192.5, 15223, 15253.5, 15284, 15314.5,
15345.5, 15375, 15404.5, 15435, 15465.5, 15496, 15526.5, 15557.5, 15588,
15618.5, 15649, 15679.5, 15710.5, 15740, 15769.5, 15800, 15830.5, 15861,
15891.5, 15922.5, 15953, 15983.5, 16014, 16044.5, 16075.5, 16105,
16134.5, 16165, 16195.5, 16226, 16256.5, 16287.5, 16318, 16348.5, 16379,
16409.5, 16440.5, 16470, 16499.5, 16530, 16560.5, 16591, 16621.5,
16652.5, 16683, 16713.5, 16744, 16774.5, 16805.5, 16835, 16864.5, 16895,
16925.5, 16956, 16986.5, 17017.5, 17048, 17078.5, 17109, 17139.5,
17170.5, 17200, 17229.5, 17260, 17290.5, 17321, 17351.5, 17382.5, 17413,
17443.5, 17474, 17504.5, 17535.5, 17565, 17594.5, 17625, 17655.5, 17686,
17716.5, 17747.5, 17778, 17808.5, 17839, 17869.5, 17900.5, 17930,
17959.5, 17990, 18020.5, 18051, 18081.5, 18112.5, 18143, 18173.5, 18204,
18234.5 ;
lev = 4.99993801116943, 15.0002880096436, 25.0017585754395, 
35.0054130554199, 45.0133171081543, 55.0294990539551, 65.0618057250977,
75.1255111694336, 85.2503662109375, 95.4942855834961, 105.969917297363,
116.896240234375, 128.697906494141, 142.195251464844, 158.96061706543,
181.962783813477, 216.647918701172, 272.476715087891, 364.303039550781,
511.534759521484, 732.200866699219, 1033.21728515625, 1405.69750976562,
1830.88500976562, 2289.76782226562, 2768.24243164062, 3257.47875976562,
3752.44213867188, 4250.4013671875, 4749.9130859375, 5250.2265625 ;
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 ;
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 ;
volume: volcello" ;
thetao:missing_value = 1.e+20f ;
thetao:_FillValue = 1.e+20f ;

Replies (6)

RE: how to switch the irrgular model output to the regular data ? - Added by Jaison-Thomas Ambadan almost 12 years ago

Hi,

try converting to NetCDF cell format before remapping, i,e,

cdo remapbil,r360x180 -setgridtype,cell -selname,thetao in.nc out.nc

OR

cdo remapbil,r360x180 -selname,thetao -setgridtype,cell in.nc out.nc

Also, it is always good to use the latest version of CDO.

Cheers,
J.

RE: how to switch the irrgular model output to the regular data ? - Added by yong sun almost 12 years ago

Hi,

thanks for your kind reply. I have a try as your method but it still doesn't work.
I have get the first time step from the original file and call it "step1.nc" (attachment)
can you help  me a test ?

Thanks so much.

Yong

step1.nc (4.25 MB) step1.nc

RE: how to switch the irrgular model output to the regular data ? - Added by Jaison-Thomas Ambadan almost 12 years ago

Hi,

Please ignore my previous post. You just need to update the CDO! that's it. The remapping works perfectly fine.

cdo remapbil,r360x180 step1.nc test.nc
Warning (cdfScanVarAttributes) : NetCDF: Variable not found - areacello
cdo remapbil: Processed 840658 values from 1 variable over 1 timestep. ( 5.83s )

Although CDO complains that "areacello" the grid cell area info is missing, it okay I guess. You can also compute the grid area/weights and repeat the same. so just update your cdo to latest version!

Cheers,
J.

RE: how to switch the irrgular model output to the regular data ? - Added by yong sun almost 12 years ago

Thanks so much for helping to finally solve this problem.

Thanks again.

Best regards

Yong

RE: how to switch the irrgular model output to the regular data ? - Added by Waheed Iqbal almost 12 years ago

Dear

I used this command in order to make same grid of my two different data sets (but both with same resolution of 0.5 degree).
After this command my final file has 1 deg resolution(twice of the original) but i don't want to make it 1 deg i want that 0.5.

Could you please suggest me what to do?

RE: how to switch the irrgular model output to the regular data ? - Added by Jaison-Thomas Ambadan almost 12 years ago

Hi,

just double the no. of grids (lon-lat), i.e.,

cdo remapbil,r720x360 ifile ofile

Cheers,
J.

    (1-6/6)