depth2pressure

 Pressure Computes pressure given the depth at some latitude
          P=Pressure(D,LAT) gives the pressure P (dbars) at a depth D (m)
          at some latitude LAT (degrees).

          This probably works best in mid-latitude oceans, if anywhere!

          Ref: Saunders, "Practical Conversion of Pressure to Depth",
              J. Phys. Oceanog., April 1981.

0001 function P80=depth2pressure(DPTH,XLAT)
0002 % Pressure Computes pressure given the depth at some latitude
0003 %          P=Pressure(D,LAT) gives the pressure P (dbars) at a depth D (m)
0004 %          at some latitude LAT (degrees).
0005 %
0006 %          This probably works best in mid-latitude oceans, if anywhere!
0007 %
0008 %          Ref: Saunders, "Practical Conversion of Pressure to Depth",
0009 %              J. Phys. Oceanog., April 1981.
0010 %
0011 
0012 %Notes: RP (WHOI) 2/Dec/91
0013 %         I copied this directly from the UNESCO algorithms.
0014 
0015 
0016 % CHECK VALUE: P80=7500.004 DBARS;FOR LAT=30 DEG., DEPTH=7321.45 METERS
0017 
0018       PLAT=abs(XLAT*pi/180.);
0019       D=sin(PLAT);
0020       C1=5.92E-3+(D.*D)*5.25E-3;
0021       P80=((1-C1)-sqrt(((1-C1).^2)-(8.84E-6*DPTH)))/4.42E-6;