Changeset 521

Timestamp:

04/16/12 18:02:17 (12 years ago)

Author:

pinsard

Message:

isolate bathy reading

Location:

trunk/src

Files:

• cm_project.pro (modified) (2 diffs)
• extract_amsua.pro (modified) (5 diffs)
• file_bathy_to_mem.pro (copied) (copied from trunk/src/extract_amsua.pro) (7 diffs)
• traite_amsuab.sh (modified) (3 diffs)

trunk/src/cm_project.pro

@@ -24 +24 @@
 ; $URL$
 ;
+; - fplod 20120416
+;
+;   * add arrays for bathy
+;
 ; - fplod 20110422T132320Z aedon.locean-ipsl.upmc.fr (Darwin)
 ;
@@ -30 +34 @@
 ;-
 COMMON project, project_env, project_id_env, project_od_env
+COMMON bathy, bate, xxe, yye

trunk/src/extract_amsua.pro

@@ -16 +16 @@
 ; decode les noms des fichiers donnés en argument dans **files_list**
 ; dans la date choisie, puis appelle le prgm de lecture
+;
+; utilise les tableaux bate, xxe, yye définis dans le common :file:`cm_project.pro` remplis par :func:`file_bathy_to_mem` 
 ;
 ; appelle :ref:`interpol_correc.pro`, qui fournit les fonctions de correction au
@@ -109 +111 @@
 ; called by :ref:`traite_amsuab.sh`
 ;
-; previous step : :func:`search_amsufiles`
+; previous step : :func:`search_amsufiles`,  :func:`file_bathy_to_mem`
 ;
 ; use :ref:`read_amsua1c.pro`, :ref:`read_amsub1c.pro`,
@@ -126 +128 @@
 ; interpolswath ne gere pas cette situation
 ;
-; comprendre les points bizarres autour de -40,30 visibles avec l'exemple
-; de 20060801
-;
 ; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
 ; crains des zeros louches!!).louche
@@ -158 +157 @@
 ;
 ; $URL$
+;
+; - fplod 20120416
+;
+;   * remove reading of bathy (see :func:`file_bathy_to_mem`)
 ;
 ; - fplod 20120323
@@ -299 +302 @@
 tbmax=350
 PRINT, 'www : debut programme extract_amsu',SYSTIME()
-; 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 lon_min=-25 &
-;lon_max=40 & lat_min=-5 & lat_max=40
-
-initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500
-domdef,lon_min-15, lon_max+15, lat_min-15, lat_max+15
-bathy=read_ncdf('z',file=file)
-bate=bathy.arr gt 0
-xxe1=reform(glamt[*,0])
-zonx=where(xxe1 ge lon_min-15 and xxe1 le lon_max+15)
-xxe=xxe1[zonx]
-yye1=reform(gphit[0,*])
-zony=where(yye1 ge lat_min-15 and yye1 le lat_max+15)
-yye=yye1[zony]
 jpie = n_elements(xxe)
 jpje = n_elements(yye)
 ;print,jpie,jpje,size(bate)
-;PRINT, 'lecture bathy terminee',SYSTIME()
+;
 ; appel au ssprgm qui interpole les fonctions de correction (inutile
 ; de l'appeler a chaque fichier....)

trunk/src/file_bathy_to_mem.pro

@@ -2 +2 @@
 ;+
 ;
-; .. _extract_amsua.pro:
+;  .. function:: file_bathy_to_mem(lon_min, lon_max, lat_min, lat_max)
 ;
-; =================
-; extract_amsua.pro
-; =================
+; .. _file_bathy_to_mem.pro:
+;
+; =====================
+; file_bathy_to_mem.pro
+; =====================
 ;
 ; DESCRIPTION
 ; ===========
 ;
-; prgm de lecture des fichiers AMSU (A seulement, a valider avant de
-; generaliser pour amsub aussi)
+; lit la carte de bathymetrie de S. Masson et remplit les tableaux du common ++
 ;
-; decode les noms des fichiers donnés en argument dans **files_list**
-; dans la date choisie, puis appelle le prgm de lecture
-;
-; appelle :ref:`interpol_correc.pro`, qui fournit les fonctions de correction au
-; nadir et les applique, en fonction du masque terre-mer issu de la
-; carte de bathymetrie de S. Masson
-;
-; applique une interpolation tenant compte de la dimension de la tache
-; au sol selon la position dans la fauchee
-;
-; ecrit les donnees pour la zone (lat/long) choisie
-;
-; ++ le fichier de sortie contient dans l'ordre mm, jj (jour de
-; l'annee), le temps (identique pour les points d'une meme fauchee), le
-; no satellite, le no de fov (1 a 40 ou 1 a 90), long, lat, tb
-;
-; Les valeurs manquantes sont codées par NaN.
-;
-; .. only:: man
+; .. only:: man or text
 ;
 ;    Figure is visible on PDF and HTML documents only.
@@ -40 +23 @@
 ;     .. graphviz::
 ;
-;        digraph extract_amsu {
+;        digraph file_bathy_to_mem {
 ;
-;           amsu_t1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/AMSU/AMSU*/L1C/yyyy/yyyy_mm_dd/*.L1C"];
-;           amsu_t2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax.dat"];
+;           bathy [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/etopo++"];
+;           common++ [shape=ellipse,fontname=Courier,label="tableaux"];
 ;
-;           extract_amsua [shape=box,
+;           file_bathy_to_mem [shape=box,
 ;           fontname=Courier,
 ;           color=blue,
-;           URL="http://forge.ipsl.jussieu.fr/varamma/browser/trunk/src/extract_amsua.pro",
-;           label="${PROJECT}/src/extract_amsua.pro"];
+;           URL="http://forge.ipsl.jussieu.fr/varamma/browser/trunk/src/file_bathy_to_mem.pro",
+;           label="${PROJECT}/src/file_bathy_to_mem.pro"];
 ;
-;           {amsu_t1} -> {extract_amsua} -> {amsu_t2};
+;           {bathy} -> {file_bathy_to_mem} -> {common++};
 ;       }
 ;
-; :param numch:
-; :param files_list:
-; :param yyyy:
-; :param mm:
-; :param resol:
 ; :param lon_min:
 ; :param lon_max:
@@ -67 +45 @@
 ; ========
 ;
-; Creating an articial list of AMSU-A files::
+; ++::
 ;
-;   files_list = strarr(2)
-;   files_list[0] = project_id_env + '/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C'
-;   files_list[1] = project_id_env + 'AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C'
-;
-; This above tip can replace the call to :func:`search_amsufiles`.
-;
-; Pour ne garder que le dernier fichier::
-;
-;   files_list = search_amsufiles(numch, yyyy, mm, dd)
-;   help, files_list
-;   files_list_last_array = strarr(1)
-;   files_list_last_array[0] = files_list[N_ELEMENTS(files_list) - 1]
-;   help, files_list_last_array
-;   files_list = files_list_last_array
-;   help, files_list
-;
-; Using AMSU-A channel 5::
-;
-;   numch = 'a5'
-;   yyyy=2006L
-;   mm=8
-;   dd=13
-;   resol=1
 ;   lon_min=-25.
 ;   lon_max=25.
 ;   lat_min=-5.
 ;   lat_max=20.
-;   files_list = search_amsufiles(numch, yyyy, mm, dd)
-;   extract_amsua, numch, files_list, yyyy, mm, dd, resol, lon_min, lon_max, lat_min, lat_max
-;
-; :file:`${PROJECT_ID}/AMSU/2006/08/a5_20060813_025w05s_025e20n.dat`
-; must have been created.
+;   result = file_bathy_to_mem(lon_min, lon_max, lat_min, lat_max)
 ;
 ; SEE ALSO
 ; ========
 ;
-; :ref:`data_amsu`
+; :ref:`data_bathy`
 ;
 ; called by :ref:`traite_amsuab.sh`
 ;
-; previous step : :func:`search_amsufiles`
-;
-; use :ref:`read_amsua1c.pro`, :ref:`read_amsub1c.pro`,
-; :func:`geolocation_to_string_idl`, :func:`mem_to_file_amsu_t2`
+; use ++
 ;
 ; next step :  ++
@@ -119 +67 @@
 ; ====
 ;
-; track why end with::
-;
-;  % Program caused arithmetic error: Floating illegal operand
-;
-; traiter les longitudes autour de 180 degres (passage de 180 a -180)
-; interpolswath ne gere pas cette situation
-;
-; comprendre les points bizarres autour de -40,30 visibles avec l'exemple
-; de 20060801
-;
-; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
-; crains des zeros louches!!).louche
-;
-; vérifier que pas NaN dans les lignes écrites
-;
-; faire un test avec un zone géographique dans le pacifique pour voir comment
-; le système réagi (devrait dire : rien à faire !)
-;
-; ajouter MHS
-;
-; eviter la creation d'un fichier vide (si pas de points dans la zone)
-;
-; tester la cohérence entre les fichiers présents dans files_list
-; en terme de date avec yyyy, mm et dd donnés en paramètre
+; make it work : % Program caused arithmetic error: Floating illegal operand dans le read_netcf à la main
 ;
 ; improve log
 ;
-; get rid of uppercase
-;
 ; check args
 ;
-; decrire la limitation avec les seuils
+; get rid of hard coded values -15, +15
+;
+; get rid of hard coded values XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500
+; thoses ones oare valid  for West Africa ONLY
 ;
 ; EVOLUTIONS
@@ -159 +85 @@
 ; $URL$
 ;
-; - fplod 20120323
+; - fplod 20120416
 ;
-;   * homogenize examples to traite_amsuab.sh
-;   * fix output name in example
-;   * add a paramater files_list to remplace usage of :file:`list_file`
-;   * replace usage of :file:`list_file` by the usage of the string array
-;     files_list
-;   * add a test if files_list no empty (see nofile label)
-;   * reorder SEE ALSO and TODO
-;
-; - lelod 20120223
-;
-;   * modif lecture du fichier bathy (ETOPO) pour limiter le temps
-;     necessaire. Introduction de coordonnees limite dans initncdf, puis
-;     domdef
-;
-;     Attention, si on change la zone, il faut changer les coordonnees!!!
-;
-; - lelod 20111215
-;
-;   * tests de fonctionnement effectues: la correction nadir est
-;     correctement implementee
-;
-;     l'interpolation de la fauchee fonctionne
-;     mais attention test sur la longitude code en dur:
-;     pbs avec les fauchees qui traversent le meridien 180deg!!!!
-;     l'interpolation ajoute des points parasites dans la zone...
-;
-;
-; - fplod 20111214T082300Z aedon.locean-ipsl.upmc.fr (Darwin)
-;
-;   * remove old style writing output file
-;   * forget idldoc
-;
-; - fplod 20111213T092550Z aedon.locean-ipsl.upmc.fr (Darwin)
-;
-;   * fix syntax error : rewrite labfile label
-;   * fix nocanal string to integer type
-;   * fix example (missing resol in call)
-;
-; - lelod 20111211
-;
-;   * reintegre les modules de correct_nadir_amsu dans extract
-;     prevu decodage de numch pour connaitre le nom instrument et le
-;     numero du canal (pas fini, au debut) --> a faire aussi dans les
-;     ssprgm appeles (interpol_correc)
-;     espere que c'est le dernier avatar de cette !!!!?XX^! de
-;     chaine de traitement AMSU
-;
-; - fplod 20111209T125312Z aedon.locean-ipsl.upmc.fr (Darwin)
-;
-;   * chgt de terminologie
-;   * ajout parametre numch
-;   * transition pour utilisation de :ref:`mem_to_file_amsu_t2`
-;
-; - lelod 20111121
-;
-;   * introduction du parametre "descending" (orbite descendante 1,
-;     montante 0) dans le fichier de sortie : a tester!
-;
-; - fplod 20110818T160435Z aedon.locean-ipsl.upmc.fr (Darwin)
-;
-;   * no more /bdd in examples
-;   * need to check if missing exist before NaN affection on idl7 cratos
-;
-; - pinsard 2011-08-18T10:34:53Z loholt1.ipsl.polytechnique.fr (Linux)
-;
-;   * correction of examples
-;   * examples with 2 files in list_filea and list_fileb
-;   * match technique from
-;     http://cerebus.as.arizona.edu/~ioannis/teaching/idl/tutorial_strings.html
-;   * revision du cas use_amsua=1 : la liste de fichiers amsu-a n'est ouverte
-;     qu'une seule fois
-;
-; - pinsard 2011-08-18T09:10:08Z loholt1.ipsl.polytechnique.fr (Linux)
-;
-;   * move bad values from -9999.99 to NaN
-;
-; - pinsard 2011-08-17T14:58:13Z loholt1.ipsl.polytechnique.fr (Linux)
-;
-;   * use geolocation_to_string_idl
-;
-; - pinsard 2011-08-10T10:46:37Z loholt1.ipsl.polytechnique.fr (Linux)
-;
-;   * remove blank in output filename
-;   * revision of articifial example with use_amsua=0
-;   * add articifial example with use_amsua=1 and try (but fail) to
-;     make it work on loholt1
-;   * remove length path dependencies using IDL function file_basename
-;   * use same date and box of examples in correct_nadir_amsu.pro
-;
-; - fplod 20110810T085644Z aedon.locean-ipsl.upmc.fr (Darwin)
-;
-;   * revision of header
-;
-; - pinsard 2011-05-26T08:44:38Z loholt1.ipsl.polytechnique.fr (Linux)
-;
-;   * strmid(fich,51) is not anymore pertinent to isolate the correct information
-;
-;     It was ok when path as /bdd/AMSU-1C/AMSUAN15/L1C/2006/2006_07_28/::
-;
-;      c='/bdd/AMSU-1C/AMSUAN15/L1C/2006/2006_07_28/NSS.AMAX.NK.D06209.S1323.E1455.B4266162.GC.L1C'
-;      b=strmid(c,51)
-;      print, b
-;      NK.D06209.S1323.E1455.B4266162.GC.L1C
-;
-;     but now path can be anything (PROJECT_ID environment variable)
-;     so the correct way is to catch the last 36 characters
-;
-;     test sur skylla: fonctionne au moins pour un fichier!
+;   * creation from extract_amsua.pro
 ;
 ;-
-PRO extract_amsua, numch, files_list, yyyy, mm, dd, resol $
+function file_bathy_to_mem $
                  , lon_min, lon_max, lat_min, lat_max
 ;
@@ -279 +98 @@
 @common
 
-; test if some files to read
-CASE size(files_list,/DIMENSION) OF
-   0L : BEGIN
-     print, 'pas de fichiers'
-     goto, nofile
-        END
-   ELSE: BEGIN
-       print, 'nb de fichers', size(files_list,/DIMENSION)
-        END
-ENDCASE
-nb_file_array = size(files_list,/DIMENSION)
-nb_file = (nb_file_array)[0]
+; Return to caller if errors
+ON_ERROR, 2
 ;
-nomcanal=strmid(numch,0,1)
-nocanal=0
-reads, strmid(numch,1,1),nocanal,format='(i1.1)'
+result = -1
 
-tbmin=100
-tbmax=350
-PRINT, 'www : debut programme extract_amsu',SYSTIME()
+PRINT, 'www : debut programme file_bathy_to_mem',SYSTIME()
 ; lecture fichier land - sea (S. Masson)
 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
@@ -306 +111 @@
 
 initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500
-domdef,lon_min-15, lon_max+15, lat_min-15, lat_max+15
+domdef,lon_min-15., lon_max+15., lat_min-15., lat_max+15.
 bathy=read_ncdf('z',file=file)
 bate=bathy.arr gt 0
 xxe1=reform(glamt[*,0])
-zonx=where(xxe1 ge lon_min-15 and xxe1 le lon_max+15)
+zonx=where(xxe1 ge lon_min-15. and xxe1 le lon_max+15.)
 xxe=xxe1[zonx]
 yye1=reform(gphit[0,*])
-zony=where(yye1 ge lat_min-15 and yye1 le lat_max+15)
+zony=where(yye1 ge lat_min-15. and yye1 le lat_max+15.)
 yye=yye1[zony]
-jpie = n_elements(xxe)
-jpje = n_elements(yye)
-;print,jpie,jpje,size(bate)
-;PRINT, 'lecture bathy terminee',SYSTIME()
-; appel au ssprgm qui interpole les fonctions de correction (inutile
-; de l'appeler a chaque fichier....)
-interpol_correc,numch,nbpixel,cor_l,cor_s,swath,track
-; ouverture des fichiers liste (annee, mois, jour, tous satellites) pour
-; chaque instrument
-; boucle sur les elements de la liste
 
-print, 'iii : traitement du jour ',yyyy,mm,dd
-a = STRARR(nb_file)
-;PRINT,'demarrage boucle sur lichiers', SYSTIME()
-FOR ifile = 0, nb_file - 1 do begin
-   filea = files_list[ifile]
-
-   COMMON amsua_header,ama_head
-   COMMON amsua_data  ,ama_scan
-   print, 'ouverture et lecture du fichier ', filea, SYSTIME()
-   openr,lu1,filea,Error=erra,/get_lun
-   read_amsua1c,filea, flag1
-   close,lu1
-   free_lun,lu1
-   if (flag1 eq 0) then goto, labfile
-   na=SIZE(reform(ama_scan.btemps[0,*,*]))
-   na_scan=na[2]                ; nb lignes (fauchees) dans l'orbite
-   na_fov=na[1]                 ; nb points dans la fauchee
-   nbpix= na_fov
-   nosat=ama_head.h_satid
-                                ; nb pixels dans la fauchee AMSUA
-   fov=indgen(nbpix)
-   n=na
-   n_scan=na_scan
-   n_fov=na_fov
-
-   ntime_1=fltarr(nbpix,n_scan)
-   fovy_1=fltarr(nbpix,n_scan)
-
-   for j=0L, nbpix-1L do begin
-      fovy_1[j,*]=j+1           ; on remplit un tableau de nscan lignes avec les valeurs de no de pixel
-   endfor
-;
-; on lit la structure et on extrait les infos: temps, longitude,
-; latitude, et les Tbs des differents canaux (15 pour AMSUA)
-
-; extraction du canal traite dans un domaine enveloppant la zone choisie
-   midpix=fix(nbpix/2)
-   amalong=REFORM(ama_scan.latlon[1,*,*]/1.E4)
-   amalati=REFORM(ama_scan.latlon[0,*,*]/1.E4)
-   amcha=REFORM(ama_scan.btemps[nocanal,*,*]/100.)
-   amzeni=REFORM(ama_scan.angles[0,*,*]/100.)
-   ttt=REFORM(ama_scan.scnlintime/3600000.)
-   jnd=where(amalong[midpix,*] gt lon_min-15 and amalong[midpix,*] lt lon_max+15 and amalati[midpix,*] gt lat_min-15 and amalati[midpix,*] lt lat_max+15 ,nzon)
-   print,"iii : nb de points du fichier dans le domaine geographique +15deg ",nzon, SYSTIME()
-   if nzon gt 1 then begin
-      amalat=amalati[*,jnd]
-      amalon=amalong[*,jnd]
-      amch=amcha[*,jnd]
-      amzen=amzeni[*,jnd]
-      tt=ttt[jnd]
-      dims = SIZE(amalon, /DIMENSIONS)
-      print,'dimension tableaux extraits',dims
-      n_scan=dims[1]
-      PRINT,'il y a des donnees dans la region '
-      amaday=REFORM(ama_scan.scnlindy)
-      amafov= fovy_1
-      jour_ama=amaday[0]        ; jour du passage
-; recherche du sens de l'orbite : descendante (1) ou montante (0)
-      lat0=amalat[*,0]
-      lat1=amalat[*,1]
-      if (mean(lat1-lat0) lt 0) then begin
-         desc=1
-      endif else begin
-         desc=0
-      endelse
-;
-; correction nadir des donnees
-         ch_nadir=fltarr(nbpix,nzon)
-         landseamask=intarr(nbpix,nzon)+2 ; valeur hors zone selectionnee
-        ; PRINT,'boucle sur les points du fichier, correction nadir', SYSTIME()
-
-         for isc=0L,nzon-1L do begin
-;recherche de la zone xxe,yye englobant la fauchee
-            utilex=where(xxe ge amalon[midpix,isc]-15 and xxe le amalon[midpix,isc]+15,nxu)
-            if nxu ne 0 then begin
-               xxutile=xxe[utilex]
-               batutilx=bate[utilex,*]
-               utiley=where(yye ge amalat[midpix,isc]-8 and yye le amalat[midpix,isc]+8,nyu)
-               if nyu ne 0 then begin
-                  yyutile=yye[utiley]
-                  batutil=batutilx[*,utiley]
-
-                  for ifo=0,nbpix-1 do begin
-; approximation: 1km= 1deg/100, et on considere le quart en lon et lat
-; de la surface du pixel pour optimiser la localisation du centre du
-; pixel attention approximation brutale des distances en fractions dedegre
-; recherche de la cote par rapport a la resolution AMSUA
-                     xind=where (xxutile ge amalon[ifo,isc]-swath[ifo]/400. and xxutile le amalon[ifo,isc]+swath[ifo]/400.,nxlandsea)
-                     yind=where (yyutile ge amalat[ifo,isc]-track[ifo]/400. and yyutile le amalat[ifo,isc]+track[ifo]/400.,nylandsea)
-                     if(nxlandsea ne 0 and nylandsea ne 0) then begin
-                        bate1=batutil[xind,*]
-                        bate2=bate1[*,yind]
-                        if (mean(bate2) le 0.5) then begin
-                           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[ifo]
-                           landseamask[ifo,isc]=1 ; terre
-                        endelse
-                     endif
-                  endfor
-               endif
-            endif
-
-         endfor
-        ; PRINT, 'fin boucle ',SYSTIME()
-         moych=fltarr(nbpix)
-
-; appel a interpolswath pour ajuster les pixels amsua sur une grille
-; reguliere dans la fauchee
-; et selection de la zone conservee
-; attention test sur la longitude code en dur:
-; pbs avec les fauchees qui traversent le meridien 180deg!!!!
-; l'interpolation ajoute des points parasites dans la zone...
-         nbgrid=0
-         cont=0L
-         chint=fltarr(1)
-         for i=0L,nzon-1L do begin
-            tb=ch_nadir[*,i]
-            lon=amalon[*,i]
-            lat=amalat[*,i]
-            mask=landseamask[*,i]
-            ind=where(tb gt tbmin,nbon)
-; test pour eviter pb d'interpolation de longitude (meridien 180)
-            if nbon eq nbpix then begin
-               cont=cont+1L
-               nbgrid1=nbgrid
-               interpolswath,tb,lat,lon,mask,resol,nbgrid,tbgrid,latgrid,longrid,maskgrid                  
-               fovgrid=indgen(n_elements(tbgrid))+1
-;print,'wwwtest: nb points fov interpole ',n_elements(tbgrid)
-; 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) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt)
-;print,'wwwtest: nb points zone d interet ',npt
-               if npt ne 0 then begin
-                  if (n_elements(chint) LE 1) then begin
-                     timeint=replicate(tt[i],npt)
-                     chint=tbgrid[zone]
-                     latint=latgrid[zone]
-                     lonint=longrid[zone]
-                     fovint=fovgrid[zone]
-                     maskint=maskgrid[zone]
-                  endif else begin
-                     chint=[chint,tbgrid[zone]]
-                     latint=[latint,latgrid[zone]]
-                     lonint=[lonint,longrid[zone]]
-                     fovint=[fovint,fovgrid[zone]]
-                     maskint=[maskint,maskgrid[zone]]
-                     timeint=[timeint,replicate(tt[i],npt)]
-                  endelse
-               endif else begin
-            ;print, 'www : aucun point dans la zone'
-               endelse
-            endif
-         endfor
-     ; PRINT,'fin interpolation et selection des points de la zone ', SYSTIME()
-
-         nn=n_elements(chint)
-
-         print,'www : nb points dans la zone en fin de traitement de l''orbite ',nn
-         if nn gt 1 then begin
-            for i=0L, nn-1L do begin
-               if (n_elements(data) EQ 0) then begin
-                  data = { desc : desc $
-                           , hour : timeint[i] $
-                           , fov : fovint[i] $
-                           , lon : lonint[i] $
-                           , lat: latint[i] $
-                           , landseamask: maskint[i] $
-                           , tb : chint[i] $
-                         }
-               endif else begin
-                  data = [data , { desc : desc $
-                                   , hour : timeint[i] $
-                                   , fov : fovint[i] $
-                                   , lon : lonint[i] $
-                                   , lat: latint[i] $
-                                   , landseamask:  maskint[i] $
-                                   , tb : chint[i] $
-                                 }]
-               endelse
-            endfor
-         endif
-;
-; fin boucle sur les fichiers lus
-   endif
-   PRINT,'www : passage au fichier suivant ', SYSTIME()
-endfor
-labfile:
-if (n_elements(data) NE 0) then begin
-    ; write outputfile
-    header1 = { yyyy : yyyy $
-             ,  mm : mm $
-             ,  dd : dd $
-             ,  numch : numch $
-             ,  nbpix : nbgrid $
-             ,  resol : resol $
-             }
-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=''
-    fileout = mem_to_file_amsu_t2(amsu_t2, testfilename)
-
-endif else begin
-    print, 'www : no data to write'
-    goto, realend
-endelse
-
-nofile:
-    print, 'www : no data to read'
-    goto, realend
-
-realend:
-    PRINT, 'www : fin programme extract_amsua ',SYSTIME()
+result = 0
+return, result
 end

trunk/src/traite_amsuab.sh

@@ -141 +141 @@
 #
 # $URL$
+#
+# - pinsard 20120416
+#
+#   * update min max yyyymmdd
 #
 # - pinsard 20120323
@@ -265 +269 @@
 hostname=$(hostname)
 #
-yyyymmddb_min=20060420
-yyyymmdde_max=20061103
+yyyymmddb_min=20000101
+yyyymmdde_max=20121231
 #
 # default
@@ -476 +480 @@
 ;
 files_list = search_amsufiles('${numch}', ${yyyy},${mm}, ${dd})
+; ++ test si files_list non vide
+;
+result = file_bathy_to_mem(${lonmin}, ${lonmax}, ${latmin}, ${latmax})
+; ++ test result ok
 ;
 ; uncomment the following lines to test with the last file of the list