Re: [idvusers] Plotting views across the prime meridian using IDD GFS

NCO (netCDF Operator) is an utility to be used for processing netcdf files. It could be change dimensions, extracts sub regions etc. The web site is http://nco.sourceforge.net/

i hope this helps you,

--ufuk


webmaster-IWC wrote:
Hi
Thanks for replying.

PINGO does not support GRIB2 products.
What is the NCO ?

Ian


-----Original Message-----
From: idvusers-bounces@xxxxxxxxxxxxxxxx
[mailto:idvusers-bounces@xxxxxxxxxxxxxxxx] On Behalf Of [BE] Ufuk Utku
Turuncoglu
Sent: 03 January 2008 13:01
To: idvusers@xxxxxxxxxxxxxxxx
Subject: Re: [idvusers] Plotting views across the prime meridian using
IDD GFS


Hi,

For GRIB files you can use PINGO (http://www.mad.zmaw.de/Pingo/pingohome.html), for example

for time subset:
grb selrec,1,100809 ../A2/EH5_OM_A2_1_TEMP2_1-146096.grb EH5_OM_A2_1_TEMP2_1-100809.grb

and for lat lon subset;
grb sellonlatbox,290,360,80,10 ./EH5_OM_A2_1_TEMP2_1-100809.grb ./EH5_OM_A2_1_TEMP2_1-100809_290-360-80-10.grb grb sellonlatbox,0,70,80,10 ./EH5_OM_A2_1_TEMP2_1-100809.grb ./EH5_OM_A2_1_TEMP2_1-100809_0-70-80-10.grb
rm EH5_OM_A2_1_TEMP2_1-100809.grb

but i don't know about GRIB2 support and CDO (Climate Data Operator).

For netcdf files NCO can modify the netcdf file as you want. The following PyNGL scripts helps to you,

import Nio
from Numeric import *
import os

# slab dimensions
os.system('rm nnrp.slab.nc')
os.system('ncks -d lon,330.0,80.0 -d lat,17.5,60.0 nnrp.nc
nnrp.slab.nc')

# open netcdf file
nc = Nio.open_file("./nnrp.slab.nc", "rw")
lon = nc.variables['lon']

# fix dimension >> monotonic
for i in xrange (0,12):
        lon[i] = lon[i][0]-360.0

# close file
nc.close()

The os.system part of the code gets subset of the file and rest of the code correct 360 problem and generates monotonic dimension.

best

Ufuk Utku Turuncoglu
Istanbul Technical University
Informatics Institute

webmaster-IWC wrote:
Hi Dan,
I have not worked out how to directly reply to your question on the forum, so I have done it this way. I have had the same problem and I know the cause and a possible solution, but the solution requires creating a regional subset of the grib or NetCDF file. The problem lies with the grib file starting from 0 - 359 and IDV not able to join the two halves of the map together. If you can create the

grib file to start from -60 to 60 then this draws this Ok.
This was not a problem as software is available for grib 1 products to

create regional subsets is available, but the NCEP will be switching to grib2 and that's where the my problem is. Where I used to used GGRIB to create sub regions, this will not work for grib2 products. So if any programmers read this, I am in desperate need of a program which will create sub regions of grib2 products complied for windows XP. Ian Gooch

No virus found in this outgoing message.
Checked by AVG Free Edition.
Version: 7.5.516 / Virus Database: 269.17.13/1207 - Release Date:
02/01/2008 11:29
_______________________________________________
idvusers mailing list
idvusers@xxxxxxxxxxxxxxxx
For list information, to unsubscribe, visit: http://www.unidata.ucar.edu/mailing_lists/

_______________________________________________
idvusers mailing list
idvusers@xxxxxxxxxxxxxxxx
For list information, to unsubscribe, visit:
http://www.unidata.ucar.edu/mailing_lists/
No virus found in this incoming message.
Checked by AVG Free Edition. Version: 7.5.516 / Virus Database: 269.17.13/1207 - Release Date:
02/01/2008 11:29
No virus found in this outgoing message.
Checked by AVG Free Edition. Version: 7.5.516 / Virus Database: 269.17.13/1207 - Release Date:
02/01/2008 11:29



  • 2008 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the idvusers archives: