[274] | 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 | !!---------------------------------------------------------------------- |
---|
[260] | 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 |
---|