/[MITgcm]/MITgcm/pkg/cfc/cfc11_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/cfc/cfc11_forcing.F

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

revision 1.2 by stephd, Tue May 3 17:02:26 2005 UTC revision 1.7 by stephd, Tue Dec 6 17:01:06 2005 UTC
# Line 7  c modified to have carbonate and biologi Line 7  c modified to have carbonate and biologi
7  c modified for cfc: swd Sep 2003  c modified for cfc: swd Sep 2003
8  C  C
9    
 #include "CPP_OPTIONS.h"  
10  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
11    
12  CStartOfInterFace  CStartOfInterFace
# Line 28  C     == GLobal variables == Line 27  C     == GLobal variables ==
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "GRID.h"  #include "GRID.h"
29  #include "CFC.h"  #include "CFC.h"
30    #include "GCHEM.h"
31    #include "PTRACERS_SIZE.h"
32    #include "PTRACERS.h"
33    
34  C     == Routine arguments ==  C     == Routine arguments ==
35        INTEGER myIter        INTEGER myIter
# Line 46  C     == Local variables == Line 48  C     == Local variables ==
48         INTEGER myYear, lastYear, thisYear         INTEGER myYear, lastYear, thisYear
49         _RL dtinc, aWght, bWght         _RL dtinc, aWght, bWght
50         _RL ACFC11north, ACFC11south         _RL ACFC11north, ACFC11south
51         _RL maxYear         INTEGER maxYear
52         _RL a1, a2         _RL a1, a2
53         _RL yNorth, ySouth         _RL yNorth, ySouth
54                
# Line 58  C     == Local variables == Line 60  C     == Local variables ==
60           ENDDO           ENDDO
61    
62  c find atmospheric CFC  c find atmospheric CFC
63         myYear=float(myIter-nIter0)*deltaTclock         myYear=float(myIter-PTRACERS_Iter0)*deltaTclock
64       &                       /(360.d0*24.d0*3600.d0)       &                       /(360.d0*24.d0*3600.d0)
65         lastYear=1+int(myYear+0.5)         lastYear=1+int(myYear+0.5)
66         thisYear=lastYear+1         thisYear=lastYear+1
# Line 74  c       &write(0,*) 'myYear = ',myYear,l Line 76  c       &write(0,*) 'myYear = ',myYear,l
76           ACFC11south = ACFC11(lastYear,2)*bWght           ACFC11south = ACFC11(lastYear,2)*bWght
77       &                + ACFC11(thisYear,2)*aWght       &                + ACFC11(thisYear,2)*aWght
78         else         else
79           ACFC11north = ACFC11(thisYear,1)           ACFC11north = ACFC11(maxYear,1)
80           ACFC11south = ACFC11(thisYear,1)           ACFC11south = ACFC11(maxYear,2)
81         endif         endif
82  c      print*,'ACFC11north,ACFC11south', ACFC11north,ACFC11south,  c      print*,'ACFC11north,ACFC11south', ACFC11north,ACFC11south,
83  c    &                                   lastYear,thisYear  c    &                                   lastYear,thisYear
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22