4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE readsulfate_preind(r_day, first, pi_sulfate) |
SUBROUTINE readsulfate_preind(dayvrai, time, first, pi_sulfate) |
8 |
|
|
9 |
! Read in /calculate pre-industrial values of sulfate |
! Read in /calculate pre-industrial values of sulfate |
10 |
|
|
32 |
! |
! |
33 |
! Input: |
! Input: |
34 |
! ------ |
! ------ |
35 |
real, intent(in):: r_day ! Day of integration |
integer, intent(in):: dayvrai |
36 |
|
! current day number, based at value 1 on January 1st of annee_ref |
37 |
|
|
38 |
|
REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour |
39 |
|
|
40 |
LOGICAL, intent(in):: first ! First timestep |
LOGICAL, intent(in):: first ! First timestep |
41 |
! (and therefore initialization necessary) |
! (and therefore initialization necessary) |
42 |
! |
! |
64 |
|
|
65 |
|
|
66 |
|
|
67 |
iday = INT(r_day) |
iday = dayvrai |
68 |
|
|
69 |
! Get the year of the run |
! Get the year of the run |
70 |
iyr = iday/360 |
iyr = iday/360 |
73 |
iday = iday-iyr*360 |
iday = iday-iyr*360 |
74 |
|
|
75 |
! 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 |
76 |
lnewday = (r_day-FLOAT(iday).LT.0.02) |
lnewday = time < 0.02 |
77 |
|
|
78 |
! --------------------------------------------- |
! --------------------------------------------- |
79 |
! All has to be done only, if a new day begins! |
! All has to be done only, if a new day begins! |