ADCP de Coque Campagne DRAKE A bord du POLARSTERN Mars
Contents
1. NMEA Port Setup NMEANavComEnable TRUE NmeaNavComPortName COM3 NmeaNavComBaudRate 9600 NmeaNavComParity NOPARITY NmeaNavComStopBits 1 NmeaNavComDataBits 8 NMEARPHComEnable FALSE NmeaRPHComPortName None NmeaRPHComBaudRate 4800 NmeaRPHComParity NOPARITY NmeaRPHComStopBits 1 NmeaRPHComDataBits 8 Serial Port for Binary Ensemble Data Output BinaryEnsembleOutputComEnable FALSE BinaryEnsembleOutputComPortName None BinaryEnsembleOutputComBaudRate 9600 BinaryEnsembleOutputComParity NOPARITY BinaryEnsembleOutputComStopBits 1 BinaryEnsembleOutputComDataBits 8 BinaryEnsembleOutputDataType 0 none l enr 2 enx 3 sta 4 lta 0 BinaryEnsembleOutputRefVelType 0 none 1 Bottom 2 Mean 0 BinaryEnsembleOutputStartBin 1 BinaryEnsembleOutputEndBin 4 BinaryEnsembleOutputMeanStartBin 1 BinaryEnsembleOutputMeanEndBin 4 BinaryEnsembleOutputLeader 0 no 1 yes TRUE BinaryEnsembleOutputBottomTrack 0 no 1 yes TRUE BinaryEnsembleOutputNavigation 0 no 1 yes TRUE BinaryEnsembleOutputVelocity 0 no 1 yes TRUE BinaryEnsembleOutputIntensity 0 no l yes TRUE BinaryEnsembleOutputCorrelation 0 no 1 yes TRUE BinaryEnsembleOutputPercentGood 0 no 1 yes TRUE BinaryEnsembleOutputStatus 0 no 1 yes TRUE Serial Port for ASCII Ensemble Data Output AsciiEnsemble2. Predictions are based on global and or regional barotropic inverse tidal solutions obtained with OTIS and available for download at http www coas oregonstate edu research po research tide KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 2 MODEL FILE FORMAT The files are big endian binary recorded in bytes on a UNIX system There is one header record which gives the grid size n grid cells in longitude and m grid cells in latitude number of constituents nc limits of the area theta_lim 2 phi_lim 2 and constituent names c_id nc The header is followed by nc records each giving the either elevation fields m ocean tide OR transports fields m 2 s for one constituent Everything is given as complex amplitudes so that the partial tide for a single constituent of frequency w is given by h t x Re h x exp i w t t0 VO tO where VO t0 is the astronomical argument for the constituent at t0 Note that with the usual conventions amplitude and phase are given by amp h phase atan Im h Re h To read with a simple FORTRAN program integer n m nc real theta_lim 2 phi_lim 2 character 4 c_id 21 complex allocatable h uv complex allocatable u v open unit 1 file h_tpxo6 2 form unformatted X status old convert big_endian CEE header read 1 n m nc theta_lim phi_lim
3. La figure suivante montre le r sultat de ce fit du premier degr pour l ensemble de la campagne pour laquelle on disposait de 90 points de calibration seulement On constate que apr s la rotation la pente de la droite passant au mieux dans le nuage de points est quasi nulle que la phase est galement quasiment annul e et que l cart type est fortement diminu Ces fichiers d angles de correction ont donc t utilis s pour effectuer une rotation des donn es de navigation C est sur ces donn es corrig es qu ont t faites la v rification et la correction de l alignement du transducteur m thode du water tracking 16 Drake 2009 raw total phase 9 fixes et fit ordre 1 avant rotation Drake 2009 rotif total phase 9 fixes et fit ordre 1 apres rotation r r r T r r r T r T r Avant rotation fit 5 5731 0 0970 latitude std 0 2248 nb points 90 Apres rotation fit 1 1538 0 0201 latitude std 0 0462 nb points 94 si 2 4 ak 2 P Pl M ala A g g a de x k x x d M da E of spa nea Frs RU ag ng eae a a a tx t 2 Keg xO ik 2 x Pi sek 4 at 40 L L L L L 1 L L 40 L L 1 L 63 62 61 60 59 57 56 55 54 63 62 61 60 59 L L L 57 56 55 54 latitude Fic 6 2 DRAKE 2009 Trac de la phase en fonction de la l
4. alignement du transducteur Mise en forme d finitive Introduction des donn es de navigation dans la base ADCP Contr le de la qualit des donn es Construction de diff rentes grilles spatiales et ou temporelles permettant de tracer les graphiques de pr sentation et de cr er des fichiers ascii pour la communication Chapitre 3 Chargements des fichiers bruts Le logiciel VMDAS permet l acquisition de tout un ensemble de fichiers dont seuls certains ayant d j subi un pr traitement seront utilis s ENR ENS ENX Fichiers totaux monopings de tous les ensembles donn es brutes Beam Coordinates NIR Fichiers des trames NMEA de navigation LTA STA Fichiers moyenn s par le logiciel VMDAS partir des fichiers monopings LTA sur la longue dur e choisie dans le fichier de configuration 20 mn STA sur la courte dur e 1mn donn es pr trait es Earth Coordinates Nous utiliserons les fichiers STA dans la suite Ces fichiers sont cr s chaque arr t remise en route du syst me d acquisition LOG Fichier Log enregistrant les commandes RDI ex cut es chaque remise en route du logiciel d acquisition Un exemple est donn en ANNEXE III VMO Copie du fichier configuration cr e chaque remise en route du logiciel d acquisition Le logiciel CODAS permet le chargement en deux temps partir des fichiers STA ou LTA quivalents pour le traitement CODAS La premi re op ra
5. 00 to 2009 04 04 04 45 00 DRAKE 3 2009 04 04 04 45 00 to 2009 04 08 22 22 42 Plusieurs outils permettent de v rifier le contenu de la base de donn es et d obtenir des informa tions utiles La commande Ist block permet d obtenir le nombre la longueur nombre de profils et les dates des extr mit s des blocs Database summary for drake Number of blocks 63 0 n 400 drake001 blk 2009 03 22 17 54 42 to 2009 03 23 00 33 39 1 n 400 drake002 blk 2009 03 23 00 34 40 to 2009 03 23 07 13 42 La commande showdb permet d obtenir la structure pr cise de la base de donn es Les informa tions classiques sont contenues dans les structures 35 Configuration 1 37 et 38 Ancillary 1 et 2 L exemple ci dessous donne la liste des variables de la base num name units v ty ac ac0 acl 0 DEPTH m 3 1 17192 160 7 AMP SOUND SCAT none 1 2 0 0 8 U m s 3 2 1 1 9 V m s 3 2 2 2 32 PROFILE COMMENTS none 2 2 3 3 33 BLOCK COMMENTS none 2 1 0 0 34 PROFILE FLAGS none 1 2 4 4 35 CONFIGURATION 1 none 11 1 17352 76 36 CONFIGURATION 2 none 11 1 0 0 37 ANCILLARY 1 none 11 2 5 5 38 ANCILLARY 2 none 11 2 6 6 39 ACCESS VARIABLES none 11 2 7 T 54 W m s 3 2 8 8 55 ERROR VEL m s 3 2 9 9 56 PERCENT GOOD none 1 2 10 10 57 PERCENT 3 BEAM none 1 2 11 11 58 SPECTRAL WIDTH none 1 2 12 12 59 U STD DEV m s 3 2 13 13 60 V STD DEV m s 3 2 14 14 61 W STD DEV m s 3 2 15 15 62 EV STD DEV m s 3 2 16 16 63 AMP STD DEV none 0 2 17 17 64
6. 2009 03 22 13 04 10 863 CF CF 11110 Set Ctrl Flags e p b s gt 2009 03 22 13 04 10 881 CF 01110 2009 03 22 13 04 10 899 CF 01110 gt 2009 03 22 13 04 10 918 gt Setting the sensor source 2009 03 22 13 04 10 918 EZ 2009 03 22 13 04 11 027 EZ EZ 10111010 Sensor Source c d h p r s t u gt 2009 03 22 13 04 11 027 EZ10222010 2009 03 22 13 04 11 063 EZ10222010 gt 2009 03 22 13 04 11 063 gt Disabling ADCP internal heading or tilt corrections 2009 03 22 13 04 11 063 EAO 2009 03 22 13 04 11 081 EAO gt 2009 03 22 13 04 11 100 EVO 2009 03 22 13 04 11 118 EVO gt 2009 03 22 13 04 11 118 EIO 2009 03 22 13 04 11 154 EIO gt 2009 03 22 13 04 11 154 EJO 2009 03 22 13 04 11 172 EJO gt 2009 03 22 13 04 11 172 gt Setting beam radial output mode 2009 03 22 13 04 11 190 EX00000 2009 03 22 13 04 11 209 EX00000 gt 2009 03 22 13 04 12 136 CS 2009 03 22 13 39 06 345 NMEA Nav communication timeout 2009 03 22 13 39 16 567 NMEA Nav communication timeout 2009 03 22 13 40 24 003 NMEA Nav communication timeout 35 36 KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK Proc_gyro doc FAR xxx Notes on running Proc_gyro m to create a model of gyrocompass behaivor STARE KAKA KKKAX xxx For additional notes on gyrocompass RAA xxx modelling see Gyro doc PARA KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK This routine is set up to create a mo
7. Die Ausgabe von Heading Roll und Pitch erfolgt ber Synchrosignal und digital Das MINS Interface Fa WERUM wandelt die Daten die auf 3 Datenbussen mit unterschiedlicher update rate ankommen und speist sie mit 1 Hz in PODAS ein MINS liefert Broad Cast Bus RS422A 200 Hz update rate Baudrate 307 2k HDLC Protocol Heading Roll Pitch Heading Roll Pitch Rates Status Broad Cast Bus II RS422A 5 Hz update rate Baudrate 307 2k HDLC Protocol Velocity Vx Vy Vz Speedlog velocity Position Status Course Bus RS422A Es gibt 2 verschiedene M glichkeiten der GPS Empf nger die die MINS st tzen Als zugelassene Version wird die MINS mit einem Raytheon Antennenempf nger Raystar 112 Product Code M93537 12 Kan le betrieben In wissenschaftlichen Fahrtgebieten wird die MINS durch einen Geod tischen Zweifrequenzempf nger Trimble MS750 gest tzt Es kann immer nur eine MINS Daten an PODAS weiterleiten Es handelt sich um getrennte Anlagen mit jeweils eigenem System auch GPS jeweils einen Raystar und einen Trimble Empf nger 30 ANNEXE II Expert only optio Turned0ffBeam 0 PashrImuFlagUseNormalInterpretation TRU ADCP Port Setup CONFIGURATION DE 1 ADCP ns La AdcpComPortName COM1 AdcpComBaudRate 96 00 AdcpComParity NOPARITY AdcpComStopBits 1 AdcpComDataBits 8 AdcpConfigFilename C Program Files RD Instruments VmDas OS38NBDEF txt
8. Prof D but Fin Nb points Nb points Amplitude Amplitude Phase DeltaT T DRAKE EU 03 09 17 42 64 68 1 0196 0 9997 0 1221 0 0504 5 36 5 31 18243 PEN ae eet Au final apr s obtention de ces donn es de calibration soit un fichier d angles pour la correction du cap puis un angle phase et une amplitude constants pour corriger l alignement du transducteur on effectue la rotation des donn es de la base programme rotate Il faut ensuite extraire le fichier nav adcpsect partir des donn es ayant subi la rotation puis calculer la vitesse absolue de la couche de r f rence refabs fichier ref et ensuite lisser ce dernier fichier smoothr fichier sm Le logiciel CODAS permet des tests de qualit globale sur l ensemble des profils Ils permettent surtout de visualiser la diff rence entre les donn es en route et en station DRAKE 2008 eR B ARR Bish Underway DRAKE 2000 ERROR AMSA Basn a 8 std_dov U ere onae zo PRA ARGO Ah an 4 2009 5 9 16 13 z vinda DRAKE 2009 AB On Sigi Dash Underway so z so so _ 100 Wolf ol E fi E PA asol 5 160 J 5 so aft 2001 200 si mL 080100 iso 200 20 o 20 40 Cg mean AMP SOUND SCAT std dev AMP SOUND SCAT prake 2088 SIR bn Ba Underway 40 80 80 20 30 4 mean PERCENT GOOD sl dev PERCENT GOOD prake 2065 Said Gn SIA Bash Underway 2009 5 9
9. SetTiltSensorType TRUE TiltSensorType 0 internal l external 1 SetProcessingMode TRUE BandwidthType 0 Wide 1 Narrow 1 ADCPTimeBetweenEnsemblesSel 1 ADCPTimeBetweenEnsembles 1 UserNavScreening FALSE UserTransform FALSE Shiptrack Options ShipTrack1Source 0 Nav 1 BT 2 WT 3 Layer 0 ShipTrack2Source 0 Nav 1 BT 2 WT 3 Layer 1 ShipTrack1RedStickEnable TRUE ShipTrack1GreenStickEnable TRUE ShipTrack1BlueStickEnable TRUE ShipTrack2RedStickEnable TRUE ShipTrack2GreenStickEnable TRUE ShipTrack2BlueStickEnable TRUE ShipTracklRedBin 1 ShipTracklGreenBin 2 ShipTrack1BlueBin 3 ShipTrack2RedBin 1 ShipTrack2GreenBin 2 ShipTrack2BlueBin 3 ShipTrack1DisplaySelect 0 Lat Lon 1 Distance ShipTrack2DisplaySelect 0 Lat Lon 1 Distance ShipTracklWaterLayerStartBin 3 ShipTracklWaterLayerEndBin 5 ShipTrack2WaterLayerStartBin 3 ShipTrack2WaterLayerEndBin 5 ShipTrackDistanceUnit 0 S J G UserWinAdcpPath C Program Files RD Instruments WinAdcp WinAdcp exe Type 0 enr l enx 2 sta 3 lta 3 0 0 ea 33 34 ANNEXE III 1 0 log for RD Instruments VmDas Version 1 42 2009 03 22 13 04 08 936 NMEA Nav communication established 2009 03 22 13 04 10 154 Ocean Surveyor Broadband Narrowband ADCP Teledyne RD Instruments c 1997 2007 All rights reserved Firmware Version 23 16 gt 2009 03 22 13 04 10 154 TS 090322 130410 2009 03 22 13 04 10 209
10. TS 090322 130410 gt 2009 03 22 13 04 10 227 gt 2009 03 22 13 04 10 227 CR1 2009 03 22 13 04 10 299 CRI Parameters set to FACTORY defaults gt 2009 03 22 13 04 10 299 TE00000000 2009 03 22 13 04 10 336 TE00000000 gt 2009 03 22 13 04 10 336 gt Setting Bottom Track parameters 2009 03 22 13 04 10 336 BPO 2009 03 22 13 04 10 372 BPO gt 2009 03 22 13 04 10 372 TP000000 2009 03 22 13 04 10 409 TP000000 gt 2009 03 22 13 04 1 091 WPO 2009 03 22 13 04 10 427 WPO oO D gt 2009 03 22 13 04 10 427 NPI 2009 03 22 13 04 10 463 NPI gt 2009 03 22 13 04 10 463 gt Determining Surveyor mode 2009 03 22 13 04 10 463 WP 2009 03 22 13 04 10 554 WP WP 000 Number of Pings gt 2009 03 22 13 04 10 554 NP 2009 03 22 13 04 10 645 NP NP 001 Number of Pings gt 2009 03 22 13 04 10 645 gt Setting profiling parameters 2009 03 22 13 04 10 645 NN80 2009 03 22 13 04 10 681 NN80 gt 2009 03 22 13 04 10 681 NS400 2009 03 22 13 04 10 699 NS400 gt 2009 03 22 13 04 10 718 NF400 2009 03 22 13 04 10 736 NF400 gt 2009 03 22 13 04 10 409 gt Setting Surveyor low resolution mode 2009 03 22 13 04 10 736 ED110 2009 03 22 13 04 10 772 ED110 gt 2009 03 22 13 04 10 772 gt Setting manual ping mode 2009 03 22 13 04 10 772 CF
11. c_id l nc allocate h nc n m do ic 1 nc read l h ic enddo allocate uv 2 n m u nc n m v nc n m 41 open unit 1 file u_tpx06 2 form unformatted status old convert big_endian ede skip header read 1 do ic 1 nc read 1 uv u ic uv 1 v ic uv 2 enddo you now have elevations m s in array h and transports in arrays u and v for nc constituents on a nx m grid more precisely gt u 1 gives the complex amplitude of zonal transport m 2 s to the East call this U gt u 2 gives meridional transport m 2 s to the North call this V KKAKKAKKAAKKAAKAAKAAKAAKAAKAAKAAKAKAAKAKAAKAAKAAKAAKAAKAAKAAKAAKAAKAAAKAKAKAKAKAKAKAAKAKAAKAAAAKAA 3 GRID There are n divisions in longitude m in latitude The elevations and transports are given on a C grid so the locations at which U and V and h are given are not exactly the same The layout of the grid looks like this U 2 1 h 2 1 U 2 2 h 2 2 v 2 1 V 2 2 U 1 1 h 1 1 U 2 1 h 2 1 lt theta_lim 1 van V 2 1 lt theta_lim 1 dy 2 gt phi_lim 1 phi_lim 1 dx 2 The lats and lons in degrees for individual grid nodes are For h nodes lat i j theta_lim 1 j 5 dy lon i j phi _lim 1 i 5 dx For U nodes lat i j theta_lim 1 j 5 dy lon i j phi_lim 1 i 1 dx For V nodes lat i j theta_lim 1
12. des sorties du mod le rg polar tide_models Model_TPXO62_plot_l Model_TPXO6 2 90 60 304 serioso en EN o o 3 0 Latitude 30 180 Longitude 21 La r solution du mod le global est de 1 4 Les donn es assimil es proviennent surtout de l al tim trie radar de T P entre 66 N et 66 S Le mod le inclut aussi des patches pour plusieurs zones c ti res en particulier pour l Antarctique CADA00 10 Circum Antarctic inverse model et Ross_Inv_2002 Ross Sea height based inverse model Le mode d emploi est d taill en Annexe V Pour valuer l effet des mar es par rapport aux mesures ADCP nous avons fourni au mod le deux fichiers de positions horaires en route et en station La figure ci dessous illustre cet effet tr s faible en mer ouverte mais qui se renforce l approche des terres et particuli rement le long de la Terre de Feu dans le d troit de Lemaire et King George Island Drake 09 Maree en Route Drake 09 Maree en Stations umax 14 02 vmax 31 88 umean 0 08 vmean 0 22 mas PAS ak 82 47 Umeak DISS mean 6 46 s NU m N bonus so ji p 1 Z ei H g aH Li SE F P 2 4A Zl 60 A 4 7 N Tw dE on su don aw M La no e ei a Fic 6 5 Tests de l influence de la mar e en route gauche et en station droite En station la vitesse moyenn
13. est pas le cas ici On a aussi souvent utilis dans d autres campagnes la proc dure de correction Proc gyro de E Firing qui cr e un mod le des erreurs de gyrocompas plus volu que la formule brute ci dessus utilisant deux fichiers cr s dans le logiciel CODAS soit le fichier de navigation liss suffixe sm et le fichier des caps enregistr s dans la base ADCP suffixe hdg Le mode d emploi de cette proc dure est donn en Annexe IV Le fichier d angles cr par cette proc dure est liss et invers pour tre utilis dans les programmes de rotation de CODAS au moyen d une adaptation du logiciel Corryge de G Eldin Mais ici le cap provenant d une centrale inertielle gyro laser acc l rom tres etc ces proc dures ne semblent plus applicables Nous l avaons toutefois test e et les r sultats sont compl tement simi laires ceux de la m thode suivante du fit en latitude Par contre comme dans le cas des campagnes TIP2000 et INTERPOLE 99 et les campagne DRAKE 2006 et 2008 le grand nombre de stations CTD effectu es permet une autre approche pour valuer les corrections angulaires dans la m thode du water tracking du logiciel CODAS plus d taill e au praragraphe 5 3 suivant le nombre de points extraits par le programme de d tection des acc l rations d c l rations timslip est th oriquement gal au minimum au double du nombre des stations et ou des ralentissements et ou des chang
14. fichiers de contours con associ s des fichiers de statistiques sta soit des fichiers de vecteurs diff rents niveaux vec Ces derniers peuvent tre trac s au moyen du programme vector de CODAS Les figures des pages 24 27 montrent les trac s des vecteurs de courant pour les 2 radiales sur les donn es apr s calibration compl te La construction de grilles temporelles se fait galement en deux temps d abord avec le programme timegrid dans lequel on entre un pas de temps r gulier si on veut des moyennes temporelles ou variable par exemple le d but et la fin des stations CTD ce qui permettra d obtenir des profils moyenn s sur la dur e de chaque station Si l on utilise une grille avec la dur e des ensembles comme pas de temps et un programme adcpsect ne faisant aucune moyenne sur les donn es on obtient tous les ensembles profils directement en ASCII la base CODAS est en binaire Un formattage post rieur permet de pr senter les profils sous forme de matrices facilement exploitables sous Matlab Les sections pr sent es pages 24 27 sont trac es sur ces donn es 6 6 Evaluation de l erreur Erreur sur la vitesse relative erreurs sur l amplitude La premiere source d erreur est comme nous l avons vu plus haut l erreur sur la vitesse du son Nous avons vu que l erreur sur la salinit dans les fichiers de configuration donc sur la vitesse du son entrainait une erreur de 2 au maximum e
15. in ODAS typically run in the cal rotate directory This creates n ascii file with ensemble time mean gyrocompass heading and last gyrocompass heading given the extension hdg Note that this file must be manually edited to remove the comments at the beginning before Matlab can read it These input files are first processed by the splitset m script to break the data at large gaps As a general rule the model output will not be sensitive to small changes in the input data Thus a complete ADCP calibration does not need to be done before generating the files for creating the model In most cases where gyrocompass errors are relatively small 2 3 degrees at most the smoothr and lst_hdg files can be generated before ADCP editing and calibration have been completed 37 Model Parameters At the top of the Proc gyro m script there is a section of model parameters which must be edited 1 Compensation terms Two arrays must be defined to specify the gyrocompass compensation used with a particular data set Gyrocompass compass compensation is set based upon latitude and north south velocity Latitude compensation is specified with an array called clat with the same size as the lat array loaded from the sm input file As with the input lat values clat should be in degrees Several possibl xamples are shown below a No latitude compe
16. none botrk hd misalign 0 deg botrk scale factor 1 none pit misalign 0 deg rol misalign 0 deg unused1 le 38 none last temp le 38 C last heading le 38 deg last pitch le 38 deg last roll le 38 deg mn pitch 0 1 deg mn roll 0 7 deg std temp 1e 38 C std heading le 38 deg std pitch le 38 deg std roll le 38 deg ocean depth 32767 m max amp bin 14 none last good bin 11 none unused2 32767 none unused3 32767 none unused4 32767 none TAB 3 4 Contenu de la variable 38 ANCILLARY 2 Chapitre 4 Obtention et v rification de la navigation Une commande du logiciel CODAS permet d obtenir partir de la base ADCP charg e pr c demment un fichier des caps moyens suivis pendant les ensembles ADCP suffixe hdg On utilise les deux premi res colonnes de ce fichier la date CODAS en premi re colonne et le cap moyen suivi pendant la dur e de l ensemble 1 minute en deuxi me colonne Ce fichier servira par la suite de r f rence pour obtenir les dates des ensembles converties en unit s sp cifiques au logiciel CODAS jours d cimaux plus nombre de secondes dans le jour La calibration effectu e ult rieurement par water tracking permettra de d terminer quel type de mesure du cap est transmise ADCP A ce stade du traitement on ne dispose que des deux fichiers positions gpsl en principe heure de d but des ensembles et gps2 heure de fin des ensembles Ces fichiers p
17. taient correctement transmises et n avons donc pas fait tous les contr les effectu s en 2006 Par contre tant donn la moins bonne qualit vidente des donn es pendant les transits en route comme en 2006 nous avons pris la d cision de privil gier les donn es acquises en station en extrayant un profil moyen pour chaque station moyenne de tous les profils valid s effectu s pendant la dur e de la station 6 2 Erreur sur le cap On a ensuite tent d valuer les erreurs sur ce cap gyro MINS en cr ant un fichier de corrections angulaires qui servirait la fin de la phase de calibration effectuer une rotation des valeurs de la base de donn es Pendant la phase de tests de calibration on effectue cette rotation programme rotnav sur le fichier des vitesses du navire par rapport la couche de r f rence fichier nav pour ne pas modifier la base elle m me 15 Dans le cas d un gyrocompas classique nous avions d j utilis une formule de correction en fonction du cap de la vitesse du navire et du cosinus de la latitude Angle de correction Vnavire Cos Cap Vangulaireterre x Cos Latitude 6 1 L inconv nient de cette formule est qu elle pr suppose que le gyrocompas n est pas du tout compens ce qui n est sans doute pas le cas Les r sultats donneraient donc une estimation maximum de l erreur sur le gyrocompas D autre part il faut avoir constamment la vitesse du bateau ce qui n
18. 0267 82136 Full No drv 10267 by 1 10267 82136 Full No gamma 1 by 1 1 8 Full No have_ash 1 by 1 1 8 Full No hd 10267 by 1 10267 82136 Full No lat 10267 by 1 10267 82136 Full No lon 10267 by 1 10267 82136 Full No phi 10267 by 1 10267 82136 Full No phi_ss 10267 by 1 10267 82136 Full No phi_ssf 10267 by 1 10267 82136 Full No rv 10267 by 1 10267 82136 Full No rcv 10267 by 1 10267 82136 Full No relat 10267 by 1 10267 82136 Full No rlat 10267 by 1 10267 82136 Full No time 10267 by 1 10267 82136 Full No ws 1 by 1 L 8 Full No The results are stored as shown below phi gyrocompass error based on differential model phi_ss gyrocompass error based on steady state model phi_ssf filtered steady state model output Input Execute Output 39 kn32pOl mat output of splitset m Matlab Proc_gyro which data file to load no extension kn32p01 Loading data Defining constants and system compensation Running dynamic modelling routine Setting initial conditions Solving array of differential equations Running steady state modelling routine Saving Results kn32p_gyro mat 40 ANNEXE V OSU TIDAL PREDICTION Software OTPS KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 1 INTRODUCTION OTPS accomplish 2 tasks extracting harmonic constants from barotropic tidal solutions in OTIS format see below at given locations predicting tides at given times and locations
19. 1 03 2009 09 04 2009 80 W 70 W 60 W 50 W 50 S F x 50 S 55 S 55 S 60 S 60 S 65 S 65 S 80 W 70 W 60 W 50 W ie ON ubany Dallmann pt Lah Punta Arenas Ba X RV Polarstern en 9x ANT XXV 4 ANTO Punta Arenas Punta Arenas LA FIL ited Werner Bain March 21st till April 9th 2009 Total distance covered 2433nm DIISLI Meemerhaven Chapitre 2 Acquisition des donn es La navigation est aquise par le biais d une centrale inertielle MINS Marine Inertial Navigation System du constructeur Raytheon dont les caract ristiques sont donn es en ANNEXE I Les donn es ADCP ont t acquises au moyen du logiciel RDI VMDAS version 1 42 L ADCP de coque du POLARSTERN Profileur acoustique effet Doppler est du type Ocean Surveyor Broadband Narrowband RDI 150 KHz Firmware version 23 14 Le d tail de la configuration est donn en ANNEXE II L appareil est configur en mode Narrow Band pour une meilleure p n tration port e maximum de 500m La configuration choisie est une longueur d ensemble de 1 minutes en fait 60 0003 s pour les fichiers Short Term Averaging et de 20 minutes pour les fichiers Long term Averaging respectivement LTA et STA pour le logiciel VMDAS La longueur de bin est de 4m la configuration nominale par d faut pour la fr quence de 150 KHz est de 8m L immersion r elle du transducteur est de 11 m La figure suivante donne le d t
20. 16 14 i 200 f b Wir 0 50 200 250 20 0 60 o so 100 o 20 20 40 mean SPECTRAL WOTH sid dev SPECTRAL WIDTH mean PERCENT 3 BEAM std dev PERCENT 3 BEAM Fic 6 4 Tests de qualit globale 6 5 Mise en forme d finitive de la base de donn es Tout d abord a t cr e une base des donn es dit es donc avec des profils propres mais avec des donn es de navigation brutes Puis nous avons effectu la calibration Donn es juste dit es suppression de tous les points ADCP aberrants Rotation par un fichier d angles calcul s par le fit en latitude sur les points issus du water tracking puis par une constante r siduelle correction de l alignement du transducteur pour les trois parties Le dernier stade de l laboration de la base est l int gration des donn es de navigation calibrees fichier sm dans la base programme putnav Il est alors possible de construire des grilles de diff rents formats suivant un maillage spatial ou temporel qui serviront de base pour la conversion des donn es en ASCII sous forme de fichiers formattables pour tre lus sous Matlab et pour les trac s 19 La construction des grilles spatiales se fait en deux temps d abord avec le programme Ilgrid qui demande la d finition de la grille en latitude et longitude puis le r sultat un tableau de temps sert de contr le pour le programme adcpsect soit pour cr er les
21. 61 60 59 58 latitude FIG 7 2 Section en latitude pour la radiale DRAKE 2009 Nord Sud 25 80 60 40 20 100 80 60 40 20 40 80 100 26 DRAKE 2009 3 au 8 avril 2009 Editees fit en latitude rot cte 52 S Layer 30a50m 54 S 55 S 57 S 58 s Lg 4 60 S 61 S 63 S l 70 W 65 W 60 W 55 W 0 100 Speed cm s FIG 7 3 Vecteurs de courant dans la couche 30 50 m pour la radiale DRAKE 2009 Sud Nord DRAKEO9 Sud Nord Vitesse zonale U cm s 50 100 150 r Profondeur m 200 250 61 60 59 58 57 56 55 latitude DRAKEO9 Sud Nord Vitesse meridienne V cm s 1007 Profondeur m 150 r 200 250 61 60 59 58 57 56 55 latitude Fic 7 4 Section en longitude pour la radiale DRAKE 2009 Sud Nord 27 80 60 40 20 100 80 60 40 20 40 80 100 28 29 ANNEXE I SYSTEME DE NAVIGATION MINS Marine Inertial Navigation System Raytheon Marine GmbH ehemals Ansch tz Zeyestra e16 24 24106 Kiel P O Box 1166 24106 Kiel Genauigkeit nach Spezifikation maximale Abweichung Heading 0 279 Roll 0 057 Stampf 0 057 im Betrieb 1 mil Rad im Bogenma bei 3 wird der Navigationsmode zugeschaltet Die Genauigkeit nach Pr fprotokoll liegt bei 0 025 und besser
22. OutputComEnable FALSE AsciiEnsembleOutputComPortName None AsciiEnsembleOutputComBaudRate 9600 AsciiEnsembleOutputComParity NOPARITY AsciiEnsembleOutputComStopBits 1 AsciiEnsembleOutputComDataBits 8 AsciiEnsembleOutputDataType 0 none l enr 2 enx 3 sta 4 lta 0 AsciiEnsembleOutputRefVelType O none 1 Bottom 2 Mean 0 AsciiEnsembleOutputStartBin 1 AsciiEnsembleOutputEndBin 4 AsciiEnsembleOutputStoreToDisk 0 no 1 yes FALSE AsciiEnsembleOutMeanStartBin 1 AsciiEnsembleOutputMeanEndBin 4 AsciiEnsembleOutputLeader 0 no 1 yes TRUE AsciiEnsembl AsciiEnsemb AsciiEnsemb AsciiEnsemb AsciiEnsemb AsciiEnsemb AsciiEnsemb SpeedLogCom Speed Log Speed Log Speed Log Speed Log Speed Log SpeedLogW Log Output ComPortName None ComBaudRate 9600 ComParity NOP ComStopBits 1 ComDataBits 8 SpeedLogDataSource STA SpeedLogW SpeedLogW Source WP StartBin 3 EndBin 5 IP Port BinaryEnsemb for BinaryEnsemb IP Port Asciil AsciiEnsemb Ensembl for TP Port SpeedLogNet T Binary eOutputNetEnable FALSE PortNumber 5433 ASCII Ensemble Data Output leOutputNetE eOutputI for Speed Log Enable FALSE eOutputIP ARITY T Ensemble Data Output 2 L nable FALSE Output SpeedLogHostName IPAddress 5434 Fake Data Options Ad
23. RAW DOPPLER none 3 2 18 18 65 RAW AMP none 1 2 19 19 66 RAW SPECTRAL WIDTH none 11 2 20 20 67 BEAM STATS none 1 2 21 21 68 NAVIGATION none 11 2 22 22 69 BOTTOM TRACK none 11 2 23 23 75 USER BUFFER none 10 2 24 24 76 ADCP CTD none 11 2 25 25 TAB 3 1 Liste des variables de la base adcpdb La commande showdb permet aussi de connaitre le contenu d une variable On a pu ainsi v rifier que la dur e d ensemble tait de 192 997 s pour le premier bloc de donn es et de 60 0003 s ensuite Le premier bloc ne sera donc pas conserv name value unit avg interval 60 0003 s compensation 1 none num bins 80 none tr depth 11 m bin length 4 m pls length 4 m blank length 4 m ping interval 0 s bot track 0 none pgs ensemble 1 none ens threshold 0 none ev threshold 0 mm s hd oflset 45 deg pit offset 0 deg rol offset 0 deg unusedl 0 none unused2 0 none unused3 0 none freq transmit 150000 Hz top ref bin 0 none bot ref bin 0 none unused4 0 none heading bias 0 deg TAB 3 2 Contenu de la variable 35 CONFIGURATION 1 name value unit tr temp 10 26 C snd spd used 1490 m s best snd spd le 38 m s mn heading 18 77 deg pgs sample 19 none unassignedl 32767 none unassigned2 32767 none unassigned3 32767 none unassigned4 32767 none unassigned5 32767 none TAB 3 3 Contenu de la variable 37 ANCILLARY 1 name value unit watrk hd misalign 0 deg watrk scale factor 1
24. U ADCP de Coque Campagne DRAKE A bord du POLARSTERN Mars Avril 2009 Rapport Interne LOCEAN Mai 2009 Laboratoire D Oc anographie et du Climat Exp rimentation et Approche Num rique UNITE MIXTE DE RECHERCHE 7159 CNRS IRD UNIVERSITE PIERRE amp MARIE CURIE INSTITUT PIERRE SIMON LAPLACE AWI Chapitre 1 Chapitre 2 Chapitre 3 Chapitre 4 Chapitre 5 Chapitre 6 Chapitre 7 ANNEXE I NNEXE II NNEXE III NNEXE IV A A A A NNEXE V SOMMAIRE La campagne Acquisition des donn es Chargement des fichiers bruts Obtention et v rification de la navigation Edition des donn es ADCP Calibration des donn es Section 1 Erreur sur la vitesse de r f rence Section 2 Erreur sur le cap Section 3 Alignement du transducteur Section 4 R sultats de la calibration Section 5 Mise en forme d finitive de la base de donn es Section 6 Evaluation de l erreur Section 7 Obtention des profils moyenn s en station Section 8 Evaluation de l influence de la mar e Section 9 References Section 10 Remerciements Pr sentation rapide des donn es Syst me de navigation Configuration de l ADCP Fichier de contr le du logiciel VMDAS Proc dure de calibration du gyro compas E Firing Mod le de Mar e OTPS 11 13 13 14 16 17 18 19 20 20 21 22 23 29 30 34 36 40 Chapitre 1 La campagne La figure montre le trajet effectu ANT XXV 4 2
25. ail de la fen tre o est install PADCP Elle est remplie d eau de mer ADCP Transducer I ri TN E ET Seawater i N 4 V Ai N okustisches Fenster z Acoustic window Polycarbonate material Thickness 79 mm Fic 2 1 Fen tre ADCP Le transducteur est situ 20 m l avant et 3 m babord de la centrale inertielle Un programme de la routine CODAS de traitement des donn es corags permet d valuer et de corriger l influence de la distance du transducteur au systme d acquisition de la navigation La figure ci dessous illustre ces positions Port Starboard F1G 2 2 Positions de ADCP et de la centrale inertielle Le traitement proprement dit se d roule en 5 phases d crites ci dessous auxquelles on rajoute l valuation des erreurs Chargement des fichiers bruts Cette phase comporte aussi une premi re limination des donn es non enregistrables erreurs de chronologie mauvaises dur es d ensembles etc Obtention et V rification de la navigation Il est pr f rable d avoir une acquisition des donn es de navigation ind pendante de l acquisition ADCP Edition des donn es ADCP Cette phase tr s importante permet d liminer toutes les mesures mauvaises ou douteuses et de rep rer le fond si n cessaire Calibration des donn es En deux temps au moyen des donn es de navigation correction des donn es de cap puis des d fauts d
26. aling FALSE VelocityScaleFactorForBTVelocities unitless 1 VelocityScaleFactorForProfileAndWTVelocities unitless 1 EnableTiltAlignmentErrorCorrection TRUE TiltAlignmentHeadingCorr deg 0 TiltAlignmentPitchCorr deg 0 TiltAlignmentRollCorr deg 0 Averaging Options AvgMethod 0 time l dist 0 FirstAvgTime 60 SecondAvgTime 1200 FirstAvgDistance 10 SecondAvgDistance 1000 EnableRefLayerAvg TRUE RefLayerStartBin 3 RefLayerEndBin 10 Reference Velocity Options RefVelSelect 0 none 1 BT 2 WT 3 LYR 4 NDP 5 NAP 6 NSPD 1 VelRefLayerStartBin 3 VelRefLayerEndBin 5 RefVelUnitVel O mm s 1 m s 2 knots 3 ft s 1 RefVelUnitDepth 0 m 1 cm 2 ft 0 User Exit Options Qu Fh t HA HM UserWinAdcpEnable FALSE UserWinAdcpUpdateInterval sec 10 UserWinAdcpFilel UserAdcpScreening FALSE ADCP Setup Options etProfileParameters TRU NumberOfBins 60 BinSize meters 8 BlankDistance meters 8 TransducerDepth meters 11 SetBTEnable 0 SendBPCmd 1 Don tSendBPCmd TRU ADCPSetupMethod 0 Options 1 CommandFile 0 BtmTrkEnable 0 SendBPO 1 SendBP1 1 MaxRange meters 500 SetHdgSensorType TRUE HdgSensorType 0 internal l external 1
27. articuli rement l ing nieur lectronicien bord du POLARSTERN sans qui l acquisition de ces donn es n aurait pas pu se faire Chapitre 7 Pr sentation rapide des donn es Nous pr sentons les trac s des vecteurs de courant moyenn s sur une grille spatiale de 1 20 me de degr pour la premi re couche soit de 30 50 m pour chacune des radiales donn es en route et en station m lang es Les trac s pr sentent les donn es avec une correction du cap par le fit en latitude expliqu page 15 suivie d une correction par une constante Sur ces donn es on pr sente galement les sections en fonction de la latitude pour les composantes U et V de la vitesse de l ensemble des profils pour les radiales Nord Sud et Sud Nord en travers du Passage de Drake 23 24 DRAKE 2009 22 mars 3 avril 2009 Editees fit en latitude rot cte 52 S Layer 30a50m 53 S 54 S 55 S 56 S 578 58 S 59 S 60 S 62 S 63 S l 70 W 65 W 60 W 55 W 0 100 Speed cm s FIG 7 1 Vecteurs de courant dans la couche 30 50 m pour la radiale DRAKE 2009 Nord Sud Profondeur m Profondeur m DRAKEO9 Nord Sud Vitesse zonale U cm s 50 100 150 F 200 F 250 61 60 59 58 latitude DRAKEO9 Nord Sud Vitesse meridienne V cm s 50 1007 150 r 200 250
28. atitude avant et apr s calibration 6 3 Alignement du transducteur On peut utiliser deux m thodes comparer le d placement du bateau par rapport au fond donn par la navigation positions GPS au d placement par rapport au fond d duit de la d tection du fond par VADCP bottom tracking ou bien comparer l acc l ration par rapport l eau mesur e au moyen de ADCP avec l acc l ration calcul e au moyen des points de positions water tracking Ces deux m thodes donnent g n ralement des r sultats identiques 0 5 pr s pour la determination de orientation du transducteur selon l axe longitudinal de la coque Dans les deux cas il faut galement examiner la qualit des positions satellites validation du fichier navigation et des donn es de cap Les enregistrements de bottom tracking n ont pas t acquis Pour la calibration par water tracking que nous avons utilis e pour corriger les d fauts d ali gnement du transducteur on utilise le fichier nav extrait de la base ADCP dit e et sur lequel on a effectu la rotation des angles de correction du cap On d tecte ensuite les acc l rations du navire partir des vitesses ADCP puis on les compare aux acc l rations issues des donn es GPS en faisant glisser les deux fichiers vitesses de la couche de r f rence et vitesses GPS de mani re faire correspondre les temps le mieux possible programme timslip Les acc l rati
29. au sal e 35 PSU eau douce 0 PSU sur la vitesse relative Par exemple pour une vitesse de bateau de 10 noeuds soit environ 5 m s l erreur est alors de 10 cm s pour la composante longitudinale dans l axe du bateau L incertitude sur la temp rature du transducteur entraine ici une erreur de 0 2 qu on consid rera comme n gligeable Les bulles peuvent aussi r duire la vitesse du son Nous n avons pas trop d informations sur l in fluence de la forme de la coque du POLARSTERN brise glace coque ronde mais vu les fichiers de configuration utilis s Narrow Band bins de Am nous pensons que les bulles doivent tre bien pr sentes D autre part les eaux polaires sont souvent pauvres en particules ce qui explique aussi cette configuration particuli re meilleure p n tration Erreur sur la vitesse absolue erreurs sur la phase En supposant que la vitesse du bateau est bien connue pr cision des points GPS les erreurs sur la composante transversale proviennent des erreurs sur la mesure du cap et sur le positionnement horizontal et vertical du transducteur La m thode de la couche de r f rence permet de corriger les erreurs ventuelles sur les points GPS Ces erreurs taient fr quentes avec l ancien syst me Transit tr s rares avec le syst me GPS m me d grad La m thode est toutefois conserv e pour ses outils de visualisation et de contr le Elle a t fort utile dans le cas de DRAKE La distan
30. avire mesures ADCP par rapport une couche de r f rence programme adcpsect fichier cr nomm nav qui est pour chaque ensemble la moyenne des vitesses entre les bins 10 et 25 soit 40m et 100m La couche de r f rence id ale serait une couche de vitesse constante alors que les vitesses dans l oc an r el varient Il faut donc choisir une couche de r f rence suffisamment fine pour que les varia tions soient le plus lentes possible mais assez paisse pour contenir le maximum de donn es le trac du PGOOD ou pourcentage de bonnes donn es donne une indication utile On rappelle que la vitesse du navire par rapport la couche de r f rence est l inverse de la vitesse ADCP par rapport la couche de r f rence Les donn es de position permettent d tablir la vitesse du navire dans le r f rentiel terrestre Si on admet que l dition des donn es ADCP est correcte suppression de tous les bins et profils mauvais ou simplement douteux lorsqu on soustrait de la vitesse du navire par rapport la couche de r f rence la vitesse absolue du navire programme refabs fichier cr ref on obtient alors les vitesses absolues de la couche de r f rence et les anomalies proviennent d erreurs dans le fichier des positions ou d erreurs d horloge vitesse distance temps Historiquement l poque du systeme Transit les erreurs de positionnement distance taient pr pond rantes Actuellement avec la qual
31. ce entre l antenne et le transducteur ADCP peut galement influer Cette derni re erreur a t corrig e au moyen de la proc dure corags Mais ces erreurs sont tr s faibles en regard des suivantes L erreur sur le positionnement horizontal angle par rapport Vaxe longitudinal du navire peut tre valu e et corrig e au moyen de la m thode du water tracking Pollard et Read condition d avoir un nombre suffisant de points de calibration id alement au moins 64 pour avoir une pr cision de 0 2 sur la calibration L erreur sur le positionnement vertical roulis tangage assiette est normalement corrig e en configuration Janus les axes du transducteur sont en X par rapport l axe du navire 20 L erreur sur le cap a t valu e approximativement au moyen du fit en latitude Une erreur de 0 6 entraine une erreur de 1 de la vitesse du bateau sur la composante transverse soit 5 cm s pour une vitesse de 10 noeuds En r sum compte tenu des difficult s d appr ciation des erreurs sur le cap issu de la centrale inertielle une incertitude de 10cm s sur le module de la vitesse parait une estimation raisonnable Toutefois les directions fortiori dans ces r gions de forts courants sont r alistes Comme il l est rappel dans la plupart des publications techniques sur les mesures ADCP de coque ces mesures ne semblent pas utilisables pour le calcul de transports 6 7 Obtention des profils moyenn s e
32. cpFakeData AdcpFakeDa FakeDataTimeB Enable FALSE a L taFilename Si mAdcp enr tweenEns D EnableDual FileRecordPath C FileRecordBackupPath D NMEAFakeDatal NMEAFakeDataFilename SimNav nmr File Name Enable FALSE mbles 2 L Components RecordDir FALSE ADCPData ARK21_1b BI BI Dep ErrScreen VertScree FishScree TAmplitude Correlati WI WI ErrorVelT FishThres rAmpScreen CorScreen ErrScreen VertScree FishScree TAmplitude rAmpScreen CorScreen Enable FALSE Enable FALSE L Enable FALSE ploymentName ARK21_1b loymentNumber 1 MaximumFileSize 10 Bottom Track Data Screening Options nable FALSE L a L L L n n T o R Enable FALSE ay L L hreshold VerticalVelThreshold 1000 hold 50 PctGoodThreshold 25 Water Track Data Screening Options IPctGoodScreenEnable T RUE hreshold 30 nable FALS nable FALS ET ET EE ES PI El nable FALS n n T nThreshold 220 1000 Enable FALSI Enable FALSI CJ DI IPctGoodScreenEnable T hreshold 30 PPortNumber 5433 eOutputBottomTrack 0 no l yes eOutputNavigation 0 no 1 yes T eOutputVelocity 0 no l yes TRU eOutputIntensity 0 no l yes TR leOutputCorr
33. del of gyrocompass errors based upon input data from an ADCP CODAS database The model is based upon theory from Gyroscopic Theory Design and Instrumentation by Wrigley Hollister and Denhard Specifically this model is for a Sperry gyrocompass as described in Chapter 10 Equations 10 36 to 10 39 This script employs a 3rd order Runge Kutta solution of the second order differential equation modelling gyrocompass error eq 10 38b Note this can be changed to a fifth order Runge Kutta solution by switching the call in Gyro_dyn m from RK3 m to RK5 m However given noise levels in input data the higher order solution does not appear significantly better In addition to the solution to the differential model the script creates a steady state solution which makes no attempt to track compass response time or oscillations A filtered version of the steady state response is also generated which typically shows some improvement over the steady state solution Input Data Two input files are needed for the time period to be modelled The first is output by the smoothr routine in CODAS navigation processing run in the nav directory as an ascii flat file with the extension sm This file contains ensemble times and the best estimates of ship s u ship s v reference layer u r T a ference layer v latitude and longitude data he second output file is generated by the lst_hdg routine
34. e de courant de mar e est de 0 23 cm s avec un maximum de 34 8 cm s En route la vitesse moyenne est de 1 38 cm s avec un maximum de 111 3 cm s 6 9 R f rences Firing E Report from the WOCE NOAA Workshop on ADCP measurements held in Austin Texas March 1 2 1988 U S WOCE Planning Report No 13 U S Planning Office for WOCE College Station TX Joyce T M On in situ calibration of shipboard ADCPs J Atmos Oceanic Technol 6 169 172 1989 Pollard R and Read J A method for calibrating shipmounted acoustic Doppler profilers and the limitations of gyro compasses J Atmos Oceanic Technol 6 859 865 1989 RD Instruments Acoustic Doppler Current Profilers Principles of Operation A Practical Primer 1989 RD Instruments Businesspark AVenue San Diego CA 92131 Fischer J and Visbeck M Deep velocity profiling with self contained ADCPs J Atmos Ocea nic Technol 10 5 764 773 1993 22 Firing E and Bahr F Caldwell P Ranada J Zhu W Processing ADCP Data with the CO DAS software System Version 3 1 unpublished manuscript 1995 JIMAR University of Ha waii 1000 Pope Road Honolulu Hawaii 96822 Eldin G Correction des erreurs gyro compas par GPS System version 1 2 unpublished manus cript ao t 2000 Egbert G D and Erofeeva S Y Efficient inverse modeling of barotropic ocean tides J Atmos Oceanic Technol 19 2 183 204 2002 6 10 Remerciements Nous remercions tout p
35. el_ into subdirectory DATA of the OTPS directory KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK Any questions contact Gary Egbert 541 737 2947 egbert coas oregonstate edu or Lana Erofeeva 541 737 3964 serofeev coas oregonstate edu
36. elation 0 no l yes T leOutputPercentGood 0 no l yes leOutputStatus 0 no 1 yes TRU Serial Port for Speed Enable FALSE R E FG zz N N O PEP 31 32 WTCorrelationThreshold 180 WTErrorVelThreshold 1000 WTVerticalVelThreshold 1000 WTFishThreshold 50 WTPctGoodThreshold 25 Profile Data Screening Options RAmpScreenEnable FALS RCorScreenEnable FALS El bi T L ErrScreenEnable FALSE VertScreenEnable FALSE FishScreenEnable FALSE PctGoodScreenEnable TRUE MarkBadBelowBottom FALSE R R R R R RAmplitudeThreshold 30 R R R R R CorrelationThreshold 180 ErrorVelThreshold 1000 VerticalVelThreshold 1000 FishThreshold 50 PctGoodThreshold 25 Transformation Options formToEarth TRUE llow3Beam TRUE nMap TRUE eamAngleSrc 0 auto l man 0 nualBeamAngle 30 adingSource 0 adcp 1 navHDT 2 navHDG 3 navPRDID 4 manual 3 nualHeading 0 S a Ei zwu wr X LU youne H a e a iltSource 0 adcp 1 nav 2 man 1 anualPitch 0 a e o z z nualRoll 0 nsorConfigSrc 0 PRfixed 1 Pfixed 2 auto ncavitySource 0 convex l concave 2 auto UpDownSource 0 dn 1 up 2 auto 2 EnableHeadingCorrections TRUE SinCorrectionAmplitudeCoefficient 0 SinCorrectionPhaseCoefficient 0 MagneticOffsetEV 0 AlignmentOffsetEA 45 EnableVelocitySc
37. ements de direction Ce programme timslip permet ainsi d valuer l erreur sur le cap dans le cas o l on dispose de nombreux points de calibration campagnes CTD avec nombreuses stations et peu de transits cap constant en valuant l influence de la latitude sur l angle entre l axe du transducteur et l axe longitudinal de la coque misalignment d termin par le programme timslip phase Rappelons que cet angle devrait id alement tre une constante On trace la phase extraite au moyen du programme timslip d crit ci dessus en fonction de la latitude et on effectue un fit polynomial sur le nuage de points On privil gie la variation en fonction de la latitude et non en fonction de la vitesse du navire voir formule 5 1 car l exp rience a montr que c est le facteur le plus important Si l erreur sur la mesure du cap tait nulle la phase resterait constante quelle que soit la latitude Dans le cas contraire le polyn me obtenu permet d valuer la variation du gyrocompas en fonction de la latitude en calculant un fichier de diff rences d angles variant avec la latitude partir du fichier de navigation Ce fichier d angles est construit en prenant comme constante du polyn me la diff rence entre la constante calcul e par le fit et la phase constante angle du transducteur obtenue au moyen de la calibration par water tracking adcpcal Le terme du premier degr est celui du polyn me
38. es croix et les cercles marquent tous les bins rejet s D une fa on g n rale les profils sont mauvais en route libre et meilleurs pendant les stations CTD Les bins marqu s comme mauvais les profils supprim s les bins correspondant au fond sont en registr s respectivement dans 3 fichiers ASCII badbin badprf bottom Ces fichiers seront consult s ensuite pour la validation de la base de donn es Chapitre 6 Calibration des donn es Le principe g n ral de la calibration est de passer de vitesses de courant relatives au bateau des vitesses absolues dans un r f rentiel terrestre Pour obtenir ces vitesses de courant absolues il faut soustraire aux mesures relatives la vitesse du navire par rapport au r f rentiel terrestre cette vitesse tant donn e par la navigation Les erreurs proviennent de trois sources positions horloges Mauvaises valuations de la vitesse cap Variation du gyrocompas ou du cap donn par toute autre source transducteur D faut ventuel d alignement du transducteur ADCP par rapport l axe du bateau Ce dernier est d termin par deux m thodes bottom tracking et ou water tracking 6 1 Erreur sur la vitesse de r f rence Apr s les premi res v rifications effectu es lors du traitement des fichiers bruts de navigation le contr le final de la qualit n cessite l extraction partir de la base ADCP valid e chapitre dition d un fichier des vitesses du n
39. he model m2 s2 n2 k2 kl1 ol pl ql 7 oce geo extract ocean geocentric HC for elevations only Geocentric tide is appropriate for comparison to or correction of altimetry data May leave blank if not z on line 3 KKAKKAKKAAKKAAKAAKAAKAAKAAKAKAKAAKAAKAKKAAKAKAAKAAAKAAKAAKAAKAAKAKAKAAAKAKAKAKAKAAKAKAAKAKAAKAKAAAKA 5 COMPILING and RUNNING OTPS To compile make extract_HC make predict_tide If you are NOT on a Unix system you have to provide options to convert from big_endian binary when compiling convert big_endian assume byterecl Other possibility is to insert the option convert big_endian in the open statement in the source codes Your specific OS Fortran compiler might not have these options In this case please consult your local computer support on how to read big_endian binary files on your system or contact us on your particular purpose Run AFTER editing setup inp extract_HC lt setup inp have lat lon file ready predict_tide lt setup inp have lat lon time file ready KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KK 6 ADDING new model Download compressed and tarred model from http www coas oregonstate edu research po research tide a into the OTPS directory uncompress and tar xvf DONE b into any destination uncompress and tar xvf cd DATA edit Model_ control file for exact path to the model files cp file Mod
40. hld 28209 d2w_threshld 368 d2uv_threshld 718 Ces valeurs sont tr s lev es on peut s attendre des mesures de mauvaise qualit cause des mauvaises conditions de mer rencontr es La visualisation des profils se fait aussi au moyen de Matlab Les diff rentes erreurs r pertori es sont marqu es par des signes de couleurs et de formes diff rentes L op rateur peut accepter les sup pressions de donn es automatiques en supprimer plus ou bien d cider d annuler le marquage de valeurs erronn es Cette derni re possibilit est tr s utile dans le cas de couches de diffusion plancton qui font apparaitre des pics d amplitude du signal confondus dans le traitement automatique avec les pics d amplitude produits par la r flexion sur le fond Toutefois la forme et la hauteur de ces pics ainsi que la position g ographique permettent de lever l ambiguit La figure ci dessous montre l aspect des profils d amplitude pour DRAKE 2009 la sortie du d troit de Le Maire afin de visualiser l aspect du fond C ADRAKE DRAKE_09 DATA adcp_coq adcpdb_noedit drake AMPLITUDE offset 10 T SOD 20 F Bin Number A o T 80 2 1 1 0 100 200 500 600 Lon 65 0941 to 65 1201 Time 82 7949 to 82 8296 Lat 54 8762 to 55 0025 B P 7 150 to 7 200 2009 5 9 14 43 Fic 5 2 Profils d amplitude du signal ADCP Les pics de la figure decroissance entre les bins 25 et 75 correspondent au fond L
41. imal day time base is used R earth s radius is meters 6 37e6 meters 38 ws natural frequency of system equal to 2 pi 84 4 1 min if the compass is assumed to be Schuler tuned Schuler tuning is used to minimize errors and can be assumed to be present on all gyrocompasses Set to 2 pi 84 4 1440 1 day to match decimal day time base Running the model Once the beginning section of the Proc_gyro m script has been edited the model is run from within Matlab by simply typing Proc_gyro The routine prompts for the input file name containing the lat lon U V and hd heading data Thus the script can be run on multiple input files using the same compensation without editing the script each time Note that the model takes some time to run To process one days worth of data at 5 minute ensembles about 280 samples requires about 20 25 minutes on a Sun workstation Output File The model is saved in a Matlab format file using the name base from the input file and the addition _gyro mat The output file contains the following variables Name Size Elements Bytes Density Complex R 1 by 1 1 8 Full No U 10267 by 1 10267 82136 Full No V 10267 by 1 10267 82136 Full No Wie 1 by 1 1 8 Full No cV 10267 by 1 10267 82136 Full No clat 10267 by 1 10267 82136 Full No ddrv 10267 by 1 10267 82136 Full No dphi 10267 by 1 1
42. it des points GPS en date et position la seule source d erreur importante ne peut provenir que de la date ADCP Ce contr le permet de retrouver les erreurs d horloge du syst me d acguisition de ADCP 13 14 Pour un contr le graphique au moyen d un programme Matlab callrefp m les donn es sont liss es programme smoothr utilisant une fen tre de Blackman sortie sm Un exemple est donn ci dessous Absolute Refers prada soy IY _ corags2 ref 10 a 25 92 92 2 92 4 92 6 92 8 93 93 2 93 4 93 6 93 8 94 2009 Days Longitude 60 5 Latitude 2 61 5 92 92 2 92 4 926 92 8 93 93 2 93 4 93 6 93 8 94 2009 Days 2009 5 9 15 0 FIG 6 1 Exemple de trac de contr le de la couche de r f rence Alors que les trac s taient catastrophiques lors de la campagne de 2006 ici malgr un bruit certain et des valeurs parfois lev es de la vitesse de la couche de r f rence l aspect des donn es est globalement raisonnable Toutefois une variation brutale comme entre les jours 92 9 et 93 n est explicable que par un saut d horloge al atoire et reste tr s difficile corriger Contrairement 2006 nous avons choisi une couche de r f rence plus paisse et plus profonde De plus lors de l dition le bin de surface tait syst matiquement marqu comme mauvais Il a t pratiquement limin partout Nous avons consid r que les donn es de cap
43. itions des ensembles ADCP Fichier gps1 corrig Chapitre 5 Edition des donn es ADCP Une premi re v rification est celle de la vitesse du son qui peut tre fig e une vitesse constante ou varier suivant la temp rature du transducteur la salinit a une valeur constante dans le fichier de configuration En effet une erreur sur la vitesse du son entraine une erreur sur l amplitude du courant mesur par l ADCP Rappelons que la constante de proportionnalit entre le d calage Doppler et la vitesse de d placement de la masse d eau implique la vitesse du son au niveau du transducteur seulement On s assure tout d abord que la temp rature du transducteur est correcte en la comparant la temp rature du thermosalinographe et que la vitesse du son utilis e par ADCP est galement correcte L cart de 1 degr C observ ici entrainerait une difference de 0 2 sur la vitesse du son ce qui se v rifie sur le trac ci dessous Un programme Matlab permet de tracer les deux temp ratures ainsi que la vitesse du son La figure montre les trac s de la temp rature du transducteur rouge du thermosalinographe magenta de la vitesse du son utilis e par ADCP bleu recalcul e avec une salinit de 35 PSU une immersion de 11m et la temp rature du transducteur cyan et de la vitesse du son calcul e avec la salinit et la temp rature du thermosalinographe vert L immersion de 11m est la position r elle du tran
44. j 1 dy lon i j phi_lim 1 i 5 dx Here dx dy is resolution in longitude and latitude correspondingly Most of OTIS solutions are given on C grids uniform in lats lons but some of them i e Arctic solution are given on C grids uniform in kilometers The OTPS is applicable for such solutions also 42 rr 4 Thi KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK SETUP file setup inp s is EXAMPLE of the file containing your input information The lat Z m2 AP oce I tmp 1 0 correct for minor constituents output file ASCII original version of file is included in OTPS and given below 1 tidal model control file _lon_time 2 latitude longitude time file 3 z U V u v s2 n2 k2 kl ol pl ql 4 tidal constituents to include 5 AP RI 6 oce geo LT LL Comments on lines 1 7 di Tidal model control file ASCII supplied with the model starting from Model_ contains AT LEAST 3 lines elevation model file name transport model file name bathymetry grid file name name of function to convert x y to lon lat Leave line 1 blank if are not sure which model to use In this case all models from model list are checked File model list should contain all your DATA Model_ files and is done as ls 1 DATA Model_ gt model list Please have all Model files ready in subdirector
45. merait l existence d un probl me d horloge lors de l acquisition d calage entre l horloge de la centrale MINS et celle acquise par ADCP qui proviendrait du GPS Trimble 1 Mais nous n avons pas trouv de variation relativement r guli re Il y aurait donc des sauts Pour la phase l cart type admissible est de l ordre de 1 5 degr et pour l amplitude il est de 2 Les valeurs de la calibration sont donc raisonnables 1 2 degr et 3 Les valeurs de l amplitude brute tr s peu sup rieures 1 confirment que la fen tre de logement du transducteur est bien remplie d eau sal e mais que la temp rature r elle est sans doute inf rieure celle mesur e par le transducteur voir Fig 5 1 Le tableau suivant r capitule les dates de d but et de fin le nombre de points de calibration par water tracking qui permettent d obtenir la correction de l alignement du transducteur ainsi que les valeurs de la phase de l amplitude et de l erreur d horloge obtenues par cette calibration La derni re colonne indique le nombre de profils ADCP obtenus Dans la colonne Nb de points le premier chiffre est le nombre de points retenu le deuxi me chiffre le nombre de points de d part 18 Dans les colonnes Amplitude Phase et Delta T le premier nombre est la valeur avant la calibration le deuxi me nombre est la valeur apr s la calibration les nombres de la deuxi me ligne sont les carts type associ s
46. n station Un programme du logiciel CODAS adcpsect permet d obtenir les statistiques sur les profils moyenn s dans une gamme de temps voir la construction des grilles temporelles On a donc construit un fichier de gammes de temps avec les dates de d but et de fin de chaque station CTD pour les trois radiales Un exemple de ce fichier est donn ci dessous pour la radiale DRAKE 1 2009 03 22 12 16 00 to 2009 03 22 12 27 00 2009 03 22 13 16 00 to 2009 03 22 13 24 00 2009 03 31 21 33 00 to 2009 04 01 00 17 00 En sortie apr s reformattage les profils moyenn s en station se pr sentent sous la forme d un en t te pour chaque profil avec les dates de d but et de fin de chaque station la position latitude et longitude d cimales le nombre de profils moyenn s pendant la dur e de la station et dans le corps du profil la profondeur les composantes zonale et m ridienne de la vitesse du courant avec pour chacune l cart type sur le calcul de la moyenne le minimum et le maximum Ce fichier peut tre lu sous Matlab 6 8 Evaluation de l influence de la mar e Nous avons utilis le mod le TPXO de l Oregon State University TOPEX POseidon Global Inverse Solution dans sa version TPXO6 2 mod le global de mar es adapt aux traces de TO PEX Poseidon et Jason Le mod le n cessite le chargement du logiciel de Pr diction de Mar es de POSU OTPS pour OSU Tidal Prediction Software La figure donne une id e
47. ng e et horloge ADCP tait celle de la centrale MINS Un contr le rapide au moyen des fichiers de navigation bruts NIR a montr qu il n y avait pas de d calage syst matique entre les heures navigation trames GPGGA et ADCP trames PADCP 10 Heures des diff rents fichiers cap positions hdg 80 746319 2009 03 22 17 54 42 gps1 80 7456019 2009 03 22 17 53 40 ler enregistrement gps1 80 7463310 2009 03 22 17 54 43 2 me enregistrement gps2 80 7463194 2009 03 22 17 54 42 ler enregistrement Apr s correction au moyen de la proc dure corags pour tenir compte de la distance entre le transducteur ADCP et le syst me d acquisition de la navigation le fichier de type gps servira de base pour tous les calculs de calibration ult rieurs water tracking La figure suivante donne la correction effectu e en latitude et longitude par la proc dure corags sur le fichier gps1 30 so 20 E E z ok g k sg sg 30 30 30 9 95 so s 95 FIG 4 1 Correction en latitude et longitude due la distance entre l ADCP et le syst me de navigation La carte suivante montre le trac des donn es de navigation d apr s le fichier gps1 corrig au moyen de la proc dure corags DRAKE 2009 Route du navire Fichier gps1 S Boo 3009 VA 60 62 3000 70 W 65 W 60 W 55 W FIG 4 2 Pos
48. nsation clat zeros size lat b Constant latitude compensation clat ones size lat mean lat c Latitude compensation updated once per day tt fix time clat zeros size lat for ii min tt max tt k find tt ii clat k ones size k mean lat k end clear tt ii k North south velocity compensation is specified with an array called cV with the same size as the north south V array loaded from the input sm file As with input U and V values cV should be in m s Several examples are shown below a No velocity compensation cV zeros size V b Constant velocity compensation set at cruising speed cV 5 14 cos hd pi 180 10 knot hr 5 14 m s c Velocity compensation set high when underway zero otherwise cV zeros size V k find sqrt U 2 V 2 gt 2 underway sp gt 2 cV k 5 14 cos hd k pi 180 k find sqrt U 2 V 2 lt 2 on station sp lt 2 cV k zeros size k 2 System constants Several constants must be set which affect model behaivor Of the four constants shown gamma is the only one which would normally be modified The constants their meanings and typical values are shown below gamma eccentric offset a physical parameter of the compass itself expressed in radians Typical values 0 015 0 025 Wie earth s inertial angular velocity 7 29e 5 1 sec or 7 29e 5 86400 1 day since a dec
49. ons se produisent lors des arriv es et d parts de station CTD ou carottage lors des changements de direction et bien s r lors des ralentissements et acc l rations du navire sur sa route directe Ce programme effectue galement un filtrage des valeurs erron es On a choisi une fen tre de glissement de 9 fixes soit 9 minutes pour des ensembles de 1 mn La longueur de la fen tre doit tre sup rieure la dur e de l acc l ration ou de la d c l ration On d termine ainsi le d calage de temps entre l horloge du PC d acquisition des donn es ADCP soit le temps de la base ADCP et le temps GPS ainsi que l angle du transducteur par rapport l axe du navire phase et l influence de cet angle sur le module de la vitesse amplitude Le crit re de choix de la fen tre tait de garder le plus grand nombre de points de calibration avec des valeurs raisonnables d carts type On admet 0 6 1 2 degr s pour la phase et 1 1 5 pour Vamplitude La difference de temps nav pc doit rester de l ordre de quelques secondes Un programme Matlab adcpcal fait les moyennes et statistiques sur les valeurs donn es par le programme timslip La valeur moyenne obtenue en angle et en amplitude permet de corriger le d faut d alignement du transducteur qui est normalement une constante 17 Le tableau ci dessous donne les valeurs de tol rance maximum pour les param tres de la calibration prises en compte par le p
50. ositions comprennent pour chaque ligne un temps exprim en unit s sp cifiques CODAS une latitude et une longitude degr s centi mes Le syst me d acquisition centralis du bord permet l obtention d un fichier PODAS POlarstern Data Acquisition System A priori et selon les informations transmises par l ing nieur lectronicien en 2006 toutes les donn es de navigation transmises ADCP proviennent de la centrale inertielle MINS La v rification des donn es de navigation sera effectu e au moment du water tracking Les diff rentes sources de positionnement r pertori es dans le fichier PODAS sont les suivantes MINS scientific navigation platform Centrale inertielle NACOS ship navigation system Pas de d tails Ashtec Z 12 DGPS GPS diff rentiel Trimble 1 et 2 GPS 2 GPS classiques Leica MX400 GPS GPS classique De m me que pour le cap les probl mes d horloge ventuels ne seront visibles clairement qu au moment de la calibration par water tracking Contrairement la campagne pr c dente de 2006 les trois fichiers de positions dont on dispose hdg gps1 et gps2 sont pratiquement la m me date une seconde pres On a suppos que comme pr c demment l horloge ADCP donc hdg tait l horloge du GPS Trimble 1 Si c est le cas le d calage entre cette horloge et celle de la centrale MINS qui avait t observ en 2006 a t corrig Mais la configuration a peut tre t cha
51. rogramme adcpcal delta u min 100 00 max 100 00 delta v min 100 00 max 100 00 clip amp 0 1 clip ph 3 0 clip dt 60 clip var 0 100 La figure ci dessous illustre un exemple du r sultat du programme adcpcal DRAKEO _ raw corags1 _ 9 fixes unedited npts 104 20 DRAKEO9 raw corags1 _ 9 fixes edited npts 64 2 x 3 S amp oi 2 D 6 DWR gt amp BEGG ow gen def Berio De o phase 20 amplitude o b oa gt a i 100 amd 2 O a sel N 50 8 ine k 4 E Of te Taat beb litt PAL by o E bz 84 36 38 30 oz 94 96 38 50 i L i i i i i days 82 B4 B B W 2 M B 15 20 r r r Aa 4 a 4 amplitude A t o a 10 0 m mm 10 0 10 20 30 40 50 60 70 time nav pc o 2009 5 9 15 17 2009 5 9 15 16 Fic 6 3 D termination de la phase de l amplitude et d un d calage d horloge Sur le graphique de droite le trac time nav pc montre effectivement qu il n y a pas de tendance de variation privil gi e de cette diff rence de date mais que les sauts sont al atoires 6 4 R sultats de la calibration Le r sultat de l dition et de la calibration est le suivant Le d calage d horloge est ramen 5 secondes en moyenne avec toutefois un cart type lev 13 9 ce qui confir
52. sducteur On a galement port sur la figure 3 points rouge correspondant aux vitesses du son calcul es d apr s les mesures de 3 stations CTD DRAKE 2009 Thermosalino quille et transducteur ADCP T T T T T T T 1530 20F 41510 15H 4 1490 o 910 1470 8 vitesse du son 5 1450 or 41430 5 L L L L L L L 1 L 1410 80 82 84 86 94 96 98 100 90 92 jours decimaux 2008 Fic 5 1 Temp rature transducteur ADCP rouge thermosalinographe magenta et vitesse du son La calibration par water tracking effectu e dans la suite du traitement confirmera que la fen tre tait bien remplie d eau sal e mais que la temp rature donn e par le transducteur est sans doute trop lev e En effet si la vitesse du son est trop faible la vitesse ADCP est aussi trop faible et la m thode du water tracking donnera un facteur d amplitude A gt 1 La calibration permet de corriger cette erreur sur l amplitude 11 12 On proc de ensuite l edition des profils ADCP En ce qui concerne les donn es des profils le premier stade est le passage de programmes de calculs de statistiques qui permettent de d terminer les seuils d erreurs maxima pour plusieurs param tres vitesse composantes u et v amplitude erreur verticale Ces calculs se font sous Matlab Les valeurs utilis es pour la campagne DRAKE sont les suivantes une seule dition pour les trois parties ev_threshld 150 w_var_thres
53. tion se fait sous Matlab programmes load_sta m ou load_sta_manual m Le programme d code les fichiers STA et cr e plusieurs fichiers en sortie cmd Fichiers ex cutables de commande qui seront utilis s dans la proc dure suivante pour la fin du chargement bin Fichiers de donn es ADCP binaires gpsl gps2 Fichiers ASCII de navigation gps1 en d but d ensemble moyenn gps2 en fin La deuxi me op ration est le lancement de l ex cution des fichiers de commande cmd pour charger la base de donn es proprement dite A ce stade on dispose d un r pertoire nomm g n ralement adcpdb pour ADCP Data Base dans lequel les donn es binaire ADCP sont maintenant organis es en blocs dont la longueur et le contenu sont d finis par l utilisateur au moyen d un fichier sp cial de d finition vmadcp def plac dans le r pertoire de chargement Pour la campagne DRAKE 2009 on a choisi de faire tout le traitement en une seule partie cause du petit nombre suppos de points de calibration peu de stations Pour les trac s de pr sentation on s parerera la campagne en trois parties une section Nord Sud depuis la sortie du d troit de Magellan jusqu 61 Sud puis les stations pr s de King george Island et dans le d troit de Bransfield et ensuite la section de retour en travers du Passage de Drake de 61 Sud Varr t de l appareil DRAKE 1 2009 03 22 17 54 42 to 2009 04 01 00 17 00 DRAKE 2 2009 04 01 00 17
54. y DATA Please show FULL path to the model files in the Model_ file unless they are located in DATA download default 4th line in Model_ file is used ONLY for models calculated on uniform grid in km Converting functions are provided with current version of OTPS Please upgrade OTPS if a function is missing latitude longitude and time file is an ASCII file consisting eighter of 2 OR 8 columns latitude degrees longitude degrees yyyy mm dd hh mm ss Common sign convention lat gt 0 degrees North lon gt 0 degrees East lat lt 0 degrees South lon lt 0 degrees West yyyy mm dd hh mm ss year month day hour seconds GMT These 6 columns are needed for tide predictions only You may leave them empty if only extracting HC Extract HC predict tide for z U V u v elevation m WE transport m 2 s SN transport m 2 s WE velocity cm s SN velocity cm s NOTE Changed Nov 2004 now for any of u v U V all components that is U V u v are calculated by predict_tide Constituent names should be in LOWER case and separated by comma Leave the line blank if ALL model constituents are included AP RI output amplitude and phase GMT OR real imaginary parts 43 Only used when extracting HC May leave blank if predicting tide 6 Only needed for tide prediction 1 0 Do Not correct for the minor constituents defined in weights h Only works if the 8 tidal constituents fields are in t
Download Pdf Manuals
Related Search