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

Annotation of /MITgcm/pkg/cfc/cfc_atmos.F

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

Revision 1.1 - (hide annotations) (download)
Mon Oct 6 20:03:06 2003 UTC (20 years, 7 months ago) by stephd
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint52n_post, checkpoint52j_post, checkpoint53d_post, checkpoint54a_pre, checkpoint51o_pre, checkpoint55c_post, checkpoint54e_post, checkpoint52e_post, checkpoint51n_pre, checkpoint54a_post, checkpoint53c_post, checkpoint55d_pre, checkpoint51l_post, checkpoint51q_post, checkpoint51j_post, checkpoint55j_post, hrcube_1, branch-netcdf, checkpoint52d_pre, checkpoint52l_post, checkpoint55h_post, checkpoint51r_post, checkpoint52k_post, checkpoint52b_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint54d_post, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint54, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint51o_post, checkpoint51p_post, checkpoint52a_pre, checkpoint55f_post, checkpoint55i_post, checkpoint56, checkpoint51i_post, checkpoint53, checkpoint52, checkpoint51f_post, checkpoint52d_post, checkpoint52a_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint52c_post, checkpoint51h_pre, checkpoint51l_pre, ecco_c52_e35, hrcube5, checkpoint55e_post, checkpoint52i_post, checkpoint52j_pre, checkpoint53f_post, checkpoint55a_post, checkpoint51t_post, checkpoint53d_pre, checkpoint54c_post, checkpoint51n_post, checkpoint51i_pre, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_2, hrcube_3, checkpoint51m_post, checkpoint51s_post, checkpoint55d_post
Branch point for: branch-nonh, tg2-branch, checkpoint51n_branch, netcdf-sm0 
add CFC's for use with gchem package
1 stephd 1.1
2     #include "CPP_OPTIONS.h"
3     #include "GCHEM_OPTIONS.h"
4    
5     CStartOfInterface
6     SUBROUTINE CFC_ATMOS( myThid )
7     C /==========================================================\
8     C | SUBROUTINE CFC_ATMOS |
9     C | o read in timeseries of atmoshperic CFC
10     C |==========================================================|
11     C \==========================================================/
12     IMPLICIT NONE
13    
14     C === Global variables ===
15     #include "SIZE.h"
16     #include "EEPARAMS.h"
17     #include "PARAMS.h"
18     #include "GRID.h"
19     #include "DYNVARS.h"
20     #include "CFC.h"
21    
22     C == Routine arguments ==
23     C myThid - Number of this instance of INI_TR1
24     INTEGER myThid
25     integer i,j,bi,bj
26     CEndOfInterface
27    
28     C == Local variables ==
29     integer nTimePnts, Iunit, it
30     _RL year(100)
31    
32     #ifdef ALLOW_CFC
33    
34     _BEGIN_MASTER( mythid )
35     DO bj = myByLo(myThid), myByHi(myThid)
36     DO bi = myBxLo(myThid), myBxHi(myThid)
37     DO j=1-OLy,sNy+OLy
38     DO i=1-OLx,sNx+OLx
39     AtmosCFC11(i,j,bi,bj)=0.d0
40     AtmosCFC12(i,j,bi,bj)=0.d0
41     ENDDO
42     ENDDO
43     ENDDO
44     ENDDO
45     c
46     C read in CFC atmospheric timeseries data
47     cfc_yearbeg = 31
48     cfc_yearend = 98
49     nTimePnts=cfc_yearend-cfc_yearbeg
50     Iunit=28
51     OPEN(Iunit,FILE='cfc1112.atm',STATUS='old')
52     C skip 6 descriptor lines
53     DO i =1,6
54     READ(Iunit,*)
55     ENDDO
56     C Read in CFC11 and CFC12, N and S Hemisphere time histories
57     DO it = 1,nTimePnts
58     READ(Iunit,*) year(it),ACFC11(it,1),ACFC12(it,1),
59     & ACFC11(it,2),ACFC12(it,2)
60     write(0,*) 'year,acfc11_n,acfc12_n,acfc11_s,acfc12_s =',
61     & ACFC11(it,1),ACFC12(it,1),
62     & ACFC11(it,2),ACFC12(it,2)
63    
64     ENDDO
65     _END_MASTER(myThid)
66     c
67     #endif
68    
69    
70    
71     RETURN
72     END
  ViewVC Help
Powered by ViewVC 1.1.22