| 1 | !!====================================================================== |
|---|
| 2 | !! *** trcctl.cfc.h90 *** |
|---|
| 3 | !! TOP : Control of CFC chemical model |
|---|
| 4 | !!====================================================================== |
|---|
| 5 | !!---------------------------------------------------------------------- |
|---|
| 6 | !! History : 1.0 ! 2005-10 (C. Ethe) assign a parameter to name individual tracers |
|---|
| 7 | !!---------------------------------------------------------------------- |
|---|
| 8 | !!---------------------------------------------------------------------- |
|---|
| 9 | !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005) |
|---|
| 10 | !! $Id:$ |
|---|
| 11 | !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) |
|---|
| 12 | !!---------------------------------------------------------------------- |
|---|
| 13 | |
|---|
| 14 | |
|---|
| 15 | IF(lwp) THEN |
|---|
| 16 | WRITE(numout,*) ' CFC Model ' |
|---|
| 17 | WRITE(numout,*) ' ' |
|---|
| 18 | ENDIF |
|---|
| 19 | |
|---|
| 20 | ! Check number of tracers |
|---|
| 21 | ! ----------------------- |
|---|
| 22 | IF( jptra > 2) 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 1 or 2 in par_sms_cfc.h90 ' |
|---|
| 28 | END IF |
|---|
| 29 | STOP 'TRC_CTL' |
|---|
| 30 | END IF |
|---|
| 31 | |
|---|
| 32 | ! Check tracer names |
|---|
| 33 | ! ------------------ |
|---|
| 34 | IF( jptra == 1 ) THEN |
|---|
| 35 | IF ( jp11 == 1 ) THEN |
|---|
| 36 | IF ( ctrcnm(jp11) /= 'CFC11') THEN |
|---|
| 37 | ctrcnm(jp11) = 'CFC11' |
|---|
| 38 | ctrcnl(jp11) = 'Chlorofuorocarbone 11 concentration' |
|---|
| 39 | ENDIF |
|---|
| 40 | ENDIF |
|---|
| 41 | IF( jp12 == 1 ) THEN |
|---|
| 42 | IF ( ctrcnm(jp12) /= 'CFC12') THEN |
|---|
| 43 | ctrcnm(jp12) = 'CFC12' |
|---|
| 44 | ctrcnl(jp12) = 'Chlorofuorocarbone 12 concentration' |
|---|
| 45 | ENDIF |
|---|
| 46 | ENDIF |
|---|
| 47 | ENDIF |
|---|
| 48 | |
|---|
| 49 | IF( jptra == 2 ) THEN |
|---|
| 50 | IF ( ctrcnm(jp11) /= 'CFC11' .OR. ctrcnm(jp12) /= 'CFC12' ) THEN |
|---|
| 51 | ctrcnm(jp11) = 'CFC11' |
|---|
| 52 | ctrcnl(jp11) = 'Chlorofuorocarbone 11 concentration' |
|---|
| 53 | ctrcnm(jp12) = 'CFC12' |
|---|
| 54 | ctrcnl(jp12) = 'Chlorofuorocarbone 12 concentration' |
|---|
| 55 | ENDIF |
|---|
| 56 | ENDIF |
|---|
| 57 | |
|---|
| 58 | IF(lwp) THEN |
|---|
| 59 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 60 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 61 | WRITE (numout,*) ' we force tracer names' |
|---|
| 62 | DO jn = 1, jptra |
|---|
| 63 | WRITE(numout,*) ' tracer nb: ',jn,' name = ',ctrcnm(jn), ctrcnl(jn) |
|---|
| 64 | END DO |
|---|
| 65 | WRITE(numout,*) ' ' |
|---|
| 66 | ENDIF |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | ! Check tracer units |
|---|
| 70 | ! ------------------ |
|---|
| 71 | DO jn = 1, jptra |
|---|
| 72 | IF( ctrcun(jn) /= 'mole/m3' ) THEN |
|---|
| 73 | ctrcun(jn) = 'mole/m3' |
|---|
| 74 | IF(lwp) THEN |
|---|
| 75 | WRITE (numout,*) ' ===>>>> : w a r n i n g ' |
|---|
| 76 | WRITE (numout,*) ' ======= ============= ' |
|---|
| 77 | WRITE (numout,*) ' we force tracer unit' |
|---|
| 78 | WRITE(numout,*) ' tracer ',ctrcnm(jn), 'UNIT= ',ctrcun(jn) |
|---|
| 79 | WRITE(numout,*) ' ' |
|---|
| 80 | ENDIF |
|---|
| 81 | ENDIF |
|---|
| 82 | END DO |
|---|
