1 | !!---------------------------------------------------------------------- |
---|
2 | !! TOP 1.0 , LOCEAN-IPSL (2005) |
---|
3 | !! $Id$ |
---|
4 | !! This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
5 | !!---------------------------------------------------------------------- |
---|
6 | !!! |
---|
7 | !!! Modifications: |
---|
8 | !!! -------------- |
---|
9 | !!! 05-10 (C. Ethe): |
---|
10 | !!! assign a parameter to name individual tracers |
---|
11 | !!! |
---|
12 | |
---|
13 | IF(lwp) THEN |
---|
14 | WRITE(numout,*) ' CFC Model ' |
---|
15 | WRITE(numout,*) ' ' |
---|
16 | ENDIF |
---|
17 | |
---|
18 | ! Check number of tracers |
---|
19 | ! ----------------------- |
---|
20 | |
---|
21 | IF (jptra .GT. 2) THEN |
---|
22 | IF (lwp) THEN |
---|
23 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
24 | WRITE (numout,*) ' ======= ============= ' |
---|
25 | WRITE (numout,*) & |
---|
26 | & ' STOP, change jptra to 1 or 2 in par_sms_cfc.h90 ' |
---|
27 | END IF |
---|
28 | STOP 'TRC_CTL' |
---|
29 | END IF |
---|
30 | |
---|
31 | ! Check tracer names |
---|
32 | ! ------------------ |
---|
33 | IF ( jptra == 1 ) THEN |
---|
34 | IF ( jp11 == 1 ) THEN |
---|
35 | IF ( ctrcnm(jp11) /= 'CFC11') THEN |
---|
36 | ctrcnm(jp11)='CFC11' |
---|
37 | ctrcnl(jp11)='Chlorofuorocarbone 11 concentration' |
---|
38 | ENDIF |
---|
39 | ENDIF |
---|
40 | IF ( jp12 == 1 ) THEN |
---|
41 | IF ( ctrcnm(jp12) /= 'CFC12') THEN |
---|
42 | ctrcnm(jp12)='CFC12' |
---|
43 | ctrcnl(jp12)='Chlorofuorocarbone 12 concentration' |
---|
44 | ENDIF |
---|
45 | ENDIF |
---|
46 | ENDIF |
---|
47 | |
---|
48 | IF ( jptra == 2 ) THEN |
---|
49 | IF ( ctrcnm(jp11) /= 'CFC11' .OR. ctrcnm(jp12) /= 'CFC12' ) THEN |
---|
50 | ctrcnm(jp11)='CFC11' |
---|
51 | ctrcnl(jp11)='Chlorofuorocarbone 11 concentration' |
---|
52 | ctrcnm(jp12)='CFC12' |
---|
53 | ctrcnl(jp12)='Chlorofuorocarbone 12 concentration' |
---|
54 | ENDIF |
---|
55 | ENDIF |
---|
56 | |
---|
57 | IF (lwp) THEN |
---|
58 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
59 | WRITE (numout,*) ' ======= ============= ' |
---|
60 | WRITE (numout,*) ' we force tracer names' |
---|
61 | DO jn=1,jptra |
---|
62 | WRITE(numout,*) ' tracer nb: ',jn,' name = ',ctrcnm(jn), ctrcnl(jn) |
---|
63 | END DO |
---|
64 | WRITE(numout,*) ' ' |
---|
65 | ENDIF |
---|
66 | |
---|
67 | |
---|
68 | ! Check tracer units |
---|
69 | |
---|
70 | DO jn=1,jptra |
---|
71 | IF (ctrcun(jn) /= 'mole/m3') THEN |
---|
72 | ctrcun(jn)='mole/m3' |
---|
73 | IF (lwp) THEN |
---|
74 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
75 | WRITE (numout,*) ' ======= ============= ' |
---|
76 | WRITE (numout,*) ' we force tracer unit' |
---|
77 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ',ctrcun(jn) |
---|
78 | WRITE(numout,*) ' ' |
---|
79 | ENDIF |
---|
80 | ENDIF |
---|
81 | END DO |
---|