[MITgcm-devel] SFLUX

[Jean-Michel Campin]

Hi Dimitris,

I looked to the history of this SFLUX diagnostics:
It seems to me that you coded the first version:
MITgcm/pkg/diagnostics/Attic/diagnostics_fill_surf_flux.F
And later on, I moved this file to model/src/diags_oceanic_surf_flux.F
but did not change this "SFLUX" diagnostics.
Do you remember why you defined SFLUX this way ?
I am not 100% sure, but this older version (before your latest
change) should match the evolution of the salt content
(and we can agree that it should be like this, in the same way
as TFLUX should match the evolution of the heat content).

Now, it would be better (and also easier to understand)
to make the distinction
IF (temp_EvPrRn.EQ.UNSET_RL) THEN / ELSE / ENDIF
and
IF (salt_EvPrRn.EQ.UNSET_RL) THEN / ELSE / ENDIF

Not directly related to SFLUX:
The diagnostics that I would like to add (with filling calls in
diags_oceanic_surf_flux.F) are for Qnet, EmPmR, saltFlux
(in addition to Qsw, tRelax & sRelax, taux & tauy) to have
the complete set of oceanic forcing diagnostics
(and after I could remove the pkg/thsice diagnostics:
SIflx2oc, SIfrw2oc, SIsaltFx).

Jean-Michel

On Thu, Nov 09, 2006 at 10:10:06PM -0800, Dimitris Menemenlis wrote:
> JM, is pkg/diagnostics SFLUX literally supposed to be the "net surface salt 
> flux", that is, EmPmR is actually supposed to have zero net contribution to 
> pkg/diagnostics SFLUX as opposed to pkg/exf SFLUX, which is defined to be 
> equivalent to EmPmR?
> 
> If yes, then since Tue Nov 7 00:38:22 2006 UTC (3 days, 5 hours ago) I have 
> screwed up diags_oceanic_surf_flux.F and I need to change it back ASAP.
> 
> I changed diags_oceanic_surf_flux.F because I noticed that SFLUX was 
> exactly identical to SRELAX for:
> #ifdef NONLIN_FRSURF
>          IF ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
>      &        .AND. useRealFreshWaterFlux ) THEN
> and for: salt_EvPrRn=0.
> 
> Now I realize that you probably "did" mean SFLUX to be identical to SRELAX 
> for above conditions.  Sorry!
> 
> But this also means that pkg/diagnostics SFLUX is not a very useful 
> quantity for most oceanographers and that we should add a separate 
> freshwater flux output to pkg/diagnostics.  No?
> 
> Please advise.
> 
> Dimitris
> 
> -- 
> Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
> Jet Propulsion Lab, California Institute of Technology
> MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
> tel: 818-354-1656;  fax: 818-393-6720
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel