pro construct_levitus1998_sss_ogcm_v2d

; Purpose
; -------
; Constructs text files containing annual-mean and monthly Levitus 1998 SSSs,
; in the format read by the CSIRO OGCM
;
; History
; -------
; 2004 Sep 14	Steve Phipps	Original version
; 2004 Oct 11	Steve Phipps	Modified for v2d SSSs

; Get Levitus 1998 SSSs
ncid = ncdf_open("levitus1998_sss_ogcm_ann_mid_v2d.nc")
annid = ncdf_varid(ncid, "sssann")
monid = ncdf_varid(ncid, "sss")
ncdf_varget, ncid, annid, sssann
ncdf_varget, ncid, monid, sss
ncdf_close, ncid

; Write annual means to output file
openw, 1, "sss.dat.levitus1998_v2d_ann"
for i = 0, 11 do begin
  header = $
    '                                                                      '
  strput, header, "Annual-mean Levitus 1998 v2d SSSs"
  printf, 1, format='(a70)', header
  printf, 1, format='(8e10.3)', sssann
endfor
close, 1

; Write monthly values to output file
openw, 1, "sss.dat.levitus1998_v2d_mid"
for i = 0, 11 do begin
  header = $
    '                                                                      '
  strput, header, "Levitus 1998 v2d SSSs, month " + strtrim(string(i+1)), 0
  printf, 1, format='(a70)', header
  printf, 1, format='(8e10.3)', sss(*, *, i)
endfor
close, 1

end