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