;
;
@common
;---------------------------------------------------------
; bornes de la grille
;----------------------------------------------------------
computegrid, 0.5, -89.5, 1, 1, 360, 180
;
iodir = '/usr1/antiope/smasson/DATA/'
;
;
jpk = 33
gdept = ncdf_lec('Levitus98_1m_01_12_S.nc', iodir = iodir, var = 'deptht')
gdepw = .5*(gdept+shift(gdept, 1))
gdepw[0] = 0
e3t = shift(gdepw, -1)-gdepw
e3t[jpk-1] = gdept[jpk-1]-gdepw[jpk-1]
e3w = gdept-shift(gdept, 1)
e3w[0] = 0
;
mer = ncdf_lec('Levitus98_1m_01_12_S.nc', iodir = iodir, var = 'vosaline', count = [jpi, jpj, jpk, 1])
mer = where(mer GT -50)
tmask = bytarr(jpi, jpj, jpk)
tmask[temporary(mer)] = 1b
;
jpt = 1
time = 0
;------------------------------------------------------
; lecture de la grille
;------------------------------------------------------
key_shift = 0
key_stride = [1, 1, 1]
;-------------------------------------------------------------
domdef
;
key_periodique = 1
triangles = triangule()
;----------------------------------------------------------
; bornes pour la lecture des donnees
;----------------------------------------------------------
jpidta = jpiglo
jpjdta = jpjglo
jpkdta = jpkglo
ixmindta = 0
ixmaxdta = jpidta-1
iymindta = 0
iymaxdta = jpjdta-1
izmindta = 0
izmaxdta = jpkdta-1