[MITgcm-support] Depth at location and hFac?

[Jean-Michel Campin]

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