--- trunk/phylmd/readsulfate.f 2015/02/13 18:22:38 129 +++ trunk/phylmd/readsulfate.f 2015/02/24 15:43:51 130 @@ -4,7 +4,7 @@ contains - SUBROUTINE readsulfate(r_day, first, sulfate) + SUBROUTINE readsulfate(dayvrai, time, first, sulfate) ! From LMDZ4/libf/phylmd/readsulfate.F, version 1.2 2005/05/19 ! 08:27:15 fairhead @@ -33,7 +33,11 @@ ! Input: - real, intent(in):: r_day ! Day of integration + integer, intent(in):: dayvrai + ! current day number, based at value 1 on January 1st of annee_ref + + REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour + LOGICAL, intent(in):: first ! First timestep ! (and therefore initialization necessary) @@ -65,7 +69,7 @@ !-------------------------------------------------------------------- - iday = INT(r_day) + iday = dayvrai ! Get the year of the run iyr = iday/360 @@ -74,7 +78,7 @@ iday = iday-iyr*360 ! 0.02 is about 0.5/24, namly less than half an hour - lnewday = (r_day-FLOAT(iday).LT.0.02) + lnewday = time < 0.02 ! All has to be done only, if a new day begins!