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