Re: Re : Re: [GMT-HELP] Re : Re: [GMT-HELP] Netcdf file help

So, it is possible to extract a 5D layer (temp, lev, lat, lon, value) in GMT?


On May 22, 2014, at 9:38 AM, fabien.retif <[email protected]> wrote:

> Yes it's ok. grd2xyz extracts a 2D layer (lon, lat ,value). If the GMT command starts by grd*, you have to give a 2D grid layer (grd) .
> So you can select a layer 2D from a netcdf file with N-dim by netcdf.nc?myVar[0,0,...,0]
> If you want to extract a specific point (1D), you can use grdtrack. You have to put the coordinate of your point in a ascii file (one coordinate (lon lat) by line) and to give your grd file with -Gmynetcdf.nc?myVar[0,0]
> I hope it's clear.
> .
> Fabien
> 
> ---- Le jeu., 22 mai 2014 16:55:56 +0200 Sumant Jha<[email protected]> a écrit ---- 
> 
> That helped.
> I used the following command:
> grd2xyz geoschem_nestedNA_10_20041011.nc?IJ_AVG_S__BCPI[time,lev,lat,lon] > IJ_AVG_S__BCPI_10_20041011_1.xyz
> 
> and got the following output: 
> 220  70 0.00774827273563
> 220.666666667 70 0.00774506758898
> 221.333333333 70 0.00774536235258
> 222  70 0.00741489697248
> 
> This output file has just three columns, whereas in my grd2xyz command to gmt, i asked for time, lev, lat, lon, as well. Did I do something wrong?
> 
> 
> On May 22, 2014, at 8:41 AM, fabien.retif <[email protected]> wrote:
> 
>> grd2xyz geoschem_nestedNA_10_20041011.nc?IJ_AVG_S__BCPI[time,lev] > IJ_AVG_S__BCPI_10_20041011.xyz
>> 
>> netcdf.nc?varName[time_index,lev_index,.....,Ndim_index]
> 
> Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first. To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu To unsubscribe, send the message "signoff gmt-help" to [email protected] Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.
> 
> 
> 
> _.-~-._.-~-._.-~-._.-~-._.-~-._ 
> PhD. Fabien Rétif
> Géosciences Montpellier
> UMR 5243 - CC 60
> Universite Montpellier II
> Bat 22 - Bureau 468
> Place Eugène Bataillon
> 34095 Montpellier Cedex 5
> France
> Tel : +33 (0)4.67.14.35.85
> _.-~-._.-~-._.-~-._.-~-._.-~-._ 
> Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first. To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu To unsubscribe, send the message "signoff gmt-help" to [email protected] Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.


Mailing list for GMT discussions of all kinds.  If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
To unsubscribe, send the message "signoff gmt-help" to [email protected]
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.