[350] | 1 | ! $Id$ |
---|
| 2 | !!! |
---|
| 3 | !!! Modifications: |
---|
| 4 | !!! -------------- |
---|
| 5 | !!! 05-10 (C. Ethe): |
---|
| 6 | !!! assign a parameter to name individual tracers |
---|
| 7 | !!! |
---|
[338] | 8 | |
---|
| 9 | IF(lwp) THEN |
---|
| 10 | WRITE(numout,*) ' CFC Model ' |
---|
| 11 | WRITE(numout,*) ' ' |
---|
| 12 | ENDIF |
---|
| 13 | |
---|
[350] | 14 | ! Check number of tracers |
---|
| 15 | ! ----------------------- |
---|
[338] | 16 | |
---|
| 17 | IF (jptra .GT. 2) 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 1 or 2 in par_sms_cfc.h90 ' |
---|
| 23 | END IF |
---|
| 24 | STOP 'TRC_CTL' |
---|
| 25 | END IF |
---|
| 26 | |
---|
[350] | 27 | ! Check tracer names |
---|
| 28 | ! ------------------ |
---|
[338] | 29 | IF ( jptra == 1 ) THEN |
---|
| 30 | IF ( jp11 == 1 ) THEN |
---|
[350] | 31 | IF ( ctrcnm(jp11) /= 'CFC11') THEN |
---|
| 32 | ctrcnm(jp11)='CFC11' |
---|
[338] | 33 | ctrcnl(jp11)='Carbonate 11 concentration' |
---|
| 34 | ENDIF |
---|
| 35 | ENDIF |
---|
| 36 | IF ( jp12 == 1 ) THEN |
---|
[350] | 37 | IF ( ctrcnm(jp12) /= 'CFC12') THEN |
---|
| 38 | ctrcnm(jp12)='CFC12' |
---|
[338] | 39 | ctrcnl(jp12)='Carbonate 12 concentration' |
---|
| 40 | ENDIF |
---|
| 41 | ENDIF |
---|
| 42 | ENDIF |
---|
| 43 | |
---|
| 44 | IF ( jptra == 2 ) THEN |
---|
[350] | 45 | IF ( ctrcnm(jp11) /= 'CFC11' .OR. ctrcnm(jp12) /= 'CFC12' ) THEN |
---|
| 46 | ctrcnm(jp11)='CFC11' |
---|
[338] | 47 | ctrcnl(jp11)='Carbonate 11 concentration' |
---|
[350] | 48 | ctrcnm(jp12)='CFC12' |
---|
[338] | 49 | ctrcnl(jp12)='Carbonate 12 concentration' |
---|
| 50 | ENDIF |
---|
| 51 | ENDIF |
---|
| 52 | |
---|
| 53 | IF (lwp) THEN |
---|
| 54 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
| 55 | WRITE (numout,*) ' ======= ============= ' |
---|
| 56 | WRITE (numout,*) ' we force tracer names' |
---|
| 57 | DO jn=1,jptra |
---|
| 58 | WRITE(numout,*) ' tracer nb: ',jn,' name = ',ctrcnm(jn), ctrcnl(jn) |
---|
| 59 | END DO |
---|
| 60 | WRITE(numout,*) ' ' |
---|
| 61 | ENDIF |
---|
| 62 | |
---|
| 63 | |
---|
[350] | 64 | ! Check tracer units |
---|
[338] | 65 | |
---|
| 66 | DO jn=1,jptra |
---|
| 67 | IF (ctrcun(jn) /= 'mole/m3') THEN |
---|
| 68 | ctrcun(jn)='mole/m3' |
---|
| 69 | IF (lwp) THEN |
---|
| 70 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
---|
| 71 | WRITE (numout,*) ' ======= ============= ' |
---|
| 72 | WRITE (numout,*) ' we force tracer unit' |
---|
| 73 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ',ctrcun(jn) |
---|
| 74 | WRITE(numout,*) ' ' |
---|
| 75 | ENDIF |
---|
| 76 | ENDIF |
---|
| 77 | END DO |
---|