[MITgcm-devel] density

[Alistair Adcroft]

Sounds like a nice idea but WHICH density would you make a 3D array? 
In-situ density is useless for PV and neutral physics. Sigma-theta is 
most widely used but d/dz sigma-theta is meaningless (except for PV?). I 
guess you could do it all with two 3D arrays, one with sigma-theta 
referenced to k+1/2 and the other referenced to k-1/2.

-A.

Baylor Fox-Kemper wrote:

> I'd be very much in favor.  It would make the PV calculations a lot  
> easier too.  However, with nonlinear EOS, one would still need to  
> calculate density again to figure out N^2.  But, it would be nicer to  
> have that as a special case that gets called only with nonlinear EOS...
>   -Baylor
>
>
>
> On Mar 10, 2006, at 9:21 AM, Patrick Heimbach wrote:
>
>>
>> Hi there,
>>
>> the density calculation and diagnostics
>> in do_oceanic_physics is a mess (I think).
>> Can't we just make density a 3-d array
>> once and for all?
>>
>> It would also save me a lot of fooling
>> with TAF to keep diagnostics that
>> would otherwise be thrown out
>> (e.g. simple call to diags_rho requires
>> some stupid recomp.).
>>
>> -Patrick
>>
>> -- 
>> --------------------------------------------------------
>> Patrick Heimbach   Massachusetts Institute of Technology
>> FON: +1/617/253-5259                  EAPS, Room 54-1518
>> FAX: +1/617/253-4464             77 Massachusetts Avenue
>> mailto:heimbach at mit.edu               Cambridge MA 02139
>> http://www.mit.edu/~heimbach/                        USA
>>
>> _______________________________________________
>> 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