| 1 | !!---------------------------------------------------------------------- |
|---|
| 2 | !! TOP 1.0, LOCEAN-IPSL (2005) |
|---|
| 3 | !! $Header$ |
|---|
| 4 | !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
|---|
| 5 | !!---------------------------------------------------------------------- |
|---|
| 6 | !!! |
|---|
| 7 | !!! Modifications: |
|---|
| 8 | !!! -------------- |
|---|
| 9 | !!! 00-12 (E. Kestenare): |
|---|
| 10 | !!! assign a parameter to name individual tracers |
|---|
| 11 | !!! 01-02 (E. Kestenare): |
|---|
| 12 | !!! introduce jpno3 instead of jpnut |
|---|
| 13 | !!! |
|---|
| 14 | |
|---|
| 15 | #if defined key_trc_npzd |
|---|
| 16 | |
|---|
| 17 | IF(lwp) WRITE(numout,*) ' use NPZD biological model ' |
|---|
| 18 | IF(lwp) WRITE(numout,*) ' ' |
|---|
| 19 | |
|---|
| 20 | ! Check number of tracers |
|---|
| 21 | ! ----------------------- |
|---|
| 22 | IF (jptra /= 4) THEN |
|---|
| 23 | IF (lwp) THEN |
|---|
| 24 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 25 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 26 | WRITE (numout,*) & |
|---|
| 27 | & ' STOP, change jptra to 4 in ' & |
|---|
| 28 | & ,'parameter.passivetrc.npzd.h ' |
|---|
| 29 | END IF |
|---|
| 30 | STOP 'TRCCTL' |
|---|
| 31 | END IF |
|---|
| 32 | |
|---|
| 33 | ! Check tracer names |
|---|
| 34 | ! ------------------ |
|---|
| 35 | IF ( ctrcnm(jpdet) /= 'DET' .OR. ctrcnm(jpzoo) /= 'ZOO' & |
|---|
| 36 | & .OR. ctrcnm(jpphy) /= 'PHY' .OR. ctrcnm(jpno3 /= 'NUT' & |
|---|
| 37 | & .OR. ctrcnl(jpdet) /= 'DETRITUS' & |
|---|
| 38 | & .OR. ctrcnl(jpzoo) /= 'ZOOPLANKTON' & |
|---|
| 39 | & .OR. ctrcnl(jpphy) /= 'PHYTOPLANKTON' & |
|---|
| 40 | & .OR. ctrcnl(jpno3) /= 'NUTRIENTS') THEN |
|---|
| 41 | |
|---|
| 42 | ctrcnm(jpdet)='DET' |
|---|
| 43 | ctrcnl(jpdet)='DETRITUS' |
|---|
| 44 | ctrcnm(jpzoo)='ZOO' |
|---|
| 45 | ctrcnl(jpzoo)='ZOOPLANKTON' |
|---|
| 46 | ctrcnm(jpphy)='PHY' |
|---|
| 47 | ctrcnl(jpphy)='PHYTOPLANKTON' |
|---|
| 48 | ctrcnm(jpno3)='NUT' |
|---|
| 49 | ctrcnl(jpno3)='NUTRIENTS' |
|---|
| 50 | |
|---|
| 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 | ! Check tracer units |
|---|
| 63 | |
|---|
| 64 | DO jn=1,jptra |
|---|
| 65 | IF (ctrcun(jn) /= 'mmole/m3') THEN |
|---|
| 66 | ctrcun(jn)='mmole/m3' |
|---|
| 67 | IF (lwp) THEN |
|---|
| 68 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 69 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 70 | WRITE (numout,*) ' we force tracer unit' |
|---|
| 71 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ' ,ctrcun(jn) |
|---|
| 72 | WRITE(numout,*) ' ' |
|---|
| 73 | ENDIF |
|---|
| 74 | ENDIF |
|---|
| 75 | END DO |
|---|
| 76 | #endif |
|---|
