1 | ; |
---|
2 | ; |
---|
3 | @common |
---|
4 | ;------------------------------------------------------------ |
---|
5 | ;---------------------------------------------------------- |
---|
6 | ; bornes de la grille |
---|
7 | ;---------------------------------------------------------- |
---|
8 | jpi = 1020 |
---|
9 | jpj = 1286 |
---|
10 | jpk = 1 |
---|
11 | gdept = 0 |
---|
12 | gdepw = 0 |
---|
13 | e3t = 10 |
---|
14 | e3w = 10 |
---|
15 | tempsun = systime(1) |
---|
16 | iodir = '/usr1/antiope/smasson/DATA/MERCATOR/' |
---|
17 | cdfid = ncdf_open(iodir+'test_mercator.nc') |
---|
18 | ncdf_varget, cdfid, 'glamt', glamt |
---|
19 | ncdf_varget, cdfid, 'gphit', gphit |
---|
20 | ncdf_varget, cdfid, 'e1t' , e1t |
---|
21 | ncdf_varget, cdfid, 'e2t', e2t |
---|
22 | ncdf_varget, cdfid, 'glamf', glamf |
---|
23 | ncdf_varget, cdfid, 'gphif', gphif |
---|
24 | ncdf_varget, cdfid, 'e1f', e1f |
---|
25 | ncdf_varget, cdfid, 'e2f', e2f |
---|
26 | ncdf_varget, cdfid, 'bat', bat |
---|
27 | ncdf_varget, cdfid, 'tmask', tmask |
---|
28 | ncdf_varget, cdfid, 'sss', sss |
---|
29 | ncdf_close, cdfid |
---|
30 | print, 'temps lecture fichier NetCdf 1ere fois', systime(1)-tempsun |
---|
31 | tempsun = systime(1) |
---|
32 | cdfid = ncdf_open(iodir+'test_mercator.nc') |
---|
33 | ncdf_varget, cdfid, 'glamt', glamt |
---|
34 | ncdf_varget, cdfid, 'gphit', gphit |
---|
35 | ncdf_varget, cdfid, 'e1t' , e1t |
---|
36 | ncdf_varget, cdfid, 'e2t', e2t |
---|
37 | ncdf_varget, cdfid, 'glamf', glamf |
---|
38 | ncdf_varget, cdfid, 'gphif', gphif |
---|
39 | ncdf_varget, cdfid, 'e1f', e1f |
---|
40 | ncdf_varget, cdfid, 'e2f', e2f |
---|
41 | ncdf_varget, cdfid, 'bat', bat |
---|
42 | ncdf_varget, cdfid, 'tmask', tmask |
---|
43 | ncdf_varget, cdfid, 'sss', sss |
---|
44 | ncdf_close, cdfid |
---|
45 | print, 'temps lecture fichier NetCdf 2eme fois', systime(1)-tempsun |
---|
46 | glamu = glamt & glamv = glamt |
---|
47 | gphiu = gphit & gphiv = gphit |
---|
48 | e1u = e1t & e1v = e1t |
---|
49 | e2v = e1t & e2v = e2t |
---|
50 | |
---|
51 | jpt = 1 |
---|
52 | time = 0 |
---|
53 | |
---|
54 | ;------------------------------------------------------------- |
---|
55 | ;------------------------------------------------------------- |
---|
56 | domdef |
---|
57 | triangles = triangule() |
---|
58 | ;---------------------------------------------------------- |
---|
59 | ; environnement graphique (window et ps) |
---|
60 | ;---------------------------------------------------------- |
---|
61 | ;---------------------------------------------------------- |
---|
62 | sauveps = 0 |
---|
63 | ;---------------------------------------------------------- |
---|
64 | ; ajout perso |
---|
65 | ;---------------------------------------------------------- |
---|
66 | key_performance = 1 |
---|
67 | plt, sss, 33, int = .2 |
---|