source: issm/oecreview/Archive/14312-15392/ISSM-14381-14382.diff

../trunk-jpl/src/pro/plot_cielo.pro

@@ -1 @@
-print,'Cielo ..'
-openr,1,'Antarctica_CIELO_1km.bin'
-vel=fltarr(5601,5601)
-readu,1,vel & close,1
-
-print,'Read MOA ..'
-read_jpeg,'moa1000_r1_hp1.bamber_white_edit.jpg',im0,order=1
-
-vel=abs(vel)<3000.
-
-val=float(im0)/max(im0)
-indzero=where(vel EQ 0)
-ind=where(vel EQ 0) & ind0=where(vel LT 1.5) & vel(ind0)=1.5
-h=bytscl(alog(vel))/255.*360. & h(ind)=0. & vel(ind)=0.
-s=(0.5+vel/125)/1.5 & s=s<1 & s(ind)=0.
-
-color_convert,h,s,val,r,g,b,/hsv_rgb
-
-a=bytarr(3,5601,5601)
-a(0,*,*)=r
-a(1,*,*)=g
-a(2,*,*)=b
-ss=size(a) & six=ss(2) & siy=ss(3)
-
-re0=bytarr(5601,5601) & re0(*,*)=a(0,*,*)
-gr0=bytarr(5601,5601) & gr0(*,*)=a(1,*,*)
-bl0=bytarr(5601,5601) & bl0(*,*)=a(2,*,*)
-
-write_jpeg,'Antarctic_CIELO_1km_no_label'+'.jpg',a,true=1,order=1
-
-titi=bytarr(5601,5601)
-
-name=['/u/adelie1/eric/FOUND10196/Foundation_basin_index.save',$
-'/u/adelie1/eric/FOUND10196/Mollereisstrom_basin_index.save',$
-'/u/adelie1/eric/FOUND10196/Institute_basin_index.save',$
-'/u/adelie0/eric/Antarctic_GIS/Ian/IceStreamAW_basin_index.save',$
-'/u/adelie0/eric/Antarctic_GIS/Ian/IceStreamC_basin_index.save',$
-'/u/adelie0/eric/Antarctic_GIS/Ian/IceStreamD_basin_index.save',$
-'/u/adelie0/eric/Antarctic_GIS/Ian/IceStreamE_basin_index.save',$
-'/u/adelie0/eric/Antarctic_GIS/Ian/IceStreamF_basin_index.save',$
-'/u/adelie2/eric/RSAT9822-10165SLE/Slessor_basin_index.save',$
-'/u/adelie2/eric/RSAT9822-10165SLE/Bailey_basin_index.save',$
-'/u/adelie7/eric/Jutulstraumen94/Jutul_basin_index.save',$
-'/u/pib3b/eric/Getz/basin_getz_index.save',$
-'/u/adelie7/eric/Land_237/Land_basin_index.save',$
-'/u/adelie10/eric/Getz_437/basin_Hull_index.save',$
-'/u/pib3b/eric/Georges_1994/Track_009/basin_George9_index.save',$  ;Rym
-'/u/adelie9b/eric/RSAT9852-10195REC/rec_basin_index_ext_new.save',$
-'/u/pib1/eric/RSAT9806-10149Byrd/byrd_basin_index.save',$
-'/u/enderby3/eric/Nimrod_1997/basin_nimrod_index.save',$
-'/u/adelie7/eric/PIG22625/THW_basin_index.save',$
-'/u/adelie7/eric/PIG22625/PIG_basin_index.save',$
-'/u/adelie7/eric/PIG22625/KOH_SMI_POP_basin_index.save',$
-'/u/adelie7/eric/PIG22625/ISR_basin_index.save',$
-'/u/adelie7/eric/PIG22625/HAY_basin_index.save',$
-'/u/adelie45/adelie4/eric/Evans96_94/Evans_basin_index.save',$
-'/u/adelie45/adelie4/eric/Rutford96_t65/RUT_basin_index.save',$
-'/u/adelie45/adelie4/eric/Rutford96_t65/CAR_basin_index.save',$
-'/u/adelie45/adelie4/eric/Lambert_318/lambert_basin_index.save',$
-'/u/adelie13c/eric/ScottDenman/basin_DenmanScott_index.save',$
-'/u/adelie13b/eric/Totten24868/Totten_basin_index.save',$
-'/u/adelie45/eric/David23763-24765/david_basin_indexbis.save',$
-'/u/enderby2/eric/Mertz/Mertz_basin_index.save',$
-'/u/enderby2/eric/Mertz/Ninnis_basin_index.save',$
-;'/u/pib1/eric/Stancomb_450/stancomb_basin_index.save',$
-'/u/pib1/eric/Shirase_229/shirase_basin_index.save',$
-'/u/adelie10/eric/Cook/basin_Cook_index.save.new',$
-'/u/adelie9b/eric/RSAT25731_Rayner2000/Rayner_basin_index.save',$
-'/u/adelie1/eric/RSAT9806-10149MUL/Mulock_basin_index.save',$
-'/u/adelie13b/eric/Sabrina_431/Moscow_basin_index.save',$
-;'/u/pib3a/eric/RL_450/Vest_basin_index.save',$
-'/u/adelie1/eric/RSAT9805_10148BEARD/Beardmore_basin_index.save',$
-'/u/adelie2/eric/Frost_402/frost_basin_index_new.save',$
-'/u/adelie13b/eric/Dibble/Dibble_basin_index.save',$
-'/u/enderby2/eric/Phil/basin_Phil_new_index.save',$
-'/u/adelie1/eric/VenableIS/Venable_basin_index.save',$
-'/u/enderby0/eric/Track_482/basin_Abbot92_10_index.save',$
-'/u/enderby0/eric/Track_482/basin_Abbot92_24_index.save',$
-'/u/enderby0/eric/Track_482/basin_Eights_index.save',$
-'/u/adelie2/eric/Eltanin92/Eltanin_basin_index.save',$
-'/u/adelie1/eric/FOUND10196/SUF_basin_index.save',$
-'/u/pib3b/eric/Ragnhild/basin_Ragnhild_index.save',$
-'/u/pib3b/eric/Ragnhild/basin_Astrid_index.save',$
-'/u/enderby2/eric/Scylla25569/basin_american_index.save',$
-'/u/enderby2/eric/Scylla25569/basin_scylla_index.save',$
-'/u/pib3b/eric/Budd_231/basin_Buddetal_index.save',$
-'/u/enderby2/eric/Scylla25269/Part0/basin_robert_index.save',$
-'/u/amundsen1/eric/GeorgeVI_1996_66/basin_lidke_index.save',$
-'/u/amundsen1/eric/GeorgeVI_1996_66/basin_berg_index.save',$
-'/u/enderby3/eric/RL/basin_Ekstrom_index.save',$
-'/u/enderby3/eric/RL/basin_Jelbart_index.save',$
-'/u/enderby3/eric/RL/basin_Riiser_index.save',$
-'/u/enderby3/eric/RL/basin_Stancomb_index.save',$
-'/u/enderby3/eric/RL/basin_Veststraumen_index.save',$
-'/u/enderby3/eric/RL/basin_Coats_index.save',$
-;'/u/adelie12c/eric/Scylla/basin_track166_index.save',$
-'/u/enderby2/eric/Track166/basin_Eng_index.save',$ ;Eng
-'/u/enderby3/eric/RSAT25469-25812-26155Rennick_2000/basin_index_Rennick.save',$
-'/u/enderby3/eric/RSAT25362-25705-26048Ronne_2000/basin_index_Lassiter.save',$
-'/u/pib3b/eric/Georges_1994/Track_019/basin_SmithEwing_index.save',$
-'/u/enderby3/eric/Sulzberger_424/basin_Sulzberger_index.save']
-
-xmin_bamber=-560*5e03
-ymax_bamber=+560*5e03
-for i=0,n_elements(name)-1 do begin
-restore,name(i)
-xmin=xmin0 & ymax=ymax0
-pix=(xfill*1e3+xmin-xmin_bamber)/1000. & lin=(ymax_bamber-(ymax-yfill*1e3))/1000.
-for k=-1,1 do begin
-for l=-1,1 do begin
-titi(pix+k,lin+l)=255
-endfor
-endfor
-endfor
-
-xmin0=-560.*5e3 & ymax0=+560*5e3
-
-namebis=['/u/adelie7/eric/Wordie_381/basin_index_fleming.save',$
-'/u/adelie9b/eric/AP_018/basin_index_Crane.save',$
-'/u/adelie9b/eric/AP_018/basin_index_Jorum.save',$
-'/u/adelie9b/eric/AP_061/basin_index_flask.save',$
-'/u/adelie9b/eric/AP_061/basin_index_leppard.save',$
-;'/u/adelie9b/eric/AP_233/basin_index_mobil.new.save',$
-'/u/adelie9b/eric/AP_476/basin_index_drygalski.save',$
-'/u/adelie9b/eric/AP_476/basin_index_evans.save']
-for i=0,n_elements(namebis)-1 do begin
-restore,namebis(i)
-pix=(xfill*1e2+xmin-xmin0)/1000. & lin=(ymax0-(ymax-yfill*1e2))/1000.
-for k=-1,1 do begin
-for l=-1,1 do begin
-titi(pix+k,lin+l)=255
-endfor
-endfor
-endfor
-
-ind=where(titi ne 0)
-
-help,re0,gr0,bl0,a
-
-re0(ind)=0 & gr0(ind)=0 & bl0(ind)=0
-
-a(0,*,*)=re0 & a(1,*,*)=gr0 & a(2,*,*)=bl0
-
-write_jpeg,'Antarctic_CIELO_1km'+'.jpg',a,true=1,order=1
-
-end

../trunk-jpl/src/pro/bytscl.m

@@ -1 @@
-function value=bytscl(value)
-%INPUT function value=bytscl(value)
-%Equivalent of the bytscl idl routine.
-
-Min=min(value(:));
-Max=max(value(:));
-Top=255;
-
-pow=1;
-value=(Top + 0.9999)*(value - Min).^(pow)/(Max - Min)^(pow);