Re: [ferret_users] LEVITUS Ocean Data/FERRET shows differing results

[Olivier Marti <olivier.marti@xxxxxxxxxxxx>]

Le 24 mars 2013 à 13:59, Star Physicist a écrit :

Dear FERRET Users,

I have LEVITUS Ocean data at 16 depth layers over entire globe. I was trying to calculate time series averaged over entire depth and entire global scale.
I tried the following two ways but didn't understand why the results of both two ways are not exactly overlapping?
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1st Way
use tanom_1955-2004_missing_value.nc
let RedSea=t_an[k=1:16@ave]
let anom1=RedSea[x=30E:45E@ave,y=10N:30N@ave]
plot/l=1:200 anom1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!2nd Way
let RedSea2=t_an[x=30E:45E@ave,y=10N:30N@ave]
let anom2=RedSea2[k=1:16@ave]
plot/l=1:200/over anom2
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

I may guess that it is a mask problem. At some grid points, the sea bottom is shallower than the level K=16.

OIivier

-- 

Dr. Olivier Marti  - mailto:olivier.marti@xxxxxxxxxxxx

Laboratoire des Sciences du Climat et de l'Environnement

Institut Pierre Simon Laplace

Laboratoire CEA-CNRS-UVSQ  -  UMR 8212

Tel : +33 1 69 08 77 27 - Mob : +33 6 45 36 43 74

Fax : +33 1 69 08 30 73

<http://www.lsce.ipsl.fr/Pisp/56/olivier.marti.html>