[MITgcm-support] Depth at location and hFac?

[David Hebert]

Ahh, actual depth is stored in R_low.  Excellent!

Yea, my intent is to use the model depth for computations on the fly, so 
it looks like using -R_low will make life much easier.

Thanks to you both for your answers.

David


On 08/07/09 13:21, Jean-Michel Campin wrote:
> Hi David,
>
> Yes. And if you need the actual depth in the fortran code,
> for an oceanic set-up in z coordinate, you can get it simply from
> H = - R_low
>
> Jean-Michel
>
> On Fri, Aug 07, 2009 at 02:08:46PM -0400, David Ferreira wrote:
>    
>> Hi David,
>> I think that the actual depth used by the model is an output, in
>> Depth.*.data(meta) or in grid.*.nc.
>> david
>>
>>
>> David Hebert wrote:
>>      
>>> Given the use of partial cells and discretization it seem possible
>>> (perhaps likely) that the actual depth of a water column will differ
>>> from the values specified in bathyFile. If I want to obtain the actual
>>> depth used in the model do I need to loop over the rF array and check
>>> if it is smaller than R_low? (smaller since rF and R_low are negative
>>> values).  Then apply a hFac to the cell. So something like this...
>>>
>>> do k=2,nr
>>>      if ( (rF(k-1) .GT.  R_low(i,j,bi,bj) .AND.
>>>             rF(k) .LE. R_low(i,j,bi,bj) ) then
>>>
>>>           H = rF(k) + (dF(k) * hFacC(i,j,k,bi,bj))
>>>
>>>     end if
>>>
>>> end do
>>>
>>> where H is the depth at that particular location. Not sure if it
>>> should be dF(k) or dF(k-1). Same on rF.
>>>
>>> Any suggestions/comments are appreciated.
>>>
>>> thanks,
>>>
>>> David
>>>
>>> _______________________________________________
>>> 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
>>      
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>    

-- 
David A. Hebert
NRC Postdoc
Naval Research Lab
Stennis Space Center, MS 39529
david.hebert.ctr at nrlssc.navy.mil
Phone: (228) 688-5846
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20090807/49d42661/attachment.htm>