[MITgcm-support] xmitgcm announcement

[Raphaël Dussin]

Hi Martin,

xmitgcm cannot open a dataset when all the variables are not aligned along the time axis. In other words, all variables in a dataset need to be outputted at the same frequency. If you’re saving variables with different frequencies, you need to specify which variables to read. You can create multiple dataset for each output frequency.

I am off in France this week but I will take a look at it when I get a chance.

Thank you for your feedback! 

Raphael 

> On 5 Sep 2018, at 13:18, Martin Losch <Martin.Losch at awi.de> wrote:
> 
> Hi Raphael,
> 
> sorry, I meant to send this to MITgcm-support.
> 
> Also, I was probably to quick to complain, I apologize.
> I could identify at least one problem: the diagnostics package produces a file cheapaml at iteration 11 and iceDiag at 12. I acknowledge that this is hard to catch in a generic reading software. But when I remove the cheapaml (and iceDiag) files and retry I still get:
> 
> OSError: Could not find the expected file prefixes ['W', 'Eta', 'ice_Tice1', 'ice_iceH', 'ice_Qice2', 'ice_snowH', 'S', 'ice_Tsrf', 'ice_Qice1', 'T', 'ice_snowAge', 'ice_fract', 'V', 'VICE', 'U', 'ice_Tice2', 'UICE'] at iternum 12. (Instead found ['ice_Tsrf', 'ice_Qice1', 'T', 'V', 'ice_Tice2', 'Eta', 'Qnet', 'ice_frwAtm', 'Qsw', 'EmPmR', 'ice_flxAtm', 'W', 'PH', 'ice_Tice1', 'ice_Qice2', 'FU', 'ice_snowH', 'VICE', 'U', 'PHL', 'ice_iceH', 'S', 'FV', 'ice_snowAge', 'ice_fract', 'UICE'])
> 
> To me this looks like the same list of files but in different order. ignore_unknown_vars=True does not work.
> 
> Martin
> 
> 
>> On 09/05/2018 12:17 PM, Martin Losch wrote:
>> Hi Raphael,
>> I just install the new xmitgcm.
>> A naive question that may be asked in the wrong place:
>> I would expect that with xmitgcm I can do this:
>> d=open_mdsdataset("someDir")
>> I tried this in a verification experiment (in fact it's one of the "verification_other/offline_cheapaml/tr_run.dyn", that you can access from the MITgcm_contrib area) and I get this error:
>> OSError: Could not find the expected file prefixes ['W', 'Eta', 'ice_Tice1', 'ice_iceH', 'ice_Qice2', 'ice_snowH', 'S', 'ice_Tsrf', 'ice_Qice1', 'T', 'ice_snowAge', 'ice_fract', 'V', 'VICE', 'U', 'ice_Tice2', 'UICE'] at iternum 11. (Instead found ['cheapAML'])
>> There is no iternum 11 (only 12). When I add this list (which is the "default" output of the model without the diagnostics fields), I get the same error and only with ignore_unknown_vars=True, there is no error, but only T,S,U,V,W,Eta are in the data set, everything with "ice_" or "ICE" is ignored.
>> I need to specify a prefix to read the diagnostics file (called iceDiag in this case). That works fine, but is that expected?
>> Also d.T.values does not return the values of variable T but only the enire dataset description, maybe because .T is interpreted as transpose?
>> Maybe this has all been answered elsewhere, please point me there.
>> Martin
>>> On 09/04/2018 11:04 AM, Raphael Dussin wrote:
>>> Hi all,
>>> 
>>> A quick message to the python users (or those thinking about
>>> making the switch) : the new release of xmitgcm (v0.2.2) is here.
>>> 
>>> xmitgcm allows you to read your mds files and build a xarray dataset
>>> from them that you can use in your analysis workflow.
>>> 
>>> Lots of improvements and bug fixes and even easier to install
>>> with anaconda. All the info here:
>>> 
>>> https://xmitgcm.readthedocs.io/en/latest/index.html
>>> 
>>> Best,
>>> Raphael
>>> 
>>> *-------------------------------------------------------------------------------------------* 
>>> 
>>> /Raphael Dussin/
>>> 
>>> /C'est pas l'homme qui prend la mer, c'est la mer qui prend l'homme.../
>>> 
>>> 
>>> 
>>> _______________________________________________
>>> MITgcm-support mailing list
>>> MITgcm-support at mitgcm.org
>>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>>> 
> 
> -- 
> Martin Losch
> Alfred Wegener Institute for Polar and Marine Research
> Postfach 120161, 27515 Bremerhaven, Germany;
> Tel./Fax: ++49(0471)4831-1872/1797
>