| 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 01-02 (E. Kestenare): |
|---|
| 8 | CCC introduce jpno3 instead of jpnut |
|---|
| 9 | |
|---|
| 10 | CCC |
|---|
| 11 | c |
|---|
| 12 | #if defined key_trc_npzd |
|---|
| 13 | IF(lwp) THEN |
|---|
| 14 | WRITE(numout,*) ' use NPZD biological model ' |
|---|
| 15 | WRITE(numout,*) ' ' |
|---|
| 16 | ENDIF |
|---|
| 17 | c |
|---|
| 18 | c Check number of tracers |
|---|
| 19 | c ----------------------- |
|---|
| 20 | IF (jptra.ne.4) THEN |
|---|
| 21 | IF (lwp) THEN |
|---|
| 22 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 23 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 24 | WRITE (numout,*) |
|---|
| 25 | $ ' STOP, change jptra to 4 in ' |
|---|
| 26 | $ ,'parameter.passivetrc.npzd.h ' |
|---|
| 27 | END IF |
|---|
| 28 | STOP 'TRCCTL' |
|---|
| 29 | END IF |
|---|
| 30 | c |
|---|
| 31 | c Check tracer names |
|---|
| 32 | c ------------------ |
|---|
| 33 | IF ((ctrcnm(jpdet).ne.'DET').OR.(ctrcnm(jpzoo).ne.'ZOO') |
|---|
| 34 | $ .or.(ctrcnm(jpphy).ne.'PHY').or.(ctrcnm(jpno3).ne.'NUT') |
|---|
| 35 | $ .or.(ctrcnl(jpdet).ne.'DETRITUS') |
|---|
| 36 | $ .or.(ctrcnl(jpzoo).ne.'ZOOPLANKTON') |
|---|
| 37 | $ .or.(ctrcnl(jpphy).ne.'PHYTOPLANKTON') |
|---|
| 38 | $ .or.(ctrcnl(jpno3).ne.'NUTRIENTS')) THEN |
|---|
| 39 | ctrcnm(jpdet)='DET' |
|---|
| 40 | ctrcnl(jpdet)='DETRITUS' |
|---|
| 41 | ctrcnm(jpzoo)='ZOO' |
|---|
| 42 | ctrcnl(jpzoo)='ZOOPLANKTON' |
|---|
| 43 | ctrcnm(jpphy)='PHY' |
|---|
| 44 | ctrcnl(jpphy)='PHYTOPLANKTON' |
|---|
| 45 | ctrcnm(jpno3)='NUT' |
|---|
| 46 | ctrcnl(jpno3)='NUTRIENTS' |
|---|
| 47 | |
|---|
| 48 | IF (lwp) THEN |
|---|
| 49 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 50 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 51 | WRITE (numout,*) ' we force tracer names' |
|---|
| 52 | DO jn=1,jptra |
|---|
| 53 | WRITE(numout,*) ' tracer nb: ',jn,' name = ',ctrcnm(jn) |
|---|
| 54 | $ ,ctrcnl(jn) |
|---|
| 55 | END DO |
|---|
| 56 | WRITE(numout,*) ' ' |
|---|
| 57 | ENDIF |
|---|
| 58 | ENDIF |
|---|
| 59 | c Check tracer units |
|---|
| 60 | DO jn=1,jptra |
|---|
| 61 | IF (ctrcun(jn).ne.'mmole/m3') THEN |
|---|
| 62 | ctrcun(jn)='mmole/m3' |
|---|
| 63 | IF (lwp) THEN |
|---|
| 64 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 65 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 66 | WRITE (numout,*) ' we force tracer unit' |
|---|
| 67 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ' |
|---|
| 68 | $ ,ctrcun(jn) |
|---|
| 69 | WRITE(numout,*) ' ' |
|---|
| 70 | ENDIF |
|---|
| 71 | ENDIF |
|---|
| 72 | END DO |
|---|
| 73 | #endif |
|---|
