Changeset 473 for trunk

Timestamp:

12/16/11 17:54:42 (12 years ago)

Author:

pinsard

Message:

complete header and try to fix amsu tools

Location:

trunk/src

Files:

• cresamsu.pro (modified) (7 diffs, 1 prop)
• idl_amsu_netcdf.pro (modified) (11 diffs)
• test_correc.pro (modified) (5 diffs, 1 prop)

trunk/src/cresamsu.pro

@@ -1 +1 @@
 ; docformat = 'rst'
 ;+
-; .. _cresamsu_cdf.pro:
-;
-; ======================
-; cresamsu_cdf.pro
-; ======================
+; .. _cresamsu.pro:
+;
+; ============
+; cresamsu.pro
+; ============
+;
+; ===========
+; DESCRIPTION
+; ===========
 ;
 ; programme d'interpolation avec algorithme de Cressmann
-; adaptation aux donnees Tb AMSU lues avecla nouvelle chaîne de lecture
+;
+; adaptation aux donnees Tb AMSU lues avec la nouvelle chaîne de lecture
 ; extract_amsuab.pro
 ; correct_nadir.pro
+;
 ; but: faire des fichiers "grilles" pour analyse statistique avec
 ; interpolation spatiale et temporelle
 ;
-; Uses:
-; parametres:
-;  use_amsua = 0 ou 1 (amsub ou les deux)
-;  listcanal= tableau des canaux a traiter (caracteres) ex ['a5', 'b5']
-;  cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lon_min,lon_max, lat_min, lat_max,listcanal
-; exemple
-;IDL> listcanal=['a5']
-;IDL> use_amsua=1
-;IDL> yyyyb=2006
-;IDL> yyyye=2006
-;IDL> mmb=08
-;IDL> mme=8
-;IDL> ddb=1
-;IDL> dde=3
-;IDL> lon_min=-20
-;IDL> lon_max=40
-;IDL> lat_min=0
-;IDL> lat_max=30
-;  cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lon_min,lon_max, lat_min, lat_max,listcanal
-
-; :Restrictions:
-;
-; :Pre:
-;
-; :Post:
-;
-; :Also:
-;
-;  :ref:`traite_amsuab.sh`
-;
-; :Todo:
+; EXAMPLES
+; ========
+;
+; ::
+;
+;    IDL> numch='a5'
+;    IDL> yyyyb=2006L
+;    IDL> yyyye=2006L
+;    IDL> mmb=8
+;    IDL> mme=8
+;    IDL> ddb=1
+;    IDL> dde=3
+;    IDL> lon_min=-20.
+;    IDL> lon_max=40.
+;    IDL> lat_min=0.
+;    IDL> lat_max=30.
+;    IDL> cresamsu, numch, yyyyb, mmb, ddb $
+;                        , yyyye, mme, dde $
+;                        , lon_min, lon_max, lat_min, lat_max
+;
+; SEE ALSO
+; ========
+;
+; :ref:`traite_amsuab.sh`
+;
+; :func:`file_amsu_t2_to_mem`
+;
+; TODO
+; ====
+;
+; check last line for uncomment endfor (not sure of the use of this line)
 ;
 ; lelod 2011/12/15
@@ -58 +64 @@
 ; past (temporelle, en jour decimal)
 ; et enfin l'heure d'interpolation heurinterp (en heures)
-;
 ; 
 ; coding rules
@@ -64 +69 @@
 ; new filename convention
 ;
-; :Examples:
-;
-; ::
-;
-; :History:
-;
-; 20121215 lelod
-; reprise du programme apres refonte de la chaine (extract_amsua):
-; lecture nouveau format de fichier (structure)
-; elimination des tentatives precedentes
-; ajout du mask dans l'appel a idl_amsu_ncdf
-; introduire en parametres les choix pour le fichier ncdf: grille,
-; rayons de Cressamann, selection de la portion du jour utilisee dans l'interpolation
-;
-; 2011/10/24 lelod
-; prise en compte des modifications de lecture et correction nadir
-; (correct_nadir.pro):
-; possibilite de traiter amsua et amsub ou amsub seul
-; utilisation de pixelsize.pro pour boucher les trous entres lignes de
-; pixels en bordure de fauchee de facon correcte
-;
-; 
-;  $Id: cresamsu_cdf.pro 435 2011-11-24 22:18:41Z lelod $
-;  $URL: svn+ssh://lelod@forge.ipsl.jussieu.fr/ipsl/forge/projets/varamma/svn/trunk/src/cresamsu_cdf.pro $
+; get rid of idldoc (fplod)
+;
+; EVOLUTIONS
+; ==========
+;
+;  $Id$
+;
+;  $URL$
+;
+; - fplod 20111216T164041Z cratos (Linux)
+;
+;   * fix file_amsu_t2_to_mem call
+;   * start to forgot idldoc
+;   * fix for loop syntax error
+;
+; - lelod 20121215 
+;
+;   * reprise du programme apres refonte de la chaine (extract_amsua):
+;     lecture nouveau format de fichier (structure)
+;     elimination des tentatives precedentes
+;     ajout du mask dans l'appel a idl_amsu_ncdf
+;     introduire en parametres les choix pour le fichier ncdf: grille,
+;     rayons de Cressamann, selection de la portion du jour utilisee dans l'interpolation
+;
+; - lelod 20111024
+;
+;   * prise en compte des modifications de lecture et correction nadir
+;     (correct_nadir.pro):
+;     possibilite de traiter amsua et amsub ou amsub seul
+;     utilisation de pixelsize.pro pour boucher les trous entres lignes de
+;     pixels en bordure de fauchee de facon correcte
+;
 ;
 ; - fplod 20110505T123230Z cratos.locean-ipsl.upmc.fr (Linux)
@@ -97 +110 @@
 ;-
 
-pro cresamsu_cdf,numch,yyyyb,mmb,ddb,yyyye,mme,dde,lon_min, lon_max, lat_min, lat_max
+pro cresamsu,numch,yyyyb,mmb,ddb,yyyye,mme,dde,lon_min, lon_max, lat_min, lat_max
 
 ; les parametres de l'interpolation
@@ -207 +220 @@
          for dd=ddb,dde do begin
 ; ouverture fichier journalier
-            result=file_amsu_t2_to_mem, yyyy, mm, dd, numch $
+            result=file_amsu_t2_to_mem( yyyy, mm, dd, numch $
                                  , lon_min, lon_max, lat_min, lat_max $
-                                 , testfilename
+                                 , testfilename)
             desc=result.data[0]
             hour=result.data[1]
@@ -223 +236 @@
 ; boucle sur les points du fichier
 
-            do i=0,nn do begin
+            for i=0,nn do begin
 
 ; recup de tb1,jour,hour,lon,lat
@@ -303 +316 @@
                                 ;plt,xxlon,yylat,toto
 idl_amsu_netcdf, numch, tb1fin, yylat, xxlon, jjours, mask
-endfor
+;+++ c'est la fin de quelle boucle
+;+++endfor
 
 

trunk/src/idl_amsu_netcdf.pro

@@ -27 +27 @@
 ;           lats [shape=tab, label="latitude 1d array"]
 ;           timesin [shape=tab, label="time 1d array"]
+;           maskin [shape=tab, label="landseamask 3d array lon, lat, time"]
 ;
 ;           amsu_t5 [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/AMSU/numch_yyyymmdd_yyyymmdd_geomin_geomax.nc"];
@@ -36 +37 @@
 ;           label="${PROJECT}/src/idl_amsu_netcdf.pro"];
 ;
-;           {datain lons lats timesin} -> {idl_amsu_netcdf} -> {amsu_t5}
+;           {datain lons lats timesin maskin} -> {idl_amsu_netcdf} -> {amsu_t5}
 ;
 ;          }
@@ -50 +51 @@
 ;   IDL> timesin = indgen(12) + julday(8,1,2006L)
 ;   IDL> datain = findgen(size(lons,/dim), size(lats,/dim), size(timesin,/dim))
-;   IDL> idl_amsu_netcdf, numch , datain, lats, lons, timesin
+;   IDL> maskin = findgen(size(lons,/dim), size(lats,/dim), size(timesin,/dim))
+;   IDL> idl_amsu_netcdf, numch , datain, lats, lons, timesin, maskin
 ;
 ; SEE ALSO
@@ -64 +66 @@
 ; ====
 ;
+; type of maskin and landseamask ? why fltarr
+; 
+; dimension on maksin : why time ?
+;
 ; check arg : time in julian days. check for noon !!
 ;
@@ -73 +79 @@
 ;
 ; handle missing values properly
+; 
+; handle tb threshold properly
 ;
 ; check time in netcdf
@@ -86 +94 @@
 ;
 ; $URL$
+;
+; - fplod 20111216T160934Z cratos (Linux)
+;
+;   * add maskin in the example and graph
+;   * add a print to see filename on log
+;   * fix landseamask unit
+;
+; - lelod 20111215
+;
+;   * add maskin in input parameters
 ;
 ; - fplod 20111130T111832Z cratos (Linux)
@@ -120 +138 @@
 ;
 ; build output filename
-prefix=''
-; nom du canal traité par interpolation de Cresmann
+; avec le nom et le numéro du canal traité par interpolation de Cresmann
 prefix = numch +'_'
 
@@ -161 +178 @@
 ; --------------------------------------
 ;
+print, 'ouverture pour ecriture de ', cdfout
 id = ncdf_create(cdfout,/clobber)
 ;++ncdf_control, id, /nofill
@@ -193 +211 @@
 ncdf_attput, id, varid[3], 'unit', 'degree_kelvin'
 ncdf_attput, id, varid[4], 'short_name', 'mask'
-ncdf_attput, id, varid[3], 'unit', 'number (0 ocean, 1 land)'
+ncdf_attput, id, varid[4], 'unit', 'number (0 ocean, 1 land)'
 ;
 ; remplissage de l'entete
@@ -201 +219 @@
 ncdf_varput, id, varid[2], jours
 ;
-; remplissage de la variable
+; remplissage des variables
 ncdf_varput, id, varid[3], moyenne_tb
 ncdf_varput, id, varid[4], landseamask
@@ -208 +226 @@
 ncdf_close, id
 ;
-return
 end

trunk/src/test_correc.pro

@@ +1 @@
+;+
+;
+; .. _test_correc.pro
+;
+; ===============
+; test_correc.pro
+; ===============
+;
+; DESCRIPTION
+; ===========
+;
+; permet de comparer les fichiers pour chaque canal
+;
+; EXAMPLES
+; ========
+;
+;   IDL> numch = 'a5'
+;   IDL> test_correc, numch
+;
+; SEE ALSO
+; ========
+;
+; :ref:`data_amsu`
+;
+; EVOLUTIONS
+; ==========
+;
+; $Id$
+;
+; $URL$
+;
+; - fplod 20111216T160346Z cratos (Linux)
+;
+;   * complete header
+;
 ;-
+
 PRO test_correc,numch
 ;
@@ -118 +154 @@
          t26,t27,t28,t29,t30]
    endfor
-   cor_l2006=a 
+   cor_l2006=a
    free_lun, lun1
    openr, lun1,cor_seab, /get_lun
@@ -135 +171 @@
          t26,t27,t28,t29,t30]
    endfor
-   cor_s2006=a 
+   cor_s2006=a
 
    cor_landb=path+'AMSU/CORR_LAND_AMSUB.DAT'
@@ -155 +191 @@
          t26,t27,t28,t29,t30]
    endfor
-   cor_lold=a 
+   cor_lold=a
 
    free_lun, lun1
@@ -173 +209 @@
          t26,t27,t28,t29,t30]
    endfor
-   cor_sold=a 
+   cor_sold=a
 endif
 window,10