[MITgcm-devel] something to report

[Manfredi Manizza]

Jean-Michel,

yes, your suggestion was the only modification I applied.
Pretty sure about it.

Cheers
Manfredi




On Mon, 21 Oct 2013, Jean-Michel Campin wrote:

> Hi Manfredi,
>
> On Mon, Oct 21, 2013 at 03:41:35PM -0400, Manfredi Manizza wrote:
>>
>> Jean-Michel,
>>
>> Just to clarify, is the nc file that is a bit larger
>> than the file in binray format whihc  I think is the default..
>
>>> and with some blank tiles ?
> Is it the case ?
>
>> I tried to impelment your sugggestion,
>> It blows the model up becuase it looks for a 6th tile
>> that is not there given that the domain is set for 5 tiles only, see
>> below from log file.
>
> Are you sure that you only changed W2_useE2ioLayOut
> from TRUE (= the default) to FALSE ?
> Because I've just checked and W2_useE2ioLayOut has no effect on loading
> the grid-files so I am little bit surprised (if reading grid-files
> was working before with W2_useE2ioLayOut=T it should work the same way now
> with W2_useE2ioLayOut=F).
>
> And just one more thing: not sure how to modify/adapt Gael's script
> for mnc + rdmnc.m output if you are using blank-tiles.
>
> Cheers,
> Jean-Michel
>
>> I will try to see if I use the grid nc files generated by the run
>> what happens to the matlab program that converts the output into
>> lat/lon.
>> It may work, not sure yet...
>>
>> forrtl: severe (29): file not found, unit 9, file
>> /nobackupp8/mmanizza/llc90/run_nocost/tile006.mitgrid
>> Image              PC                Routine            Line        Source
>> libintlc.so.5      00002AAAAC49CF6A  Unknown               Unknown
>> Unknown
>> libintlc.so.5      00002AAAAC49BAE5  Unknown               Unknown
>> Unknown
>> libifcore.so.5     00002AAAAB62DE12  Unknown               Unknown
>> Unknown
>> libifcore.so.5     00002AAAAB5A35CB  Unknown               Unknown
>> Unknown
>> libifcore.so.5     00002AAAAB5A2CDE  Unknown               Unknown
>> Unknown
>> libifcore.so.5     00002AAAAB5C0848  Unknown               Unknown
>> Unknown
>> mitgcmuv           00000000005A2B00  mds_facef_read_rs        1081
>> mdsio_facef_read.f
>> mitgcmuv           0000000000739A5C  ini_curvilinear_g        3482
>> ini_curvilinear_grid.f
>> mitgcmuv           0000000000741804  ini_grid_                3330
>> ini_grid.f
>> mitgcmuv           000000000075D1F1  initialise_fixed_        1904
>> initialise_fixed.f
>> mitgcmuv           00000000007721F2  the_model_main_          3119
>> the_model_main.f
>> mitgcmuv           00000000006DCF82  MAIN__                   4010  main.f
>> mitgcmuv           0000000000409ABC  Unknown               Unknown
>> Unknown
>> libc.so.6          00002AAAACA1CBC6  Unknown               Unknown
>> Unknown
>> mitgcmuv           00000000004099B9  Unknown               Unknown
>> Unknown
>>
>>
>>
>> Manfredi
>>
>>
>>
>>
>> On Sun, 20 Oct 2013, Jean-Michel Campin wrote:
>>
>>> Hi Manfredi,
>>>
>>> I don't know what is your "data.exch2" parameter file,
>>> but I assume that you have:
>>> W2_useE2ioLayOut=T (this is the default)
>>> and:
>>> W2_mapIO   = 1 (this is not the default, which is -1)
>>> and with some blank tiles ?
>>>
>>> Now, the option "W2_useE2ioLayOut" has not been implemented
>>> in pkg/mnc, and this could explain that you have different
>>> size between mnc + rdmnc.m and mdsio + rdmds.m
>>>
>>> I think that if you set:
>>> W2_useE2ioLayOut=F
>>> the plain binary output file (pkg/mdsio + rdmds.m) should have the
>>> same size the current output that you get with mnc + rdmnc.m
>>>
>>> Cheers,
>>> Jean-Michel
>>>
>>> On Sat, Oct 19, 2013 at 12:21:41PM -0400, Manfredi Manizza wrote:
>>>>
>>>> Hi,
>>>>
>>>> I am reporting a feture I found comparing the output
>>>> of the ECCO lat-lon-cap global set-up provided by Patrick to me.
>>>>
>>>> The dimensions of the files are different depending on the format.
>>>> If in netcdf :
>>>>
>>>> S =
>>>>
>>>>         attributes: [1x1 struct]
>>>>    iters_from_file: [12x1 double]
>>>>        si_AREAtave: [90x1082x12 double]
>>>>
>>>> If in binary :
>>>>
>>>>
>>>>> more pickup_seaice.0000008772.meta
>>>> nDims = [   2 ];
>>>> dimList = [
>>>>    90,    1,   90,
>>>>  1170,    1, 1170
>>>> ];
>>>> dataprec = [ 'float64' ];
>>>> nrecords = [     6 ];
>>>> timeStepNumber = [       8772 ];
>>>> timeInterval = [  3.157920000000E+07 ];
>>>> nFlds = [    6 ];
>>>> fldList = {
>>>> 'siTICE  ' 'siAREA  ' 'siHEFF  ' 'siHSNOW ' 'siUICE  ' 'siVICE  '
>>>> };
>>>>
>>>>
>>>> The grid quantities that I have , when loaded, have the dimensions
>>>> of 90 by 1170.
>>>>
>>>> In fact i can use the functions of gcmfaces only using the binary
>>>> files...
>>>>
>>>> I am not sure if there is something I am missing or if
>>>> there is something wrong in the code. I just felt I needed to
>>>> send this out.
>>>>
>>>> Thanks in advance for your attention.
>>>>
>>>> Cheers
>>>> Manfredi
>>>>
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> 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
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>