New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Users/SetupNewConfiguration/OPABAT (diff) – NEMO

Changes between Version 6 and Version 7 of Users/SetupNewConfiguration/OPABAT


Ignore:
Timestamp:
2017-04-12T22:54:12+02:00 (7 years ago)
Author:
nicolasmartin
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Users/SetupNewConfiguration/OPABAT

    v6 v7  
    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 
     3Last edition: '''[[Wikinfo(changed_ts)]]''' by '''[[Wikinfo(changed_by)]]''' 
     4 
     5OPABAT 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 
     7However, 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 
     9The 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 
     13You can download the tools here: 
     14[[AttachmentNum(1)]] 
     15 
     16And find some documentation here (in French): [[AttachmentNum(2)]] [[AttachmentNum(3)]] 
     17 
     18An example of IDL procedures can be found at the bottom of the page. 
     19 
     20== Automatic contruction of the files 
     21 
     22The 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 
     24Up 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{{{ 
     28PARAMETER( 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{{{ 
     35PARAMETER( JPIBAT = 10800, JPJBAT = 6336) 
     36}}} 
     37 
     38Nota bene : If the field of the model covers areas in high latitudes, it is necessary to resort simultaneously to the two files. 
     39 
     40Two 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 
     44This 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 
     46To 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 
     48Example of parameters used for the construction of the bathymetry of configuration ORCA : 
    7149{{{ 
    7250CC MODEL DATA SIZE 
    73  
    7451CC --------------- 
    75  
    7652CC JPI : FIRST HORIZONTAL DIMENSION 
    77  
    7853CC JPJ : SECOND 
    79  
    8054CC JPK : NUMBER OF VERTICAL LEVELS 
    81  
    8255CC 
    83  
    8456    INTEGER jpiglo,jpjglo,jpidta,jpjdta,jpi,jpj,jpk 
    85  
    8657    PARAMETER(JPIGLO = 182,JPJGLO = 149) 
    87  
    8858C 
    8959}}} 
    9060 
    91 Second execution, bathymetry of finer resolution, in the case of ORCA, the precision was increased by a factor 12 in both direction of space thus a point of bathymetry is compared with a whole of 144 values :[[BR]] 
     61Second execution, bathymetry of finer resolution, in the case of ORCA, the precision was increased by a factor 12 in both direction of space thus a point of bathymetry is compared with a whole of 144 values : 
    9262 
    9363{{{ 
    9464   INTEGER NRES 
    95  
    9665     PARAMETER(NRES=12) 
    97  
    9866C si haute resolution 
    99  
    10067    PARAMETER(JPIGLO = 182*NRES,JPJGLO = 149*NRES) 
    101  
    10268    PARAMETER(JPI = jpiglo,JPJ = jpjglo ) 
    103  
    10469     PARAMETER(JPIDTA = jpiglo,JPJDTA = jpjglo ). 
    105  
    10670}}} 
    10771 
    10872The whole of the programs is called OPABAT :[[BR]] 
    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 
     83A makefile allows to obtain a binary. 
     84 
     85A file `opabat.h` allows to initialize the size of arrays. 
     86 
     87The namelist file allows to initialize the context of the experiment and the choice of the variables. 
    12588 
    12689Nota Bene :[[BR]] 
    127  
    12890The files of bathymetry have as a unit the levels of model, the file of outline is in whole value varying from 0 (ground) with jpk, ocean floor.[[BR]] 
    12991The file of control is to express in floating value between 0. and float(jpk), to see program batwri.[[BR]] 
    13092 
    13193== Manual file improvement == 
    132  
    13394 
    13495The realization of a bathymetry file in an automatic way is not fully satisfactory for the reasons which we indicated in introduction. It is thus necessary to check the coherence of the results visually, to correct and sometimes modify after validation by a simulation with the model. The manual correction is not easy and several methods were tried in the past. The systematic checking is tiresome and it must be limited to the principal thresholds as soon as the field comprises more than 10000 points of grid. The use of IDL allows us :[[BR]] 
     
    148109 
    149110 
    150 == Example == 
    151  
    152  
    153  
    154 {{{ 
     111== Example 
     112 
     113{{{#!idl 
    155114rhodes-rrup448: idl 
    156115 
     
    172131== Initialization part == 
    173132 
    174  
    175 {{{ 
     133{{{#!idl 
    176134IDL> @batini 
    177135 
     
    192150}}} 
    193151 
    194  
    195 == Launching reading procedures == 
    196  
    197  
    198 {{{ 
    199  
     152== Launching reading procedures 
     153 
     154{{{ 
    200155IDL> batexp 
    201  
    202156% Compiled module: BATEXP. 
    203  
    204157Lecture des fichiers, noms par défaut : global.bin, global.bat 
    205  
    206158-> -> -> -> -> -> voulez-vous : 
    207159lire les bathymetries initiales , y or n(implicit)? 
    208  
    209160:y 
    210  
    211161lecture des bathymetries initiales 
    212  
    213162% Compiled module: BATLEC. 
    214  
    215163bathymetrie haute resolution nom = global.bin, y(implicit) or n? 
    216  
    217164:y 
    218  
    219165haute resolution lue 
    220  
    221166% Compiled module: BATLEC2. 
    222  
    223167bathymetrie basse resolution nom = global.bat, y(implicit) or n? 
    224  
    225168: n 
    226  
    227169------>> rentrez un nom de fichier ? 
    228170: bathy_ebauche 
    229  
    230171OK lecture bathy_ebauche 
    231172basse resolution lue 
    232173}}} 
    233174 
    234 Possible treatments:[[BR]] 
    235  
    236 For more information on the channels, refer to the appendix A3[[BR]] 
    237  
     175Possible treatments: 
     176 
     177For more information on the channels, refer to the appendix A3 
    238178 
    239179{{{ 
    240180-> -> -> -> -> -> voulez-vous : 
    241181traiter les mares et chenaux, y or n(implicit)? 
    242  
    243182:n 
    244183}}} 
    245184 
    246 Selection of a under-field by clicking the mouse on a zone[[BR]] 
    247  
    248  
     185Selection of a under-field by clicking the mouse on a zone 
    249186 
    250187{{{ 
    251188-> -> -> -> -> -> voulez-vous : 
    252189positionner le zoom , y or n(implicit)? 
    253  
    254 : y 
    255  
     190: y 
    256191% Compiled module: BATCLIC. 
    257  
    258192-> -> -> -> positionner et cliquer la souris 
    259193}}} 
    260194 
    261  
    262  
    263195{{{ 
    264196-> -> -> -> -> -> le zoom vous convient-il, y(implicit) or n? 
     
    270202}}} 
    271203 
    272  
    273 Posting if desired the file of bathymetry and file of control with a higher resolution:[[BR]] 
    274  
    275  
    276  
    277 {{{ 
    278 -> -> -> -> -> -> voulez-vous : 
    279  
     204Posting if desired the file of bathymetry and file of control with a higher resolution: 
     205 
     206{{{ 
     207-> -> -> -> -> -> voulez-vous : 
    280208affichez le zoom de la basse resolution, y or n(implicit)? 
    281  
    282209:y 
    283  
    284210affichage 
    285211}}} 
    286212 
    287  
    288  
    289 {{{ 
    290 -> -> -> -> -> -> voulez-vous : 
    291  
     213{{{ 
     214-> -> -> -> -> -> voulez-vous : 
    292215afficher le zoom de la basse resolution avec chenaux, y or n(implicit)? 
    293216}}} 
    294217 
    295  
    296 For more information on the channels please refer to appendix A2:[[BR]] 
    297  
    298  
     218For more information on the channels please refer to appendix A2: 
    299219 
    300220{{{ 
    301221n 
    302  
    303 -> -> -> -> -> -> voulez-vous : 
    304  
     222-> -> -> -> -> -> voulez-vous : 
    305223affichez le zoom de la haute resolution, y or n(implicit)? 
    306  
    307224:y 
    308  
    309225affichage 
    310226}}} 
    311227 
    312  
    313 {{{ 
    314  
     228{{{ 
    315229Fin de la procedure BATEXP 
    316230}}} 
    317231 
    318  
    319  
    320 == Launching the working procedure == 
    321  
     232== Launching the working procedure 
    322233 
    323234{{{ 
    324235IDL> batmod 
    325  
    326236% Compiled module: BATMOD. 
    327  
    328237% LOADCT: Loading table Arc en Ciel / Phase 
    329  
    330238entrez le niveau? 
    331239}}} 
     
    333241One selected a level expressed in model level[[BR]] 
    334242 
    335  
    336243{{{ 
    33724420 
    338  
    339245on va visualiser le niveau 20 
    340  
    341246On peut changer de niveau en cliquant un bouton de la souris 
    342  
    343247-> -> -> -> -> -> positionner et cliquer la souris 
    344  
    345248bouton gauche -> informations 
    346  
    347249milieu -> modifier 
    348  
    349250droite -> sortir 
    350251}}} 
    351252 
    352253Action with the mouse 
    353  
    354254{{{ 
    355255-> -> -> -> changement de niveau de 0 a jpk ou n? 
    356  
    357256: 25 
    358257}}} 
    359258 
    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  
     259The 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 
     261Possibly, we can modify the levels posted. 
     262 
     263Two 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 
     266Action with the mouse 
     267 
     268{{{ 
     269pt ( 124, 120) BR= 25.5000 
     270Haute resolution min 24.6479 moy 25.9379 max 26.6315 
     271}}} 
     272 
     273Action with the mouse 
     274 
     275{{{ 
     276pt ( 125, 120) BR= 25.5000 
     277Haute 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{{{ 
     283bouton gauche --> monter (combler) new jk=jk-1 
     284milieu --> descendre (creuser) new jk=jk+1 
     285droite --> sortir 
     286}}} 
    370287 
    371288Action with the mouse[[BR]] 
    372289 
    373  
    374 {{{ 
    375 pt ( 124, 120) BR= 25.5000 
    376  
    377 Haute resolution min 24.6479 moy 25.9379 max 26.6315 
    378 }}} 
    379  
    380 Action with the mouse[[BR]] 
    381  
    382 {{{ 
    383 pt ( 125, 120) BR= 25.5000 
    384  
    385 Haute resolution min 25.8304 moy 26.6145 max 27.2046 
    386 }}} 
    387  
    388 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) :[[BR]] 
    389  
    390 {{{ 
    391 bouton gauche --> monter (combler) new jk=jk-1 
    392  
    393 milieu --> descendre (creuser) new jk=jk+1 
    394  
    395 droite --> sortir 
    396 }}} 
    397  
    398 Action with the mouse[[BR]] 
    399  
    400 {{pt ( 123, 120) BR= 24.5000 
    401  
     290{{{ 
     291pt ( 123, 120) BR= 24.5000 
    402292Haute resolution min 23.5273 moy 24.9224 max 26.0350{ 
    403  
    404293}}} 
    405294 
     
    408297{{{ 
    409298--> descendu de 24.5000 a 25.5000 
    410  
    411299pt ( 123, 120) BR= 25.5000 
    412  
    413300Haute resolution min 23.5273 moy 24.9224 max 26.0350 
    414301}}} 
    415302 
    416 Action with the mouse[[BR]] 
    417  
     303Action with the mouse 
    418304 
    419305{{{ 
     
    421307}}} 
    422308 
    423 Exit of the procedure and save the modifications in the under-field[[BR]] 
    424  
    425 Action with the mouse[[BR]] 
     309Exit of the procedure and save the modifications in the under-field 
     310 
     311Action with the mouse 
    426312 
    427313 
    428314{{{ 
    429315-> -> -> -> changement de niveau de 0 a jpk ou n? 
    430  
    431316: n 
    432  
    433317sortir 
    434318}}} 
    435319 
    436 == Safeguard == 
    437  
     320== Safeguard 
    438321 
    439322{{{ 
    440323Voulez-vous mettre a jour le tableau BATHY 
    441  
    442324en prenant en compte les corrections faite sur ce ZOOM ? 
    443  
    444325y or n ? 
    445  
    446326:y 
    447  
    448327sauvegarde 
    449  
    450328% Compiled module: BATSAV. 
    451329}}} 
    452330 
    453  
    454 We can number the islands before saving. (further informations, see appendix A3)[[BR]] 
     331We can number the islands before saving. (further informations, see appendix A3) 
    455332 
    456333{{{ 
    457334-> -> -> -> remplissage des iles, y ou n(implicit)? 
    458  
    459335Sauvegarde de bathy 
    460  
    461336------>> rentrez un nom de fichier ? 
    462  
    463337: nouvelle_bathy 
    464  
    465338procedure terminee: 
    466  
    467339Vous pouvez faire appel a d autres procedures : 
    468  
    469340@batini : Reinitialisation des dimensions et tableaux communs 
    470  
    471341batexp : lecture des fichiers, suppression des mares 
    472  
    473342et construction d un zoom 
    474  
    475343batmod : travail sur le zoom de la bathymetrie, niveau par niveau 
    476  
    477344batsav : remise des iles a -n et sauvegarde d un fichier bathy 
    478  
    479 -> -> -> -> -> -> voulez-vous : 
    480  
     345-> -> -> -> -> -> voulez-vous : 
    481346affichez le zoom de la nouvelle bathy, y or n(implicit)? 
    482  
    483347: n 
    484  
    485348Vous pouvez faire appel a d autres procedures … 
    486  
    487349IDL> 
    488350}}} 
     
    490352A call of the procedures can be reiterated as many times as wanted to make posting and correction of under-field then leave IDL. 
    491353 
    492  
    493  
    494 == Appendix == 
    495  
    496 == A1: Mouth of the rivers == 
    497  
     354== Appendix 
     355 
     356== A1: Mouth of the rivers 
    498357 
    499358An additional procedure RUNOFF is under development for, starting from a first file of outline and layout of the coasts, to correct the localization and the distribution (on one or more points of grid) the mouth of the rivers. 
    500359 
    501  
    502  
    503 == A2: Channels removal == 
    504  
     360== A2: Channels removal 
    505361 
    506362The outline file of bathymetry created automatically using FORTRAN programs involves the appearance in the file of ponds (major isolated points) or of channels (deep corridors). Initially a subroutine made it possible to fill them. It seemed more relevant to defer this algorithm in procedure BATEXP of IDL to make it possible to the user to operate a choice on the way of correcting them. 
     
    508364{{{ 
    509365-> -> -> -> -> -> voulez-vous : 
    510  
    511366traiter les mares et cheneaux, y or n(implicit)? 
    512  
    513 : y 
    514  
    515 -> -> -> -> -> -> voulez-vous : 
    516  
     367: y 
     368-> -> -> -> -> -> voulez-vous : 
    517369supprimer les points ocean quasi-isoles , y or n(implicit) ? 
    518  
    519 : y 
    520  
     370: y 
    521371affichage des points isoles, passage numero : 1 
    522  
    523372nombre de points isoles trouves : 5 
    524  
    525373continuer la suppression ? , y or n(implicit)? 
    526  
    527 : y 
    528  
     374: y 
    529375affichage des points isoles, passage numero : 2 
    530  
    531 : y 
    532  
     376: y 
    533377nombre de points isoles trouves : 2 
    534  
    535378continuer la suppression ? , y or n(implicit)? 
    536  
    537379: n 
    538  
    539380affichage des points isoles supprimes 
    540  
    541381nombre de points retouches EN TOUT sur 3 iterations : 8 
    542  
    543382enregistrer ces suppressions sur bathy? , y(implicit) or n? 
    544  
    545 : y 
    546  
     383: y 
    547384suppression enregistree 
    548  
    549 -> -> -> -> -> -> voulez-vous : 
    550  
     385-> -> -> -> -> -> voulez-vous : 
    551386elargissement des cheneaux, y or n(implicit) ? 
    552  
    553 : y 
    554  
     387: y 
    555388affichage des cheneaux, passage numero : 1 
    556  
    557389probleme chenal 140 102 entre continent/iles ! 18 0 0 
    558  
    559390nombre de points cheneaux trouves N-S : 0 
    560  
    561391nombre de points cheneaux trouves E-W : 1 
    562  
    563392continuer la correction ? , y or n(implicit)? 
    564  
    565393affichage des points cheneaux corriges 
    566  
    567394nombre de points N-S EN TOUT sur 1 iterations : 0 
    568  
    569395nombre de points E-W EN TOUT sur 1 iterations : 1 
    570  
    571396enregistrer ces modifications sur bathy? , y(implicit) or n? 
    572  
    573397modification enregistree 
    574398}}} 
    575399 
    576  
    577 == A3: Islands classification == 
     400== A3: Islands classification 
    578401 
    579402When the field has islands, for example fourteen for model ORCA, the corresponding areas are initially put at 0 by the programs and procedures, corresponding to continent. In the file of bathymetry of model OPA, these islands are numbered with negative values. Once bathymetry to finalize, one can thus carry out this operation starting from procedure BATSAV. This algorithm rises from programs worked out by Anne-Marie Tréguier (LPO). 
    580403 
    581  
    582404{{{ 
    583405IDL> batsav 
    584  
    585406-> -> -> -> remplissage des iles, y ou n(implicit)? 
    586  
    587 : y 
    588  
     407: y 
    589408Remplissage des iles 
    590  
    591409-> -> -> -> -> -> positionner et cliquer la souris 
    592  
    593410si probleme : droite --> sortir 
    594  
    595411-1 est valeur par defaut de l ile, OK, y or n ? 
    596  
    597  
    598  
    599 }}} 
    600  
    601  
     412}}} 
     413 
     414