[MITgcm-support] xmitgcm announcement

[Martin Losch]

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