- Timestamp:
- 03/21/07 18:18:44 (17 years ago)
- Location:
- trunk/SRC
- Files:
-
- 26 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Documentation/idldoc_assistant_output/Obsolete/jourdsmois.html
r234 r236 62 62 63 63 64 used <a href=". //../../Calendar/daysinmonth.html">daysinmonth</a> instead64 used <a href="../Calendar/daysinmonth.html">daysinmonth</a> instead 65 65 66 66 -
trunk/SRC/Grid/checkperio.pro
r232 r236 47 47 48 48 IF file_test(file) EQ 0 THEN BEGIN 49 print, 'file '+file+' not found'49 ras = report( 'file '+file+' not found') 50 50 RETURN 51 51 ENDIF … … 78 78 part2 = reverse(shift(part2, -1)) 79 79 IF array_equal(part1, part2) NE 1 THEN BEGIN 80 print, vinq.name+' bad y periodicity (1)'80 ras = report( vinq.name+' bad y periodicity (1)') 81 81 IF keyword_set(write) THEN $ 82 82 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1], offset = [1, jpj-1] … … 88 88 part2 = reverse(part2) 89 89 IF array_equal(part1, part2) NE 1 THEN BEGIN 90 print, vinq.name+' bad y periodicity (1)'90 ras = report( vinq.name+' bad y periodicity (1)') 91 91 IF keyword_set(write) THEN $ 92 92 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1], offset = [1, jpj-1] 93 93 ENDIF 94 94 END 95 ELSE: print, vinq.name+' nothing to check'95 ELSE:ras = report( vinq.name+' nothing to check') 96 96 ENDCASE 97 97 ; y periodicity (2) … … 102 102 part2 = reverse(part2) 103 103 IF array_equal(part1, part2) NE 1 THEN BEGIN 104 print, vinq.name+' bad y periodicity (2)'104 ras = report( vinq.name+' bad y periodicity (2)') 105 105 IF keyword_set(write) THEN $ 106 106 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1], offset = [(jpi-1)/2+2, jpj-2] … … 109 109 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 110 110 END 111 ELSE: print, vinq.name+' nothing to check'111 ELSE:ras = report( vinq.name+' nothing to check') 112 112 ENDCASE 113 113 ; x periodicity (1) … … 115 115 ncdf_varget, cdfid, vid, part2, count = [1, jpj], offset = [jpi-2, 0] 116 116 IF array_equal(part1, part2) NE 1 THEN BEGIN 117 print, vinq.name+' bad x periodicity (1)'117 ras = report( vinq.name+' bad x periodicity (1)') 118 118 IF keyword_set(write) THEN $ 119 119 ncdf_varput, cdfid, vid, part2, count = [1, jpj], offset = [0, 0] … … 123 123 ncdf_varget, cdfid, vid, part2, count = [1, jpj], offset = [1, 0] 124 124 IF array_equal(part1, part2) NE 1 THEN BEGIN 125 print, vinq.name+' bad x periodicity (2)'125 ras = report( vinq.name+' bad x periodicity (2)') 126 126 IF keyword_set(write) THEN $ 127 127 ncdf_varput, cdfid, vid, part2, count = [1, jpj], offset = [jpi-1, 0] … … 141 141 part2 = reverse(shift(part2, -1, 0, 0), 1) 142 142 IF array_equal(part1, part2) NE 1 THEN BEGIN 143 print, vinq.name+' bad y periodicity (1)'143 ras = report( vinq.name+' bad y periodicity (1)') 144 144 IF keyword_set(write) THEN $ 145 145 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk], offset = [1, jpj-1, 0] … … 152 152 part2 = reverse(part2, 1) 153 153 IF array_equal(part1, part2) NE 1 THEN BEGIN 154 print, vinq.name+' bad y periodicity (1)'154 ras = report( vinq.name+' bad y periodicity (1)') 155 155 IF keyword_set(write) THEN $ 156 156 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk], offset = [1, jpj-1, 0] 157 157 ENDIF 158 158 END 159 ELSE: print, vinq.name+' nothing to check'159 ELSE:ras = report( vinq.name+' nothing to check') 160 160 ENDCASE 161 161 ; y periodicity (2) … … 166 166 part2 = reverse(part2, 1) 167 167 IF array_equal(part1, part2) NE 1 THEN BEGIN 168 print, vinq.name+' bad y periodicity (2)'168 ras = report( vinq.name+' bad y periodicity (2)') 169 169 IF keyword_set(write) THEN $ 170 170 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1, jpk], offset = [(jpi-1)/2+2, jpj-2, 0] … … 173 173 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 174 174 END 175 ELSE: print, vinq.name+' nothing to check'175 ELSE:ras = report( vinq.name+' nothing to check') 176 176 ENDCASE 177 177 ; x periodicity (1) … … 179 179 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk], offset = [jpi-2, 0, 0] 180 180 IF array_equal(part1, part2) NE 1 THEN BEGIN 181 print, vinq.name+' bad x periodicity (1)'181 ras = report( vinq.name+' bad x periodicity (1)') 182 182 IF keyword_set(write) THEN $ 183 183 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk], offset = [0, 0, 0] … … 187 187 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk], offset = [1, 0, 0] 188 188 IF array_equal(part1, part2) NE 1 THEN BEGIN 189 print, vinq.name+' bad x periodicity (2)'189 ras = report( vinq.name+' bad x periodicity (2)') 190 190 IF keyword_set(write) THEN $ 191 191 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk], offset = [jpi-1, 0, 0] … … 206 206 part2 = reverse(shift(part2, -1, 0, 0, 0), 1) 207 207 IF array_equal(part1, part2) NE 1 THEN BEGIN 208 print, vinq.name+' bad y periodicity (1)'208 ras = report( vinq.name+' bad y periodicity (1)') 209 209 IF keyword_set(write) THEN $ 210 210 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk, jpt], offset = [1, jpj-1, 0, 0] … … 217 217 part2 = reverse(part2, 1) 218 218 IF array_equal(part1, part2) NE 1 THEN BEGIN 219 print, vinq.name+' bad y periodicity (1)'219 ras = report( vinq.name+' bad y periodicity (1)') 220 220 IF keyword_set(write) THEN $ 221 221 ncdf_varput, cdfid, vid, part2, count = [jpi-2, 1, jpk, jpt], offset = [1, jpj-1, 0, 0] 222 222 ENDIF 223 223 END 224 ELSE: print, vinq.name+' nothing to check'224 ELSE:ras = report( vinq.name+' nothing to check') 225 225 ENDCASE 226 226 ; y periodicity (2) … … 231 231 part2 = reverse(part2, 1) 232 232 IF array_equal(part1, part2) NE 1 THEN BEGIN 233 print, vinq.name+' bad y periodicity (2)'233 ras = report( vinq.name+' bad y periodicity (2)') 234 234 IF keyword_set(write) THEN $ 235 235 ncdf_varput, cdfid, vid, part2, count = [(jpi-1)/2-2+1, 1, jpk, jpt], offset = [(jpi-1)/2+2, jpj-2, 0, 0] … … 238 238 array_equal([722, 511], vdimsz[0:1]):BEGIN ; ORCA05 239 239 END 240 ELSE: print, vinq.name+' nothing to check'240 ELSE:ras = report( vinq.name+' nothing to check') 241 241 ENDCASE 242 242 ; x periodicity (1) … … 244 244 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [jpi-2, 0, 0, 0] 245 245 IF array_equal(part1, part2) NE 1 THEN BEGIN 246 print, vinq.name+' bad x periodicity (1)'246 ras = report( vinq.name+' bad x periodicity (1)') 247 247 IF keyword_set(write) THEN $ 248 248 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [0, 0, 0, 0] … … 252 252 ncdf_varget, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [1, 0, 0, 0] 253 253 IF array_equal(part1, part2) NE 1 THEN BEGIN 254 print, vinq.name+' bad x periodicity (2)'254 ras = report( vinq.name+' bad x periodicity (2)') 255 255 IF keyword_set(write) THEN $ 256 256 ncdf_varput, cdfid, vid, part2, count = [1, jpj, jpk, jpt], offset = [jpi-1, 0, 0, 0] 257 257 ENDIF 258 258 END 259 ELSE: print, vinq.name+' nothing to check'259 ELSE:ras = report( vinq.name+' nothing to check') 260 260 ENDCASE 261 261 ENDIF -
trunk/SRC/Grid/computegrid.pro
r232 r236 330 330 1:nz = (size(zaxis))[1] 331 331 ELSE:BEGIN 332 print, 'not coded'332 ras = report( 'not coded') 333 333 stop 334 334 END -
trunk/SRC/Grid/micromeshmask.pro
r231 r236 64 64 filein_msk = (findfile(filein+'.mask.nc'))[0] 65 65 IF filein_hgr EQ '' OR filein_zgr EQ '' OR filein_msk EQ '' THEN BEGIN 66 print, 'meshmask file(s) not found...'67 print, filein+' does not exist'68 print, filein+'.mesh_hgr.nc does not exist'69 print, filein+'.mesh_zgr.nc does not exist'70 print, filein+'.mask.nc does not exist'66 ras = report( 'meshmask file(s) not found...') 67 ras = report( filein+' does not exist') 68 ras = report( filein+'.mesh_hgr.nc does not exist') 69 ras = report( filein+'.mesh_zgr.nc does not exist') 70 ras = report( filein+'.mask.nc does not exist') 71 71 return 72 72 ENDIF -
trunk/SRC/Grid/smallmeshmask.pro
r231 r236 72 72 filein_msk = (findfile(filein+'.mask.nc'))[0] 73 73 IF filein_hgr EQ '' OR filein_zgr EQ '' OR filein_msk EQ '' THEN BEGIN 74 print, 'meshmask file(s) not found...'75 print, filein+' does not exist'76 print, filein+'.mesh_hgr.nc does not exist'77 print, filein+'.mesh_zgr.nc does not exist'78 print, filein+'.mask.nc does not exist'74 ras = report( 'meshmask file(s) not found...') 75 ras = report( filein+' does not exist') 76 ras = report( filein+'.mesh_hgr.nc does not exist') 77 ras = report( filein+'.mesh_zgr.nc does not exist') 78 ras = report( filein+'.mask.nc does not exist') 79 79 return 80 80 ENDIF -
trunk/SRC/Interpolation/angle.pro
r232 r236 111 111 fileoce = findfile(fileoce, count = okfile) 112 112 IF okfile NE 1 THEN BEGIN 113 print, 'the file '+fileoce+' is not found... we stop'113 ras = report('the file '+fileoce+' is not found... we stop') 114 114 stop 115 115 ENDIF -
trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro
r232 r236 72 72 IF n_elements(amsk) EQ 1 AND amsk[0] EQ -1 THEN amsk = replicate(1b, jpia, jpja) 73 73 IF n_elements(omsk) NE jpio*jpjo THEN BEGIN 74 print, 'input grid mask do not have the good size'74 ras = report('input grid mask do not have the good size') 75 75 stop 76 76 ENDIF 77 77 IF n_elements(amsk) NE jpia*jpja THEN BEGIN 78 print, 'output grid mask do not have the good size'78 ras = report('output grid mask do not have the good size') 79 79 stop 80 80 ENDIF -
trunk/SRC/Interpolation/extrapolate.pro
r232 r236 75 75 IF n_elements(maskinput) EQ 1 AND maskinput[0] EQ -1 THEN maskinput = replicate(1b, nx, ny) 76 76 IF n_elements(maskinput) NE nx*ny THEN BEGIN 77 print, 'input grid mask do not have the good size'77 ras = report('input grid mask do not have the good size') 78 78 return, -1 79 79 ENDIF -
trunk/SRC/Interpolation/extrapsmooth.pro
r232 r236 52 52 sz = size(reform(in)) 53 53 IF sz[0] NE 2 THEN BEGIN 54 print, 'Input arrays must have 2 dimensions'54 ras = report('Input arrays must have 2 dimensions') 55 55 return, -1 56 56 ENDIF … … 59 59 IF n_elements(mskin) EQ 1 AND mskin[0] EQ -1 THEN mskin = replicate(1b, nx, ny) 60 60 IF n_elements(mskin) NE nx*ny THEN BEGIN 61 print, 'input grid mask do not have the good size'61 ras = report('input grid mask do not have the good size') 62 62 return, -1 63 63 ENDIF -
trunk/SRC/Interpolation/fromirr.pro
r231 r236 110 110 'bilinear':compute_fromirr_bilinear_weigaddr, alon, alat, mskin, olon, olat, mskout, weig, addr 111 111 ELSE:BEGIN 112 print, ' unknown interpolation method... we stop'112 ras = report(' unknown interpolation method... we stop') 113 113 stop 114 114 ENDELSE -
trunk/SRC/Interpolation/fromreg.pro
r231 r236 113 113 'imoms3': compute_fromreg_imoms3_weigaddr, alon, alat, olon, olat, weig, addr, NONORTHERNLINE = nonorthernline, NOSOUTHERNLINE = nosouthernline 114 114 ELSE:BEGIN 115 print, ' unknown interpolation method... we stop'115 ras = report(' unknown interpolation method... we stop') 116 116 stop 117 117 ENDELSE -
trunk/SRC/Interpolation/get_gridparams.pro
r231 r236 87 87 ; get longitude and latitude 88 88 IF file_test(in1) EQ 0 THEN BEGIN 89 print, 'file ' + in1 + ' does not exist'89 ras = report('file ' + in1 + ' does not exist') 90 90 stop 91 91 ENDIF … … 102 102 END 103 103 ELSE:BEGIN 104 print, 'Bad nimber of input parameters'104 ras = report('Bad number of input parameters') 105 105 stop 106 106 end … … 140 140 1:BEGIN 141 141 IF array_equal(lon, lon[*, 0] # replicate(1, jpj)) NE 1 THEN BEGIN 142 print, 'Longitudes are not the same for all latitudes, impossible to extract a 1D array of the longitudes'142 ras = report('Longitudes are not the same for all latitudes, impossible to extract a 1D array of the longitudes') 143 143 stop 144 144 ENDIF … … 161 161 1:BEGIN 162 162 IF array_equal(lat, replicate(1, jpi) # lat[0, *]) NE 1 THEN BEGIN 163 print, 'Latitudes are not the same for all longitudes, impossible to extract a 1D array of the latitudes'163 ras = report('Latitudes are not the same for all longitudes, impossible to extract a 1D array of the latitudes') 164 164 stop 165 165 ENDIF … … 182 182 1:BEGIN 183 183 IF array_equal(lon, lon[*, 0] # replicate(1, jpj)) NE 1 THEN BEGIN 184 print, 'Longitudes are not the same for all latitudes, impossible to extract a 1D array of the longitudes'184 ras = report('Longitudes are not the same for all latitudes, impossible to extract a 1D array of the longitudes') 185 185 stop 186 186 ENDIF 187 187 lon = lon[*, 0] 188 188 IF array_equal(lat, replicate(1, jpi) # reform(lat[0, *])) NE 1 THEN BEGIN 189 print, 'Latitudes are not the same for all longitudes, impossible to extract a 1D array of the latitudes'189 ras = report('Latitudes are not the same for all longitudes, impossible to extract a 1D array of the latitudes') 190 190 stop 191 191 ENDIF … … 199 199 ; lon and lat are not 1D and/or 2D arrays 200 200 ;------- 201 ELSE:stop 201 ELSE: BEGIN 202 ras = report('Longitudes and latitudes are not 1D and/or 2D arrays') 203 stop 204 END 202 205 ENDCASE 203 206 ; -
trunk/SRC/Interpolation/inquad.pro
r232 r236 261 261 ; points out of the cells 262 262 IF (where(chtest EQ 0))[0] NE -1 THEN BEGIN 263 IF NOT keyword_set(noprint) THEN print, 'Points number '+strjoin(strtrim(where(chtest EQ 0), 1), ', ')+' are out of the grid' 263 IF NOT keyword_set(noprint) THEN BEGIN 264 ras = report('Points number '+strjoin(strtrim(where(chtest EQ 0), 1), ', ')+' are out of the grid') 265 ENDIF 264 266 stop 265 267 ENDIF 266 268 ; points in more than one cell 267 269 IF (where(chtest GT 1))[0] NE -1 THEN BEGIN 268 IF NOT keyword_set(noprint) THEN print, 'Points number '+strjoin(strtrim(where(chtest GT 1), 1), ', ')+' are in more than one cell' 270 IF NOT keyword_set(noprint) THEN BEGIN 271 ras = report('Points number '+strjoin(strtrim(where(chtest GT 1), 1), ', ')+' are in more than one cell') 272 ENDIF 269 273 stop 270 274 ENDIF … … 277 281 CASE 1 OF 278 282 found[0] EQ -1:BEGIN 279 IF NOT keyword_set(noprint) THEN print, 'The point is out of the cells' 283 IF NOT keyword_set(noprint) THEN BEGIN 284 ras = report('The point is out of the cells') 285 ENDIF 280 286 IF keyword_set(onsphere) THEN BEGIN 281 287 ; restore old !map... … … 289 295 END 290 296 n_elements(found) GT ntofind:BEGIN 291 IF NOT keyword_set(noprint) THEN print, 'The point is in more than one cell' 297 IF NOT keyword_set(noprint) THEN BEGIN 298 ras = report('The point is in more than one cell') 299 ENDIF 292 300 END 293 301 ELSE: -
trunk/SRC/Picture/imdisp.pro
r232 r236 197 197 ; @file_comments 198 198 ; Display an image on the current graphics device. 199 ; IMDISP is an advanced replacement for TV and TVSCL. 199 ; IMDISP is an advanced replacement for <prodil>TV</proidl> and 200 ; <proidl>TVSCL</proidl>. 200 201 ; 201 202 ; - Supports WIN, MAC, X, CGM, PCL, PRINTER, PS, and Z graphics devices, -
trunk/SRC/ReadWrite/read_grads.pro
r232 r236 97 97 if keyword_set(timestep) then BEGIN 98 98 if date1 GT date2 then begin 99 print, 'date2 must be larger than date1'99 ras = report( 'date2 must be larger than date1') 100 100 return, -1 101 101 endif … … 106 106 jdate2 = time[0] > date2jul(date2, /grads) 107 107 if jdate1 GT jdate2 then begin 108 print, 'date2 must be larger than date1'108 ras = report('date2 must be larger than date1') 109 109 return, -1 110 110 endif … … 114 114 ENDELSE 115 115 if t2 LT t1 then begin 116 print, 'There is no date between date1 and date2'116 ras = report('There is no date between date1 and date2') 117 117 return, -1 118 118 endif … … 126 126 varid = varid[0] 127 127 if varid EQ -1 then begin 128 print, var+' not found in the variable list of '+filename128 ras = report(var+' not found in the variable list of '+filename) 129 129 return, -1 130 130 ENDIF … … 159 159 , swap_endian = swapbytes 160 160 if err ne 0 then begin 161 print,!err_string161 ras = report(!err_string) 162 162 return, -1 163 163 endif … … 179 179 infof2read=fstat(unit) 180 180 if infof2read.size NE filesize then begin 181 print, 'According to '+filename+' the file size must be '+strtrim(filesize, 1)+' instead of '+strtrim(infof2read.size, 1)182 print, 'jpi: '+strtrim(jpi, 2)183 print, 'jpj: '+strtrim(jpj, 2)184 print, 'jpt: '+strtrim(jpt, 2)185 print, 'format size in byte: '+strtrim(fmtsz, 2)186 print, '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 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)) 187 187 return, -1 188 188 endif -
trunk/SRC/ReadWrite/read_oasis.pro
r232 r236 45 45 , error=err 46 46 if err ne 0 then begin 47 print,!err_string47 ras = report(!err_string) 48 48 return, -1 49 49 endif … … 57 57 readu, unit 58 58 if EOF(unit) then begin 59 print, varname+' not found in '+filename59 ras = report(varname+' not found in '+filename) 60 60 return, -1 61 61 endif -
trunk/SRC/ReadWrite/readbat.pro
r232 r236 32 32 iname_file = findfile(filename) 33 33 if iname_file[0] EQ '' then begin 34 print, 'Bad file name'34 ras = report( 'Bad file name') 35 35 return, -1 36 36 ENDIF ELSE iname_file = iname_file[0] -
trunk/SRC/ReadWrite/readoldopadistcoast.pro
r232 r236 80 80 iname_file = findfile(filename) 81 81 if iname_file[0] EQ '' then begin 82 print, 'Bad file name'82 ras = report( 'Bad file name') 83 83 return, -1 84 84 ENDIF ELSE iname_file = iname_file[0] … … 103 103 size = numrec*reclen-toomuch 104 104 if size NE filepamameters.size then begin 105 print, 'The size of the file is not the expected one!'106 print, 'Check your file or the values of ibloc, jpiglo,'107 print, 'jpjglo, jpk, jpbyt, numrec in this program'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') 108 108 return, -1 109 109 endif … … 116 116 readu, numcost, iimlu, ijmlu, ikmlu 117 117 if iimlu NE jpiglo then begin 118 print, 'iimlu = '+strtrim(iimlu, 1)+' differs from jpiglo ='+strtrim(jpiglo, 1)118 ras = report( 'iimlu = '+strtrim(iimlu, 1)+' differs from jpiglo ='+strtrim(jpiglo, 1)) 119 119 return, -1 120 120 endif 121 121 if ijmlu NE jpjglo then begin 122 print, 'ijmlu = '+strtrim(ijmlu, 1)+' differs from jpjglo ='+strtrim(jpjglo, 1)122 ras = report( 'ijmlu = '+strtrim(ijmlu, 1)+' differs from jpjglo ='+strtrim(jpjglo, 1)) 123 123 return, -1 124 124 endif 125 125 if ikmlu NE jpk then begin 126 print, 'ikmlu = '+strtrim(ikmlu, 1)+' differs from jpk ='+strtrim(jpk, 1)126 ras = report( 'ikmlu = '+strtrim(ikmlu, 1)+' differs from jpk ='+strtrim(jpk, 1)) 127 127 return, -1 128 128 endif -
trunk/SRC/ReadWrite/readoldoparestart.pro
r231 r236 145 145 iname_file = findfile(filename) 146 146 if iname_file[0] EQ '' then begin 147 print, 'Bad file name'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 print, 'The size of the file is not the expected one!'168 print, 'Check your file or the values of ibloc, jpiglo,'169 print, 'jpjglo, jpk, jpbyt, numrec in this program'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') 170 170 return 171 171 endif -
trunk/SRC/ReadWrite/scanctl.pro
r232 r236 54 54 spawn, '\grep -i "^DTYPE" '+filename, notgood 55 55 if keyword_set(notgood) then begin 56 print, 'This program is not adapted to data type station or grib. Sorry...'56 ras = report( 'This program is not adapted to data type station or grib. Sorry...') 57 57 stop 58 58 endif -
trunk/SRC/ReadWrite/scanoasis.pro
r232 r236 28 28 , error=err 29 29 if err ne 0 then begin 30 print,!err_string30 ras = report(!err_string) 31 31 return 32 32 endif … … 35 35 WHILE NOT EOF(unit) DO BEGIN 36 36 readu, unit, char8 37 print, char837 ras = report(char8) 38 38 readu, unit 39 39 ENDWHILE -
trunk/SRC/ReadWrite/write_oasis.pro
r231 r236 39 39 , error=err, APPEND = append 40 40 if err ne 0 then begin 41 print,!err_string41 ras = report(!err_string) 42 42 return 43 43 endif -
trunk/SRC/ReadWrite/writebat.pro
r231 r236 26 26 ; basic checks 27 27 IF n_params() NE 2 THEN BEGIN 28 print, 'bad number of arguments in the call of writebat'28 ras = report( 'bad number of arguments in the call of writebat') 29 29 return 30 30 ENDIF 31 31 IF size(filename, /type) NE 7 THEN BEGIN 32 print, 'the filename should be a string'32 ras = report( 'the filename should be a string') 33 33 return 34 34 ENDIF 35 35 sbat = size(zbat) 36 36 IF sbat[0] NE 2 THEN BEGIN 37 print, 'bathymetry array should be 2d array'37 ras = report( 'bathymetry array should be 2d array') 38 38 return 39 39 ENDIF -
trunk/SRC/ToBeReviewed/INIT/initncdf.pro
r232 r236 90 90 ncfile = isafile(FILENAME = ncfilein, IODIRECTORY = iodir, _extra = ex) 91 91 if size(ncfile, /type) NE 7 then BEGIN 92 print, 'initncdf cancelled'92 ras = report( 'initncdf cancelled') 93 93 return 94 94 endif … … 139 139 zvarid = (where(namevar EQ 'nav_lev' or namevar EQ zaxisname OR namevar EQ 'level' OR namevar EQ 'lev' OR strmid(namevar, 0, 5) EQ 'depth'))[0] 140 140 if zvarid EQ -1 AND inside.ndims GT 3 then begin 141 print, 'initncdf: the zaxis was not found..., check the the use of ZAXISNAME keyword if you whant to find one...'141 ras = report( 'initncdf: the zaxis was not found..., check the the use of ZAXISNAME keyword if you whant to find one...') 142 142 ; stop 143 143 endif … … 230 230 tmask = finite(tmask) 231 231 IF min(tmask) EQ 1 THEN BEGIN 232 print, 'missing or nan values not found...'232 ras = report( 'missing or nan values not found...') 233 233 tmask = -1 234 234 ENDIF -
trunk/SRC/ToBeReviewed/LECTURE/GRIB/read_grib.pro
r232 r236 41 41 openr, num, gribfile, /GET_LUN, ERROR = err, /SWAP_IF_LITTLE_ENDIAN 42 42 if err ne 0 then begin 43 print, !err_string43 ras = report(!err_string) 44 44 return, -1 45 45 ENDIF … … 64 64 goodvar = where(codes EQ varcode) 65 65 IF goodvar[0] EQ -1 THEN BEGIN 66 print, 'no var code '+strtrim(varcode, 2)+' in the file'66 ras = report( 'no var code '+strtrim(varcode, 2)+' in the file') 67 67 return, -1 68 68 ENDIF … … 73 73 gooddate = where(dates GE date1 AND dates LE date2) 74 74 IF gooddate[0] EQ -1 THEN BEGIN 75 print, 'no dates between '+strtrim(date1, 2)+' and '+strtrim(date2, 2)+' in the file'75 ras = report( 'no dates between '+strtrim(date1, 2)+' and '+strtrim(date2, 2)+' in the file') 76 76 return, -1 77 77 ENDIF -
trunk/SRC/ToBeReviewed/PLOTS/DESSINE/tvplus.pro
r232 r236 2 2 ; 3 3 ; @file_comments 4 ; Enhanced version of tvscl4 ; Enhanced version of <proidl>tvscl</proidl> 5 5 ; 6 6 ; @categories … … 62 62 ; 63 63 ; @keyword _EXTRA 64 ; Used to pass keywords to <proidl>TV</proidl>, 65 ; <proidl>PLOT</proidl>, 64 ; Used to pass keywords to <proidl>TV</proidl>, 65 ; <proidl>PLOT</proidl>, 66 66 ; <pro>COLORBAR</pro> 67 67 ;
Note: See TracChangeset
for help on using the changeset viewer.