/[MITgcm]/MITgcm/verification/global_with_CFC11/code1x1/ini_tr1.F
ViewVC logotype

Annotation of /MITgcm/verification/global_with_CFC11/code1x1/ini_tr1.F

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


Revision 1.1.2.2 - (hide annotations) (download)
Thu Aug 25 18:21:17 2005 UTC (15 years, 10 months ago) by dimitri
Branch: release1_50yr
Changes since 1.1.2.1: +8 -5 lines
updating verification/global_with_CFC11/code1x1/external_forcing_tr.F
and verification/global_with_CFC11/code1x1/ini_tr1.F
for release1_50yr

1 dimitri 1.1.2.2 C $Header: /u/gcmpack/MITgcm/verification/global_with_CFC11/code1x1/Attic/ini_tr1.F,v 1.1.2.1 2005/08/25 16:22:17 dimitri Exp $
2     C $Name: $
3 dimitri 1.1.2.1
4     #include "CPP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: INI_TR1
8     C !INTERFACE:
9     SUBROUTINE INI_TR1( myThid )
10     C !DESCRIPTION: \bv
11     C *==========================================================*
12     C | SUBROUTINE INI_TR1
13     C | o Set initial tracer 1 distribution.
14     C *==========================================================*
15     C | Passive tracers 1-N can be initialised so as to dye fluid.
16     C | This routine is hardcoded for OCMIP CFC-11 experiment.
17     C *==========================================================*
18     C \ev
19    
20     C !USES:
21     IMPLICIT NONE
22     C === Global variables ===
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26     #include "GRID.h"
27     #include "DYNVARS.h"
28     #ifdef ALLOW_PASSIVE_TRACER
29     #include "TR1.h"
30     #endif
31    
32     C !INPUT/OUTPUT PARAMETERS:
33     C == Routine arguments ==
34     C myThid - Number of this instance of INI_TR1
35     INTEGER myThid
36    
37     #ifdef ALLOW_PASSIVE_TRACER
38    
39     C !LOCAL VARIABLES:
40     C == Local variables ==
41     C iC, jC - Center of domain
42     C iD, jD - Disitance from domain center.
43     C rad - Radius of initial patch
44     C rD - Radial displacement of point I,J
45     C iG, jG - Global coordinate index
46     C bi,bj - Loop counters
47     C I,J,K
48     INTEGER iC, jC, iD, jD
49     INTEGER iG, jG
50     INTEGER bi, bj
51     INTEGER I, J, K, localWarnings
52     _RL rad, rD
53     CHARACTER*(MAX_LEN_MBUF) msgBuf
54     CEOP
55    
56     _BARRIER
57    
58     C-- Initialise tracer to inline.
59     DO bj = myByLo(myThid), myByHi(myThid)
60     DO bi = myBxLo(myThid), myBxHi(myThid)
61     DO K=1,Nr
62     DO J=1,sNy
63     DO I=1,sNx
64     tr1(I,J,K,bi,bj) = 0. _d 0
65     ENDDO
66     ENDDO
67     ENDDO
68     ENDDO
69     ENDDO
70     C Set initial tendency terms
71     localWarnings=0
72     DO bj = myByLo(myThid), myByHi(myThid)
73     DO bi = myBxLo(myThid), myBxHi(myThid)
74     DO K=1,Nr
75     DO J=1,sNy
76     DO I=1,sNx
77     gtr1 (I,J,K,bi,bj) = 0. _d 0
78     gtr1NM1(I,J,K,bi,bj) = 0. _d 0
79     ENDDO
80     ENDDO
81     ENDDO
82     ENDDO
83     ENDDO
84     C
85     _EXCH_XYZ_R8(tr1 , myThid )
86     _EXCH_XYZ_R8(gtr1 , myThid )
87     _EXCH_XYZ_R8(gtr1NM1 , myThid )
88    
89     CALL PLOT_FIELD_XYZRL( tr1, 'Initial TR1' ,
90     & Nr, 1, myThid )
91    
92     C-- Input file names
93     FiceFile = 'FICE'
94     XkwFile = 'XKW'
95     PatmFile = 'PATM'
96    
97     C-- Compute numbers need to determine pCFC
98     CALL read_cfcatm
99     DO bj = myByLo(myThid), myByHi(myThid)
100     DO bi = myBxLo(myThid), myBxHi(myThid)
101     DO J=1,sNy
102     DO I=1,sNx
103     IF ( yC(I,J,bi,bj) .GE. 10. ) THEN
104     pCFCw1(I,J,bi,bj) = 1.
105     pCFCw2(I,J,bi,bj) = 0.
106     ELSEIF ( yC(I,J,bi,bj) .LE. -10. ) THEN
107     pCFCw1(I,J,bi,bj) = 0.
108     pCFCw2(I,J,bi,bj) = 1.
109     ELSE
110     pCFCw1(I,J,bi,bj) = ( yC(I,J,bi,bj) + 10. ) / 20.
111     pCFCw2(I,J,bi,bj) = ( 10. - yC(I,J,bi,bj) ) / 20.
112     ENDIF
113     ENDDO
114     ENDDO
115     ENDDO
116     ENDDO
117     _EXCH_XY_R8(pCFCw1 , myThid )
118     _EXCH_XY_R8(pCFCw2 , myThid )
119    
120     #endif /* ALLOW_PASSIVE_TRACER */
121    
122     RETURN
123     END
124    
125     C====================================================================
126    
127     c_ ---------------------------------------------------------------------
128     c_ RCS lines preceded by "c_ "
129     c_ ---------------------------------------------------------------------
130     c_
131 dimitri 1.1.2.2 c_ $Source: /u/gcmpack/MITgcm/verification/global_with_CFC11/code1x1/Attic/ini_tr1.F,v $
132 dimitri 1.1.2.1 c_ $Revision: 1.1.2.1 $
133 dimitri 1.1.2.2 c_ $Date: 2005/08/25 16:22:17 $ ; $State: Exp $
134 dimitri 1.1.2.1 c_ $Author: dimitri $ ; $Locker: $
135     c_
136     c_ ---------------------------------------------------------------------
137     c_ $Log: ini_tr1.F,v $
138 dimitri 1.1.2.2 c_ Revision 1.1.2.1 2005/08/25 16:22:17 dimitri
139     c_ adding ecco1x1 verification/global_with_CFC11 experiment
140     c_
141 dimitri 1.1.2.1 c_ Revision 1.1.2.1 2003/05/03 04:27:34 dimitri
142     c_ Branched off release1_p15
143     c_ Modified pkg/exf/exf_clim_param.h, exf_clim_readparms.F,
144 dimitri 1.1.2.2 c_ exf_set_climsss.F, and exf_set_climsst.F to add exf_interp
145 dimitri 1.1.2.1 c_ capability to SST and SSS
146     c_ Added verification/global_with_CFC11/code50yr and input50yr
147     c_
148     c_ Revision 1.1.2.1 2003/05/01 18:18:29 dimitri
149     c_ release1_p15
150     c_ o Added CFC-11 diagnostics, see
151     c_ verification/global_with_CFC11/README
152     c_
153     c_ Revision 1.2 1998/07/20 15:33:32 jomce
154     c_ Used new data set from Walker et al. (sent by J. Bullister 17/7/1998)
155     c_ from mid 1931 to mid 1997. Format of cfc1112.atm has changed. This
156     c_ file modified accordingly, and old atm record moved to old/cfc1112.atm.old
157     c_
158     c_ Revision 1.1 1998/07/07 15:22:00 orr
159     c_ Initial revision
160     c_
161     c_ ---------------------------------------------------------------------
162     c_
163     SUBROUTINE read_cfcatm
164    
165     C --------------------------------------------------------------------
166     C Reads temporal history of atmospheric CFC-11 and CFC-12 (both in pptv)
167     C
168     C year = year of measurement (mid-year)
169     C p11(year,nt) = pcfc11 in northern (1) and southern (2) hemisphere
170     C p12(year,nt) = pcfc12 in northern (1) and southern (2) hemisphere
171     C --------------------------------------------------------------------
172    
173     c Jean-Claude Dutay, LSCE/CEA-CNRS, Saclay, France, 10 May 1998
174     c - modfied by J. Orr (9 June 1998)
175    
176     IMPLICIT NONE
177     #include "SIZE.h"
178     #include "EEPARAMS.h"
179     #include "TR1.h"
180    
181     INTEGER i, n, iu
182     PARAMETER (iu=10)
183    
184     CHARACTER*80 filen
185     C
186     C OPEN FILE
187     C ---------
188     filen='cfc1112.atm'
189     OPEN(UNIT=iu,FILE=filen,STATUS='old')
190     c
191     c Skip over 1st six descriptor lines
192     c ----------------------------------
193     DO i=1,6
194     READ(iu,99)
195     END do
196    
197     C
198     C READ FILE
199     C ---------
200     DO n = nyrbeg, nyrend
201     READ(iu,*)
202     + CFCyear(n), CFCp11(n,1), CFCp12(n,1)
203     + , CFCp11(n,2), CFCp12(n,2)
204     WRITE(6,100)
205     + CFCyear(n), CFCp11(n,1), CFCp12(n,1)
206     + , CFCp11(n,2), CFCp12(n,2)
207     END DO
208    
209     99 FORMAT(1x)
210     100 FORMAT(f7.2, 4f8.2)
211    
212     RETURN
213     END

  ViewVC Help
Powered by ViewVC 1.1.22