1 | PRO mesh_from_file, model, file_name, ncdf_db, var_name |
---|
2 | @common |
---|
3 | @com_eg |
---|
4 | ; |
---|
5 | ; init grid, sf from file_name |
---|
6 | ; masks made in data_read |
---|
7 | ; |
---|
8 | |
---|
9 | IF strpos(ncdf_db, ':') GE 1 THEN directory = (str_sep(ncdf_db, ':'))[1] $ |
---|
10 | ELSE directory = ncdf_db |
---|
11 | |
---|
12 | s_file = directory+file_name |
---|
13 | |
---|
14 | print,' Model inits for ', model, ' from file', s_file |
---|
15 | |
---|
16 | sm_file = hom_idl+'grids/grids_'+model+'.nc' |
---|
17 | res = find(sm_file) |
---|
18 | |
---|
19 | IF res NE 'NOT FOUND' THEN BEGIN |
---|
20 | |
---|
21 | initncdf, s_file, GLAMBOUNDARY = glamboundary_box |
---|
22 | |
---|
23 | IF debug_w THEN print, ' Found mask from ',sm_file |
---|
24 | |
---|
25 | tmask = byte(read_ncdf('sftlf', 0, 0, /timestep, file = sm_file, /nostruct)) |
---|
26 | |
---|
27 | idx = where(tmask EQ valmask) |
---|
28 | |
---|
29 | IF idx(0) NE -1 THEN tmask(idx) = 0. |
---|
30 | idx = where(tmask LE 50.) |
---|
31 | tmask(idx) = 0. |
---|
32 | tmask = tmask < 1 |
---|
33 | tmask = 1-tmask |
---|
34 | triangles=triangule() |
---|
35 | |
---|
36 | ENDIF ELSE BEGIN |
---|
37 | initncdf, s_file, USEASMASK = var_name, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = 'depth' |
---|
38 | ENDELSE |
---|
39 | |
---|
40 | IF debug_w THEN print, 'vargrid = ',vargrid |
---|
41 | |
---|
42 | key_offset = [0, 0, 0] |
---|
43 | ; |
---|
44 | ; indice i pour grille j moyenne zonale |
---|
45 | ; |
---|
46 | diaznl_idx = 1 |
---|
47 | |
---|
48 | return |
---|
49 | END |
---|
50 | |
---|