Changeset 461 for trunk

Timestamp:

12/13/11 18:52:07 (13 years ago)

Author:

lelod

Message:

plein de corrections du jour

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (15 diffs)
• interpol_correc.pro (modified) (2 diffs)
• interpolswath.pro (modified) (1 diff)

trunk/src/extract_amsua.pro

@@ -77 +77 @@
 ;   IDL> yyyy=2006
 ;   IDL> mm=8
-;   IDL> dd=1
+;   IDL> dd=13
 ;   IDL> resol=1
 ;   IDL> lon_min=-60.
@@ -178 +178 @@
 ; :Todo:
 ;
-; make it work : resete un n_elements bizarre que IDL ne détecte pas à la
-; compilation
 ;
 ; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
@@ -207 +205 @@
 ; decrire la limitation avec les seuils
 ;
-; est-ce que la zon géographique ne peut pas ere déduite des min et max lus ce
-; qui eviterait des parametres dont on avait réussi à se débarasser
 ;
 ;-
@@ -218 +214 @@
 
 nomcanal=strmid(numch,0,1)
-nocanal=0L
+nocanal=0
 reads, strmid(numch,1,1),nocanal,format='(i1.1)'
 
@@ -227 +223 @@
 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
 domdef,lon_min, lon_max, lat_min, lat_max
-
 bate = ncdf_lec(file, var = 'z') GT 0
 xxe = float(ncdf_lec(file, var = 'lon'))
@@ -234 +229 @@
 jpje = n_elements(yye)
 
-; extraction zone donnees AMSU
-domdef,lon_min,lon_max,lat_min,lat_max
-bate = ncdf_lec(file, var = 'z') GT 0
-xxe = float(ncdf_lec(file, var = 'lon'))
-yye = float(ncdf_lec(file, var = 'lat'))
-jpie = n_elements;module lecture fichier a introduire ici
 
 ; inutile maintenant: ouverture du fichier de sortie
-look = 'filename'
-scale = 1.
-geomin = geolocation_to_string_idl(lon_min, lat_min, look, scale)
-geomax = geolocation_to_string_idl(lon_max, lat_max, look, scale)
-fileout=project_id_env+ 'AMSU/' $
-        + string(yyyy,format='(I4.4)') + '/' $
-        + string(mm,format='(I2.2)')   + '/' $
-        + numch + '_' $
-        + string(yyyy,format='(I4.4)') $
-        + string(mm,format='(I2.2)') $
-        + string(dd,format='(I2.2)') + '_' $
-        + geomin + '_' $
-        + geomax $
-        + '.dat'
-print, 'ouverture pour ecriture de ', fileout
-openw, lun_obs, fileout, /get_lun
+;look = 'filename'
+;scale = 1.
+;geomin = geolocation_to_string_idl(lon_min, lat_min, look, scale)
+;geomax = geolocation_to_string_idl(lon_max, lat_max, look, scale)
+;fileout=project_id_env+ 'AMSU/' $
+;        + string(yyyy,format='(I4.4)') + '/' $
+;        + string(mm,format='(I2.2)')   + '/' $
+;        + numch + '_' $
+;        + string(yyyy,format='(I4.4)') $
+;        + string(mm,format='(I2.2)') $
+;        + string(dd,format='(I2.2)') + '_' $
+;        + geomin + '_' $
+;        + geomax $
+;        + '.dat'
+;print, 'ouverture pour ecriture de ', fileout
+;openw, lun_obs, fileout, /get_lun
 
 
 ; appel au ssprgm qui interpole les fonctions de correction (inutile
 ; de l'appeler a chaque fichier....)
-if nomcanal eq 'a' then nbpixel=30
-if nomcanal eq 'b' then nbpixel=90
-interpol_correc,nomcanal,nbpixel,cor_l,cor_s,swath,track
+
+interpol_correc,numch,nbpixel,cor_l,cor_s,swath,track
 
 
@@ -282 +270 @@
 filea=''
 nlist=1
-ilist=1
+ilist=0
    ; while (not eof(lun_a)) do begin
 while ilist le nlist do begin
@@ -342 +330 @@
 ; correction nadir des donnees
    ch_nadir=fltarr(nbpix,n_scan)
-   lanseamask=intarr(nbpix,n_scan)
+   landseamask=intarr(nbpix,n_scan)
    for isc=0L,n_scan-1L do begin
       for ifo=0,nbpix-1 do begin
@@ -355 +343 @@
             bate2=bate1[*,yind]
             if (mean(bate2) le 0.5) then begin
-               ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[numcanal,ifo]
+               ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[ifo]
                landseamask[ifo,isc]=0 ; mer
             endif else begin
-               ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[numcanal,ifo]
+               ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[ifo]
                landseamask[ifo,isc]=1 ; terre
             endelse
@@ -390 +378 @@
 ; selection des taches au sol situees dans la zone d'interet
       zone=where((longrid ge lon_min) and (longrid le lon_max) $
-                 and (latgrid ge lat_min) and (latgrid le lat_max),npt)
-;print, 'npt', npt
+                 and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin), npt)
       if npt ne 0 then begin
-         time[zone]=tt[i]
- ; alerte: comment initialiser les tableaux qu'on concatene????
          if (n_elements(chint) EQ 0) then begin
+            timeint=replicate(tt[i],npt)
             chint=tbgrid[zone]
             latint=latgrid[zone]
@@ -401 +387 @@
             fovint=fovgrid[zone]
             maskint=maskgrid[zone]
-            timeint= time[zone]
          endif else begin
             chint=[chint,tbgrid[zone]]
@@ -408 +393 @@
             fovint=[fovint,fovgrid[zone]]
             maskint=[maskint,maskgrid[zone]]
-            timeint=[timeint,time[zone]]
+            timeint=[timeint,replicate(tt[i],npt)]
          endelse
       endif else begin
@@ -416 +401 @@
    endfor
 
-
-   for i=0L, npt-1L do begin
-      if (( chint[i] lt tbmin) or (chint[i] gt tbmax) ) then begin
-         ch_adj[i]=!values.f_nan
-      endif
-   endfor
-
-
-   print,'nb points dans la zone',npt
-
-   for i=0L, npt-1L do begin
+   nn=n_elements(chint)
+ ;  for i=0L, npt-1L do begin
+ ;     if (( chint[i] lt tbmin) or (chint[i] gt tbmax) ) then begin
+ ;        chint[i]=!values.f_nan
+ ;     endif
+ ;  endfor
+
+
+   print,'nb points dans la zone',nn
+
+   for i=0L, nn-1L do begin
 ;format: nosat,mois, day, lon, lat, zen,fov, time, ama_ch1:15,amb_ch1:5
 ; print, 'mois::::', mm,'-----', nab-i
-      printf,lun_obs, mm, jour_ama, desc, nosat, fovint[i], landseamask, lonint[i], latint[i], timeint[i],$
-             chint[i],format='(6(i3,1x),4(f8.4,1x))'
+;      printf,lun_obs, mm, jour_ama, desc, nosat, fovint[i], maskint[i], lonint[i], latint[i], timeint[i],$
+;             chint[i],format='(6(i3,1x),4(f8.4,1x))'
       if (n_elements(data) EQ 0) then begin
          data = { desc : desc $
@@ -460 +445 @@
     ; write outputfile
    ; ajouter resol dans le header
-    header = { yyyy : yyyy $
+    header1 = { yyyy : yyyy $
              ,  mm : mm $
              ,  dd : dd $
              ,  numch : numch $
+             ,  nbpix : nbgrid $
+             ,  resol : resol $
              }
-    amsu_t2={header: header, data: data}
+header2= { lon_min: lon_min $
+                      , lon_max: lon_max $
+                      , lat_min: lat_min $
+                      , lat_max: lat_max $
+                     }
+
+    amsu_t2={header1: header1, header2: header2, data: data}
     help, amsu_t2,/structure
     testfilename='./new_amsu_t2.dat'
-    fileout = mem_to_file_amsu_t2(yyyy, mm, dd, numch $
-                               , lon_min, lon_max, lat_min $
-                               , lat_max, amsu_t2, testfilename)
+    fileout = mem_to_file_amsu_t2(amsu_t2, testfilename)
 
 endif else begin
@@ -476 +467 @@
 endelse
 
+window,1
+plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,psym=2
 end

trunk/src/interpol_correc.pro

@@ -107 +107 @@
 ;
 ;-
-PRO interpol_correc,canal,nbpix,cor_l,cor_s,swath,track
+PRO interpol_correc,numch,nbpix,cor_l,cor_s,swath,track
 ;
 @cm_project
-; mettre au carre!!!!!!!!!
-nomcanal=canal(0)
-numcanal=str2num(canal,1)
-
+nomcanal=strmid(numch,0,1)
+numcanal=0
+reads, strmid(numch,1,1),numcanal,format='(i1.1)'
+
+if nomcanal eq 'a' then nbpixel=30
+if nomcanal eq 'b' then nbpixel=90
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
 if nomcanal eq 'a' then begin
@@ -128 +130 @@
 ndfova=nfova/2
 ndfovb=nfovb/2
+print,'nb points',numch,ndfova,ndfovb
 fovb=fltarr(ndfovb)
 fova=fltarr(ndfova)
 alongb=fltarr(ndfovb)
 alongab=alongbtot[ndfovb:nfovb-1] ; verifier l'intervalle!!!!!
-pixb=pixbtot[ndfovb:nfovb]
-alonga=alongatot[ndfova:nfova]
-pixa=pixatot[ndfova:nfova]
+pixb=pixbtot[ndfovb:nfovb-1]
+alonga=alongatot[ndfova:nfova-1]
+pixa=pixatot[ndfova:nfova-1]
 ;ifov nadir : diametre (km)
 ifova=48.05 & ifovb=16

trunk/src/interpolswath.pro

@@ -69 +69 @@
 ;print,'pixatot'
 ;print,pixatot
-window,1
-ploterr,ifov,fova,pixatot
+;window,1
+;ploterr,ifov,fova,pixatot
 
 distmax=-min(fova)+max(fova)+max(pixa)