Changeset 519 for trunk

Timestamp:

04/16/12 11:18:26 (12 years ago)

Author:

lelod

Message:

ajout parametres cresmann

Location:

trunk/src

Files:

• cresamsu.pro (modified) (8 diffs)
• idl_amsu_netcdf.pro (modified) (4 diffs)

trunk/src/cresamsu.pro

@@ -35 +35 @@
 ;    lat_min=-5.
 ;    lat_max=20.
+;    temp='m'
+;    rtemp=0.5
+;    rhdeg=1.
 ;    cresamsu, numch, yyyyb, mmb, ddb $
 ;                        , yyyye, mme, dde $
-;                        , lon_min, lon_max, lat_min, lat_max
+;                        , lon_min, lon_max, lat_min, lat_max,temp,rtemp,rhdeg
 ;
 ; SEE ALSO
@@ -74 +77 @@
 ; ==========
 ;
+; - lelod 20120410
+;  + choix du rayon de Cressmann en degres et pas en km
+;  + ajout d'un parametre temp (= "m" ou "e") pour choix de la plage
+;  + ajour en parametres des rayons de Cressmann rtemp et rhdeg
+;  horaire d'interpolation temporelle
+;
 ;  $Id$
 ;
@@ -121 +130 @@
 ;-
 
-pro cresamsu, 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,temp,rtemp,rhdeg
 
 ; les parametres de l'interpolation
@@ -146 +155 @@
 ; unites utilisees pour le temps: heure decimale
 pasthor=pastemp*24.
-hdeb=0.
-hfin=12. ; plage horaire pour l'interpolation des donnees
+
+if (temp eq 'e') then begin
+    hdeb=12.
+    hfin=24.          ; interpolation des donnees sur orbites apres midi
+endif
+if (temp eq 'm') then begin
+    hdeb=0.
+    hfin=12.          ; interpolation des donnees sur orbites du matin
+endif
 
 ; rayons de cressmann
 
-rtemp=0.5 ;en jours rayon de Cressmann (0.25, 0.5, 1, 2,... - nb entiers!)
+;rtemp=0.5 ;en jours rayon de Cressmann (0.25, 0.5, 1, 2,... - nb entiers!)
 rtemphor=rtemp*24. ; en heures
 print,'rayon temporel en heures',rtemphor
@@ -159 +175 @@
 rterre=6400.
 ;rhor=120. ; en gros deux fois la resol de AMSUA
-rhdeg=1.   ; en degres
+;rhdeg=1.   ; en degres
 ;rhdeg=rhor*180./(rterre*!pi)    ;rayon de Cressmann en degres
 rhor=rhdeg*(rterre*!pi)/180.
@@ -231 +247 @@
 nbfile = 0L
 ; lecture des fichiers
+ddmb=ddb
+ddme=dde
 for yyyy=yyyyb,yyyye do begin
    for mm=mmb,mme do begin
@@ -238 +256 @@
          if mm eq mmb then ddmb=ddb
          if mm eq mme then ddme=dde
-      endif
+     endif
+      print,"mois ",mm,"jours ",ddmb,"a ",ddme
       for dd=ddmb,ddme do begin
          print, 'iii : traitement du jour ',yyyy,mm,dd
@@ -345 +364 @@
     ;toto=reform(tb1fin[*,*,0])
                                 ;plt,xxlon,yylat,toto
-    idl_amsu_netcdf, numch, tb1fin, yylat, xxlon, jjours
+    idl_amsu_netcdf, numch, temp,rtemp,rhdeg, tb1fin, yylat, xxlon, jjours
 ENDIF  ELSE BEGIN
     print, 'www : no files at all'

trunk/src/idl_amsu_netcdf.pro

@@ -86 +86 @@
 ; $URL$
 ;
+; - lelod 20120410
+;  ajout parametre temp pour choix des orbites interpolees  (m pour
+;  morning, e pour evening)
+;
 ; - fplod 20120410
 ;
@@ -129 +133 @@
 ;
 ;-
-pro idl_amsu_netcdf, numch, datain, lats, lons, timesin
+pro idl_amsu_netcdf, numch, temp,rtemp,rhdeg, datain, lats, lons, timesin
 ;
 compile_opt idl2, strictarrsubs
@@ -140 +144 @@
 ;
 ; build output filename
-; avec le nom et le numéro du canal traité par interpolation de Cresmann
-prefix = numch +'_'
+; avec le nom et le numéro du canal traité par interpolation de
+; Cresmann
+rt=string(rtemp,format='(f3.1)')
+rh=string(rhdeg,format='(f3.1)')
+prefix = numch + temp + 'rt'+rt+'rh'+rh +'_'
 
 look = 'filename'
@@ -151 +158 @@
 timemax = max(timesin)
 yyyymmdde = jul2date(timemax)
-cdfout = project_id_env+ 'AMSU/'
+cdfout = project_id_env+ 'AMSU/' $
         + prefix+string(yyyymmddb,format='(I8.8)') + '_' $
         + string(yyyymmdde,format='(I8.8)') + '_' $