Changeset 608

Timestamp:

11/30/12 17:58:33 (12 years ago)

Author:

lelod

Message:

correction bug

File:

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

trunk/src/file_bathy_to_mem.pro

@@ -122 +122 @@
 fileout=project_id_env+'/MASK/ETOPO1_interp.nc'
 
-;  lonmin=0 & latmin=0
-;  lonmax=20 & latmax=20
+  lonmin=0 & latmin=0
+  lonmax=20 & latmax=20
 ; marge autour du domaine cible pour restriction des recherches dans l'orbite
 marge=16.
@@ -134 +134 @@
 initncdf,file,zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=xmin,XMAXMESH=xmax,YMINMESH=ymin,YMAXMESH=ymax
 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
+;bate=bathy.arr gt 0
 
 xxe1=reform(glamt[*,0])
@@ -149 +148 @@
 ;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
+xxf=indgen(max(size(xxe2))/2)*2
+yyf=indgen(max(size(yye2))/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)
+xxe=indgen(max(size(xxe2))/2)*2*(max(xxe2)-min(xxe2))/max(size(xxe2))+min(xxe2)
+yye=indgen(max(size(yye2))/2)*2*(max(yye2)-min(yye2))/max(size(yye2))+min(yye2)
 bate=result gt 0.