getting the positions of global n,m,l indices

[ipelupessy]

what is the recommended way to get grid center positions from n,m,l from the c++ ocean class?

[erik808]

Not sure if I understand, do you mean returning an array index based n,m,l,i,j,k? That would be the FIND_ROW2 macro.

[ipelupessy]

i mean, what is the position in lon/lat and depth for a given grid index..

[erik808]

I'm afraid a conversion like that is not in Ocean right now so you would have to infer it using member functions in TRIOS_Domain:

        double Hdim() const {return hdim;}
        //@}

        //@{ get physical bound of the global domain
        double Xmin() const {return xmin;}
        double Xmax() const {return xmax;}
        double Ymin() const {return ymin;}
        double Ymax() const {return ymax;}
        double Zmin() const {return zmin;}
        double Zmax() const {return zmax;}
        //@}

For a coordinate i it would be something like (xmin + i*dx)*180/pi .

[ipelupessy]

that's ok..there is also no way to retrieve the depth levels? (these are not evenly spaced)

[erik808]

Sorry, no straightforward access to depth levels. But they follow from hdim, k, dz and qz:

do k = l,1,-1
    write(6,*) "    ", k, " ", dz*dfzT(k)*hdim
 enddo

dfzT gives the stretching for T locations: dfzT(k) = dfdz(ze(k),qz)
and

REAL FUNCTION dfdz(z,qz)
  implicit none
  ! INPUT/OUTPUT
  real z,qz
  ! LOCAL
  real tth,ch
  ch = cosh (qz * (z+1))
  tth = tanh (qz)
  IF (qz.gt.1.0) THEN
     dfdz = qz  / (tth * ch * ch)
  ELSE
     !      dfdz = 1.0
     dfdz = 1.0 + (1.-qz)*(1.-2.*z)
  ENDIF
  RETURN
END FUNCTION dfdz

Here only the (qz.gt.1.0) case is relevant.

[ipelupessy]

ok, thanks( I already had that) - I think it would be nice to have direct (read) access; just to be sure of the values and also to have the bitwise exact values used in the code..its not super high priority at the moment though

[Sbte]

I think it would still be nice to have an interface for this implemented in i-emic though, instead of coding it inside of osume. I'm not sure if you have your own i-emic fork that you're working on.

[ipelupessy]

we don't use a seperate fork...& I agree!!