Changeset 204 for trunk/src/d2m_to_q2m_erai.pro
- Timestamp:
- 03/26/14 11:58:56 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/d2m_to_q2m_erai.pro
r203 r204 121 121 ; % Program caused arithmetic error: Floating illegal operand 122 122 ; 123 ; related to occur ence of NaNf in ncdump output ? do no yet123 ; related to occurrence of NaNf in ncdump output ? do no yet 124 124 ; 125 125 ; pk also see these error messages ... 126 126 ; 127 ; check for NaN value in P return by reading Netcdfmsl file128 ; solution is may be to apply *mskout to q2m like 127 ; check for NaN value in P return by reading netCDF msl file 128 ; solution is may be to apply *mskout to q2m like in interp* tools 129 129 ; not yet done to keep close as possible to pk processing 130 130 ; … … 198 198 nparam = N_PARAMS() 199 199 IF (nparam NE 2) THEN BEGIN 200 ras = report(['Incorrect number of arguments.' $201 202 203 200 ras = report(['Incorrect number of arguments.' $ 201 + '!C' $ 202 + 'Usage : ' + usage]) 203 return, result 204 204 ENDIF 205 205 ; 206 206 ; test if ${PROJECT_OD} defined 207 207 CASE project_od_env OF 208 '' : BEGIN209 210 211 return, result212 213 ELSE: BEGIN214 215 216 217 208 '' : BEGIN 209 msg = 'eee : ${PROJECT_OD} is not defined' 210 ras = report(msg) 211 return, result 212 END 213 ELSE : BEGIN 214 msg = 'iii : ${PROJECT_OD} is ' + project_od_env 215 ras = report(msg) 216 END 217 ENDCASE 218 218 ; 219 219 ; check if output data will be possible … … 222 222 ; existence and protection for reading 223 223 IF (FILE_TEST(iodirout, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 224 msg = 'eee : the directory' + iodirout + ' is not accessible.'225 ras = report(msg)226 return, result224 msg = 'eee : the directory' + iodirout + ' is not accessible.' 225 ras = report(msg) 226 return, result 227 227 ENDIF 228 228 ; … … 242 242 fullfilename_d2m = isafile(iodirout + filename_d2m, NEW=0, /MUST_EXIST) 243 243 IF fullfilename_d2m[0] EQ '' THEN BEGIN 244 msg = 'eee : the file ' + fullfilename_d2m + ' was not found.'245 ras = report(msg)246 return, result244 msg = 'eee : the file ' + fullfilename_d2m + ' was not found.' 245 ras = report(msg) 246 return, result 247 247 ENDIF 248 248 ; … … 255 255 fullfilename_t2m = isafile(iodirout + filename_t2m, NEW=0, /MUST_EXIST) 256 256 IF fullfilename_t2m[0] EQ '' THEN BEGIN 257 msg = 'eee : the file ' + fullfilename_t2m + ' was not found.'258 ras = report(msg)259 return, result257 msg = 'eee : the file ' + fullfilename_t2m + ' was not found.' 258 ras = report(msg) 259 return, result 260 260 ENDIF 261 261 ; … … 268 268 fullfilename_msl = isafile(iodirout + filename_msl, NEW=0, /MUST_EXIST) 269 269 IF fullfilename_msl[0] EQ '' THEN BEGIN 270 msg = 'eee : the file ' + fullfilename_msl + ' was not found.'271 ras = report(msg)272 return, result270 msg = 'eee : the file ' + fullfilename_msl + ' was not found.' 271 ras = report(msg) 272 return, result 273 273 ENDIF 274 274 ; … … 278 278 ; in order to avoid unexpected overwritten 279 279 IF (FILE_TEST(fullfilename_out) EQ 1) THEN BEGIN 280 msg = 'eee : the file ' + fullfilename_out + ' already exists.'281 ras = report(msg)282 return, result280 msg = 'eee : the file ' + fullfilename_out + ' already exists.' 281 ras = report(msg) 282 return, result 283 283 ENDIF 284 284 ; … … 351 351 ; 352 352 ncfields = 'q2m[longitude,latitude,*time]=q2m:q2m_attr; ' $ 353 354 355 356 353 + 'longitude[]=xlon:lon_attr; ' $ 354 + 'latitude[]=ylat:lat_attr; ' $ 355 + 'time[]=timein:time_attr ' $ 356 + ' @ globattr' 357 357 ; 358 358 @ncdf_quickwrite
Note: See TracChangeset
for help on using the changeset viewer.