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

Which way I should proceed then? How could I overcome this issue? Actually there is a function in CDO for vertical averaging over all grid points called "cdo vertmean ifile ofile" but I don't know whether it can handle bottom shallower sea level. Any suggestion?

Thanks,

Mubashar

Subject: Re: [ferret_users] LEVITUS Ocean Data/FERRET shows differing results
From: olivier.marti@xxxxxxxxxxxx
Date: Sun, 24 Mar 2013 17:50:27 +0100
CC: couragetoknowgc@xxxxxxxxxxx
To: ferret_users@xxxxxxxx

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>