1 | C $Id$ |
---|
2 | CCC |
---|
3 | CCC Modifications: |
---|
4 | CCC -------------- |
---|
5 | CCC 00-12 (E. Kestenare): |
---|
6 | CCC assign a parameter to name individual tracers |
---|
7 | CCC |
---|
8 | c |
---|
9 | #if defined key_trc_lobster1 |
---|
10 | IF(lwp) THEN |
---|
11 | WRITE(numout,*) ' use LOBSTER1 biological model ' |
---|
12 | WRITE(numout,*) ' ' |
---|
13 | ENDIF |
---|
14 | c |
---|
15 | c Check number of tracers |
---|
16 | c ----------------------- |
---|
17 | IF (jptra .ne. 6) 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 6 in ' |
---|
23 | $ ,'parameter.passivetrc.lobster1.h ' |
---|
24 | END IF |
---|
25 | STOP 'TRCCTL' |
---|
26 | END IF |
---|
27 | c |
---|
28 | c Check tracer names |
---|
29 | c ------------------ |
---|
30 | IF ( (ctrcnm(jpdet).ne.'DET').or.(ctrcnm(jpzoo).ne.'ZOO') |
---|
31 | $ .or.(ctrcnm(jpphy).ne.'PHY').or.(ctrcnm(jpno3).ne.'NO3') |
---|
32 | $ .or.(ctrcnm(jpnh4).ne.'NH4').or.(ctrcnm(jpdom).ne.'DOM') |
---|
33 | $ .or.(ctrcnl(jpdet).ne.'Detritus') |
---|
34 | $ .or.(ctrcnl(jpzoo).ne.'Zooplankton concentration') |
---|
35 | $ .or.(ctrcnl(jpphy).ne.'Phytoplankton concentration') |
---|
36 | $ .or.(ctrcnl(jpno3).ne.'Nitrate concentration') |
---|
37 | $ .or.(ctrcnl(jpnh4).ne.'Ammonium concentration') |
---|
38 | $ .or.(ctrcnl(jpdom).ne.'Dissolved organic matter') ) THEN |
---|
39 | ctrcnm(jpdet)='DET' |
---|
40 | ctrcnl(jpdet)='Detritus' |
---|
41 | ctrcnm(jpzoo)='ZOO' |
---|
42 | ctrcnl(jpzoo)='Zooplankton concentration' |
---|
43 | ctrcnm(jpphy)='PHY' |
---|
44 | ctrcnl(jpphy)='Phytoplankton concentration' |
---|
45 | ctrcnm(jpno3)='NO3' |
---|
46 | ctrcnl(jpno3)='Nitrate concentration' |
---|
47 | ctrcnm(jpnh4)='NH4' |
---|
48 | ctrcnl(jpnh4)='Ammonium concentration' |
---|
49 | ctrcnm(jpdom)='DOM' |
---|
50 | ctrcnl(jpdom)='Dissolved organic matter' |
---|
51 | IF (lwp) THEN |
---|
52 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
53 | WRITE (numout,*) ' ======= ============= ' |
---|
54 | WRITE (numout,*) ' we force tracer names' |
---|
55 | DO jn=1,jptra |
---|
56 | WRITE(numout,*) ' tracer nb: ',jn,' name = ',ctrcnm(jn) |
---|
57 | $ ,ctrcnl(jn) |
---|
58 | END DO |
---|
59 | WRITE(numout,*) ' ' |
---|
60 | ENDIF |
---|
61 | ENDIF |
---|
62 | c Check tracer units |
---|
63 | DO jn=1,jptra |
---|
64 | IF (ctrcun(jn).ne.'mmole-N/m3') THEN |
---|
65 | ctrcun(jn)='mmole-N/m3' |
---|
66 | IF (lwp) THEN |
---|
67 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
68 | WRITE (numout,*) ' ======= ============= ' |
---|
69 | WRITE (numout,*) ' we force tracer unit' |
---|
70 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ' |
---|
71 | $ ,ctrcun(jn) |
---|
72 | WRITE(numout,*) ' ' |
---|
73 | ENDIF |
---|
74 | ENDIF |
---|
75 | END DO |
---|
76 | #endif |
---|