[MITgcm-devel] Seaice wind-stress diagnostics

[Martin Losch]

Hi Jean-Michel,

thanks for doing this. I guess there used to be a reason for these diagnostics, back when exf was not entirely required for running seaice, but that’s really before I ever had a look at the code. 
As far as I remember there were very subtle differences in the heat and fw-flux diagnostics to oceFWflx and oceQnet etc, but I cannot remember them, so they can probably also be retired.

Martin

> On 9. May 2017, at 05:25, Jean-Michel Campin <jmc at mit.edu> wrote:
> 
> Hi Martin,
> 
> I added diagnostics of wind-stress over seaice and also the
> seaice fraction weighted wind stress over ocean+seaice.
> If you (or any one) found someting curious (not unlikely since I did 
> not check too much) or has some suggestion, don't hesitate.
> 
> I also remove SIfu & SIfv diags (duplicates), and also SIuwind,SIvwind
> that I found confusing (since the corresponding variables belong to 
> pkg/exf); for now, kept SIempmr, SIqnet, SIqsw.
> 
> Note that there is a good reason for having diagnostics of ocean net surface 
> fluxes within pkg/thsice since this pkg is used in atmospheric config and
> in this case, the main model surface-fluxes diag are not active.
> In the case of pkg/seaice, the reason is not so clear ...
> 
> Cheers,
> Jean-Michel
> 
> On Fri, Mar 10, 2017 at 02:49:59PM +0100, Martin Losch wrote:
>> Hi Jean-Michel,
>> 
>> you are absolutely right SIfu/fv are duplicates and we don???t need them (if fu and fv are diagnosed ***after*** calling seaice_model, which seems to be the case.). Similar: SIempmr, SIqnet, SIqsw (-oceFWfx, -oceQnet)
>> 
>> I don???t think that we have a diagnostics for sea-ice-atmosphere stress. 
>> The atmosphere sea ice stress is calculated in seaice_get_dynforcing and stored in taux/tauy, which are fields local to pkg/seaice and can be easily included in seaice_model? Or do it in seaice_get_dynforcing?
>> 
>> Martin
>> 
>>> On 9 Mar 2017, at 17:02, Jean-Michel Campin <jmc at mit.edu> wrote:
>>> 
>>> Hi Martin and others,
>>> 
>>> I was looking for a diagnostics of wind-stress over seaice,
>>> but I don't think we currently have one (although there are many
>>> diagnostics available within pkg/seaice, and I sometime get lost there).
>>> 
>>> Initially, after reading the "description" (from: available_diagnostics.log), 
>>> I thought that these 2 (SIfu,SIfv) would be what I was looking for:
>>> 'SEAICE zonal surface wind stress, >0 increases uVel '
>>> 'SEAICE merid. surface wind stress, >0 increases vVel'
>>> but there are simply an exact copy  of model main diagnostics: oceTAUX & oceTAUY
>>> (exactly from the same FFIELDS.h arrays fu & fv), i.e., the ocean surface stress 
>>> and not the seaice wind-stress.
>>> 
>>> It's not very clear to me why we need these duplicate diagnostics,
>>> but if we want to keep SIfu,SIfv, would suggest to change the description 
>>> to avoid confusion.
>>> 
>>> Also I think it would be useful to have a seaice wind-stress diagnostics
>>> (if it does not already exist ?).
>>> 
>>> Cheers,
>>> 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