Re: Potential Vorticity and spiciness

[Mick Spillane <spillane@xxxxxxxxxxxxx>]

Hi Paulo,
	The following is a script that computes "spice".  Its based on
a program that we use at PMEL but while I think the reference is
Flament, I can't seem to find it.  I also have put it in the form of
a ferret external function which I will put in the contributed list
soon.
Good luck,
Mick

can mode verify
! evaluates "spice" for specified potential temperature and salinity
! which are contained in the input variables "potemp" and "salt"
!
! note that potential temperature is available through the ferret
!      function THETA_FO(salt,temp,p,pref)
!
! example:
!       let potemp = {5,5,5,5}
!       let salt = {25,30,35,40}
!       go spice
!       list spice returns -8.044, -3.618, 0.418, 4.047
!
let b00=0          ; let b01=0.77442    ; let b02=-0.00585
let b03=-9.84E-4   ; let b04=-2.06E-4
let b10=5.1655E-2  ; let b11=2.034E-3   ; let b12=-2.742E-4
let b13=-8.5E-6    ; let b14=1.36E-5
let b20=6.64783E-3 ; let b21=-2.4681E-4 ; let b22=-1.428E-5
let b23=3.337E-5   ; let b24=7.894E-6
let b30=-5.4023E-5 ; let b31=7.326E-6   ; let b32=7.0036E-6
let b33=-3.0412E-6 ; let b34=-1.0853E-6
let b40=3.949E-7   ; let b41=03.029E-8  ; let b42=-3.8209E-7
let b43=1.0012E-7  ; let b44=4.7133E-8
let b50=-6.36E-10  ; let b51=-1.309E-9  ; let b52=6.048E-9
let b53=-1.1409E-9 ; let b54=-6.676E-10
let dsalt=salt-35
let sp0=b00+dsalt*(b01+dsalt*(b02+dsalt*(b03+b04*dsalt)))
let sp1=b10+dsalt*(b11+dsalt*(b12+dsalt*(b13+b14*dsalt)))
let sp2=b20+dsalt*(b21+dsalt*(b22+dsalt*(b23+b24*dsalt)))
let sp3=b30+dsalt*(b31+dsalt*(b32+dsalt*(b33+b34*dsalt)))
let sp4=b40+dsalt*(b41+dsalt*(b42+dsalt*(b43+b44*dsalt)))
let sp5=b50+dsalt*(b51+dsalt*(b52+dsalt*(b53+b54*dsalt)))
let
spice=sp0+potemp*(sp1+potemp*(sp2+potemp*(sp3+potemp*(sp4+potemp*sp5))))
set mode/last verify


|--****--****-*---*---***--***--|____spillane@pmel.noaa.gov____|
|-*__---*-----*--*-*--*--*-*--*-|_SCIENCE APPLICATIONS SUPPORT_|
|--***--*-----*-*---*-***--***--|____EPIC/Ferret/PlotPlus______|
|-----*-*-----*-*****-*----*----|__Room 2070 Bldg#3 NOAA/PMEL__|
|-****---****-*-*---*-*----*----|____Phone_:_(206)526-6780_____|