| 1 | ! $Id$ |
|---|
| 2 | !!! |
|---|
| 3 | !!! Modifications: |
|---|
| 4 | !!! -------------- |
|---|
| 5 | !!! 05-10 (C. Ethe): |
|---|
| 6 | !!! assign a parameter to name individual tracers |
|---|
| 7 | !!! |
|---|
| 8 | |
|---|
| 9 | IF(lwp) THEN |
|---|
| 10 | WRITE(numout,*) ' CFC Model ' |
|---|
| 11 | WRITE(numout,*) ' ' |
|---|
| 12 | ENDIF |
|---|
| 13 | |
|---|
| 14 | ! Check number of tracers |
|---|
| 15 | ! ----------------------- |
|---|
| 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 | |
|---|
| 27 | ! Check tracer names |
|---|
| 28 | ! ------------------ |
|---|
| 29 | IF ( jptra == 1 ) THEN |
|---|
| 30 | IF ( jp11 == 1 ) THEN |
|---|
| 31 | IF ( ctrcnm(jp11) /= 'CFC11') THEN |
|---|
| 32 | ctrcnm(jp11)='CFC11' |
|---|
| 33 | ctrcnl(jp11)='Chlorofuorocarbone 11 concentration' |
|---|
| 34 | ENDIF |
|---|
| 35 | ENDIF |
|---|
| 36 | IF ( jp12 == 1 ) THEN |
|---|
| 37 | IF ( ctrcnm(jp12) /= 'CFC12') THEN |
|---|
| 38 | ctrcnm(jp12)='CFC12' |
|---|
| 39 | ctrcnl(jp12)='Chlorofuorocarbone 12 concentration' |
|---|
| 40 | ENDIF |
|---|
| 41 | ENDIF |
|---|
| 42 | ENDIF |
|---|
| 43 | |
|---|
| 44 | IF ( jptra == 2 ) THEN |
|---|
| 45 | IF ( ctrcnm(jp11) /= 'CFC11' .OR. ctrcnm(jp12) /= 'CFC12' ) THEN |
|---|
| 46 | ctrcnm(jp11)='CFC11' |
|---|
| 47 | ctrcnl(jp11)='Chlorofuorocarbone 11 concentration' |
|---|
| 48 | ctrcnm(jp12)='CFC12' |
|---|
| 49 | ctrcnl(jp12)='Chlorofuorocarbone 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 | |
|---|
| 64 | ! Check tracer units |
|---|
| 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 |
|---|
