Changeset 607

Timestamp:

11/30/12 17:38:00 (12 years ago)

Author:

lelod

Message:

test interpolation bathy

File:

• trunk/src/file_bathy_to_mem.pro (modified) (2 diffs)

trunk/src/file_bathy_to_mem.pro

@@ -83 +83 @@
 ; ==========
 ;
+; - lelod 20121127
+;  introduction d'une interpolation bilineaire de la topographie
+;  pas reussi a faire marcher file_interp, donc creation du tableau
+;  interpole à l'appel (tres rapide) sur la zone choisie
+;  la zone est codee comme latmin, lonmin -marge; latmax, lonmax +marge
+;  marge code en dur = 16
+
 ; $Id$
 ;
@@ -113 +120 @@
 ; lecture fichier land - sea (S. Masson)
 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
-;attention: zone prise en compte limitee par le choix des nos en x et
-;y dans initncdf. La region traitee ne doit pas exceder lonmin=-25 &
-;lonmax=40 & latmin=-5 & latmax=40
-xmin=(lonmin-16+180)*60.
-ymin=(latmin-16+90.)*60.
-xmax=(lonmax+16+180)*60.
-ymax=(latmax+16+90.)*60.
+fileout=project_id_env+'/MASK/ETOPO1_interp.nc'
+
+;  lonmin=0 & latmin=0
+;  lonmax=20 & latmax=20
+; marge autour du domaine cible pour restriction des recherches dans l'orbite
+marge=16.
+xmin=(lonmin-marge+180)*60.
+ymin=(latmin-marge+90.)*60.
+xmax=(lonmax+marge+180)*60.
+ymax=(latmax+marge+90.)*60.
 
 ;initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500
 initncdf,file,zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=xmin,XMAXMESH=xmax,YMINMESH=ymin,YMAXMESH=ymax
-domdef,lonmin-15., lonmax+15., latmin-15., latmax+15.
+domdef,lonmin-marge, lonmax+marge, latmin-marge, latmax+marge
+;file_interp,file,fileout,[2,2], method='bilinear'
+
 bathy=read_ncdf('z',file=file)
 bate=bathy.arr gt 0
+
 xxe1=reform(glamt[*,0])
-zonx=where(xxe1 ge lonmin-15. and xxe1 le lonmax+15.)
-xxe=xxe1[zonx]
+zonx=where(xxe1 ge lonmin-marge and xxe1 le lonmax+marge)
+xxe2=xxe1[zonx]
 yye1=reform(gphit[0,*])
-zony=where(yye1 ge latmin-15. and yye1 le latmax+15.)
-yye=yye1[zony]
-;
+zony=where(yye1 ge latmin-marge and yye1 le latmax+marge)
+yye2=yye1[zony]
+;xxf=fltarr(max(size(zonx))/2,max(size(zony))/2)
+;yyf=fltarr(max(size(zonx))/2,max(size(zony))/2)
+;for j=0,max(size(zony))/2-1 do xxf(*,j)=indgen(max(size(xxe))/2)*2*(max(xxe)-min(xxe))/max(size(xxe))+min(xxe)
+;for i=0,max(size(zonx))/2-1 do yyf(i,*)=indgen(max(size(yye))/2)*2*(max(yye)-min(yye))/max(size(yye))+min(yye)
+xxf=indgen(max(size(xxe))/2)*2
+yyf=indgen(max(size(yye))/2)*2
+Result = bilinear( bathy.arr,xxf, yyf) 
+xxe=indgen(max(size(xxe))/2)*2*(max(xxe)-min(xxe))/max(size(xxe))+min(xxe)
+yye=indgen(max(size(yye))/2)*2*(max(yye)-min(yye))/max(size(yye))+min(yye)
+bate=result gt 0.
+
 IF key_performance EQ 1 THEN BEGIN
     msg = report(['ppp : ' + routine + ' : durée totale ' $