[MITgcm-devel] problem with gluemnc ?

[Jean-Michel Campin]

Hi Martin,

Thanks for your advise.
I will check-in this fix, + do something about the number at the 
beginning (since there might be some folks like me who got their 
netcdf pkg from fc11 and don't want to re-install netcdf 4.0.1).
And then I will add what Ryan gave me to apply gluemnc to 
grid.*.nc files (although it seems to drop the 1.D vector
like drF, drC ... a bit annoying).
Cheers,
Jean-Michel

On Thu, Aug 13, 2009 at 08:52:17AM +0200, Martin Losch wrote:
> Hi Jean-Michel,
>
> I can confirm that with netcdf 3.6.2 it works for my but for netcdf 4.0 
> it doesn't (on the same linux computer).
>
> Actually on the netcdf site, the bug with numbers at the beginning of  
> the file name is already documented along with a fix (basically upgrade 
> to 4.0.1)
> <http://www.unidata.ucar.edu/software/netcdf/docs/known_problems.html>
> I haven't seen a bug report for the other problem (multiple files on the 
> command line), but in glumnc this workaround could work:
> Instead of
>>       withY=$(ncdump -h $countxx.$somevar.$somepre.*.nc | grep "Y =")
>>       withYp1=$(ncdump -h $countxx.$somevar.$somepre.*.nc | grep "Yp1 
>> =")
> use this
>>        filelist=(`ls $countxx.$somevar.$somepre.*.nc`)
>>        withY=$(ncdump -h ${filelist[0]} | grep "Y =")
>>        withYp1=$(ncdump -h ${filelist[0]} | grep "Yp1 =")
>
> Martin
>
> On Aug 12, 2009, at 11:38 PM, Jean-Michel Campin wrote:
>
>>
>> Looks like it's related to the version of ncdump/netcdf:
>> It does not work with my version:
>>> netcdf library version "4.0" of Feb 26 2009 03:56:01 $
>> but on beagle it works with:
>>> netcdf library version "3.6.2" of Oct 26 2007 12:30:01 $
>>
>> Jean-Michel
>>
>> On Wed, Aug 12, 2009 at 12:24:52PM -0400, Jean-Michel Campin wrote:
>>> Hi,
>>>
>>> I tried to use gluemnc for the output of exp. ideal_2D_oce
>>> and also global_ocean.90x40x15 , statring with "state.*.t*.nc" files.
>>>
>>> I did not get the right single netcdf file at the end,
>>> but when I comment out the 2 lines 277,278 in gluemnc:
>>>>     withY=$(ncdump -h $countxx.$somevar.$somepre.*.nc | grep "Y =")
>>>>     withYp1=$(ncdump -h $countxx.$somevar.$somepre.*.nc | grep  
>>>> "Yp1 =")
>>> then it seems to work much better.
>>>
>>> Is there something wrong there ? Or did I miss something ?
>>>
>>> Jean-Michel
>>>
>>>
>>> _______________________________________________
>>> MITgcm-devel mailing list
>>> MITgcm-devel at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel