[2] | 1 | ;********************************************************************* |
---|
| 2 | PRO readmesh, filename, filetype, _extra = ex |
---|
| 3 | ; |
---|
| 4 | case filetype of |
---|
| 5 | 'OPA C-Grid.Binary IEEE: Meshmask':BEGIN |
---|
| 6 | meshlec, filename, /pasblabla, _extra = ex |
---|
| 7 | END |
---|
| 8 | 'OPA C-Grid.Net Cdf: Meshmask':BEGIN |
---|
| 9 | ncdf_meshlec, filename, _extra = ex |
---|
| 10 | END |
---|
| 11 | 'Regular 2D.Binary IEEE: mask': |
---|
| 12 | 'Regular 2D.Net Cdf: mask': |
---|
| 13 | 'Regular 3D.Binary IEEE: mask+vertical axis': |
---|
| 14 | 'Regular 3D.Net Cdf: mask+vertical axis': |
---|
| 15 | 'Gaussian 2D.T30 truncature.Binary IEEE: mask': |
---|
| 16 | 'Gaussian 2D.T30 truncature.Net Cdf: mask': |
---|
| 17 | 'Gaussian 2D.T42 truncature.Binary IEEE: mask': |
---|
| 18 | 'Gaussian 2D.T42 truncature.Net Cdf: mask': |
---|
| 19 | 'Gaussian 3D.T30 truncature.Binary IEEE: mask': |
---|
| 20 | 'Gaussian 3D.T30 truncature.Net Cdf: mask': |
---|
| 21 | 'Gaussian 3D.T42 truncature.Binary IEEE: mask': |
---|
| 22 | 'Gaussian 3D.T42 truncature.Net Cdf: mask': |
---|
| 23 | ELSE : |
---|
| 24 | ENDCASE |
---|
| 25 | ; |
---|
| 26 | return |
---|
| 27 | end |
---|
| 28 | ;********************************************************************* |
---|
| 29 | FUNCTION getgridparameter, top |
---|
| 30 | @common |
---|
| 31 | |
---|
| 32 | widget_control, widget_info(top, find_by_uname = 'xmesh'), get_value = answer |
---|
| 33 | jpiglo = long(answer[0]) |
---|
| 34 | jpi = long(answer[1]) |
---|
| 35 | ixminmesh = long(answer[2]) |
---|
| 36 | ixmaxmesh = long(answer[3]) |
---|
| 37 | widget_control, widget_info(top, find_by_uname = 'ymesh'), get_value = answer |
---|
| 38 | jpjglo = long(answer[0]) |
---|
| 39 | jpj = long(answer[1]) |
---|
| 40 | iyminmesh = long(answer[2]) |
---|
| 41 | iymaxmesh = long(answer[3]) |
---|
| 42 | widget_control, widget_info(top, find_by_uname = 'zmesh'), get_value = answer |
---|
| 43 | jpkglo = long(answer[0]) |
---|
| 44 | jpk = long(answer[1]) |
---|
| 45 | izminmesh = long(answer[2]) |
---|
| 46 | izmaxmesh = long(answer[3]) |
---|
| 47 | widget_control, widget_info(top, find_by_uname = 'key_shift') $ |
---|
| 48 | , get_value = answer |
---|
| 49 | key_shift = long(answer[0]) |
---|
| 50 | widget_control, widget_info(top, find_by_uname = 'key_periodique') $ |
---|
| 51 | , get_value = answer |
---|
| 52 | key_periodique = long(answer[0]) |
---|
| 53 | widget_control, widget_info(top, find_by_uname = 'triangulation') $ |
---|
| 54 | , get_value = answer |
---|
| 55 | triangulation = answer[0] EQ 'y' |
---|
| 56 | widget_control, widget_info(top, find_by_uname = 'Glam Boundary') $ |
---|
| 57 | , get_value = answer |
---|
| 58 | rien = execute('boundary = float('+answer[0]+')') |
---|
| 59 | res = {jpiglo:jpiglo, jpjglo:jpjglo, jpkglo:jpkglo, jpi:jpi, jpj:jpj, jpk:jpk $ |
---|
| 60 | , ixminmesh:ixminmesh, ixmaxmesh:ixmaxmesh $ |
---|
| 61 | , iyminmesh:iyminmesh, iymaxmesh:iymaxmesh $ |
---|
| 62 | , izminmesh:izminmesh, izmaxmesh:izmaxmesh $ |
---|
| 63 | , key_shift:key_shift, key_periodique:key_periodique $ |
---|
| 64 | , triangulation:triangulation, boundary:boundary} |
---|
| 65 | |
---|
| 66 | return, res |
---|
| 67 | end |
---|
| 68 | ;********************************************************************* |
---|
| 69 | pro showgridparameter, basetop, EDITABLE = editable, _EXTRA = ex |
---|
| 70 | ; |
---|
| 71 | @common |
---|
| 72 | ;------------------------------------------------------------ |
---|
| 73 | ; |
---|
| 74 | widget_control, basetop, update = 0 |
---|
| 75 | base=widget_base(basetop, /COLUMN, /align_center, _EXTRA = ex) |
---|
| 76 | nothing = cw_fourparameter(base, jpiglo, jpi, ixminmesh, ixmaxmesh, 'x', 'mesh', EDITABLE = editable, CLEAR = editable, uname = 'xmesh') |
---|
| 77 | nothing = cw_fourparameter(base, jpjglo, jpj, iyminmesh, iymaxmesh, 'y', 'mesh', EDITABLE = editable, CLEAR = editable, uname = 'ymesh') |
---|
| 78 | nothing = cw_fourparameter(base, jpkglo, jpk, izminmesh, izmaxmesh, 'z', 'mesh', EDITABLE = editable, CLEAR = editable, uname = 'zmesh') |
---|
| 79 | nothing = widget_label(base, value = 'autres parametres ...') |
---|
| 80 | basea=widget_base(base, /row, /align_center) |
---|
| 81 | nothing = widget_label(basea, value = 'key_shift') |
---|
| 82 | if NOT keyword_set(key_shift) then key_shift = 0 |
---|
| 83 | nothing = widget_text(basea, value = strtrim(key_shift*(1-keyword_set(editable)),1), uname = 'key_shift', xsize = 4, EDITABLE = editable) |
---|
| 84 | nothing = widget_label(basea, value = 'key_periodique') |
---|
| 85 | if NOT keyword_set(key_periodique) then key_periodique = 0 |
---|
| 86 | nothing = widget_text(basea, value = strtrim(key_periodique*(1-keyword_set(editable)),1), uname = 'key_periodique', xsize = 4, EDITABLE = editable) |
---|
| 87 | baseb=widget_base(base, /row, /align_center) |
---|
| 88 | nothing = widget_label(baseb, value = 'use a triangulation (y/n) ?') |
---|
| 89 | nothing = widget_text(baseb, value = 'y', uname = 'triangulation', xsize = 4, EDITABLE = editable) |
---|
| 90 | if NOT keyword_set(editable) then begin |
---|
| 91 | ming = min(glamt, max = maxg) |
---|
| 92 | value = tostr([floor(ming), ceil(maxg)]) |
---|
| 93 | ENDIF ELSE value = tostr([0, 0]) |
---|
| 94 | nothing = widget_label(baseb, value = 'Glam Boundary') |
---|
| 95 | nothing = widget_text(baseb, value = value, uname = 'Glam Boundary', xsize = 10, EDITABLE = editable) |
---|
| 96 | nothing = widget_button(base, value = 'Ok, Get the Mesh parameters!', uvalue = {name:'ok finish'}) |
---|
| 97 | widget_control, basetop, update = 1 |
---|
| 98 | ; |
---|
| 99 | return |
---|
| 100 | END |
---|
| 101 | ;********************************************************************* |
---|
| 102 | FUNCTION give_file_type_choice, type = type |
---|
| 103 | if keyword_set(type) then begin |
---|
| 104 | case type of |
---|
| 105 | 'nc':desc = [ '1\File or Grid type: ' , $ |
---|
| 106 | '0\OPA C-Grid.Net Cdf: Meshmask' , $ |
---|
| 107 | '0\Regular 2D.Net Cdf: mask' , $ |
---|
| 108 | '0\Regular 3D.Net Cdf: mask+vertical axis' , $ |
---|
| 109 | '1\Gaussian 2D' , $ |
---|
| 110 | '0\T30 truncature.Net Cdf: mask' , $ |
---|
| 111 | '2\T42 truncature.Net Cdf: mask' , $ |
---|
| 112 | '1\Gaussian 3D' , $ |
---|
| 113 | '0\T30 truncature.Net Cdf: mask+vertical axis' , $ |
---|
| 114 | '2\T42 truncature.Net Cdf: mask+vertical axis' , $ |
---|
| 115 | '2\No file' ] |
---|
| 116 | endcase |
---|
| 117 | ENDIF ELSE BEGIN |
---|
| 118 | desc = [ '1\File or Grid type: ' , $ |
---|
| 119 | '0\OPA C-Grid.Binary IEEE: Meshmask' , $ |
---|
| 120 | '0\Regular 2D.Binary IEEE: mask' , $ |
---|
| 121 | '0\Regular 3D.Binary IEEE: mask+vertical axis' , $ |
---|
| 122 | '1\Gaussian 2D' , $ |
---|
| 123 | '0\T30 truncature.Binary IEEE: mask' , $ |
---|
| 124 | '2\T42 truncature.Binary IEEE: mask' , $ |
---|
| 125 | '1\Gaussian 3D' , $ |
---|
| 126 | '0\T30 truncature.Binary IEEE: mask+vertical axis' , $ |
---|
| 127 | '2\T42 truncature.Binary IEEE: mask+vertical axis' , $ |
---|
| 128 | '2\No file' ] |
---|
| 129 | ; desc = [ '1\File or Grid type: ' , $ |
---|
| 130 | ; '1\OPA C-Grid' , $ |
---|
| 131 | ; '0\Binary IEEE: Meshmask' , $ |
---|
| 132 | ; '2\Net Cdf: Meshmask' , $ |
---|
| 133 | ; '1\Regular 2D' , $ |
---|
| 134 | ; '0\Binary IEEE: mask' , $ |
---|
| 135 | ; '2\Net Cdf: mask' , $ |
---|
| 136 | ; '1\Regular 3D' , $ |
---|
| 137 | ; '0\Binary IEEE: mask+vertical axis' , $ |
---|
| 138 | ; '2\Net Cdf: mask+vertical axis' , $ |
---|
| 139 | ; '1\Gaussian 2D' , $ |
---|
| 140 | ; '1\T30 truncature' , $ |
---|
| 141 | ; '0\Binary IEEE: mask' , $ |
---|
| 142 | ; '2\Net Cdf: mask' , $ |
---|
| 143 | ; '1\T42 truncature' , $ |
---|
| 144 | ; '0\Binary IEEE: mask' , $ |
---|
| 145 | ; '2\Net Cdf: mask' , $ |
---|
| 146 | ; '2\' , $ |
---|
| 147 | ; '1\Gaussian 3D' , $ |
---|
| 148 | ; '1\T30 truncature' , $ |
---|
| 149 | ; '0\Binary IEEE: mask+vertical axis' , $ |
---|
| 150 | ; '2\Net Cdf: mask+vertical axis' , $ |
---|
| 151 | ; '1\T42 truncature' , $ |
---|
| 152 | ; '0\Binary IEEE: mask+vertical axis' , $ |
---|
| 153 | ; '2\Net Cdf: mask+vertical axis' , $ |
---|
| 154 | ; '2\' , $ |
---|
| 155 | ; '2\No file' ] |
---|
| 156 | ENDELSE |
---|
| 157 | return, desc |
---|
| 158 | end |
---|
| 159 | ;********************************************************************* |
---|
| 160 | PRO whichgrid_event, event |
---|
| 161 | @common |
---|
| 162 | widget_control, event.id, get_uvalue = eventuvalue |
---|
| 163 | IF chkstru(eventuvalue,'name') EQ 0 THEN return |
---|
| 164 | case eventuvalue.name OF |
---|
| 165 | 'Cancel':BEGIN |
---|
| 166 | widget_control, event.handler, get_uvalue = messenger |
---|
| 167 | *messenger = -1 |
---|
| 168 | widget_control,event.handler,/destroy |
---|
| 169 | END |
---|
| 170 | 'filename':BEGIN |
---|
| 171 | ; on detruit, si il existe deja, ce que l''on va ajouter au widget |
---|
| 172 | typebaseid = widget_info(event.handler, find_by_uname = 'typebase') |
---|
| 173 | if typebaseid NE 0 then widget_control,typebaseid,/destroy |
---|
| 174 | showgridparameterid = widget_info(event.handler, find_by_uname = 'showgridparameter') |
---|
| 175 | if showgridparameterid NE 0 then widget_control,showgridparameterid,/destroy |
---|
| 176 | ; |
---|
| 177 | widget_control, event.id, get_value = filename |
---|
| 178 | filename = filename[0] |
---|
| 179 | filename = isafile(filename = filename) |
---|
| 180 | if size(filename, /type) NE 7 then return |
---|
| 181 | if rstrpos(filename, '.pro') EQ strlen(filename)-4 then begin |
---|
| 182 | createpro, '@'+strmid(filename, 0, strlen(filename)-4) $ |
---|
| 183 | , filename = isadirectory(io = homedir, title = 'Bad definition of Homedir') $ |
---|
| 184 | +'for_createpro.pro' |
---|
| 185 | showgridparameter, event.handler, group_leader = event.handler,/frame, uname = 'showgridparameter' |
---|
| 186 | ENDIF ELSE BEGIN |
---|
| 187 | basetype=widget_base(event.handler, /row, /align_center, group_leader = event.handler, uname = 'typebase') |
---|
| 188 | type = '' |
---|
| 189 | if rstrpos(filename, '.nc') EQ strlen(filename)-3 then type = 'nc' |
---|
| 190 | nothing = cw_pdmenu(basetype, give_file_type_choice(type = type), /return_full_name, uname = 'file type choice', uvalue = {name:'file type choice'}) |
---|
| 191 | nothing = widget_text(basetype, value = '', uname = 'file type', uvalue = {name:'file type'}) |
---|
| 192 | ENDELSE |
---|
| 193 | END |
---|
| 194 | 'browse filename':BEGIN |
---|
| 195 | ; on detruit, si il existe deja, ce que l''on va ajouter au widget |
---|
| 196 | typebaseid = widget_info(event.handler, find_by_uname = 'typebase') |
---|
| 197 | if typebaseid NE 0 then widget_control,typebaseid,/destroy |
---|
| 198 | showgridparameterid = widget_info(event.handler, find_by_uname = 'showgridparameter') |
---|
| 199 | if showgridparameterid NE 0 then widget_control,showgridparameterid,/destroy |
---|
| 200 | ; |
---|
| 201 | filename = isafile() |
---|
| 202 | if size(filename, /type) NE 7 then return |
---|
| 203 | widget_control, widget_info(event.handler, find_by_uname = 'Filename'), set_value = filename |
---|
| 204 | if rstrpos(filename, '.pro') EQ strlen(filename)-4 then begin |
---|
| 205 | createpro, '@'+strmid(filename, 0, strlen(filename)-4) $ |
---|
| 206 | , filename = isadirectory(io = homedir, title = 'Bad definition of Homedir') $ |
---|
| 207 | +'for_createpro.pro' |
---|
| 208 | showgridparameter, event.handler, group_leader = event.handler,/frame, uname = 'showgridparameter' |
---|
| 209 | ENDIF ELSE BEGIN |
---|
| 210 | basetype=widget_base(event.handler, /row, /align_center, group_leader = event.handler, uname = 'typebase') |
---|
| 211 | type = '' |
---|
| 212 | if rstrpos(filename, '.nc') EQ strlen(filename)-3 then type = 'nc' |
---|
| 213 | nothing = cw_pdmenu(basetype, give_file_type_choice(type = type), /return_full_name, uname = 'file type choice', uvalue = {name:'file type choice'}) |
---|
| 214 | nothing = widget_text(basetype, value = '', uname = 'file type', uvalue = {name:'file type'}) |
---|
| 215 | ENDELSE |
---|
| 216 | END |
---|
| 217 | 'file type choice':BEGIN |
---|
| 218 | widget_control, widget_info(event.handler,find_by_uname='Filename'), get_value = filename |
---|
| 219 | filename = filename[0] |
---|
| 220 | filetype = event.value |
---|
| 221 | filetype = strmid(filetype, strpos(filetype, '.')+1) |
---|
| 222 | widget_control, widget_info(event.handler, find_by_uname = 'file type'), set_value = filetype |
---|
| 223 | case filetype of |
---|
| 224 | 'OPA C-Grid.Binary IEEE: Meshmask':BEGIN |
---|
| 225 | meshlec, filename, /pasblabla, /getdimensions |
---|
| 226 | showgridparameter, event.handler, /editable, group_leader = event.handler,/frame, uname = 'showgridparameter' |
---|
| 227 | END |
---|
| 228 | 'OPA C-Grid.Net Cdf: Meshmask':BEGIN |
---|
| 229 | ncdf_meshlec, filename, /getdimensions |
---|
| 230 | showgridparameter, event.handler, /editable, group_leader = event.handler,/frame, uname = 'showgridparameter' |
---|
| 231 | END |
---|
| 232 | 'Regular 2D.Binary IEEE: mask': |
---|
| 233 | 'Regular 2D.Net Cdf: mask': |
---|
| 234 | 'Regular 3D.Binary IEEE: mask+vertical axis': |
---|
| 235 | 'Regular 3D.Net Cdf: mask+vertical axis': |
---|
| 236 | 'Gaussian 2D.T30 truncature.Binary IEEE: mask': |
---|
| 237 | 'Gaussian 2D.T30 truncature.Net Cdf: mask': |
---|
| 238 | 'Gaussian 2D.T42 truncature.Binary IEEE: mask': |
---|
| 239 | 'Gaussian 2D.T42 truncature.Net Cdf: mask': |
---|
| 240 | 'Gaussian 3D.T30 truncature.Binary IEEE: mask': |
---|
| 241 | 'Gaussian 3D.T30 truncature.Net Cdf: mask': |
---|
| 242 | 'Gaussian 3D.T42 truncature.Binary IEEE: mask': |
---|
| 243 | 'Gaussian 3D.T42 truncature.Net Cdf: mask': |
---|
| 244 | ELSE : |
---|
| 245 | endcase |
---|
| 246 | |
---|
| 247 | END |
---|
| 248 | 'file type':BEGIN |
---|
| 249 | help, event, /struct |
---|
| 250 | END |
---|
| 251 | 'ok finish':BEGIN |
---|
| 252 | widget_control, widget_info(event.handler,find_by_uname='Filename'),get_value=filename |
---|
| 253 | filename = filename[0] |
---|
| 254 | typeid = widget_info(event.handler,find_by_uname='file type') |
---|
| 255 | if typeid NE 0 then begin |
---|
| 256 | widget_control, typeid, get_value=filetype |
---|
| 257 | filetype = filetype[0] |
---|
| 258 | ENDIF ELSE filetype = 'batch file' |
---|
| 259 | widget_control, event.handler, get_uvalue = messenger |
---|
| 260 | *messenger = create_struct('filename', filename, 'filetype', filetype, getgridparameter(event.handler)) |
---|
| 261 | ; |
---|
| 262 | readmesh, filename, filetype[0], glamboundary = (*messenger).boundary |
---|
| 263 | ; |
---|
| 264 | widget_control,event.handler,/destroy |
---|
| 265 | END |
---|
| 266 | ELSE: |
---|
| 267 | endcase |
---|
| 268 | return |
---|
| 269 | END |
---|
| 270 | ;********************************************************************* |
---|
| 271 | FUNCTION whichgrid, name, IODIRECTORY = iodirectory, PARENT = parent, _EXTRA = ex |
---|
| 272 | @common |
---|
| 273 | ; |
---|
| 274 | ; |
---|
| 275 | if n_elements(name) NE 0 then begin |
---|
| 276 | filename = isafile(filename = name) |
---|
| 277 | if size(filename, /type) NE 7 then return, -1 |
---|
| 278 | ENDIF ELSE filaname = 'no file' |
---|
| 279 | ; |
---|
| 280 | ; pour recuperer les reponses possees lors de l''utilisation de ce |
---|
| 281 | ; widget on cree un pointeur que l''on place dans la uvalue. Comme ca |
---|
| 282 | ; une fois que le widget est detruit dans la procedure ...event.pro, |
---|
| 283 | ; la variable surlaquelle pointait le pointeur (contenue ds la uvalue |
---|
| 284 | ; du widget) n''est pas detruite est on peut recuperer le resultat! |
---|
| 285 | ; |
---|
| 286 | messenger = ptr_new(/allocate_heap) |
---|
| 287 | ; |
---|
| 288 | if keyword_set(parent) then BEGIN |
---|
| 289 | base=widget_base(parent, /COLUMN, title = 'whichgrid', /align_center, uvalue = messenger, _EXTRA = ex) |
---|
| 290 | ENDIF ELSE BEGIN |
---|
| 291 | base=widget_base(/COLUMN, title = 'whichgrid', /align_center, uvalue = messenger, _EXTRA = ex) |
---|
| 292 | nothing = widget_button(base, value = 'Cancel', uvalue = {name:'Cancel'}) |
---|
| 293 | ENDELSE |
---|
| 294 | ; |
---|
| 295 | basemane=widget_base(base, /row, /align_center) |
---|
| 296 | nothing = widget_label(basemane, value = 'Mesh Filename or IDL batch file') |
---|
| 297 | nothing = widget_text(basemane, value = filename, uname = 'Filename', uvalue = {name:'filename'}, /editable) |
---|
| 298 | nothing = widget_button(basemane, value = 'Browse', uvalue = {name:'browse filename'}) |
---|
| 299 | ; |
---|
| 300 | if filename NE 'no file' then begin |
---|
| 301 | if rstrpos(filename, '.pro') EQ strlen(filename)-4 then begin |
---|
| 302 | createpro, '@'+strmid(filename, 0, strlen(filename)-4) $ |
---|
| 303 | , filename = isadirectory(io = homedir, title = 'Bad definition of Homedir') $ |
---|
| 304 | +'for_createpro.pro' |
---|
| 305 | showgridparameter, base, group_leader = base,/frame, uname = 'showgridparameter' |
---|
| 306 | ENDIF ELSE BEGIN |
---|
| 307 | basetype=widget_base(base, /row, /align_center, group_leader = base, uname = 'typebase') |
---|
| 308 | type = '' |
---|
| 309 | if rstrpos(filename, '.nc') EQ strlen(filename)-3 then type = 'nc' |
---|
| 310 | nothing = cw_pdmenu(basetype, give_file_type_choice(type = type), /return_full_name, uname = 'file type choice', uvalue = {name:'file type choice'}) |
---|
| 311 | nothing = widget_text(basetype, value = '', uname = 'file type', uvalue = {name:'file type'}) |
---|
| 312 | ENDELSE |
---|
| 313 | |
---|
| 314 | endif |
---|
| 315 | ; |
---|
| 316 | widget_control,base,/realize |
---|
| 317 | xmanager,'whichgrid', base,event_handler = 'whichgrid_event', no_block = 0 |
---|
| 318 | res = *messenger |
---|
| 319 | ptr_free, messenger |
---|
| 320 | |
---|
| 321 | |
---|
| 322 | return, res |
---|
| 323 | end |
---|