Changeset 451 for trunk

Timestamp:

12/09/11 15:23:43 (13 years ago)

Author:

pinsard

Message:

transition to use writing tool for AMSU

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (17 diffs)
• mem_to_file_amsu_t2.pro (modified) (3 diffs)

trunk/src/extract_amsua.pro

@@ -15 +15 @@
 ; ecrit les donnees pour la zone (lat/long) choisie
 ;
-; le fichier de sortie contient dans l'ordre mm, jj (jour de
+; ++ 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,
-; tb1,... tbn
+; no satellite, le no de fov (1 a 40 ou 1 a 90), long, lat, tb
 ;
 ; Les valeurs manquantes sont codées par NaN.
@@ -34 +33 @@
 ;           amsu_lista [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/list_filea"];
 ;           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/amsua_yyyymmdd_geomin_geomax.dat"];
+;           amsu_t2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax.dat"];
 ;
 ;           extract_amsua [shape=box,
@@ -42 +41 @@
 ;           label="${PROJECT}/src/extract_amsuab.pro"];
 ;
-;           {amsu_lista amsu_listb amsu_t1} -> {extract_amsuab} -> {amsu_t2};
+;           {amsu_lista amsu_t1} -> {extract_amsuab} -> {amsu_t2};
 ;       }
 ;
@@ -48 +47 @@
 ;
 ; :Params:
+;   numch
 ;   yyyy
 ;   mm
@@ -57 +57 @@
 ; :Examples:
 ;
-;  Creating an articial list of AMSU-B and AMSU-A files (You may need to create ${PROJECT_ID}/AMSU/2006/08/ directory) ::
+;  Creating an articial list of AMSU-A files (You may need to create ${PROJECT_ID}/AMSU/2006/08/ directory) ::
 ;
 ;   $ mkdir -p ${PROJECT_ID}/AMSU/2006/08
-;   $ echo "${PROJECT_ID}/AMSU/AMSUBN15/L1C/2006/2006_08_01/NSS.AMBX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_fileb
-;   $ echo "${PROJECT_ID}/AMSU/AMSUBN15/L1C/2006/2006_08_01/NSS.AMBX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_fileb
 ;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_filea
 ;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_filea
 ;
-; Using AMSU-A:
-;
+; Using AMSU-A channel 5::
+;
+;   IDL> numch = 'a5'
 ;   IDL> yyyy=2006
 ;   IDL> mm=8
@@ -74 +73 @@
 ;   IDL> lat_min=-30.
 ;   IDL> lat_max=45.
-;   IDL> extract_amsub, yyyy, mm, dd, lon_min, lon_max, lat_min, lat_max
-;
-;   ${PROJECT_ID}/AMSU/2006/08/amsua_20060801_060w30s_050e45n.dat must have been created.
+;   IDL> extract_amsua, numch, yyyy, mm, dd, lon_min, lon_max, lat_min, lat_max
+;
+;   ${PROJECT_ID}/AMSU/2006/08/a5_20060801_060w30s_050e45n.dat must have been created.
 ;
 ; :History:
@@ -87 +86 @@
 ; amelioree
 ;
+; - 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!
+;
+;   * 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 
+;   * need to check if missing exist before NaN affection on idl7 cratos
 ;
 ; - pinsard 2011-08-18T10:34:53Z loholt1.ipsl.polytechnique.fr (Linux)
@@ -149 +156 @@
 ; :Todo:
 ;
+; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
+; crains des zesors louches!!).louchse
+;
+; 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 !)
+;
 ; changer le nom des fichiers de listes pour éviter des conflits lors de runs
 ; simultanés
@@ -156 +171 @@
 ; eviter la creation d'un fichier vide (si pas de points dans la zone)
 ;
-; tester la cohérence entre les fichiers présents dans list_filea et list_fileb
+; tester la cohérence entre les fichiers présents dans list_filea
 ; en terme de date avec yyyy, mm et dd donnés en paramètre
 ;
@@ -163 +178 @@
 ; get rid of uppercase
 ;
-;
 ;-
-PRO extract_amsua, yyyy, mm, dd,lon_min, lon_max, lat_min, lat_max
+PRO extract_amsua, numch, yyyy, mm, dd,lon_min, lon_max, lat_min, lat_max
 ;
 compile_opt idl2, strictarrsubs
@@ -177 +191 @@
 
 
-prefix='amsua_'
 look = 'filename'
 scale = 1.
@@ -185 +198 @@
         + string(yyyy,format='(I4.4)') + '/' $
         + string(mm,format='(I2.2)')   + '/' $
-        + prefix+string(yyyy,format='(I4.4)') $
+        + numch + '_' $
+        + string(yyyy,format='(I4.4)') $
         + string(mm,format='(I2.2)') $
         + string(dd,format='(I2.2)') + '_' $
@@ -199 +213 @@
 ; boucle sur les elements de la liste
 
-; lecture de la liste des fichiers AMSU-A 
+; lecture de la liste des fichiers AMSU-A
     ;list_filea = project_id_env+'list_filea'
      list_filea='list_filea'
@@ -313 +327 @@
   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
   if npt ne 0 then begin
      time[zone]=tt[i]
@@ -320 +335 @@
      fovint=[fovint,fovgrid[zone]]
      timeint=[timeint,time[zone]]
-  endif
+  endif else begin
+     print, 'www : aucun point dans la zone'
+  endelse
 
  endfor
 
+ if (npt ne 0) then begin
  ; partie a remplacer par module ecriture standardise FP
 ;
@@ -336 +354 @@
 print,'nb points dans la zone',nn
 
-      for i=0L, nn-1L do begin     
+      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))'
+         if (n_elements(data) EQ 0) then begin
+             data = { desc : desc $
+                    , hour : timeint[i] $
+                    , fov : fovint[i] $
+                    , lon : lonint[i] $
+                    , lat: latint[i] $
+                    , landseamask: 9$
+                    , tb : chint[i] $
+                    }
+         endif else begin
+             data = [data , { desc : desc $
+                          , hour : timeint[i] $
+                          , fov : fovint[i] $
+                          , lon : lonint[i] $
+                          , lat: latint[i] $
+                          , landseamask: 9$
+                          , tb : chint[i] $
+                       }]
+         endelse
      endfor
-      
-
+  endif else begin
+     print, 'www : 2 aucun point dans la zone'
+  endelse
 
 ; fin boucle sur les fichiers lus
@@ -350 +388 @@
 print, 'fermeture de ', fileout
 close, /all
+
+if (n_elements(data) NE 0) then begin
+    ; write outputfile
+    header = { yyyy : yyyy $
+             ,  mm : mm $
+             ,  dd : dd $
+             ,  numch : numch $
+             }
+    amsu_t2={header: header, 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)
+
+endif else begin
+    print, 'www : no data to write'
+endelse
+;
 end
-

trunk/src/mem_to_file_amsu_t2.pro

@@ -15 +15 @@
 ; amsu_t2 terminology will be applied.
 ;
-; :returns: 0 if ok or -1 if error
+; :returns: name of the output file if ok or -1 if error
 ; :rtype: integer
 ;
@@ -70 +70 @@
 ;
 ; $URL$
+;
+; - fplod 20111208T173432Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * name of output file as result if ok
 ;
 ; - fplod 20111208T105205Z cratos (Linux)
@@ -148 +152 @@
 FREE_LUN, lun
 ;
-result = 0
+result = fullfilename
 ;
 return, result