1 | | |
2 | | == OPABAT Tool == |
3 | | |
4 | | |
5 | | OPABAT is a set of Fortran and IDL procedures allowing to build a bathymetry filed on the grid of OPA, the oceanic dynamic component of the NEMO system. The bathymetry's construction is an essential step when setting up an new configuration of the sydtem : the choice of the trsehholds, of the passage's width and dorsal's height, to sum up the choice of the realism of the bathymetry at the given resolution conditions at the first order the oceanic flows computed by the model. This aspect of the problem imposes a manual intervention during the construction, in order to control the precision.[[BR]] |
6 | | |
7 | | However, this manual intervention is increasingly light with the raise in resolution, the limit of resolution equalizing the precision of the data. With OPA model, bathymetry is calculated in vertical levels at the point of mass T of the C-Arakawa grid, with for value, zero, for the continents, of 1 with jpk for the oceanic bottom.[[BR]] |
8 | | |
9 | | The development of the field of bathymetry is carried out in two times :[[BR]] |
10 | | |
11 | | - By the automatic construction of an outline through FORTRAN programs by using a bathymetric data file high resolution.[[BR]] |
12 | | - In one second step by the manual correction under IDL of the field.[[BR]] |
13 | | |
14 | | You can download the tools here:[[BR]] |
15 | | [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/Users/SetupNewConfiguration/OPABAT/IDL_OPABAT3.tar IDL_OPABAT3] |
16 | | |
17 | | And find some documentation here (in French):[[BR]] |
18 | | [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/Users/SetupNewConfiguration/OPABAT/opabat_RAPPOR.pdf opabat_RAPPORT][[BR]] |
19 | | [https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/Users/SetupNewConfiguration/OPABAT/opabat_ANNEXE.pdf opabat_ANNEXE][[BR]] |
20 | | |
21 | | |
22 | | An example of IDL procedures can be found at the bottom of the page.[[BR]] |
23 | | |
24 | | |
25 | | == Automatic contruction of the files == |
26 | | |
27 | | |
28 | | The construction schedules of the outline file of bathymetry were conceived for fields of the model with relative small dimension, 100x100 with 1000x1000 points on the horizontal one, based on the handling of large data files and on expensive interpolations in computing times. These programs were initially conceived in 1993 to build the bathymetry of the first version world ocean OPA and were based on a set of FORTRAN programs and programs of layouts NCAR to validate the results. They were re-examined into 1997/98 per Anne-Marie Tréguier (LPO) within the framework of the CLIPPER project and with the LODYC for construction, inter alia, of the bathymetry for model ORCA. The graphic part and certain developments was given up or appreciably modified by the use of IDL which adds employment opportunities, in particular the use of the mouse. The programs, developed in FORTRAN 77 for an initial application, then used in several other contexts of the years afterwards, do not present the desired finished aspect and comprise sometimes abandoned developments, obsolete comments and certainly of the errors which we will perhaps have time and courage to begin again.[[BR]] |
29 | | |
30 | | Up to now, two data files of bathymetry were used :[[BR]] |
31 | | |
32 | | - ETOP5[[BR]] |
33 | | Released by the Marine Geology and Geophysics Division/National Geophysical Dated Center, covering the sphere and 5 minutes per 5 minutes of[[BR]] |
34 | | |
35 | | {{{ |
36 | | PARAMETER( JPITP5 = 4320, JPJTP5 = 2160 ) |
37 | | }}} |
38 | | |
39 | | - Smith and Sandwel[[BR]] |
40 | | file covering the sphere of 72 S with 72 N and integrating the recent data satellitales, 2 minutes resolution per 2 minutes: [[BR]] |
41 | | |
42 | | (David T. Sandwell, Walter H. F. Smith, Stuart M. Smith, and Christopher Small,[[BR]] |
43 | | |
44 | | Institute for Geophysics and Planetary Physics, SIO, Geodynamics Laboratory, NOAA - Geological Data Center, SIO - Lamont Doherty Earth Observatory)[[BR]] |
45 | | |
46 | | (look at [http://topex.ucsd.edu/marine_topo/mar_topo.html http://topex.ucsd.edu/marine_topo/mar_topo.html] ).[[BR]] |
47 | | |
48 | | |
49 | | {{{ |
50 | | PARAMETER( JPIBAT = 10800, JPJBAT = 6336) |
51 | | }}} |
52 | | |
53 | | |
54 | | Nota bene : If the field of the model covers areas in high latitudes, it is necessary to resort simultaneously to the two files.[[BR]] |
55 | | |
56 | | Two methods of interpolation are employed:[[BR]] |
57 | | |
58 | | A bilinear interpolation using the four points of data framing the point of grid.[[BR]] |
59 | | Variable lopa_method=.false. in the namelist file.[[BR]] |
60 | | |
61 | | An average of the data of bathymetry contained in an areaof the grid of the model.[[BR]] |
62 | | Variable lopa_method=.true. in the namelist file.[[BR]] |
63 | | |
64 | | This last method is more complex to carry out when the model uses an irregular grid with a symmetry east- west. It is recommended when the model is of weak resolution. The bilinear method is necessary when the resolution approaches the resolution of the data, the averages becoming nonsignificant.[[BR]] |
65 | | |
66 | | To use IDL utility, it is necessary to build two files of outline of bathymetry, i.e. to carry out these programs twice to obtain a file with the resolution of the model, a file of highher resolution which will be used as file of control in the choices of the levels of bathymetry and will replace the use of the initial files of topography which are often difficult to handle and covering all the sphere. Typically the file of the model will be built by using an average by area, the file of control, to which the resolution is closer to the data, with a bilinear interpolation.[[BR]] |
67 | | |
68 | | Example of parameters used for the construction of the bathymetry of configuration ORCA :[[BR]] |
69 | | |
70 | | |
| 1 | == OPABAT Tool |
| 2 | |
| 3 | Last edition: '''[[Wikinfo(changed_ts)]]''' by '''[[Wikinfo(changed_by)]]''' |
| 4 | |
| 5 | OPABAT is a set of Fortran and IDL procedures allowing to build a bathymetry filed on the grid of OPA, the oceanic dynamic component of the NEMO system. The bathymetry's construction is an essential step when setting up an new configuration of the sydtem : the choice of the trsehholds, of the passage's width and dorsal's height, to sum up the choice of the realism of the bathymetry at the given resolution conditions at the first order the oceanic flows computed by the model. This aspect of the problem imposes a manual intervention during the construction, in order to control the precision. |
| 6 | |
| 7 | However, this manual intervention is increasingly light with the raise in resolution, the limit of resolution equalizing the precision of the data. With OPA model, bathymetry is calculated in vertical levels at the point of mass T of the C-Arakawa grid, with for value, zero, for the continents, of 1 with jpk for the oceanic bottom. |
| 8 | |
| 9 | The development of the field of bathymetry is carried out in two times : |
| 10 | - By the automatic construction of an outline through FORTRAN programs by using a bathymetric data file high resolution. |
| 11 | - In one second step by the manual correction under IDL of the field. |
| 12 | |
| 13 | You can download the tools here: |
| 14 | [[AttachmentNum(1)]] |
| 15 | |
| 16 | And find some documentation here (in French): [[AttachmentNum(2)]] [[AttachmentNum(3)]] |
| 17 | |
| 18 | An example of IDL procedures can be found at the bottom of the page. |
| 19 | |
| 20 | == Automatic contruction of the files |
| 21 | |
| 22 | The construction schedules of the outline file of bathymetry were conceived for fields of the model with relative small dimension, 100x100 with 1000x1000 points on the horizontal one, based on the handling of large data files and on expensive interpolations in computing times. These programs were initially conceived in 1993 to build the bathymetry of the first version world ocean OPA and were based on a set of FORTRAN programs and programs of layouts NCAR to validate the results. They were re-examined into 1997/98 per Anne-Marie Tréguier (LPO) within the framework of the CLIPPER project and with the LODYC for construction, inter alia, of the bathymetry for model ORCA. The graphic part and certain developments was given up or appreciably modified by the use of IDL which adds employment opportunities, in particular the use of the mouse. The programs, developed in FORTRAN 77 for an initial application, then used in several other contexts of the years afterwards, do not present the desired finished aspect and comprise sometimes abandoned developments, obsolete comments and certainly of the errors which we will perhaps have time and courage to begin again. |
| 23 | |
| 24 | Up to now, two data files of bathymetry were used : |
| 25 | - ETOP5 |
| 26 | Released by the Marine Geology and Geophysics !Division/National Geophysical Dated Center, covering the sphere and 5 minutes per 5 minutes of |
| 27 | {{{ |
| 28 | PARAMETER( JPITP5 = 4320, JPJTP5 = 2160 ) |
| 29 | }}} |
| 30 | - Smith and Sandwel |
| 31 | file covering the sphere of 72 S with 72 N and integrating the recent data satellitales, 2 minutes resolution per 2 minutes: |
| 32 | (David T. Sandwell, Walter H. F. Smith, Stuart M. Smith, and Christopher Small, Institute for Geophysics and Planetary Physics, SIO, Geodynamics Laboratory, NOAA - Geological Data Center, SIO - Lamont Doherty Earth Observatory) |
| 33 | (look at [http://topex.ucsd.edu/marine_topo/mar_topo.html http://topex.ucsd.edu/marine_topo/mar_topo.html] ). |
| 34 | {{{ |
| 35 | PARAMETER( JPIBAT = 10800, JPJBAT = 6336) |
| 36 | }}} |
| 37 | |
| 38 | Nota bene : If the field of the model covers areas in high latitudes, it is necessary to resort simultaneously to the two files. |
| 39 | |
| 40 | Two methods of interpolation are employed: |
| 41 | - A bilinear interpolation using the four points of data framing the point of grid. Variable `lopa_method=.false.` in the namelist file. |
| 42 | - An average of the data of bathymetry contained in an areaof the grid of the model. Variable `lopa_method=.true.` in the namelist file. |
| 43 | |
| 44 | This last method is more complex to carry out when the model uses an irregular grid with a symmetry east- west. It is recommended when the model is of weak resolution. The bilinear method is necessary when the resolution approaches the resolution of the data, the averages becoming nonsignificant. |
| 45 | |
| 46 | To use IDL utility, it is necessary to build two files of outline of bathymetry, i.e. to carry out these programs twice to obtain a file with the resolution of the model, a file of highher resolution which will be used as file of control in the choices of the levels of bathymetry and will replace the use of the initial files of topography which are often difficult to handle and covering all the sphere. Typically the file of the model will be built by using an average by area, the file of control, to which the resolution is closer to the data, with a bilinear interpolation. |
| 47 | |
| 48 | Example of parameters used for the construction of the bathymetry of configuration ORCA : |
109 | | |
110 | | -OPABAT : main program, reading of the file namelist, call of the routines.[[BR]] |
111 | | -DOMHGR : programme OPA, définition de la grille horizontale et des facteurs d'échelle.[[BR]] |
112 | | -HGRCOO : program OPA, definition of the horizontal grid and scale factors.programme OPA lecture de la grille horizontale si NTOPO=1.[[BR]] |
113 | | -DOMZGR : program OPA, definition of the vertical grid.[[BR]] |
114 | | -DOMMSK : programe OPA , masks computing.[[BR]] |
115 | | -BATLEG : reading of the data file.[[BR]] |
116 | | -BATGUE : interpolation and creation of the outline file.[[BR]] |
117 | | -BATMOD : possible corrections of the file.[[BR]] |
118 | | -BATWRI : saving the file.[[BR]] |
119 | | |
120 | | A makefile allows to obtain a binary.[[BR]] |
121 | | |
122 | | A file opabat.h allows to initialize the size of arrays.[[BR]] |
123 | | |
124 | | The namelist file allows to initialize the context of the experiment and the choice of the variables.[[BR]] |
| 73 | - OPABAT : main program, reading of the file namelist, call of the routines. |
| 74 | - DOMHGR : programme OPA, définition de la grille horizontale et des facteurs d'échelle. |
| 75 | - HGRCOO : program OPA, definition of the horizontal grid and scale factors.programme OPA lecture de la grille horizontale si NTOPO=1. |
| 76 | - DOMZGR : program OPA, definition of the vertical grid. |
| 77 | - DOMMSK : programe OPA , masks computing. |
| 78 | - BATLEG : reading of the data file. |
| 79 | - BATGUE : interpolation and creation of the outline file. |
| 80 | - BATMOD : possible corrections of the file. |
| 81 | - BATWRI : saving the file. |
| 82 | |
| 83 | A makefile allows to obtain a binary. |
| 84 | |
| 85 | A file `opabat.h` allows to initialize the size of arrays. |
| 86 | |
| 87 | The namelist file allows to initialize the context of the experiment and the choice of the variables. |
360 | | The under-field is posted with a pallet of colors corresponding to the two levels directly above (hot colors) or below (cold colors) the level selected (yellow color), in isolignes black contours of the file of control :[[BR]] |
361 | | |
362 | | |
363 | | Possibly, we can modify the levels posted.[[BR]] |
364 | | |
365 | | Two modes can be selected via the mouse :[[BR]] |
366 | | |
367 | | |
368 | | - An indicative mode, while clicking on the mouse, we post the level of an area and the maximum and minimum values of the file of control in this area.[[BR]] |
369 | | |
| 259 | The under-field is posted with a pallet of colors corresponding to the two levels directly above (hot colors) or below (cold colors) the level selected (yellow color), in isolignes black contours of the file of control : |
| 260 | |
| 261 | Possibly, we can modify the levels posted. |
| 262 | |
| 263 | Two modes can be selected via the mouse : |
| 264 | - An indicative mode, while clicking on the mouse, we post the level of an area and the maximum and minimum values of the file of control in this area. |
| 265 | |
| 266 | Action with the mouse |
| 267 | |
| 268 | {{{ |
| 269 | pt ( 124, 120) BR= 25.5000 |
| 270 | Haute resolution min 24.6479 moy 25.9379 max 26.6315 |
| 271 | }}} |
| 272 | |
| 273 | Action with the mouse |
| 274 | |
| 275 | {{{ |
| 276 | pt ( 125, 120) BR= 25.5000 |
| 277 | Haute resolution min 25.8304 moy 26.6145 max 27.2046 |
| 278 | }}} |
| 279 | |
| 280 | - A modification mode, where, with the mouse, we increase or decrease by a model level, example: passage from 24 to 25 for the item (123, 120) : |
| 281 | |
| 282 | {{{ |
| 283 | bouton gauche --> monter (combler) new jk=jk-1 |
| 284 | milieu --> descendre (creuser) new jk=jk+1 |
| 285 | droite --> sortir |
| 286 | }}} |