- Timestamp:
- 03/28/07 14:17:24 (17 years ago)
- Location:
- trunk/SRC
- Files:
-
- 24 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Grid/computegrid.pro
r238 r240 608 608 ; check that xtest is now sorted in the increasing order 609 609 IF array_equal(sort(xtest), lindgen(jpi)) EQ 0 THEN BEGIN 610 print, 'WARNING: we cannot sort the xaxis with a simple shift...'611 print, 'we force key_periodic = 0 and key_shift = 0'612 print, 'only horizontal plot may be ok...'610 ras = report(['WARNING: we cannot sort the xaxis with a simple shift...', $ 611 'we force key_periodic = 0 and key_shift = 0', $ 612 'only horizontal plot may be ok...']) 613 613 key_periodic = 0 614 614 xnotsorted = 1 -
trunk/SRC/Grid/micromeshmask.pro
r238 r240 69 69 filein_msk = (findfile(filein+'.mask.nc'))[0] 70 70 IF filein_hgr EQ '' OR filein_zgr EQ '' OR filein_msk EQ '' THEN BEGIN 71 ras = report( 'meshmask file(s) not found...')72 ras = report( filein+' does not exist')73 ras = report( filein+'.mesh_hgr.nc does not exist')74 ras = report( filein+'.mesh_zgr.nc does not exist')75 ras = report( filein+'.mask.nc does not exist')71 ras = report(['meshmask file(s) not found...', $ 72 filein+' does not exist', $ 73 filein+'.mesh_hgr.nc does not exist', $ 74 filein+'.mesh_zgr.nc does not exist', $ 75 filein+'.mask.nc does not exist']) 76 76 return 77 77 ENDIF -
trunk/SRC/Grid/ncdf_meshread.pro
r232 r240 101 101 tempsun = systime(1) ; for key_performance 102 102 IF keyword_set(CHECKDAT) THEN BEGIN 103 print, ' The keyword CHECKDAT has been suppressed (it could create bugs).'104 print, ' Remove it from the call of ncdf_meshread'105 print, ' Please use smallmeshmask.pro or micromeshmask.pro to create a'106 print, ' meshmask that has manageable size'103 ras = report([' The keyword CHECKDAT has been suppressed (it could create bugs).', $ 104 ' Remove it from the call of ncdf_meshread', $ 105 ' Please use smallmeshmask.pro or micromeshmask.pro to create a', $ 106 ' meshmask that has manageable size']) 107 107 return 108 108 ENDIF … … 218 218 ; 219 219 IF array_equal(sort(xaxis), lindgen(jpi)) NE 1 THEN BEGIN 220 print, 'the x axis (1st line of glamt) is not sorted in the increasing order after the automatic definition of key_shift'221 print, 'please use the keyword shift (and periodic) to suppress the automatic definition of key_shift (and key_periodic) and define by hand a more suitable value...'220 ras = report(['the x axis (1st line of glamt) is not sorted in the increasing order after the automatic definition of key_shift', $ 221 'please use the keyword shift (and periodic) to suppress the automatic definition of key_shift (and key_periodic) and define by hand a more suitable value...']) 222 222 widget_control, noticebase, bad_id = nothing, /destroy 223 223 return -
trunk/SRC/Grid/ncdf_meshroms.pro
r231 r240 196 196 ; 197 197 IF array_equal(sort(xaxis), lindgen(jpi)) NE 1 THEN BEGIN 198 print, 'the x axis (1st line of glamt) is not sorted in the increasing order after the automatic definition of key_shift'199 print, 'please use the keyword shift (and periodic) to suppress the automatic definition of key_shift (and key_periodic) and define by hand a more suitable value...'198 ras = report (['the x axis (1st line of glamt) is not sorted in the increasing order after the automatic definition of key_shift', $ 199 'please use the keyword shift (and periodic) to suppress the automatic definition of key_shift (and key_periodic) and define by hand a more suitable value...']) 200 200 widget_control, noticebase, bad_id = nothing, /destroy 201 201 return -
trunk/SRC/Grid/numbisland.pro
r238 r240 90 90 land = (where(msk EQ 0, count))[0] 91 91 WHILE count NE 0 DO BEGIN 92 IF (islnumb-9) MOD 10 EQ 0 THEN print, 'island number :'+strtrim(islnumb-9, 1) 93 92 IF (islnumb-9) MOD 10 EQ 0 THEN BEGIN 93 ras = report('island number :'+strtrim(islnumb-9, 1)) 94 ENDIF 94 95 mskneig, msk, nx, ny, land, islnumb 95 96 -
trunk/SRC/Grid/smallmeshmask.pro
r238 r240 78 78 filein_msk = (findfile(filein+'.mask.nc'))[0] 79 79 IF filein_hgr EQ '' OR filein_zgr EQ '' OR filein_msk EQ '' THEN BEGIN 80 ras = report( 'meshmask file(s) not found...')81 ras = report( filein+' does not exist')82 ras = report( filein+'.mesh_hgr.nc does not exist')83 ras = report( filein+'.mesh_zgr.nc does not exist')84 ras = report( filein+'.mask.nc does not exist')80 ras = report(['meshmask file(s) not found...', $ 81 filein+' does not exist', $ 82 filein+'.mesh_hgr.nc does not exist', $ 83 filein+'.mesh_zgr.nc does not exist', $ 84 filein+'.mask.nc does not exist']) 85 85 return 86 86 ENDIF -
trunk/SRC/Interpolation/inquad.pro
r238 r240 127 127 ; 128 128 IF good[0] EQ -1 THEN BEGIN 129 IF NOT keyword_set(noprint) THEN print, 'The point is out of the cells' 129 IF NOT keyword_set(noprint) THEN BEGIN 130 ras = report('The point is out of the cells') 131 ENDIF 130 132 ; restore the input parameters 131 133 x = temporary(xin) … … 168 170 ; 169 171 IF good2[0] EQ -1 THEN BEGIN 170 IF NOT keyword_set(noprint) THEN print, 'The point is out of the cells' 172 IF NOT keyword_set(noprint) THEN BEGIN 173 ras = report('The point is out of the cells') 174 ENDIF 171 175 ; restore the input parameters 172 176 x = temporary(xin) -
trunk/SRC/ReadWrite/read_grads.pro
r236 r240 179 179 infof2read=fstat(unit) 180 180 if infof2read.size NE filesize then begin 181 ras = report( 'According to '+filename+' the file size must be '+strtrim(filesize, 1)+' instead of '+strtrim(infof2read.size, 1))182 ras = report( 'jpi: '+strtrim(jpi, 2))183 ras = report( 'jpj: '+strtrim(jpj, 2))184 ras = report( 'jpt: '+strtrim(jpt, 2))185 ras = report( 'format size in byte: '+strtrim(fmtsz, 2))186 ras = report( 'number of xy arrays: '+strtrim(nxybloc, 2))181 ras = report(['According to '+filename+' the file size must be '+strtrim(filesize, 1)+' instead of '+strtrim(infof2read.size, 1), $ 182 'jpi: '+strtrim(jpi, 2), $ 183 'jpj: '+strtrim(jpj, 2), $ 184 'jpt: '+strtrim(jpt, 2), $ 185 'format size in byte: '+strtrim(fmtsz, 2), $ 186 'number of xy arrays: '+strtrim(nxybloc, 2)]) 187 187 return, -1 188 188 endif -
trunk/SRC/ReadWrite/readoldopadistcoast.pro
r238 r240 103 103 size = numrec*reclen-toomuch 104 104 if size NE filepamameters.size then begin 105 ras = report( 'The size of the file is not the expected one!')106 ras = report( 'Check your file or the values of ibloc, jpiglo,')107 ras = report( 'jpjglo, jpk, jpbyt, numrec in this program')105 ras = report(['The size of the file is not the expected one!', $ 106 'Check your file or the values of ibloc, jpiglo,', $ 107 'jpjglo, jpk, jpbyt, numrec in this program']) 108 108 return, -1 109 109 endif -
trunk/SRC/ReadWrite/readoldoparestart.pro
r236 r240 145 145 iname_file = findfile(filename) 146 146 if iname_file[0] EQ '' then begin 147 ras = report( 147 ras = report('Bad file name') 148 148 return 149 149 ENDIF ELSE iname_file = iname_file[0] … … 165 165 size = numrec*reclen-toomuch 166 166 if size NE filepamameters.size then begin 167 ras = report( 'The size of the file is not the expected one!')168 ras = report( 'Check your file or the values of ibloc, jpiglo,')169 ras = report( 'jpjglo, jpk, jpbyt, numrec in this program')167 ras = report(['The size of the file is not the expected one!', $ 168 'Check your file or the values of ibloc, jpiglo,', $ 169 'jpjglo, jpk, jpbyt, numrec in this program']) 170 170 return 171 171 endif -
trunk/SRC/ReadWrite/scanctl.pro
r238 r240 273 273 files = strsplit(files,/extract) 274 274 if n_elements(files) NE 2 then begin 275 ras = report( 'Bad definition of the filename. There should be 2 elements:')276 ras = report( 'DEST and 1 filename that may define many files')275 ras = report(['Bad definition of the filename. There should be 2 elements:', $ 276 'DEST and 1 filename that may define many files']) 277 277 stop 278 278 endif -
trunk/SRC/ToBeReviewed/CALCULS/grossemoyenne.pro
r237 r240 556 556 if dirt EQ 1 AND keyword_set(spatialfirst) then BEGIN 557 557 IF (reverse(size(res, /dimension)))[0] NE jpt THEN BEGIN 558 print, 'the last dimension of res is not equal to jpt: '+strtrim(jpt, 2)558 ras = report('the last dimension of res is not equal to jpt: '+strtrim(jpt, 2)) 559 559 if keyword_set(savedbox) THEN restoreboxparam, 'boxparam4grmoyenne.dat' 560 560 return, -1 -
trunk/SRC/ToBeReviewed/GRILLE/domdef.pro
r232 r240 248 248 IF (dom[0] EQ -1) THEN BEGIN 249 249 IF keyword_set(findalways) THEN BEGIN 250 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'250 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 251 251 neig1 = neighbor(lon1, lat1, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 252 252 neig2 = neighbor(lon2, lat2, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 284 284 CASE 1 OF 285 285 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 286 print, 'WARNING, empty U points box... we use the same index as T points...'286 ras = report('WARNING, empty U points box... we use the same index as T points...') 287 287 firstxu = firstxt & lastxu = lastxt & nxu = nxt 288 288 firstyu = firstyt & lastyu = lastyt & nyu = nyt 289 289 END 290 290 ELSE:BEGIN 291 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'291 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 292 292 neig1 = neighbor(lon1, lat1, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 293 293 neig2 = neighbor(lon2, lat2, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 327 327 CASE 1 OF 328 328 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 329 print, 'WARNING, empty V points box... we use the same index as T points...'329 ras = report('WARNING, empty V points box... we use the same index as T points...') 330 330 firstxv = firstxt & lastxv = lastxt & nxv = nxt 331 331 firstyv = firstyt & lastyv = lastyt & nyv = nyt 332 332 END 333 333 (where(gridtype eq 'U'))[0] NE -1:BEGIN 334 print, 'WARNING, empty V points box... we use the same index as U points...'334 ras = report('WARNING, empty V points box... we use the same index as U points...') 335 335 firstxv = firstxu & lastxv = lastxu & nxv = nxu 336 336 firstyv = firstyu & lastyv = lastyu & nyv = nyu 337 337 END 338 338 ELSE:BEGIN 339 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'339 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 340 340 neig1 = neighbor(lon1, lat1, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 341 341 neig2 = neighbor(lon2, lat2, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 375 375 CASE 1 OF 376 376 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 377 print, 'WARNING, empty F points box... we use the same index as T points...'377 ras = report('WARNING, empty F points box... we use the same index as T points...') 378 378 firstxf = firstxt & lastxf = lastxt & nxf = nxt 379 379 firstyf = firstyt & lastyf = lastyt & nyf = nyt 380 380 END 381 381 (where(gridtype eq 'U'))[0] NE -1:BEGIN 382 print, 'WARNING, empty F points box... we use the same index as U points...'382 ras = report('WARNING, empty F points box... we use the same index as U points...') 383 383 firstxf = firstxu & lastxf = lastxu & nxf = nxu 384 384 firstyf = firstyu & lastyf = lastyu & nyf = nyu 385 385 END 386 386 (where(gridtype eq 'V'))[0] NE -1:BEGIN 387 print, 'WARNING, empty F points box... we use the same index as V points...'387 ras = report('WARNING, empty F points box... we use the same index as V points...') 388 388 firstxf = firstxv & lastxf = lastxv & nxf = nxv 389 389 firstyf = firstyv & lastyf = lastyv & nyf = nyv 390 390 END 391 391 ELSE:BEGIN 392 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'392 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 393 393 neig1 = neighbor(lon1, lat1, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 394 394 neig2 = neighbor(lon2, lat2, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 497 497 IF (dom[0] EQ -1) THEN BEGIN 498 498 IF keyword_set(findalways) THEN BEGIN 499 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'499 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 500 500 neig1 = neighbor(lon1, lat1, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 501 501 neig2 = neighbor(lon2, lat2, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 528 528 CASE 1 OF 529 529 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 530 print, 'WARNING, empty U points box... we use the same index as T points...'530 ras = report( 'WARNING, empty U points box... we use the same index as T points...') 531 531 firstyu = firstyt & lastyu = lastyt & nyu = nyt 532 532 END 533 533 ELSE:BEGIN 534 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'534 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 535 535 neig1 = neighbor(lon1, lat1, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 536 536 neig2 = neighbor(lon2, lat2, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 566 566 CASE 1 OF 567 567 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 568 print, 'WARNING, empty V points box... we use the same index as T points...'568 ras = report('WARNING, empty V points box... we use the same index as T points...') 569 569 firstyv = firstyt & lastyv = lastyt & nyv = nyt 570 570 END 571 571 (where(gridtype eq 'U'))[0] NE -1:BEGIN 572 print, 'WARNING, empty V points box... we use the same index as U points...'572 ras = report('WARNING, empty V points box... we use the same index as U points...') 573 573 firstyv = firstyu & lastyv = lastyu & nyv = nyu 574 574 END 575 575 ELSE:BEGIN 576 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'576 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 577 577 neig1 = neighbor(lon1, lat1, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 578 578 neig2 = neighbor(lon2, lat2, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 608 608 CASE 1 OF 609 609 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 610 print, 'WARNING, empty F points box... we use the same index as T points...'610 ras = report('WARNING, empty F points box... we use the same index as T points...') 611 611 firstyf = firstyt & lastyf = lastyt & nyf = nyt 612 612 END 613 613 (where(gridtype eq 'U'))[0] NE -1:BEGIN 614 print, 'WARNING, empty F points box... we use the same index as U points...'614 ras = report('WARNING, empty F points box... we use the same index as U points...') 615 615 firstyf = firstyu & lastyf = lastyu & nyf = nyu 616 616 END 617 617 (where(gridtype eq 'V'))[0] NE -1:BEGIN 618 print, 'WARNING, empty F points box... we use the same index as V points...'618 ras = report('WARNING, empty F points box... we use the same index as V points...') 619 619 firstyf = firstyv & lastyf = lastyv & nyf = nyv 620 620 END 621 621 ELSE:BEGIN 622 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'622 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 623 623 neig1 = neighbor(lon1, lat1, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 624 624 neig2 = neighbor(lon2, lat2, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 665 665 IF (dom[0] EQ -1) THEN BEGIN 666 666 IF keyword_set(findalways) THEN BEGIN 667 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'667 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 668 668 neig1 = neighbor(lon1, lat1, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 669 669 neig2 = neighbor(lon2, lat2, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 696 696 CASE 1 OF 697 697 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 698 print, 'WARNING, empty U points box... we use the same index as T points...'698 ras = report('WARNING, empty U points box... we use the same index as T points...') 699 699 firstxu = firstxt & lastxu = lastxt & nxu = nxt 700 700 END 701 701 ELSE:BEGIN 702 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'702 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 703 703 neig1 = neighbor(lon1, lat1, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 704 704 neig2 = neighbor(lon2, lat2, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 734 734 CASE 1 OF 735 735 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 736 print, 'WARNING, empty V points box... we use the same index as T points...'736 ras = report('WARNING, empty V points box... we use the same index as T points...') 737 737 firstxv = firstxt & lastxv = lastxt & nxv = nxt 738 738 END 739 739 (where(gridtype eq 'U'))[0] NE -1:BEGIN 740 print, 'WARNING, empty V points box... we use the same index as U points...'740 ras = report('WARNING, empty V points box... we use the same index as U points...') 741 741 firstxv = firstxu & lastxv = lastxu & nxv = nxu 742 742 END 743 743 ELSE:BEGIN 744 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'744 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 745 745 neig1 = neighbor(lon1, lat1, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 746 746 neig2 = neighbor(lon2, lat2, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 776 776 CASE 1 OF 777 777 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 778 print, 'WARNING, empty F points box... we use the same index as T points...'778 ras = report('WARNING, empty F points box... we use the same index as T points...') 779 779 firstxf = firstxt & lastxf = lastxt & nxf = nxt 780 780 END 781 781 (where(gridtype eq 'U'))[0] NE -1:BEGIN 782 print, 'WARNING, empty F points box... we use the same index as U points...'782 ras = report('WARNING, empty F points box... we use the same index as U points...') 783 783 firstxf = firstxu & lastxf = lastxu & nxf = nxu 784 784 END 785 785 (where(gridtype eq 'V'))[0] NE -1:BEGIN 786 print, 'WARNING, empty F points box... we use the same index as V points...'786 ras = report('WARNING, empty F points box... we use the same index as V points...') 787 787 firstxf = firstxv & lastxf = lastxv & nxf = nxv 788 788 END 789 789 ELSE:BEGIN 790 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'790 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 791 791 neig1 = neighbor(lon1, lat1, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 792 792 neig2 = neighbor(lon2, lat2, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) -
trunk/SRC/ToBeReviewed/HOPE/findlineandpointtype.pro
r231 r240 85 85 if (jpi NE 128 and jpi NE 256) $ 86 86 OR (jpj NE 121 AND jpj NE 121*2 AND jpj NE 195 AND jpj NE 195*2) then begin 87 print, 'CASE NOT coded...'87 ras = report('CASE NOT coded...') 88 88 stop 89 89 ENDIF … … 106 106 END 107 107 ELSE:BEGIN 108 print, 'CASE NOT coded...'108 ras = report('CASE NOT coded...') 109 109 stop 110 110 END … … 134 134 END 135 135 ELSE:BEGIN 136 print, 'CASE NOT coded...'136 ras = report('CASE NOT coded...') 137 137 stop 138 138 END … … 151 151 1.4:vargrid = 'U' 152 152 ELSE:BEGIN 153 print, 'CASE NOT coded...'153 ras = report('CASE NOT coded...') 154 154 stop 155 155 END … … 162 162 2.8:vargrid = 'U' 163 163 ELSE:BEGIN 164 print, 'CASE NOT coded...'164 ras = report('CASE NOT coded...') 165 165 stop 166 166 END … … 168 168 END 169 169 ELSE:BEGIN 170 print, 'CASE NOT coded...'170 ras = report('CASE NOT coded...') 171 171 stop 172 172 END … … 175 175 'xz':BEGIN 176 176 if (jpi NE 128 and jpi NE 256) then begin 177 print, 'CASE NOT coded...'177 ras = report('CASE NOT coded...') 178 178 stop 179 179 ENDIF … … 188 188 test = (where(lat EQ yaxis[0]))[0] 189 189 if test EQ -1 then begin 190 print, 'CASE NOT coded...'190 ras = report('CASE NOT coded...') 191 191 stop 192 192 endif … … 206 206 test = (where(lat EQ yaxis[0]))[0] 207 207 if test EQ -1 then begin 208 print, 'CASE NOT coded...'208 ras = report('CASE NOT coded...') 209 209 stop 210 210 endif … … 230 230 'yz':BEGIN 231 231 if (jpj NE 121 AND jpj NE 195) then begin 232 print, 'CASE NOT coded...'232 ras= report('CASE NOT coded...') 233 233 stop 234 234 ENDIF … … 243 243 (where(lonTodd+xstep EQ xaxis[0]))[0] NE -1:vargrid = 'U' 244 244 ELSE:BEGIN 245 print, 'CASE NOT coded...'245 ras = report('CASE NOT coded...') 246 246 stop 247 247 END … … 256 256 (where(lonTodd-xstep EQ xaxis[0]))[0] NE -1:vargrid = 'T' 257 257 ELSE:BEGIN 258 print, 'CASE NOT coded...'258 ras = report('CASE NOT coded...') 259 259 stop 260 260 END … … 269 269 (where(lonTodd+xstep EQ xaxis[0]))[0] NE -1:vargrid = 'U' 270 270 ELSE:BEGIN 271 print, 'CASE NOT coded...'271 ras = report('CASE NOT coded...') 272 272 stop 273 273 END … … 282 282 (where(lonTodd-xstep EQ xaxis[0]))[0] NE -1:vargrid = 'T' 283 283 ELSE:BEGIN 284 print, 'CASE NOT coded...'284 ras = report('CASE NOT coded...') 285 285 stop 286 286 END … … 288 288 END 289 289 ELSE:BEGIN 290 print, 'CASE NOT coded...'290 ras = report('CASE NOT coded...') 291 291 stop 292 292 END … … 295 295 END 296 296 else:BEGIN 297 print, 'case not coded'297 ras = report('CASE NOT coded...') 298 298 stop 299 299 END -
trunk/SRC/ToBeReviewed/LECTURE/ncdf_lec.pro
r232 r240 65 65 cdfid=ncdf_open(iodir+nom) 66 66 ;------------------------------------------------------------ 67 ; Wh t does the file contain??67 ; What does the file contain?? 68 68 ;------------------------------------------------------------ 69 69 wathinside=ncdf_inquire(cdfid) 70 70 ; print,'in the file, ',iodir+nom,', there are:' 71 71 if keyword_set(dim) then begin 72 print,'nombre de dimensions: ',strtrim(wathinside.ndims,1)73 print,'numero de la dimension dont la valeur est infini: ',strtrim(wathinside.recdim,1)72 ras = report(['number of dimensions: ' + strtrim(wathinside.ndims,1), $ 73 'number of the dimension which value is infinite : ',strtrim(wathinside.recdim,1)]) 74 74 endif 75 75 if keyword_set(var) then $ 76 if size(var, /type) NE 7 then print,'nombre de variables :',strtrim(wathinside.nvars,1)76 if size(var, /type) NE 7 then ras = report('number of variables :'+strtrim(wathinside.nvars,1)) 77 77 if keyword_set(att) then begin 78 78 if strlowcase(att) ne 'global' then goto,nonglobal 79 print,'nombre de attributs globaux :',strtrim(wathinside.ngatts,1)79 ras = report('number of global attributes :' + strtrim(wathinside.ngatts,1)) 80 80 endif 81 81 ;------------------------------------------------------------ … … 88 88 name=ncdf_attname(cdfid,attiq,/global) ;attribute's name 89 89 ncdf_attget,cdfid,name,value,/global ;attribute's value 90 print,name,': ',string(value)90 ras = report(name + ': ' + string(value)) 91 91 endfor 92 92 endif … … 106 106 tailledim[dimiq]=value 107 107 if keyword_set(dim) then begin 108 print,'dimension numero ',strtrim(dimiq,1),', nom: ',nomdim[dimiq]$109 ,', valeur: ' ,strtrim(tailledim[dimiq],1)108 ras = report('dimension number ' + strtrim(dimiq,1) + ', name: ' + $ 109 nomdim[dimiq] + ' value: ' + strtrim(tailledim[dimiq],1)) 110 110 endif 111 111 endfor … … 130 130 if strlowcase(att) eq strlowcase(varcontent.name) or keyword_set(var) $ 131 131 then begin 132 print,'variable numero: ',strtrim(varid,1),', nom:',varcontent.name $ 133 ,', type:' ,varcontent.datatype,', dimensions:',nomdim[varcontent.dim] 132 ras = report('variable number: ' + strtrim(varid,1) + $ 133 ', name:' + varcontent.name + $ 134 ', type:' + varcontent.datatype + $ 135 ', dimensions:' + nomdim[varcontent.dim]) 134 136 if strlowcase(att) eq strlowcase(varcontent.name) then begin 135 137 for attiq=0,varcontent.natts-1 do begin 136 138 name=ncdf_attname(cdfid,varid,attiq) 137 139 ncdf_attget,cdfid,varid,name,value 138 print,' ',strtrim(attiq),' ',name,': ',strtrim(string(value),1) 140 ras = report(' ' + strtrim(attiq) + ' ' + $ 141 ' name : ' + strtrim(string(value),1)) 139 142 endfor 140 143 goto, sortie -
trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro
r238 r240 486 486 missing_value = strmid(missing_value, 0, strlen(missing_value)-1) 487 487 IF isnumber(string(missing_value), tmp) EQ 1 THEN missing_value = tmp ELSE BEGIN 488 print, 'Warning: missing value is not a number: ', missing_value488 ras = report('Warning: missing value is not a number: ' + missing_value) 489 489 missing_value = - 1 490 490 ENDELSE -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/bar_plot.pro
r232 r240 141 141 ; 142 142 if (n_params(d) eq 0) then begin ;Print call & return if no parameters 143 print,'bar_test,values,baselines=baselines,colors=colors,barnames=barnames,$'144 print,' title=title,xtitle=xtitle,ytitle=ytitle,baserange=baserange, $'145 print,' barwidth=barwidth,barspace=barspace,baroffset=baroffset, $'146 print,' outline=outline,overplot=overplot,background=background, $'147 print,' rotate=rotate'143 ras = report(['bar_test,values,baselines=baselines,colors=colors,barnames=barnames,$', $ 144 ' title=title,xtitle=xtitle,ytitle=ytitle,baserange=baserange, $', $ 145 ' barwidth=barwidth,barspace=barspace,baroffset=baroffset, $', $ 146 ' outline=outline,overplot=overplot,background=background, $', $ 147 ' rotate=rotate']) 148 148 return 149 149 endif -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltt.pro
r232 r240 277 277 278 278 IF type EQ 't' AND keyword_set(filter) THEN BEGIN 279 print, ' Applying a running mean filter of width '+string(filter, format = '(I3)')279 ras = report(' Applying a running mean filter of width '+string(filter, format = '(I3)')) 280 280 z2d = smooth(z2d, filter) 281 281 z2d[0:filter/2-1] = 0. -
trunk/SRC/ToBeReviewed/PLOTS/LABEL/label_gmt.pro
r231 r240 59 59 max=level_z2d[Ncontour-1]+intervalle 60 60 61 print, ' Number of contour intervals, plotting min & max ', ncontour, min, max62 print, ' ' 61 ras = report([' Number of contour intervals, plotting min & max ' + $ 62 string(ncontour) + string(min)+ string(max),' ']) 63 63 64 64 ; color index -
trunk/SRC/ToBeReviewed/POSTSCRIPT/chcolps.pro
r231 r240 162 162 strput, ligne, color, pos-18 163 163 ENDIF ELSE BEGIN 164 ; print, 'erreur ligne :', nl164 ; ras = report('erreur ligne :' + string(nl)) 165 165 dist = abs(r-red)+abs(g-green)+abs(b-blue) 166 166 ind = (where(dist EQ min(dist)))[0] 167 167 ind = ind[0] 168 ; print, 'je trouve ', long([r, g, b]) 169 ; print, 'je remplace par ', [red[ind], green[ind], blue[ind]] 168 ; ras = report( $ 169 ; ['I found ' + string(long([r, g, b])), $ 170 ; 'I replace it by ' + string([red[ind], green[ind], blue[ind]])]) 170 171 r1 = red1[ind]/255. 171 172 g1 = green1[ind]/255. … … 191 192 ENDIF 192 193 ; 193 ; Ecrit le fichier de sorti 194 ; Ecrit le fichier de sortie 194 195 ; 195 196 printf, numout, ligne, format = '(A)' -
trunk/SRC/ToBeReviewed/POSTSCRIPT/ps.pro
r238 r240 23 23 ; We have to close the journal! 24 24 ;------------------------------------------------------------ 25 if !journal EQ 0 then print, 'ERROR: The Journal is not open !!'25 if !journal EQ 0 then ras = report('ERROR: The Journal is not open !!') 26 26 journal 27 27 ;------------------------------------------------------------ -
trunk/SRC/ToBeReviewed/STRING/delchr.pro
r232 r240 40 40 41 41 if (n_params(0) lt 2) or keyword_set(hlp) then begin 42 print,' Delete all occurrences of a character from a text string.'43 print,' new = delchr(old, char)'44 print,' old = original text string. in'45 print,' char = character to delete. in'46 print,' new = resulting string. out'42 ras = report(['Delete all occurrences of a character from a text string.', $ 43 ' new = delchr(old, char)', $ 44 ' old = original text string. in' , $ 45 ' char = character to delete. in' , $ 46 ' new = resulting string. out']) 47 47 return, -1 48 48 endif -
trunk/SRC/ToBeReviewed/STRING/getfile.pro
r237 r240 58 58 59 59 if (n_params(0) lt 1) or keyword_set(hlp) then begin 60 print,' Read a text file into a string array.' 61 print,' s = getfile(f)' 62 print,' f = text file name. in' 63 print,' s = string array. out' 64 print,' Keywords:' 65 print,' ERROR=err error flag: 0=ok, 1=file not opened,' 66 print,' 2=no lines in file.' 67 print,' /QUIET means give no error message.' 68 print,' LINES=n Number of lines to read (def=all).' 69 print,' Much faster if number of lines is known.' 70 print,' Automatic for IDL 5.6 or later.' 60 ras = report([ $ 61 ' Read a text file into a string array.', $ 62 ' s = getfile(f)', $ 63 ' f = text file name. in', $ 64 ' s = string array. out', $ 65 ' Keywords:', $ 66 ' ERROR=err error flag: 0=ok, 1=file not opened,', $ 67 ' 2=no lines in file.', $ 68 ' /QUIET means give no error message.', $ 69 ' LINES=n Number of lines to read (def=all).', $ 70 ' Much faster if number of lines is known.', $ 71 ' Automatic for IDL 5.6 or later.']) 71 72 return, -1 72 73 endif -
trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/inserthistory.pro
r231 r240 42 42 index1 = where(globalcommand EQ line1) & index1 = index1[0]+1 43 43 if index1 EQ -1 then begin 44 rien = report(line1+' n ''est pas trouve dsglobalcommand')44 rien = report(line1+' not found in globalcommand') 45 45 return 46 46 endif 47 47 index2 = where(globalcommand EQ line2) & index2 = index2[0] 48 48 if index2 EQ -1 then begin 49 rien = report(line2+' n ''est pas trouve dsglobalcommand')49 rien = report(line2+' not found in globalcommand') 50 50 return 51 51 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.