[4728] | 1 | ! (C) Copyright 1989- ECMWF. |
---|
| 2 | ! (C) Copyright 1989- Meteo-France. |
---|
| 3 | ! |
---|
| 4 | ! This software is licensed under the terms of the Apache Licence Version 2.0 |
---|
| 5 | ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. |
---|
| 6 | ! |
---|
| 7 | ! In applying this licence, ECMWF does not waive the privileges and immunities |
---|
| 8 | ! granted to it by virtue of its status as an intergovernmental organisation |
---|
| 9 | ! nor does it submit to any jurisdiction. |
---|
| 10 | |
---|
| 11 | ! ------------------------------------------------------------------ |
---|
| 12 | |
---|
| 13 | ! - Time functions |
---|
| 14 | ! the descriptions are in the annex 1 of the documentation |
---|
| 15 | |
---|
| 16 | ! TIME |
---|
| 17 | |
---|
| 18 | ! NDD : extraxt dd from ccaammdd |
---|
| 19 | ! NMM : extract mm from ccaammdd |
---|
| 20 | ! NAA : extract aa from ccaammdd |
---|
| 21 | ! NCCAA : extract ccaa from ccaammdd |
---|
| 22 | ! NAMD : extract aammdd from ccaammdd |
---|
| 23 | ! NCENT : return centuary of ccaammdd |
---|
| 24 | ! NYEARC: returns year of the centuary from ccaammdd |
---|
| 25 | ! NCONSTRUCT_DATE : returns ccaammdd given centuary,year,month and day |
---|
| 26 | ! NCTH : turn seconds into hours |
---|
| 27 | ! RTIME : returns the time of the model (in seconds of course!) |
---|
| 28 | ! RINCDAY : length of a day (in seconds) |
---|
| 29 | |
---|
| 30 | INTEGER(KIND=JPIM) :: NDD,NMM,NCCAA,NAA,NAMD,NCTH,NZZAA,NZZMM,NCENT,NYEARC,& |
---|
| 31 | &NCONSTRUCT_DATE |
---|
| 32 | REAL(KIND=JPRD) :: RJUDAT,RTIME |
---|
| 33 | REAL(KIND=JPRB) :: RINCDAY |
---|
| 34 | INTEGER(KIND=JPIM) :: KGRDAT,KSEC,KAAAA,KMM,KDD,KSS |
---|
| 35 | INTEGER(KIND=JPIM) :: KCENT,KYEARC,KMONTH,KDAY |
---|
| 36 | |
---|
| 37 | NDD(KGRDAT) =MOD(KGRDAT,100) |
---|
| 38 | NMM(KGRDAT) =MOD((KGRDAT-NDD(KGRDAT))/100,100) |
---|
| 39 | NCCAA(KGRDAT)=KGRDAT/10000 |
---|
| 40 | NAA(KGRDAT)=MOD(NCCAA(KGRDAT),100) |
---|
| 41 | NAMD(KGRDAT)=MOD(KGRDAT,1000000) |
---|
| 42 | NCTH(KSEC)=KSEC/3600 |
---|
| 43 | NCENT(KGRDAT)=NCCAA(KGRDAT)/100+MIN(NAA(KGRDAT),1) |
---|
| 44 | NYEARC(KGRDAT)=NAA(KGRDAT)+100*(1-MIN(NAA(KGRDAT),1)) |
---|
| 45 | NCONSTRUCT_DATE(KCENT,KYEARC,KMONTH,KDAY)=& |
---|
| 46 | &(KCENT-1)*10**6+KYEARC*10**4+KMONTH*10**2+KDAY |
---|
| 47 | |
---|
| 48 | NZZAA(KAAAA,KMM)=KAAAA-( (1-SIGN(1,KMM-3))/2 ) |
---|
| 49 | NZZMM(KMM)=KMM+6*(1-SIGN(1,KMM-3)) |
---|
| 50 | RJUDAT(KAAAA,KMM,KDD)=1720994.5_JPRD + REAL(& |
---|
| 51 | &2-NZZAA(KAAAA,KMM)/100 + (NZZAA(KAAAA,KMM)/100)/4 & |
---|
| 52 | &+ INT(365.25_JPRD*REAL(NZZAA(KAAAA,KMM),JPRD))& |
---|
| 53 | &+ INT(30.601_JPRD*REAL(NZZMM(KMM)+1,JPRD))& |
---|
| 54 | &+ KDD,JPRD) |
---|
| 55 | RTIME(KAAAA,KMM,KDD,KSS,RINCDAY)=(RJUDAT(KAAAA,KMM,KDD)-2451545._JPRD)& |
---|
| 56 | &*RINCDAY+REAL(KSS,JPRD) |
---|
| 57 | ! ------------------------------------------------------------- |
---|