[Mitgcm-support] Re: W field

[mitgcm-support at dev.mitgcm.org]

Hi Tony/Dimitri,

 Happy New Year. 

 I'll check back at the old code. Its a while since I looked back at c34
et. al. Which flavor of run are the fields from e.g. pure forward or
adjoint assimilated or reduced kalman filter assimilated. If its the
reduced kalman filter output, is it guaranteed that the w field be
consistent with the  u,v or is there an extra wind pumping term? 

 I think Mick, Galen and Steph came across something similar with W not
equal to div_h(U) in the JPL stored fields (so I've cc'd them) but I
thought that was limited to the reduced kalman filter assimilation - M,
G or S any comments?

Chris
 
On Fri, 2003-01-03 at 20:39, Tony Lee wrote:
> Chris,
> 
> I found that our 10-day averaged vertical velocity (the Wave* file)
> generated from our code (c19) is inconsistent
> with Uave and Vave in the sense that the divergence for a given cell box
> is not zero.  This creates a problem
> in doing heat budget because the net volume flux due to the divergence
> is multiplied with temperature (say
> 30 degree of SST) to enhance that divergence problem.   Dimitris noticed
> that, in do_time_average subroutine
> (for any code before c34 I believe), the 10-day average of wvel are
> different from that of uvel and vvel.  It appears
> that they the wvel is offset from the uvel and vvel by half a time
> step.  That might be the cause of the divergence
> I found.    I just want to confirm with you that for any code before
> c34, Wave is offset from Uave and Vave.
> If that's confirm, I'll pull the averaged Wave* files from the data
> server and ask Detlef to do the same for the SIO
> server.   People who want W will have to compute it by themself using
> the continuity equation.
> 
> Thanks in advance for the clarification.
> 
> Tony
> 
> --
> Dr. Tong (Tony) Lee
> Research Scientist
> 
> MS 300-323
> Jet Propulsion Laboratory
> 4800 Oak Grove Drive, Pasadena, CA 91109
> 
> Tel: (818) 354-1401  Fax: (818) 393-6720  http://oceans-www.jpl.nasa.gov/
> 
>