Changeset 586

Timestamp:

07/17/12 22:25:21 (12 years ago)

Author:

lelod

Message:

adaptation amsub

File:

• trunk/src/interpolswath.pro (modified) (4 diffs)

trunk/src/interpolswath.pro

@@ -24 +24 @@
 ; resolution fixee au diametre du pixel central ou la moitie (resol
 ; vaut 1 ou 2)
+; type= 'a' ou 'b' pour choisir le traitement adapte
 ;
 ; on interpole les longitudes et latitudes en correspondance
@@ -53 +54 @@
 ;
 ;-
-pro interpolswath, tb, latlu, lonlu, masklu, resol, nbgrid, tbint, latgrid, longrid, mask
+pro interpolswath, tb, latlu, lonlu, masklu, type, resol, nbgrid, tbint, latgrid, longrid, mask
 
 @common
@@ -64 +65 @@
 ;
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
-na=max(size(pixatot))
-nnadir=na/2
-pixa=pixatot[nnadir:(na-1)]
+if type='a' then begin
+   na=max(size(pixatot))
+   nnadir=na/2
+   pixa=pixatot[nnadir:(na-1)]
+endif else begin
+   na=max(size(pixbtot))
+   nnadir=na/2
+   pixa=pixbtot[nnadir:(na-1)]
+endelse
 fova=fltarr(na)
 fova[nnadir]=pixa[0]/2.
@@ -134 +141 @@
 ;endfor
 
-
+; interpolation des positions
 latgrid=interpol(latlu,fova,grid)
-longrid=interpol(lonlu,fova,grid)
+; pour la longitude, traitement specifique des longitudes proches de -180
+if (mean(lonlu) ge -120.) then begin
+   longrid=interpol(lonlu,fova,grid)
+endif else begin
+   ind=where(lonlu lt 0.,neg)
+   if neg ne 0 then begin
+      lonlu(ind)=360.+lonlu(ind)
+      longrid=interpol(lonlu,fova,grid)
+   endif
+   ind=where(longrid lt 0.,neg)
+   if neg ne 0 then begin
+      longrid(ind)=longrid(ind)-360.
+      longrid=interpol(lonlu,fova,grid)
+   endif
+endelse
 
 ;print,tbint