Changeset 236 for trunk/SRC/Interpolation/get_gridparams.pro
- Timestamp:
- 03/21/07 18:18:44 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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 ;
Note: See TracChangeset
for help on using the changeset viewer.