4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE readsulfate(r_day, first, sulfate) |
SUBROUTINE readsulfate(dayvrai, time, first, sulfate) |
8 |
|
|
9 |
! From LMDZ4/libf/phylmd/readsulfate.F, version 1.2 2005/05/19 |
! From LMDZ4/libf/phylmd/readsulfate.F, version 1.2 2005/05/19 |
10 |
! 08:27:15 fairhead |
! 08:27:15 fairhead |
33 |
|
|
34 |
! Input: |
! Input: |
35 |
|
|
36 |
real, intent(in):: r_day ! Day of integration |
integer, intent(in):: dayvrai |
37 |
|
! current day number, based at value 1 on January 1st of annee_ref |
38 |
|
|
39 |
|
REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour |
40 |
|
|
41 |
LOGICAL, intent(in):: first ! First timestep |
LOGICAL, intent(in):: first ! First timestep |
42 |
! (and therefore initialization necessary) |
! (and therefore initialization necessary) |
43 |
|
|
69 |
|
|
70 |
!-------------------------------------------------------------------- |
!-------------------------------------------------------------------- |
71 |
|
|
72 |
iday = INT(r_day) |
iday = dayvrai |
73 |
|
|
74 |
! Get the year of the run |
! Get the year of the run |
75 |
iyr = iday/360 |
iyr = iday/360 |
78 |
iday = iday-iyr*360 |
iday = iday-iyr*360 |
79 |
|
|
80 |
! 0.02 is about 0.5/24, namly less than half an hour |
! 0.02 is about 0.5/24, namly less than half an hour |
81 |
lnewday = (r_day-FLOAT(iday).LT.0.02) |
lnewday = time < 0.02 |
82 |
|
|
83 |
! All has to be done only, if a new day begins! |
! All has to be done only, if a new day begins! |
84 |
|
|