/[lmdze]/trunk/Sources/phylmd/readsulfate.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/readsulfate.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/readsulfate.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/phylmd/readsulfate.f revision 130 by guez, Tue Feb 24 15:43:51 2015 UTC
# Line 4  module readsulfate_m Line 4  module readsulfate_m
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
# Line 28  contains Line 28  contains
28    
29      USE dimens_m, ONLY: iim, jjm      USE dimens_m, ONLY: iim, jjm
30      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
31        USE dynetat0_m, ONLY: annee_ref
32      use getso4fromfile_m, only: getso4fromfile      use getso4fromfile_m, only: getso4fromfile
     USE temps, ONLY: annee_ref  
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    
# Line 48  contains Line 52  contains
52      INTEGER j, iday, ny, iyr, iyr1, iyr2      INTEGER j, iday, ny, iyr, iyr1, iyr2
53      parameter (ny=jjm+1)      parameter (ny=jjm+1)
54    
55      INTEGER ismaller      CHARACTER(len=4) cyear
     !JLD      INTEGER idec1, idec2 ! The two decadal data read ini  
     CHARACTER*4 cyear  
56    
57      INTEGER im, day1, day2, im2      INTEGER im, day1, day2, im2
58      double precision so4_1(iim, jjm+1, klev, 12)      double precision so4_1(iim, jjm+1, klev, 12)
# Line 67  contains Line 69  contains
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
# Line 76  contains Line 78  contains
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    
Legend:
Removed from v.76  
changed lines
  Added in v.130
  ViewVC Help
Powered by ViewVC 1.1.21