Ignore:
Timestamp:
05/11/06 12:35:53 (18 years ago)
Author:
smasson
Message:

debug + new xxx

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ToBeReviewed/CALENDRIER/caldat.pro

    r9 r69  
    113113      year = TEMPORARY(year) - (month GT 2) 
    114114      year = year - (year LE 0) 
    115  
    116115; see if we need to do hours, minutes, seconds 
    117116      IF (nParam GT 4) THEN BEGIN 
     
    176175    END  
    177176    'noleap':BEGIN 
     177 
     178      jul = long(julian) 
     179      year = jul/365 + 1 
     180      day = jul MOD 365L 
     181; 
     182; 
     183      zero = where(day EQ 0) 
     184 ; 
     185      month = 1 + (day GT 31) + (day GT 59) + (day GT 90) + (day GT 120) $ 
     186              + (day GT 151) + (day GT 181) + (day GT 212) + (day GT 243) $ 
     187              + (day GT 273) + (day GT 304) + (day GT 334) 
     188      month = long(month) 
     189;  
     190      day = day - 31L * (day GT 31) - 28L * (day GT 59) - 31L * (day GT 90) $ 
     191              - 30L * (day GT 120) - 31L * (day GT 151) - 30L * (day GT 181) $ 
     192              - 31L * (day GT 212) - 31L * (day GT 243) - 30L * (day GT 273) $ 
     193              - 31L * (day GT 304) - 30L * (day GT 334) 
     194; 
     195      IF zero[0] NE -1 THEN BEGIN 
     196        year[zero] = year[zero]-1 
     197        month[zero] = 12L 
     198        day[zero] = 31L 
     199      ENDIF  
     200 
    178201    END  
    179202    ELSE:BEGIN 
     
    183206  ENDCASE 
    184207; 
     208  zero = where(year ge 600000L, cnt) 
     209  IF cnt NE 0 THEN year[zero] = year[zero]-654321L 
     210; 
    185211  return 
    186212 
Note: See TracChangeset for help on using the changeset viewer.