UPD_NRH2DFITS $SSW/radio/nrh/idl/2dnrh/fichier/upd_nrh2dfits.pro
[Previous] [Next]
 NAME:
	UPD_NRH2DFITS

 PURPOSE:
	Cette procedure modifie les HEADER des fichiers FITS NRH2D, afin de
	les rendre compatible avec les modifications du logiciel

 CATEGORY:
	Fichiers FITS

 CALLING SEQUENCE:
	UPD_NRH2DFITS, File

 INPUTS:
	FILE 	Nom du fichier a modifier

 KEYWORD PARAMETERS:

 OUTPUTS:

 CALLED BY:
	COMP_DATE, FOPEN_NRH2
 PROCEDURE:
	Cette procedure devra etre adaptee chaque fois qu'une modification
	des Headers aura des repercussions sur le traitement.
	Elle modifie le contenu des KEYWORDS directement dans le fichier. Il
	faut donc avoir un droit d'ecriture sur le fichier pour l'executer.

	Cette procedure remet la date dans le keyword DATE du header primaire
	dans le cas de modifications

	Elle est appelee par CW_NRHFILE, dans le cas de fichiers 2D.

 LISTE DES MISES A JOUR:
	AVRIL 98	Remplace STOKES I et STOKES V par STOKESI et STOKESV

 CALLS: ***
	BIN_DATE, COMP_DATE, FORM_DATE, FXHMODIFY [1], FXHMODIFY [2], FXPAR [1], FXPAR [2]
	HEADFITS [1], HEADFITS [2], HEADFITS [3], STR_SEP
 EXAMPLE:

 MODIFICATION HISTORY:
	Ecrit par: J Bonmartin le 07/04/98


UPDATE_TPRECS $SSW/radio/ovsa/idl/util/update_tprecs.pro
[Previous] [Next]
 NAME:
     UPDATE_TPRECS
 PURPOSE:
     Routine to generate the TPUPD and (optionally) TPCAL segment records
     from CTRCAL data, and return them for writing into the current .ARC
     file.
 CATEGORY:
     OVRO APC DATA CALIBRATION
 CALLING SEQUENCE:
     result = UPDATE_tprecs(a,nrec,header,antpol,f,afac,arms,sdaryrec,primrec,$
                   updflag,updflag2m[,/reset])
 INPUTS:
     a          The ASSOC variable associated with the .ARC file, as returned
                  from the call to OPENARC().
     nrec       The number of records in the file, as returned from the call to
                  OPENARC().
     header     The header structure of the CTRCAL scan that contains the data
                  on which the TPUPDATE is based.
     antpol     The ANTPOL array from the CONFIG segment, specifying the identity
                  of the feeds on each antenna.
     f          The frequency list, in GHz, of size NF.
     afac       An array containing the total power update factors,
                  of size (3,nant,NF), where the first index is the polarization
                  (0=RCP, 1=LCP, 2=Linear), the second is the antenna,
                  and the third is the frequencies.  The calibration factors are
                  near unity (hopefully) and are to be multiplied by the TPCAL
                  factors to update them.
     arms       An array of identical size to AFAC, giving the corresponding RMS
                  deviations derived from multiple measurements, in the same
                  units as for AFAC.
     updflag    A flag that indicates whether we are updating the primary or
                  secondary calibration for 27 m antennas.  UPDFLAG = 1 => Secondary
                                                            UPDFLAG = 2 => Primary
     updflag2m  A flag that indicates whether we are updating the primary or
                  secondary calibration for 2m antennas.  UPDFLAG2M = 1 => Secondary
                                                          UPDFLAG2M = 2 => Primary
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
     reset      A switch that sets all entries to NaN (even those not in freq list)
                  before placing data into the records
 ROUTINES CALLED:
     bin_date, tl_encode, tl_decode, dayofyr, getdata
 OUTPUTS:
     result     A flag specifying success [1] or failure [0]
     sdaryrec   The secondary calibration (TPUPD) segment, in four records.  The
                  size is [1024,4]
     primrec    The primary calibration (TPCAL) segment, in four records as for
                  SDARYREC, or -1 if no primary calibration is requested.
 COMMENTS:
     The values written to the file are in scaled, 2-byte integer form.  The
     conversion factors are different for various data sources.  Primary calibration
     conversion factors are 10 for 27 m, 100 for 2 m; Secondary calibration
     conversion factors are 1000 for both sources of antenna.  Values for which
     there is no measurement (e.g. for certain frequencies) the factors are unity
     and the rms is zero in the secondary calibration segment.
 CALLS: ***
	BIN_DATE, DAYOFYR [1], GETDATA, OVSA_GET_INDEX, TL_DECODE, TL_ENCODE, dayofyr [2]
 CALLED BY:
	CTRCALCHEK
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
     Written 17-Nov-1999 by Dale E. Gary
     04-Dec-1999  DG
       Changed name from MAKE_TPUPDATE to WRITE_TPUPDREC, and modified the routine
       to parallel the similar function routine WRITE_GAINREC.
     08-Dec-1999  DG
       Changed to update either the primary or secondary calibration depending on
       the value of UPDFLAG.
     11-Jan-2000  DG
       Eliminated hardwired directory locations.
     03-Feb-2000  DG
       Considerable rewrite, to allow the routine to handle both secondary and
       primary update to calibration.
     04-Feb-2000  DG
       Changed name from WRITE_TPUPDREC to UPDATE_TPRECS to better reflect its
       purpose.
     18-Jul-2001 GN
       Replaced get_index by ovsa_get_index to avoid conflict with Yohkoh
     08-Nov-2003  DG
       Added ANTPOL argument to function call, and made output of factors
       dependent on contents of that array, so that columns reflect the correct
       antenna polarization (not hardwired to 1I 1L 2R 2I, etc.)
     24-Aug-2005  DG
       Considerable changes to reflect R and L feeds on all antennas.  TPCAL and
       TP-UPDATE segments are now doubled in size (changed from 2 to 4 records),
       and the first two records contain R data, while last two have L data.  Now,
       PRIMREC and SDARYREC have 4 records in all cases, but for data before 2004.231,
       only the first two are to be used.
     27-Aug-2005  DG
       Added RESET switch, to set all values to NaN before updating.


UPTIME $SSW/radio/ovsa/idl/util/uptime.pro
[Previous] [Next]
 NAME:
     UPTIME
 PURPOSE:
     Given an optional date (system date used if not specified) and a
     flux limit, return a list of calibrators with flux above the flux
     limit, along with times of rising and setting.
 CATEGORY:
     OVRO APC DATA CALIBRATION
 CALLING SEQUENCE
     upt_stc = uptime(fmin[,jd=jd])
 INPUTS:
     fmin     the minimum flux (over the entire band) that a matching
                calibrator can have.  Calibrators with a smaller
                flux density at any frequency will not be returned in
                the list.
 OPTIONAL (KEYWORD) INPUT PARAMETERS:
     jd       the Julian date for which to calculate the rising and
                setting times
 ROUTINES CALLED:
     cvt_radec (internal), caldat, dayofyr, lst, precess, msec2str
 OUTPUTS:
     upt_stc  a structures containing the quantities being returned, as defined below:
         UPT_STC = {epoch: '' $ ; String containing date, format [YYYY-MM-DD]
                     info: replicate(uptime,  $ ; Variable number of copies
                           name: '', $
                             ra: 0L, $ ; RA of Epoch [0.1 mdeg]
                            dec: 0L, $ ; Dec of Epoch [0.1 mdeg]
                           rise: '', $ ; Rise time of day string format [HH:MM]
                            set: '', $ ; Set time of day string format [HH:MM]
                           flux: {L:0, C:0, X:0, U:0}}, $ ; Flux [Jy] in L,C,U and X bands
                           nsel)
                    }
 COMMENTS:
 CALLS: ***
	CALDAT, CVT_RADEC, DAYOFYR [1], GET_SOLEPH, LST, MSEC2STR, PRECESS, dayofyr [2]
 CALLED BY:
	AUTOCAL, WUPTIME, check_cals
 SIDE EFFECTS:
 RESTRICTIONS:
 MODIFICATION HISTORY:
     Written 13-Mar-1999 by Dale E. Gary
     15-Mar-1999  DG
       Significant modification of format of output structure.
     16-May-1999  DG
       Added SUN to uptime list.
     11-Jan-2000  DG
       Eliminated hardwired directory locations.
     12-Jan-2000  DG
       Include calculation for restricted HA range at low declinations,
       given [in minutes of time] by:
         restricted_ha = -51.15 - 1.00287*dec + 0.0736*dec^2, for dec < -21.5 degrees).
     10-Dec-2002  DG
       Bring in low-dec corridor limits by 2 min on each side, to avoid hitting
       the limit.


UTPLOT [4] $SSW/radio/ethz/idl/ragview/display/utplot.pro
[Previous] [Next]
 NAME:
	UTPLOT
 PURPOSE:
	Plot X vs Y with Universal time labels on bottom X axis.  X axis
	range can be as small as 5 msec or as large as 20 years.
 CALLING SEQUENCE:
	UTPLOT,X,Y,BASE_TIME( or Utstring or Xst0), $
		timerange = timerange, LABELPAR=LBL, /SAV,TICK_UNIT=TICK_UNIT,$
		MINORS=MINORS, /NOLABEL, ERROR=ERROR,$
               [& ALL KEYWORDS AVAILABLE TO PLOT]
       UTPLOT, roadmap, y
       utplot, x, y, '1-sep-91'
       utplot, x, y, '1-sep-91', timerange=['1-sep-91', '2-sep-92']
       utplot, x, y, '1-sep-91', timerange=[index(0), index(i)], xstyle=1
 INPUTS:
	X -	X array to plot in seconds relative to base time.
               (MDM) Structures allowed
	Y -	Y array to plot.
	Base_time (optional)- reference time, it's purpose is as a fiducial
		time to preserve the precision of the graphics, that is
		the range over which a plot is drawn must not be too small
		with respect to the absolute size of the start of the plot or
		roundoff errors will create strange effects.  In previous
		versions of Utplot, Base_time was called Utstring or Xst which
		are explained below:
	        UTSTRING - 
		(Optional) ASCII string containing base time for X axis. 
		Format for time is YY/MM/DD,HHMM:SS.XXX.  If this parameter 
		isn't present  	and hasn't been set yet (by a previous call to 
		UTPLOT, via routine SETUTBASE, or directly), user is 
		prompted for base time.
      		xst0 -   
		Optional. The reference time to use for converting a structure
               into a seconds array (OR) the time for the first value if
               passing a double prec. array.
		
		N.B.-
		While there are 

       timerange - Optional. This can be a two element time range to be
               plotted.  It can be a string representation or structure.

       LBL -   (Labelpar=LBL) 2 element vector selecting substring from publication format
               of ASCII time (YY/MM/DD, HH:MM:SS.XXX).  For example,
               LBL=[11,18] would select the HH:MM:SS part of the string.
       SAV -   If set, UTPLOT labels, tick positions, etc. in !X... system
               variables will remain set so that they can be used by
               subsequent plots (normally !x... structure is saved in
               temporary location before plot and restored after plot).
               To clear !x... structure, call CLEAR_UTPLOT.
       TICK_UNIT - Time (in seconds) between Major tick marks is forced to TICK_UNIT
                   Has no effect for axis longer than 62 days.
       MINORS    - Number of minor tick marks is forced to MINORS
       NOLABEL If set then UTLABEL isn't printed on plot
       XTITLE - text string for x-axis label - If the string contains
                4 asterisks ('****'), the UT time will be substituted
                for that substring
	Year - Force the year in the x axis labels
	YOHKOH - Use Yohkoh style labels, e.g. '03-May-93 18:11:30.732'
 OPTIONAL OUTPUT PARAMETERS:
	ERROR - 0/1.  1 means there was an error in plotting.
 CALLS: ***
	ANYTIM2UTPLOT, ATIME [1], ATIME [2], CHECKVAR [1], Ex2Int [1], Ex2Int [2], FCHECK
	GETUTBASE [1], GETUTBASE [2], SET_UTPLOT [1], SET_UTPLOT [2], anytim2ex [1]
	anytim2ex [2], checkvar [2], int2secarr [1], int2secarr [2], str_replace [1]
	str_replace [2]
 CALLED BY:
	ACRIM, ACRIM_EV, ACRIM_OUT, ACRIM_READER, AR_AREAS, BATSE_FLARES_WEEK
	BCS_24HR_PLOT [1], BCS_24HR_PLOT [3], BCS_SPMOVIE [1], BCS_SPMOVIE [2]
	BORN_AGAIN, BOX_LC, CHECK_SFC_PNT, COMPUTE_MONEXP_FACTORS, COMPUTE_MONEXP_POLY
	DEMODULATOR_TEST, ECLIPSE_PLOTTER, EVAL_SHERB, FE25_BSC_TEMP, FIND_FAST_DATA
	FLAREMETER, FS_GRAPH, GET_HXRBS_FITS, GIS_UTPLOT, GLITCH_TRACKER
	HESSI FASTRATE CLASS DEFINITION, HESSI PACKET CLASS DEFINITION
	HESSI ROLL_DB CLASS DEFINITION, HSI_ANY_SOHDATA, HSI_DO_QLOOK_PLOTS
	HSI_NPAK_HISTOGRAM, HSI_VCFILE_PACKET_RATE, HXARCHIVE, HXFITS, HXISLIFE, HXISTPLOT
	HXT_4CHPLOT, HXT_AUTOIMG, HXT_BACK, HXT_FANS, HXT_LCURVE, HXT_SPECTRA, HXT_SUMMER
	LCBDA, NEUPERT_PLOT, NIS_BIAS_DEMO, NOV93_ECLIPSE, NOV_ECLIPSE, OBS_PLOT, ORB_EXTRAP
	OVER_THE_LIMB, OVSA_EXPLORER formerly OVSA_PRESUB, OVSA_PRESUB, OVSA__DEFINE
	PIXEL_LC, PLOTBATSE_QL, PLOT_BSC [1], PLOT_BSC [2], PLOT_BSC_AS2, PLOT_CDS_CRAYS
	PLOT_ERR, PLOT_EXPINT, PLOT_EXP_FACTOR, PLOT_GCPARM, PLOT_GRSPCH, PLOT_GRSPCL
	PLOT_GRSPHH, PLOT_GRSPHL, PLOT_HKLZ, PLOT_HT, PLOT_HXSPC, PLOT_HXSPH, PLOT_MED
	PLOT_MONEXP_STD, PLOT_OCCS, PLOT_OCC_ONDAY, PLOT_SUBHTR, PROFILES2, Plot
	QL_DAY_PLOT, QL_ORBIT_PLOT, QUICKATT [1], QUICKATT [2], QUICKDARK [2]
	QUICKDRIFT [1], QUICKDRIFT [2], QUICKSTRAY [2], QUICK_DPE, READ_HXT_CAL, ROLL_PLOT
	SET_UTAXIS, SOLAID, SPARTAN_PLANNER, SPEX_THISTORY [1], SPEX_THISTORY [2]
	SPLINE_PLOT, SXSPC, SXT_HYST, SYNOPTIC_STATS, Sxt_goes_teem, TERM_QUICK, TERM_REVIEW
	TERM_REVIEW_GIF, TEST_ATT, TV_SSX, TV_SYNOP, UTPLOT_IMAGE, UTPLOT_IO [1]
	UTPLOT_IO [2], UTPLOT_IO [3], UTPLOT__DEFINE, WINDMETER, WL_CUBE_II, WPLOT_NONOP
	WPLOT_SUBHTR [2], W_HISTOGRAM, XPLOT_HT, auto_bck_find, auto_bck_find0, avsig_temp
	bcs_survey, box_lc_plot, box_sxthxt_fsp, cal_fig_mich, db_gbo, edac_summary
	eis_status_sts1_conversions [5], eis_status_sts2_conversions [1]
	eis_status_sts3_conversions [1], eit_proton_summary, emi_plot [1]
	fit_bsc_plot, get_utevent [1], get_utevent [2], go_teem_plots, goes_plot [1]
	goes_plot [2], goes_plot [3], goes_plot [4], goes_plot [5], goes_teem, goes_widget
	goesem_eqn [1], goest_eqn [1], gr_summary_plot, grs32_fsp, grsl_fsp, high_mag_lat
	hsi_as_plquality, hsi_obs_times, hsi_sas_plot_scpointing, hsi_test_packet_rate
	hxs_fsp, hxsgrs_fsp, hxt_fsp, hxt_survey, hxt_utplot, hxtbox_fsp, imspec [10]
	imspec [11], imspec [12], imspec [13], imspec [14], imspec [15], imspec [16]
	imspec [17], imspec [18], imspec [19], imspec [1], imspec [20], imspec [21]
	imspec [22], imspec [23], imspec [24], imspec [25], imspec [26], imspec [27]
	imspec [28], imspec [29], imspec [2], imspec [30], imspec [31], imspec [32]
	imspec [33], imspec [34], imspec [35], imspec [36], imspec [3], imspec [4]
	imspec [5], imspec [6], imspec [7], imspec [8], imspec [9], jitter_gif_xyimg
	last_lc, lcbsd, lchda, lcobs, lcur_plot [1], lcur_plot [2], lcwda, lsplot
	mon_health [1], mon_health [2], mwlt_select, new_edac_summary
	new_mon_health [1], new_mon_health [2], ovsa_interact, pl_dt_genx
	pl_monitor_center, plot_door_open, plot_evn, plot_fem, plot_fft, plot_gbl [1]
	plot_gbl [2], plot_img_cadence, plot_lcur [1], plot_lcur [2], plot_loi_mmad
	plot_loi_summary [1], plot_loi_summary [2], plot_pzt_corr, plot_shutter_perf
	plot_ssl, plot_ssw_fov, plot_sxg, plot_sxl [1], plot_sxl [2], plot_temps2 [1]
	plot_trav [1], plot_trav [2], plot_trav [3], plotbft [1], plotbft [2], plotltc [1]
	plotltc [2], plott_bda, plott_bda_dp, plott_hda, plott_hi_sum, plott_low_sum
	plott_med1_sum, plott_med2_sum, plott_wda, ploty [1], ploty [2], pmtras_analysis
	psk_correlate, rd_so_at_ftr, rd_sxs_pha [1], rd_sxs_pha [2], rd_therm_rs232
	read_yohkoh_4_spex [1], read_yohkoh_4_spex [2], ref_day_plot, res_freq
	show_obs2, show_obs3, show_obs4, sld_sls_vstrips, sld_vstrips, soon_cadence
	soon_catstat, sooncat_cadence, soup_obs2, spectro_plot2, ssw_apkpbar, stt_plot
	sxt_align, sxt_temps, sxtbox_fsp, tim2tfss, tprofiles, uvspplot, wdefroi [1]
	wdefroi [2], yo_height
 COMMON BLOCKS:
	COMMON UTCOMMON, UTBASE, UTSTART, UTEND, xst = base, start, and 
	end time for X axis in double precision variables containing
	seconds since 79/1/1, 00:00. xst is the fully referenced start
	time of the plot.
       COMMON UTSTART_TIME, PRINT_START_TIME = 0/1 ( don't/do print
       start time label on plot), ATIME_FORMAT ('YOHKOH'  or 'HXRBS')
 SIDE EFFECTS:
	X vs Y plot is produced on current graphics device.  The normal 
	X-axis is replaced by an axis with ticks and labels showing the 
	universal time.  The start time of the plot is displayed in the 
	upper right inside corner of the plot if SET_UTLABEL,0 hasn't been
	called.
 RESTRICTIONS:
	Cannot be used for times before 1-JAN-1950 or after 31-DEC-2049
	Range of X axis can be anywhere between 5 msec and 20 years.
	Keywords DATA, DEVICE, NORMAL are not available
	The utbase time is only set if either base_time
	is passed or required as input.
 PROCEDURE:
	If start or end time hasn't been set, autoscale X axis. 
	If either has been set (via routines SETUTSTART	and SETUTEND, 
	or directly), only data between times selected will be displayed; 
	i.e. X min = UTSTART - UTBASE; X max = UTEND - UTBASE.  
	Calls SET_UTPLOT using keyword LABELPAR to customize X 
	axis labels and tickmarks.  Otherwise it uses all normal plotting 
	procedures and 	the !X and !Y system variables.
	Note:   Format of time written in labels differs slightly from format
	used to pass times to routines.  Input format contains only one colon
	between minutes and seconds (makes the meaning of a partial string 
	precise) while labels include an extra colon between hours and minutes 
	(more acceptable for publication).
 MODIFICATION HISTORY:
	Written by Kim Tolbert, 4/88.
	Mod. by Richard Schwartz to Version 2 91/02/17
       Mod. by Richard Schwartz for compatibility with OPLOT. 3/26/91
	Mod. by RAS, keywords TICK_UNIT and MINORS added 91/05/02
	MOD. BY RAS TO ACCEPT ALL PLOT KEYWORDS, 91/10/03
       MOD. by AKT 11/12/92. Added error keyword argument.
       Mod. ras, yohkoh style keyword
	removed DATA, DEVICE, NORMAL keywords 14-May-93, ras
		so version 2 can compile the many keywords
	RAS, 93/6/8 -
	If X window, make sure window is open before saving any bang variables
	otherwise Linecolors will fail!
       Mod by MDM to work with YOHKOH data (structure data types)
       Mod by MDM to expand the COMMON block so that OUTPLOT will work
       Mod by MDM 21-Apr-92 to not set !quiet
       Mod by SLF 26-Apr-92 to reinstate the xtitle option
       Mod by MDM 28-Aug-92 to add keyword YEAR to print the year on the
                 tick label
       Mod by MDM  9-Apr-93 to added TIMERANGE option
       Mod by MDM 12-Apr-93 to removed ZMARGIN since older IDL versions had
                 trouble with 9-Apr addition of TIMERANGE parameter
	1-NOV-93, ras, integrate Yohkoh and SDAC software packages
	10-Nov-93, ras, accept single element arguments for x0 and y0
	14-Dec-93, MDM, - Added keyword XTHICK and YTHICK
			- Set the x-axis label to be Start Time ... by default
	4-jan-94, ras   - undeclared base_times are set to start of day
			  for fully referenced times, used anytim to make time processing
			  more resilient to the type of the x axis variable
	28-jan-94, ras  - fixed ynozero problem
	7-apr-94,  ras  - clarified some of the documentation on
			  tick_unit and labelpar
	28-jul-94, ras - added keyword inheritance and eliminated explicit reference to zkeywords
       8-Aug-94,  MDM  - Modified so that /NOLABEL is recognized with Yohkoh time
                         format
	10-aug-94, ras  - xtype is set explicitly to zero to facilitate later use of oplot
		          if utplot was used after a call to plot_oi, xtype not set properly
	16-aug-94, ras  - added zvalue keyword and removed _extra from call to oplot
			  treatment of symsize is changed to allow control by !p.symsize
	2-aug-95, ras   - symsize was not being utilized through !p.symsize
			  This is a clear IDL bug in version 4, also present in 
			  version 3.6.1.c, !p.symsize is not being used, where
			  using the keyword symsize explicitly results in the
			  expected symbol size?
			  The second fix is to have the routine explicitly check for
			  data values within the requested range and to plot only the
			  plot box without data points.  So utplot will
			  be much more resilient in the future and will not cause
			  a code-stopping error.  However, an informative
			  message is issued using "message,/continue"


UVALUE_REGION $SSW/radio/nrh/idl/2dnrh/gen/cw_region.pro
[Previous] [Next]
 NAME:
	UVALUE_REGION

 PURPOSE:
	Cette procedure reactualise le contenu de UVALUE pour le widget
		CW_REGION

 CATEGORY:
	Widgets

 CALLING SEQUENCE:
	UVALUE_REGION,Xdim,Ydim,Vdim

 INPUTS:
	XDIM	dimension X du tableau contenant l'image
	YDIM	dimension Y du tableau contenant l'image
	VDIM	valeur superieur de la dimension de l'image

 KEYWORD PARAMETERS:

 OUTPUTS:
	remise a jour du contenu de la structure STATE
		state{size_imx, size_imy, size_vis}
	remisse a jour des UVALUE de CURSOR et RESET

 CALLS: ***
	BOX_CURSOR2, CW_FIELD, CW_REGION, IND_RS, REGION_EVENT, RS_IND, STR_SEP
 CALLED BY:
	CW_REGION, REGION_EVENT
 COMMON BLOCKS:
	IDENT_REG: 
		cursor	identification du widget 'CURSOR'
		reset	identification du widget 'RESET'
		textx	identification du widget X:...
		texty	identification du widget Y:
		state	structure

 EXAMPLE:
	UVALUE_REG,state.size_imx,state.size_imy,state.size_vis
		pour reactualiser apres changement de donnees ou de 
		dimension
 MODIFICATION HISTORY:
 	Ecrit par: J Bonmartin le 03/03/98 (bonmartin@obspm.fr)