Changeset 548

Timestamp:

05/04/12 14:01:42 (12 years ago)

Author:

lelod

Message:

merge version 543 et mine

File:

• trunk/src/cresamsu.pro (modified) (14 diffs)

trunk/src/cresamsu.pro

@@ -1 +1 @@
 ; docformat = 'rst'
 ;+
-; .. _cresamsu.pro:
 ;
 ; ============
@@ -7 +6 @@
 ; ============
 ;
+; .. function: cresamsu(numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax, temp, rtemp, rhdeg)
 ; DESCRIPTION
 ; ===========
@@ -12 +12 @@
 ; programme d'interpolation avec algorithme de Cressmann
 ;
+;; grille de sortie en long / lat
+; recoit en entree les jours de l'annee jdeb et jfin
+;correspondant au debut et la fin du mois
+;
 ; adaptation aux donnees Tb AMSU lues avec la nouvelle chaîne de lecture
-; extract_amsuab.pro
+; extract_amsu.pro
 ; correct_nadir.pro
 ;
 ; but: faire des fichiers "grilles" pour analyse statistique avec
 ; interpolation spatiale et temporelle
+;
+; .. only:: man
+;
+;    Figure is visible on PDF and HTML documents only.
+;
+; .. only:: html or latex
+;
+;     .. graphviz::
+;
+;        digraph cresamsu {
+;           amsu_t2 [shape=ellipse,
+;               fontname=Courier,
+;               label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax.dat"];
+;           amsu_t5 [shape=ellipse,
+;               fontname=Courier,
+;               label="${PROJECT_OD}/AMSU/numch_temprtrtemprhrhdeg_yyyymmdd_yyyymmdd_geomin_geomax.nc"];
+;           cresamsu [shape=box,
+;           fontname=Courier,
+;           color=blue,
+;           URL="http://forge.ipsl.jussieu.fr/varamma/browser/trunk/src/cresamsu.pro",
+;           label="${PROJECT}/src/cresamsua.pro"];
+;
+;           {amsu_t2} -> {cresamsua} -> {amsu_t5};
+;       }
 ;
 ; EXAMPLES
@@ -38 +66 @@
 ;    rtemp=0.5
 ;    rhdeg=1.
-;    cresamsu, numch, yyyyb, mmb, ddb $
-;                        , yyyye, mme, dde $
-;                        , lonmin, lonmax, latmin, latmax,temp,rtemp,rhdeg
+;    result = cresamsu(numch, yyyyb, mmb, ddb $
+;            , yyyye, mme, dde $
+;            , lonmin, lonmax, latmin, latmax $
+;            , temp, rtemp, rhdeg);
 ;
 ; SEE ALSO
 ; ========
 ;
+;  Previous step :ref:`extract_amsu.pro` (:ref:`traite_amsuab.sh`)
 ; :ref:`traite_amsuab.sh`
 ;
-; :func:`file_amsu_t2_to_mem`
+; Use :func:`file_amsu_t2_to_mem
+; :func:`idl_amsu_netcdf`
 ;
 ; TODO
 ; ====
 ;
-; add temp,rtemp,rhdeg in example
 ;
 ; fix pb with fov_int2 and fov_int3
@@ -60 +90 @@
 ; aligner lal/lon de l'exemple sur les autres modules amsu
 ;
-; lelod 2011/10/24
-; actuellement codes en dur
-; passer le choix du rayon de Cressmann en parametres
-; rayon temporel rtemp (en jour decimal)
-; rayon horizontal rhor (km)
-; et la resolution
-; pash (horizontale, en degres)
-; past (temporelle, en jour decimal)
-; et enfin l'heure d'interpolation heurinterp (en heures)
-;
+; 
 ; coding rules
 ;
@@ -79 +100 @@
 ; ==========
 ;
-;  $Id$
+;  $Id: 
 ;
 ;  $URL$
 ;
-; - lelod 20120410
-;
+; - fplod 20120430
+;
+;   * pro -> function
+;   * idl_amsu_netcdf is a function
+;   * add digraph
+;
+; - lelod 20120410 puis 20120503
+;
+;   * parametres de la grille dependant des rayons de cresmann
+;   * revisite de la grille temporelle
 ;   * choix du rayon de Cressmann en degres et pas en km
 ;   * ajout d'un parametre temp (= "m" ou "e") pour choix de la plage
@@ -133 +162 @@
 ;-
 
-pro cresamsu, numch, yyyyb, mmb, ddb, yyyye,mme,dde,lonmin,lonmax,latmin,latmax,temp,rtemp,rhdeg
-
-; les parametres de l'interpolation
-
-; grille de sortie en long / lat
-; recoit en entree les jours de l'annee jdeb et jfin correspondant au debut et la fin du mois
-;=============================================================================================
+function cresamsu $
+   , numch $
+   , yyyyb $
+   , mmb $
+   , ddb $
+   , yyyye $
+   , mme $
+   , dde $
+   , lonmin $
+   , lonmax $
+   , latmin $
+   , latmax $
+   , temp $
+   , rtemp $
+   , rhdeg
 
 compile_opt idl2, strictarrsubs
@@ -150 +187 @@
 if key_performance EQ 1 THEN BEGIN
     time1 = SYSTIME(1)
+ENDIF
+;
+;
+; Return to caller if errors
+ON_ERROR, 2
+result = -1
+;
+usage = 'result = cresamsu(numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax, temp, rtemp, rhdeg)'
+;
+nparam = N_PARAMS()
+IF (nparam NE 14) THEN BEGIN
+    ras = report(['Incorrect number of arguments.' $
+          + '!C' $
+          + 'Usage : ' + usage])
+    return, result
 ENDIF
 ;
@@ -210 +262 @@
 hrdeb=jdeb*24
 hrfin=jfin*24
-print,'debut et fin', jdeb,jfin
+msg = report(['iii : ' + routine + ' : debut et fin = ' + string(jdeb) + ' ' + string(jfin)])
+
 jmax=fix((hrfin-hrdeb)/pasthor) ; valeurs de 0 a jmax
 nbjour=jmax+1
@@ -233 +286 @@
 lmax=fix((la2-la1)/pash) ; valeurs de 0 a lmax
 nblat=lmax+1
-print,'nb pas de temps',nbjour, ' nb lat',nblat,' nblon',nblon
-
+msg = report(['iii : ' + routine + ' : nb pas de temps = ' + string(nbjour)])
+msg = report(['iii : ' + routine + ' : nb lat = ' + string(nblat)])
+msg = report(['iii : ' + routine + ' : nb lon = ' + string(nblon)])
+msg = report(['iii : ' + routine + ' : rayon temporel en heures' + string(rtemphor)])
 
 
@@ -249 +304 @@
 la2=max(yylat)
 hrfin=max(jjours)*24.
-print, 'pas de temps du fichier ncdf en jours decimaux',jjours
-print,'derniers points de la grille',max(xxlon),max(yylat),hrfin
-;-----------------------------------
+msg = report(['iii : ' + routine + ' : pas de temps du fichier ncdf en jours decimaux= ' + string(jjours)])
+msg = report(['iii : ' + routine + ' : derniers points de la grille (lon) = ' + string(max(xxlon))])
+msg = report(['iii : ' + routine + ' : derniers points de la grille (lat) = ' + string(max(yylat))])
+msg = report(['iii : ' + routine + ' : derniers points de la grille (temps) = ' + string(max(hrfin))]);-----------------------------------
 t=indgen(nbjour)*pasthor+hrdeb ; en heures
 print,'t',t
@@ -372 +428 @@
                endif
             endif
-            ; fin test sur donnees valides (inutil normalement)
+            ; fin test sur donnees valides (inutile normalement)
          endfor
          ; fin lecture donnees fichier journalier
@@ -394 +450 @@
 if nbfile GT 0 THEN BEGIN
     good=where (kont1 ge 1,nb)
-    print,'points valides',nb
+    msg = report(['iii : ' + routine + ' : points valides = ' + string(nb)])
     if (nb ne 0) then tb1fin[good]=stb1[good]/z1[good]
     bad=where (kont1 eq 0,nbad)
@@ -400 +456 @@
     if (nbad ne 0) then tb1fin[bad]=!values.f_Nan
                                     ; print,tb1fin
-    print,'extremes du compteur',min(kont1),max(kont1)
-    
-    ;toto=reform(tb1fin[*,*,0])
-                                ;plt,xxlon,yylat,toto
+    msg = report(['iii : ' + routine + ' : extremes du compteur = ' + string(min(kont1)) + ' ' + string(max(kont1))])
+
     result=idl_amsu_netcdf(numch, temp,rtemp,rhdeg, tb1fin, yylat, xxlon, jjours)
+    result_type=size(result,/type)
+    IF (result_type EQ 3) THEN BEGIN
+       IF (result NE 0) THEN BEGIN
+           msg = report(['eee : ' + routine + ' : pb idl_amsu_netcdf'])
+           ;return, result
+           STOP
+       ENDIF
+    ENDIF
 ENDIF  ELSE BEGIN
-    print, 'www : no files at all'
-ENDELSE
-
+      msg = report(['www : ' + routine + 'no files at all'])
+    goto, realend
+ ENDELSE
+
+realend:
 IF key_performance EQ 1 THEN BEGIN
-     PRINT, 'ppp : durée totale cresamsu', SYSTIME(1)-time1, format='(a, 1x, F12.6)'
+    msg = report(['ppp : ' + routine $
+         + ' : durée totale = ' $
+         + string(SYSTIME(1)-time1, format='(a, 1x, F12.6)')])
 ENDIF
 
-print,'fin cresamsu'
-
+msg = report(['iii : ' + routine + ': fin'])
+        
+result = 0
+return, result  
 end