[MITgcm-support] diagnostics : ADVr_TH,DFrE_TH and WVELTH

[Krishnakumar Rajagopalan]

Hi Martin,
 
Thanks a lot for your reply. From diagnostics_main_init.F, the L in column 9 means the diagnostics are at the top of the layer. So  as you suggested, it is clear why the fluxes should be zero for first layer, but WVELTH need not be. 
 
Best regards
 
Krishnakumar

--- On Thu, 6/23/11, Martin Losch <Martin.Losch at awi.de> wrote:


From: Martin Losch <Martin.Losch at awi.de>
Subject: Re: [MITgcm-support] diagnostics : ADVr_TH,DFrE_TH and WVELTH
To: mitgcm-support at mitgcm.org
Date: Thursday, June 23, 2011, 8:03 AM


Hi,
ADVr_TH is the advective flux through the surface of the ocean. I guess it's clear that you cannot advect a property through the sea surface (the same is true for the diffusive flux DFrE_TH). Depending on your parameters the vertical velocity is W = dEta/dt, which is non-zero in the non-stationary case. As a consequence W*THETA is non-zero, which also makes sense as you transport temperature vertically with W, but you do not flux temperature throught the surface.

see DIAGNOSTICS_MAIN_INIT (or documentation?) for explanation of the different identifiers (LR, etc)

M.

On Jun 22, 2011, at 11:10 PM, Krishnakumar Rajagopalan wrote:

> Hi All,
>  
> I am reading the diagnostics ADVr_TH, DFrE_TH and WVELTH  from *.data, *.meta files in tutorial tutorial_global_oce_latlon with DELR setting and GMREDI scheme.  All the 3 diagnostics are [WM        LR].  I think M in [WM      LR] means these diagnostics are at center of the cell (thus center of the layer).  What does L and R mean?
>  
> I would like to calculate the advective and diffusive flux of pot. temp from the first layer. I am finding that ADVr_Th and DFrE_Th are always zero in the first layer. But WVELTH in the first layer takes non zero values. Why is  ADVr_TH  zero since it is a product of WVELTH and area. 
>  
> Thanks in advance
>  
> Krishnakumar
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support


_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20110623/195fda5f/attachment.htm>