<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<small><big>Hello,<br>
<br>
I have been struggling to calculate Barotropic streamfunction
and AMOC using online examples and function but I have recently
been told it is due to the fact that the data I use contains
missing values, therefore all values are set to the missing
value. <br>
<br>
How can I somehow ignore the missing values during calculation?<br>
<br>
Here is my code for the bartropic streamfunction:</big><br>
<br>
----------------------------------------------------------------------------------------------------------------------------------------------------<br>
<br>
<small><font size="4"><small><small><font color="FF0000">load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
<br>
begin<br>
<br>
;##########################################################################<br>
;##########################################################################<br>
;-------------------------Barotropic
Plots---------------------------------<br>
;--------------------------------------------------------------------------<br>
;##########################################################################<br>
;##########################################################################<br>
<br>
;----------Load in data<br>
<br>
;ICON<br>
<br>
noshupw = addfile
("con.ym.atlbox.noshift.100y.nc", "r")<br>
u = noshupw->u_vint_acc<br>
u&lon = u&lon-360 ; Fix the
longitudes to be < 360<br>
lat = noshupw->lat<br>
<br>
<br>
v = noshupw->v_acc<br>
<br>
printVarSummary(u)<br>
printVarSummary(v)<br>
<br>
dimUV= dimsizes(barai)<br>
<br>
klev = 20 <font color="FF0000">; number of levels in
total</font><br>
<br>
sf = v(:,{klev},:,:)<br>
vp = v(:,{klev},:,:)<br>
sf@long_name = "stream function"<br>
vp@long_name = "velocity potential"<br>
sf@units = "m^2/s"<br>
vp@units = "m^2/s"<br>
<br>
printVarSummary(sf)<br>
printVarSummary(vp) <br>
<br>
uv2sfvpg(u,v,sf,vp)</font></small></small></font></small><br>
</small><br>
<small><br>
----------------------------------------------------------------------------------------------------------------------------------------------------<br>
<br>
<big><big><small>and my AMOC code is largely based on the online
example moc_1.ncl<br>
<br>
<br>
<br>
thanks,<br>
Michael</small><br>
</big></big><br>
</small><br>
</body>
</html>