[2038] | 1 | MODULE trcnam_cfc |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE trcnam_cfc *** |
---|
| 4 | !! TOP : initialisation of some run parameters for CFC chemical model |
---|
| 5 | !!====================================================================== |
---|
| 6 | !! History : 2.0 ! 2007-12 (C. Ethe, G. Madec) from trcnam.cfc.h90 |
---|
| 7 | !!---------------------------------------------------------------------- |
---|
| 8 | #if defined key_cfc |
---|
| 9 | !!---------------------------------------------------------------------- |
---|
| 10 | !! 'key_cfc' CFC tracers |
---|
| 11 | !!---------------------------------------------------------------------- |
---|
| 12 | !! trc_nam_cfc : CFC model initialisation |
---|
| 13 | !!---------------------------------------------------------------------- |
---|
| 14 | USE oce_trc ! Ocean variables |
---|
| 15 | USE par_trc ! TOP parameters |
---|
| 16 | USE trc ! TOP variables |
---|
| 17 | USE trcsms_cfc ! CFC specific variable |
---|
[3294] | 18 | USE iom ! I/O manager |
---|
[2038] | 19 | |
---|
| 20 | IMPLICIT NONE |
---|
| 21 | PRIVATE |
---|
| 22 | |
---|
| 23 | PUBLIC trc_nam_cfc ! called by trcnam.F90 module |
---|
| 24 | |
---|
| 25 | !!---------------------------------------------------------------------- |
---|
[2287] | 26 | !! NEMO/TOP 3.3 , NEMO Consortium (2010) |
---|
[2281] | 27 | !! $Id$ |
---|
[2287] | 28 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
[2038] | 29 | !!---------------------------------------------------------------------- |
---|
| 30 | |
---|
| 31 | CONTAINS |
---|
| 32 | |
---|
| 33 | SUBROUTINE trc_nam_cfc |
---|
| 34 | !!------------------------------------------------------------------- |
---|
| 35 | !! *** ROUTINE trc_nam_cfc *** |
---|
| 36 | !! |
---|
| 37 | !! ** Purpose : Definition some run parameter for CFC model |
---|
| 38 | !! |
---|
| 39 | !! ** Method : Read the namcfc namelist and check the parameter |
---|
[3294] | 40 | !! values called at the first timestep (nittrc000) |
---|
[2038] | 41 | !! |
---|
| 42 | !! ** input : Namelist namcfc |
---|
| 43 | !!---------------------------------------------------------------------- |
---|
[4147] | 44 | INTEGER :: numnatc_ref = -1 ! Logical unit for reference CFC namelist |
---|
| 45 | INTEGER :: numnatc_cfg = -1 ! Logical unit for configuration CFC namelist |
---|
| 46 | INTEGER :: numonc = -1 ! Logical unit for output namelist |
---|
| 47 | INTEGER :: ios ! Local integer output status for namelist read |
---|
[8353] | 48 | INTEGER :: jl, jn, cnt |
---|
[3294] | 49 | TYPE(DIAG), DIMENSION(jp_cfc_2d) :: cfcdia2d |
---|
[2038] | 50 | !! |
---|
[8353] | 51 | NAMELIST/namcfcdate/ ndate_beg, nyear_res, clnamecfc |
---|
[3294] | 52 | NAMELIST/namcfcdia/ cfcdia2d ! additional diagnostics |
---|
[4147] | 53 | !!---------------------------------------------------------------------- |
---|
| 54 | ! ! Open namelist files |
---|
| 55 | CALL ctl_opn( numnatc_ref, 'namelist_cfc_ref' , 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) |
---|
| 56 | CALL ctl_opn( numnatc_cfg, 'namelist_cfc_cfg' , 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) |
---|
[4624] | 57 | IF(lwm) CALL ctl_opn( numonc, 'output.namelist.cfc', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) |
---|
[2038] | 58 | |
---|
[4147] | 59 | REWIND( numnatc_ref ) ! Namelist namcfcdate in reference namelist : CFC parameters |
---|
| 60 | READ ( numnatc_ref, namcfcdate, IOSTAT = ios, ERR = 901) |
---|
| 61 | 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdate in reference namelist', lwp ) |
---|
[2038] | 62 | |
---|
[4147] | 63 | REWIND( numnatc_cfg ) ! Namelist namcfcdate in configuration namelist : CFC parameters |
---|
| 64 | READ ( numnatc_cfg, namcfcdate, IOSTAT = ios, ERR = 902 ) |
---|
| 65 | 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdate in configuration namelist', lwp ) |
---|
[4624] | 66 | IF(lwm) WRITE ( numonc, namcfcdate ) |
---|
[2038] | 67 | |
---|
[8353] | 68 | jn = jp_cfc0 - 1 |
---|
| 69 | ! Variables setting |
---|
| 70 | IF( lp_cfc11 ) THEN |
---|
| 71 | jn = jn + 1 |
---|
| 72 | ctrcnm (jn) = 'CFC11' |
---|
| 73 | ctrcln (jn) = 'Chlorofluoro carbon 11 Concentration' |
---|
| 74 | ctrcun (jn) = 'umolC/L' |
---|
| 75 | ln_trc_ini(jn) = .false. |
---|
| 76 | ln_trc_wri(jn) = .true. |
---|
| 77 | ENDIF |
---|
| 78 | ! |
---|
| 79 | IF( lp_cfc12 ) THEN |
---|
| 80 | jn = jn + 1 |
---|
| 81 | ctrcnm (jn) = 'CFC12' |
---|
| 82 | ctrcln (jn) = 'Chlorofluoro carbon 12 Concentration' |
---|
| 83 | ctrcun (jn) = 'umolC/L' |
---|
| 84 | ln_trc_ini(jn) = .false. |
---|
| 85 | ln_trc_wri(jn) = .true. |
---|
| 86 | ENDIF |
---|
| 87 | ! |
---|
| 88 | IF( lp_sf6 ) THEN |
---|
| 89 | jn = jn + 1 |
---|
| 90 | ctrcnm (jn) = 'SF6' |
---|
| 91 | ctrcln (jn) = 'Sulfur hexafluoride Concentration' |
---|
| 92 | ctrcun (jn) = 'umol/L' |
---|
| 93 | ln_trc_ini(jn) = .false. |
---|
| 94 | ln_trc_wri(jn) = .true. |
---|
| 95 | ENDIF |
---|
| 96 | |
---|
[2038] | 97 | IF(lwp) THEN ! control print |
---|
[8353] | 98 | WRITE(numout,*) ' ' |
---|
| 99 | WRITE(numout,*) ' CFCs' |
---|
| 100 | WRITE(numout,*) ' ' |
---|
[2038] | 101 | WRITE(numout,*) ' trc_nam: Read namdates, namelist for CFC chemical model' |
---|
| 102 | WRITE(numout,*) ' ~~~~~~~' |
---|
| 103 | WRITE(numout,*) ' initial calendar date (aammjj) for CFC ndate_beg = ', ndate_beg |
---|
| 104 | WRITE(numout,*) ' restoring time constant (year) nyear_res = ', nyear_res |
---|
[8353] | 105 | WRITE(numout,*) ' Atmospheric CFC concentrations file clnamecfc = ', TRIM(clnamecfc) |
---|
| 106 | WRITE(numout,*) ' Compute dynamics for CFC-11 lp_cfc11 = ', lp_cfc11 |
---|
| 107 | WRITE(numout,*) ' Compute dynamics for CFC-12 lp_cfc12 = ', lp_cfc12 |
---|
| 108 | WRITE(numout,*) ' Compute dynamics for SF6 lp_sf6 = ', lp_sf6 |
---|
[2038] | 109 | ENDIF |
---|
| 110 | nyear_beg = ndate_beg / 10000 |
---|
| 111 | IF(lwp) WRITE(numout,*) ' initial year (aa) nyear_beg = ', nyear_beg |
---|
| 112 | ! |
---|
[8353] | 113 | ! check consistency between CFC namelist and par_cfc setting |
---|
| 114 | if ( jn - jp_cfc0 + 1 .ne. jp_cfc ) & |
---|
| 115 | CALL ctl_stop( 'trc_nam_cfc: Number of selected CFCs is different from total CFC number (jp_cfc) specified in par_cfc.F90' ) |
---|
| 116 | ! |
---|
[2038] | 117 | |
---|
[3294] | 118 | IF( .NOT.lk_iomput .AND. ln_diatrc ) THEN |
---|
| 119 | ! |
---|
| 120 | ! Namelist namcfcdia |
---|
| 121 | ! ------------------- |
---|
[4147] | 122 | REWIND( numnatc_ref ) ! Namelist namcfcdia in reference namelist : CFC diagnostics |
---|
| 123 | READ ( numnatc_ref, namcfcdia, IOSTAT = ios, ERR = 903) |
---|
| 124 | 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdia in reference namelist', lwp ) |
---|
[2038] | 125 | |
---|
[4147] | 126 | REWIND( numnatc_cfg ) ! Namelist namcfcdia in configuration namelist : CFC diagnostics |
---|
| 127 | READ ( numnatc_cfg, namcfcdia, IOSTAT = ios, ERR = 904 ) |
---|
| 128 | 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdia in configuration namelist', lwp ) |
---|
[4624] | 129 | IF(lwm) WRITE ( numonc, namcfcdia ) |
---|
[2038] | 130 | |
---|
| 131 | DO jl = 1, jp_cfc_2d |
---|
| 132 | jn = jp_cfc0_2d + jl - 1 |
---|
[3294] | 133 | ctrc2d(jn) = TRIM( cfcdia2d(jl)%sname ) |
---|
| 134 | ctrc2l(jn) = TRIM( cfcdia2d(jl)%lname ) |
---|
| 135 | ctrc2u(jn) = TRIM( cfcdia2d(jl)%units ) |
---|
| 136 | END DO |
---|
| 137 | |
---|
| 138 | IF(lwp) THEN ! control print |
---|
| 139 | WRITE(numout,*) |
---|
| 140 | WRITE(numout,*) ' Namelist : natadd' |
---|
| 141 | DO jl = 1, jp_cfc_2d |
---|
| 142 | jn = jp_cfc0_2d + jl - 1 |
---|
| 143 | WRITE(numout,*) ' 2d diag nb : ', jn, ' short name : ', ctrc2d(jn), & |
---|
| 144 | & ' long name : ', ctrc2l(jn), ' unit : ', ctrc2u(jn) |
---|
| 145 | END DO |
---|
[2038] | 146 | WRITE(numout,*) ' ' |
---|
[3294] | 147 | ENDIF |
---|
| 148 | ! |
---|
[2038] | 149 | ENDIF |
---|
| 150 | |
---|
[4624] | 151 | IF(lwm) CALL FLUSH ( numonc ) ! flush output namelist CFC |
---|
[4147] | 152 | |
---|
[2038] | 153 | END SUBROUTINE trc_nam_cfc |
---|
| 154 | |
---|
| 155 | #else |
---|
| 156 | !!---------------------------------------------------------------------- |
---|
| 157 | !! Dummy module : No CFC |
---|
| 158 | !!---------------------------------------------------------------------- |
---|
| 159 | CONTAINS |
---|
| 160 | SUBROUTINE trc_nam_cfc ! Empty routine |
---|
| 161 | END SUBROUTINE trc_nam_cfc |
---|
| 162 | #endif |
---|
| 163 | |
---|
| 164 | !!====================================================================== |
---|
| 165 | END MODULE trcnam_cfc |
---|