- Timestamp:
- 03/22/07 10:02:19 (17 years ago)
- Location:
- trunk
- Files:
-
- 1 deleted
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/condmag_on_orca.pro
r14 r18 24 24 ; @keyword DRAKKAR_EXP {type=string} 25 25 ; code for Drakkar experiment 26 ; only used when orcares = ORC 02526 ; only used when orcares = ORCA025 27 27 ; must be G42 ++ G70 28 28 ; … … 37 37 ; provide tools to plot output files 38 38 ; produce a NetCDF GDT or CF compliant 39 ; introduce <progeomag>../<progeomag> and <pro>...</pro> in postprocessing 40 ; of idldoc outputs 41 ; use <pro>fromreg</pro> 39 42 ; 40 43 ; @pre … … 45 48 ; for ORCA2 filename is meshmask_bab.nc 46 49 ; ++ au pif entre mesh_hgr, mesh_zgr et mask 47 ; for ORCA025 2filename is ORCA025-G42_mesh_hgr.nc50 ; for ORCA025 filename is ORCA025-G42_mesh_hgr.nc 48 51 ; 49 52 ; @post … … 83 86 ; by Sebastien Masson 84 87 ; to reproduce /usr/work/sur/fvi/OPA/geomag/cond_sed_ORCA2.nc 85 ; main differences : no yyyy parameter, no time loop 88 ; main differences : no yyyy parameter, no time loop, no mask and no missing 86 89 ; values (++ to be checked ) in data input (condmag.nc), no scale factor, 87 90 ; no OFFSET, no save of weight and addresses … … 208 211 ; 209 212 ; existence and protection of ${GEOMAG_ID} 210 IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN213 IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 211 214 msg = 'eee : the directory' + iodirin + ' is not accessible.' 212 215 PRINT, msg … … 218 221 ; 219 222 ; check if this file exists 220 fullfilename_condmag = isafile(iodirin + filename_condmag, NEW=0, /MUST_EXIST)223 fullfilename_condmag = isafile(iodirin + filename_condmag, NEW=0, /MUST_EXIST) 221 224 IF fullfilename_condmag[0] EQ '' THEN BEGIN 222 225 msg = 'eee : the file ' + fullfilename_condmag + ' was not found.' … … 244 247 ; protection 245 248 IF (FILE_TEST(fullfilename_oce[0], /READ) EQ 0) THEN BEGIN 246 msg = 'eee : the file ' + fullfile oce_condmag[0] + ' is not readable.'249 msg = 'eee : the file ' + fullfilename_oce[0] + ' is not readable.' 247 250 PRINT, msg 248 251 RETURN … … 267 270 ; 268 271 ; existence and protection 269 IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN272 IF (FILE_TEST(iodirout, /DIRECTORY, /WRITE) EQ 0) THEN BEGIN 270 273 msg = 'eee : the directory' + iodirout + ' was not found.' 271 274 PRINT, msg … … 283 286 get_gridparams, fullfilename_condmag[0], $ 284 287 condmaglonname, condmaglatname, $ 285 condmaglon, condmaglat, jpia, jpja, 1, /DOUBLE288 condmaglon, condmaglat, jpia, jpja, 1, /DOUBLE 286 289 ; 287 290 ;---- … … 293 296 get_gridparams, fullfilename_oce[0], $ 294 297 olonname, olatname, $ 295 olon, olat, jpio, jpjo, 2, /DOUBLE298 olon, olat, jpio, jpjo, 2, /DOUBLE 296 299 msg = 'iii : ' + fullfilename_oce[0] + ' opened for read' 297 300 PRINT, msg -
trunk/condmag_output.pro
r16 r18 13 13 ; 14 14 ; @param variable {in}{required}{type=structure return by NCDF_VARINQ} 15 ; { NAME:"", DATATYPE:"", NDIMS:0L, NATTS:0L, DIM:LONARR(NDIMS) } 16 ; variable of condmag.nc to be written 15 ; { NAME:"", DATATYPE:"", NDIMS:0L, NATTS:0L, DIM:LONARR(NDIMS) } 16 ; variable of condmag.nc to be written 17 17 ; variable.NAME must be 'cond_sed' or 'Br' 18 18 ; … … 55 55 ; 56 56 ; @pre 57 ; reading of condmag.nc 57 ; reading of condmag.nc 58 58 ; see condmag_on_orca.pro 59 59 ; … … 215 215 PRINT, msg 216 216 ENDELSE 217 218 217 ; 219 218 ; check jpjo … … 238 237 RETURN 239 238 ENDIF ELSE BEGIN 240 msg = 'iii : valid nb of elements olon parameter' 239 msg = 'iii : valid nb of elements olon parameter' 241 240 PRINT, msg 242 241 ENDELSE 243 242 ; 244 243 ; check values ++ 245 246 244 ; test if ${GEOMAG_OD} defined 247 245 geomag_od_env=GETENV('GEOMAG_OD') … … 255 253 msg = 'iii : ${GEOMAG_OD} is ' + geomag_od_env 256 254 PRINT, msg 257 END 255 END 258 256 ENDCASE 259 257 ; … … 308 306 varid[0] = NCDF_VARDEF(netcdf_id, 'nav_lon' , [dimidx, dimidy], /FLOAT) 309 307 NCDF_ATTPUT, netcdf_id, varid[0], 'units' , 'degrees_east' 310 NCDF_ATTPUT, netcdf_id, varid[0], 'valid_min', min(olon, max = omax), /FLOAT311 NCDF_ATTPUT, netcdf_id, varid[0], 'valid_max', omax, /FLOAT308 NCDF_ATTPUT, netcdf_id, varid[0], 'valid_min', min(olon, max = omax), /FLOAT 309 NCDF_ATTPUT, netcdf_id, varid[0], 'valid_max', omax, /FLOAT 312 310 NCDF_ATTPUT, netcdf_id, varid[0], 'long_name', 'Longitude at t-point' 313 311 ; 314 312 varid[1] = NCDF_VARDEF(netcdf_id, 'nav_lat' , [dimidx, dimidy], /FLOAT) 315 313 NCDF_ATTPUT, netcdf_id, varid[1], 'units' , 'degrees_north' 316 NCDF_ATTPUT, netcdf_id, varid[1], 'valid_min', min(olat, max = omax), /FLOAT317 NCDF_ATTPUT, netcdf_id, varid[1], 'valid_max', omax, /FLOAT314 NCDF_ATTPUT, netcdf_id, varid[1], 'valid_min', min(olat, max = omax), /FLOAT 315 NCDF_ATTPUT, netcdf_id, varid[1], 'valid_max', omax, /FLOAT 318 316 NCDF_ATTPUT, netcdf_id, varid[1], 'long_name', 'Latitude at t-point' 319 317 ; … … 325 323 ; pour min et max, il faut avoir lu la variable ... cf. plus bas ++ 326 324 ; donc pour l'instant on les met à valeur manquante 327 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_min', !VALUES.F_NAN ,/FLOAT328 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_max', !VALUES.F_NAN ,/FLOAT325 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_min', !VALUES.F_NAN, /FLOAT 326 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_max', !VALUES.F_NAN, /FLOAT 329 327 ; 330 328 ;--------------------------- … … 337 335 NCDF_VARPUT, netcdf_id, 'nav_lat', olat 338 336 NCDF_VARPUT, netcdf_id, varid[2], FLOAT(0.5) ; ++ c'est quoi cette valeur 339 340 ;---------------------------341 337 ; 342 338 ; compute min max … … 351 347 ; update min max attributes 352 348 NCDF_CONTROL, netcdf_id, /REDEF 353 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_min', minarr, /FLOAT354 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_max', maxarr, /FLOAT349 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_min', minarr, /FLOAT 350 NCDF_ATTPUT, netcdf_id, varid[3], 'valid_max', maxarr, /FLOAT 355 351 NCDF_CONTROL, netcdf_id, /ENDEF 356 352 ; 357 ;---------------------------358 353 ; close the netcdf files 359 354 NCDF_CLOSE, netcdf_id
Note: See TracChangeset
for help on using the changeset viewer.