Changeset 511

Timestamp:

04/09/12 22:51:28 (12 years ago)

Author:

lelod

Message:

correction bug boucle sur les dates

File:

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

trunk/src/cresamsu.pro

@@ -150 +150 @@
 
 ; rayons de cressmann
-rhor=120. ; en gros deux fois la resol de AMSUA
+
 rtemp=0.5 ;en jours rayon de Cressmann (0.25, 0.5, 1, 2,... - nb entiers!)
-r2=rhor*rhor
 rtemphor=rtemp*24. ; en heures
 print,'rayon temporel en heures',rtemphor
@@ -159 +158 @@
 coef=!pi/180.
 rterre=6400.
-rhdeg=rhor*180./(rterre*!pi)    ;rayon de Cressmann en degres
+;rhor=120. ; en gros deux fois la resol de AMSUA
+rhdeg=1.   ; en degres
+;rhdeg=rhor*180./(rterre*!pi)    ;rayon de Cressmann en degres
+rhor=rhdeg*(rterre*!pi)/180.
+r2=rhor*rhor
 
 ; debut et fin de la grille temporelle, a partir des fichiers d'entree
@@ -230 +233 @@
 for yyyy=yyyyb,yyyye do begin
    for mm=mmb,mme do begin
-      for dd=ddb,dde do begin
+      if (mme ne mmb) then begin
+         ddmb=1
+         ddme=31
+         if mm eq mmb then ddmb=ddb
+         if mm eq mme then ddme=dde
+      endif
+      for dd=ddmb,ddme do begin
          print, 'iii : traitement du jour ',yyyy,mm,dd
          ; lecture du fichier journalier