Changeset 17 for trunk/step1_diff.pro
- Timestamp:
- 03/21/07 14:18:07 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/step1_diff.pro
r16 r17 13 13 ; 14 14 ; @examples 15 ; test (for header and delta) 15 ; test (for header and delta) 16 16 ; to compute difference between cond_sed in cond_sed_ORCA2.nc and 17 17 ; itself and write it in ginette.nc : 18 ; rhodes\@reee522$ cd /home/rech/eee/reee522/incas/geomag/geomag_ws 19 ; rhodes\@reee522$ idl 20 ; IDL> step1_diff, 'cond_sed_ORCA2.nc', 'cond_sed_ORCA2.nc', 'ginette.nc' 18 ; $ cd ${GEOMAG} 19 ; $ idl 20 ; IDL> file1=getenv('GEOMAG') + '/cond_sed_ORCA2.nc' 21 ; IDL> file2=getenv('GEOMAG_OD') + '/cond_sed_ORCA2.nc' 22 ; IDL> step1_diff, file1, file2, 'ginette.nc' 21 23 ; values of delta must be 0 everywhere 22 ; rhodes\@reee522$ rm ginette.nc24 ; $ rm ginette.nc 23 25 ; idem with BR file 24 ; IDL> step1_diff, 'Br_ORCA2.nc', 'Br_ORCA2.nc', 'ginette.nc' 26 ; IDL> file1=getenv('GEOMAG') + '/Br_ORCA2.nc' 27 ; IDL> file2=getenv('GEOMAG_OD') + '/Br_ORCA2.nc' 28 ; IDL> step1_diff, file1, file2, 'ginette.nc' 25 29 ; 26 30 ; real life … … 28 32 ; /home/rech/eee/reee522/geomag_d/cond_sed_ORCA2.nc and write it in 29 33 ; cond_sed_ORCA2_diff.nc : 30 ; rhodes\@reee522$ cd /home/rech/eee/reee522/incas/geomag/geomag_ws31 ; rhodes\@reee522$ idl32 ; IDL> step1_diff, 'cond_sed_ORCA2.nc', $33 ; '/home/rech/eee/reee522/geomag_d/cond_sed_ORCA2.nc', $34 ; 35 ; IDL> step1_diff, 'Br_ORCA2.nc', $36 ; '/home/rech/eee/reee522/geomag_d/Br_ORCA2.nc', $37 ; 34 ; $ cd ${GEOMAG} 35 ; $ idl 36 ; IDL> file1=getenv('GEOMAG') + '/cond_sed_ORCA2.nc' 37 ; IDL> file2=getenv('GEOMAG_OD') + '/cond_sed_ORCA2.nc' 38 ; IDL> step1_diff, file1, file2, 'cond_sed_ORCA2_diff.nc' 39 ; IDL> file1=getenv('GEOMAG') + '/Br_ORCA2.nc' 40 ; IDL> file2=getenv('GEOMAG_OD') + '/Br_ORCA2.nc' 41 ; IDL> step1_diff, file1, file2, 'Br_ORCA2_diff.nc' 38 42 ; 39 43 ; to see the difference, for example … … 41 45 ; IDL> xxx,'Br_ORCA2_diff.nc' 42 46 ; select xy on plt wigdet 43 ; 47 ; 44 48 ; @history 45 49 ; reee522 2006-12-18T10:38:49Z rhodes (IRIX64) … … 53 57 ; $Id$ 54 58 ; 55 PRO step1_diff, file1, file2,file359 PRO step1_diff, file1, file2, file3 56 60 ; 57 61 ncverbose=1 … … 102 106 checknonzero=where(delta NE 0.,count) 103 107 IF count EQ 0 THEN BEGIN 104 msg = ' 108 msg = 'iii : delta is zero everywhere' 105 109 PRINT, msg 106 110 ENDIF ELSE BEGIN 107 msg = ' 111 msg = 'iii : delta is not zero ' + STRING(count) + ' times' 108 112 PRINT, msg 109 113 ENDELSE
Note: See TracChangeset
for help on using the changeset viewer.