Project

General

Profile

remapbil without gridfile?

Added by Juanma Sancho almost 12 years ago

Hello,
I am interested in using remapbil but defining the new grid specifications directly in the cdo command line without using an external 'gridfile'. This would add a point of efficiency in loops interpolating multiple files, por instance. In case of regular latlon grids interpolations I know the possibility to define the grid in the command (,r50x50) but what happens when the input file is a rotated one (like the attached) Is it mandatory to provide a gridfile defining the grid and the rotated pole coordinates? It would be nice to define everything (xfirst, xinc, etc..) in just one command without accessing for every input file to the disk.

This is the example:

cdo= Cdo()
fileIn = 'zrot.grb'
gridfile= 'gridfile'
interpArray=cdo.remapbil(gridfile,input="-copy"+" "+fileIn,options = "-f nc",returnCdf=True).variables['var61'][0,0,:,:]

Thanks in advance

gridfile (120 Bytes) gridfile file defining grid
zrot.grb (154 KB) zrot.grb data