Changeset 945 for trunk/NEMO/TOP_SRC/trcdtr.F90
- Timestamp:
- 2008-05-14T18:14:53+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/trcdtr.F90
r730 r945 1 1 MODULE trcdtr 2 !!======================================================================================= 3 !! 4 !! *** MODULE trcdtr *** 5 !! 6 !! Computes or READ initial DATA for passive tracer 7 !! 8 !!======================================================================================= 9 !! TOP 1.0, LOCEAN-IPSL (2005) 10 !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/trcdtr.F90,v 1.8 2007/10/17 14:48:56 opalod Exp $ 11 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 12 !!---------------------------------------------------------------------- 13 !!---------------------------------------------------------------------- 14 !! * Modules used 15 !! ============== 16 USE oce_trc 17 USE trc 18 USE sms 19 USE trcdta 20 USE lib_mpp 21 22 IMPLICIT NONE 23 PRIVATE 24 !! * Accessibility 25 PUBLIC trc_dtr 2 !!====================================================================== 3 !! *** MODULE trcdtr *** 4 !! TOP : computes or READ initial DATA for passive tracer 5 !!====================================================================== 6 !! History : - ! 1996-11 () original code 7 !! ! 2000-12 (O. Aumont, E. Kestenare) add for POC in sediments 8 !! 1.0 ! 2005-12 (O. Aumont, A. El Moussaoui) F90 9 !!---------------------------------------------------------------------- 10 #if defined key_top 11 !!---------------------------------------------------------------------- 12 !! 'key_top' TOP models 13 !!---------------------------------------------------------------------- 14 !! trc_dtr : computes or READ initial DATA for passive tracer 15 !!---------------------------------------------------------------------- 16 USE oce_trc 17 USE trc 18 USE sms 19 USE trcdta 20 USE lib_mpp 21 22 IMPLICIT NONE 23 PRIVATE 24 25 PUBLIC trc_dtr ! called in ??? 26 27 !!---------------------------------------------------------------------- 28 !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005) 29 !! $Header:$ 30 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 31 !!---------------------------------------------------------------------- 26 32 27 33 CONTAINS 28 34 29 #if defined key_passivetrc 30 31 SUBROUTINE trc_dtr 32 !!--------------------------------------------------------------------- 33 !! 34 !! ROUTINE trci_dtr 35 !! ****************** 36 !! PURPOSE : 37 !! --------- 38 !! computes or READ initial DATA for passive tracer 39 !! ----- 40 !! COMMON 41 !! /comdom/ : domain PARAMETER 42 !! /comcoo/ : orthogonal curvilinear coordinates 43 !! and scale factors 44 !! /comask/ : masks, bathymetry 45 !! OUTPUT : 46 !! ------ 47 !! COMMON 48 !! /cottrc/ : passive tracer field now and before 49 !! 50 !! 51 !! History: 52 !! -------- 53 !! original : 96-11 54 !! additions : 99-9 55 !! : 00-12 (O. Aumont, E. Kestenare) add for POC in sediments 56 !! add for POC in sediments 57 !! 03/05 O. Aumont and A. El Moussaoui F90 58 !!---------------------------------------------------------------------- 59 !!---------------------------------------------------------------------- 60 !! local declarations 61 !! ================== 62 INTEGER :: ji,jj,jk,jn 63 #if defined key_trc_pisces 64 REAL(wp) :: alka0,oxyg0,calc0,bioma0, & 65 silic1,po4,no3,caralk,bicarb,co3 66 #endif 67 !!--------------------------------------------------------------------- 68 !! OPA.9 69 !!--------------------------------------------------------------------- 70 !! 0. initialisations 71 !! ------------------ 72 73 IF(lwp) WRITE(numout,*) ' ' 74 IF(lwp) WRITE(numout,*) ' *** trcdtr initialisation for ' 75 IF(lwp) WRITE(numout,*) ' passive tracers' 76 IF(lwp) WRITE(numout,*) ' ' 77 35 SUBROUTINE trc_dtr 36 !!--------------------------------------------------------------------- 37 !! *** ROUTINE trc_dtr *** 38 !! 39 !! ** Purpose : computes or READ initial DATA for passive tracer 40 !! 41 !! ** Method : 42 !!--------------------------------------------------------------------- 43 INTEGER :: ji, jj, jk, jn 44 # if defined key_pisces 45 REAL(wp) :: alka0, sco2, oxyg0, calc0, bioma0 46 REAL(wp) :: silic1, po4, no3, caralk, bicarb, co3 47 # endif 48 !!--------------------------------------------------------------------- 49 50 IF(lwp) WRITE(numout,*) 51 IF(lwp) WRITE(numout,*) 'trc_dtr : initialisation of the passive tracers' 52 IF(lwp) WRITE(numout,*) '~~~~~~~' 78 53 79 54 #if defined key_cfc 80 trn(:,:,:,:)=0.0 81 #elif defined key_trc_pisces 82 83 sco2 = 2.3e-3 84 alka0 = 2.39e-3 85 oxyg0 = 1.8e-4 86 po4 = 2.165e-6/po4r 87 bioma0 = 1.e-8 88 silic1 = 91.51e-6 89 calc0 = 1.e-6 90 no3 = 30.88E-6*7.6 55 trn(:,:,:,:) = 0.0 ! CFC initialisation to zero 56 57 #elif defined key_pisces 58 ! PISCES initialisation 59 ! --------------------- 60 sco2 = 2.312e-3 61 alka0 = 2.423e-3 62 oxyg0 = 177.6e-6 63 po4 = 2.174e-6 / po4r 64 bioma0 = 1.000e-8 65 silic1 = 91.65e-6 66 no3 = 31.04e-6 * 7.6 91 67 92 68 trn(:,:,:,jpdic) = sco2 … … 94 70 trn(:,:,:,jptal) = alka0 95 71 trn(:,:,:,jpoxy) = oxyg0 96 trn(:,:,:,jpcal) = calc072 trn(:,:,:,jpcal) = bioma0 97 73 trn(:,:,:,jppo4) = po4 98 74 trn(:,:,:,jppoc) = bioma0 99 # if ! defined key_trc_kriest75 # if ! defined key_kriest 100 76 trn(:,:,:,jpgoc) = bioma0 101 trn(:,:,:,jpbfe) = bioma0 *5E-6102 # else103 trn(:,:,:,jpnum) = bioma0 /(6.*xkr_massp)104 # endif77 trn(:,:,:,jpbfe) = bioma0 * 5.e-6 78 # else 79 trn(:,:,:,jpnum) = bioma0 / ( 6. * xkr_massp ) 80 # endif 105 81 trn(:,:,:,jpsil) = silic1 106 trn(:,:,:,jpbsi) = bioma0 *0.15107 trn(:,:,:,jpdsi) = bioma0 *5.E-682 trn(:,:,:,jpbsi) = bioma0 * 0.15 83 trn(:,:,:,jpdsi) = bioma0 * 5.e-6 108 84 trn(:,:,:,jpphy) = bioma0 109 85 trn(:,:,:,jpdia) = bioma0 … … 111 87 trn(:,:,:,jpmes) = bioma0 112 88 trn(:,:,:,jpfer) = 0.6E-9 113 trn(:,:,:,jpsfe) = bioma0 *5.E-6114 trn(:,:,:,jpdfe) = bioma0 *5.E-6115 trn(:,:,:,jpnfe) = bioma0 *5.E-6116 trn(:,:,:,jpnch) = bioma0 *12./55.117 trn(:,:,:,jpdch) = bioma0 *12./55.89 trn(:,:,:,jpsfe) = bioma0 * 5.e-6 90 trn(:,:,:,jpdfe) = bioma0 * 5.e-6 91 trn(:,:,:,jpnfe) = bioma0 * 5.e-6 92 trn(:,:,:,jpnch) = bioma0 * 12. / 55. 93 trn(:,:,:,jpdch) = bioma0 * 12. / 55. 118 94 trn(:,:,:,jpno3) = no3 119 95 trn(:,:,:,jpnh4) = bioma0 120 96 121 122 !! Initialization of chemical variables of the carbon cycle 123 !! -------------------------------------------------------- 124 125 DO jk = 1,jpk 126 DO jj = 1,jpj 127 DO ji = 1,jpi 128 caralk = trn(ji,jj,jk,jptal)- & 129 borat(ji,jj,jk)/(1.+1.E-8/(rtrn+akb3(ji,jj,jk))) 130 co3 = (caralk-trn(ji,jj,jk,jpdic))*tmask(ji,jj,jk) & 131 & +(1.-tmask(ji,jj,jk))*.5e-3 132 bicarb = (2.*trn(ji,jj,jk,jpdic)-caralk) 133 hi(ji,jj,jk) = (ak23(ji,jj,jk)*bicarb/co3) & 134 *tmask(ji,jj,jk)+(1.-tmask(ji,jj,jk))*1.e-9 135 ENDDO 136 ENDDO 137 ENDDO 138 139 140 !! initialize the half saturation constant for silicate 141 !! ---------------------------------------------------- 142 143 xksi(:,:)=2.E-6 97 ! Initialization of chemical variables of the carbon cycle 98 ! -------------------------------------------------------- 99 DO jk = 1, jpk 100 DO jj = 1, jpj 101 DO ji = 1, jpi 102 caralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / ( 1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) ) ) 103 co3 = ( caralk - trn(ji,jj,jk,jpdic) ) * tmask(ji,jj,jk) & 104 & + 0.5e-3 * ( 1. - tmask(ji,jj,jk) ) 105 bicarb = (2.*trn(ji,jj,jk,jpdic)-caralk) 106 hi(ji,jj,jk) = ( ak23(ji,jj,jk) * bicarb / co3 ) * tmask(ji,jj,jk) & 107 & + 1.e-9 * ( 1. - tmask(ji,jj,jk) ) 108 END DO 109 END DO 110 END DO 111 112 ! initialize the half saturation constant for silicate 113 ! ---------------------------------------------------- 114 xksi(:,:) = 2.e-6 144 115 145 116 IF(lwp) WRITE(numout,*) 'Initialization of PISCES tracers done' 146 117 IF(lwp) WRITE(numout,*) ' ' 147 118 148 #elif defined key_trc_lobster1 && ( defined key_eel_r6 || defined key_eel_r2 ) 149 ! analytical initialisation used in Levy et al. (2001) 119 # elif defined key_lobster && ( defined key_eel_r6 || defined key_eel_r2 ) 120 ! LOBSTER initialisation for EEL 121 ! ---------------------- 122 ! here: analytical initialisation used in Levy et al. (2001) 150 123 151 DO jk =1,7152 trn(:,:,jk,jpdet)=0.016*tmask(:,:,jk)153 trn(:,:,jk,jpzoo)=0.018*tmask(:,:,jk)154 trn(:,:,jk,jpphy)=0.036*tmask(:,:,jk)155 trn(:,:,jk,jpno3)=1.e-5*tmask(:,:,jk)156 trn(:,:,jk,jpnh4)=0.0005*tmask(:,:,jk)157 trn(:,:,jk,jpdom)=0.017*tmask(:,:,jk)158 END DO 159 160 trn(:,:, 8,jpdet)=0.020*tmask(:,:,1)161 trn(:,:, 8,jpzoo)=0.027*tmask(:,:,1)162 trn(:,:, 8,jpphy)=0.041*tmask(:,:,1)163 trn(:,:, 8,jpno3)=0.00022*tmask(:,:,1)164 trn(:,:, 8,jpnh4)=0.0033*tmask(:,:,1)165 trn(:,:, 8,jpdom)=0.021*tmask(:,:,1)166 167 trn(:,:, 9,jpdet)=0.0556*tmask(:,:,1)168 trn(:,:, 9,jpzoo)=0.123*tmask(:,:,1)169 trn(:,:, 9,jpphy)=0.122*tmask(:,:,1)170 trn(:,:, 9,jpno3)=0.028*tmask(:,:,1)171 trn(:,:, 9,jpnh4)=0.024*tmask(:,:,1)172 trn(:,:, 9,jpdom)=0.06*tmask(:,:,1)173 174 trn(:,:,10,jpdet) =0.025*tmask(:,:,1)175 trn(:,:,10,jpzoo) =0.016*tmask(:,:,1)176 trn(:,:,10,jpphy) =0.029*tmask(:,:,1)177 trn(:,:,10,jpno3) =2.462*tmask(:,:,1)178 trn(:,:,10,jpnh4) =0.04*tmask(:,:,1)179 trn(:,:,10,jpdom) =0.022*tmask(:,:,1)180 181 trn(:,:,11,jpdet) =0.0057*tmask(:,:,1)182 trn(:,:,11,jpzoo) =0.0005*tmask(:,:,1)183 trn(:,:,11,jpphy) =0.0006*tmask(:,:,1)184 trn(:,:,11,jpno3) =3.336*tmask(:,:,1)185 trn(:,:,11,jpnh4) =0.005*tmask(:,:,1)186 trn(:,:,11,jpdom) =0.004*tmask(:,:,1)187 188 trn(:,:,12,jpdet) =0.002*tmask(:,:,1)189 trn(:,:,12,jpzoo) =1.e-6*tmask(:,:,1)190 trn(:,:,12,jpphy) =5.e-6*tmask(:,:,1)191 trn(:,:,12,jpno3) =4.24*tmask(:,:,1)192 trn(:,:,12,jpnh4) =0.001*tmask(:,:,1)193 trn(:,:,12,jpdom) =3.e-5*tmask(:,:,1)124 DO jk = 1, 7 125 trn(:,:,jk,jpdet) = 0.016 * tmask(:,:,jk) 126 trn(:,:,jk,jpzoo) = 0.018 * tmask(:,:,jk) 127 trn(:,:,jk,jpphy) = 0.036 * tmask(:,:,jk) 128 trn(:,:,jk,jpno3) = 1.e-5 * tmask(:,:,jk) 129 trn(:,:,jk,jpnh4) = 5.e-4 * tmask(:,:,jk) 130 trn(:,:,jk,jpdom) = 0.017 * tmask(:,:,jk) 131 END DO 132 133 trn(:,:, 8,jpdet) = 0.020 * tmask(:,:, 8) 134 trn(:,:, 8,jpzoo) = 0.027 * tmask(:,:, 8) 135 trn(:,:, 8,jpphy) = 0.041 * tmask(:,:, 8) 136 trn(:,:, 8,jpno3) = 0.00022 * tmask(:,:, 8) 137 trn(:,:, 8,jpnh4) = 0.0033 * tmask(:,:, 8) 138 trn(:,:, 8,jpdom) = 0.021 * tmask(:,:, 8) 139 140 trn(:,:, 9,jpdet) = 0.0556 * tmask(:,:, 9) 141 trn(:,:, 9,jpzoo) = 0.123 * tmask(:,:, 9) 142 trn(:,:, 9,jpphy) = 0.122 * tmask(:,:, 9) 143 trn(:,:, 9,jpno3) = 0.028 * tmask(:,:, 9) 144 trn(:,:, 9,jpnh4) = 0.024 * tmask(:,:, 9) 145 trn(:,:, 9,jpdom) = 0.06 * tmask(:,:, 9) 146 147 trn(:,:,10,jpdet) = 0.025 * tmask(:,:,10) 148 trn(:,:,10,jpzoo) = 0.016 * tmask(:,:,10) 149 trn(:,:,10,jpphy) = 0.029 * tmask(:,:,10) 150 trn(:,:,10,jpno3) = 2.462 * tmask(:,:,10) 151 trn(:,:,10,jpnh4) = 0.04 * tmask(:,:,10) 152 trn(:,:,10,jpdom) = 0.022 * tmask(:,:,10) 153 154 trn(:,:,11,jpdet) = 0.0057 * tmask(:,:,11) 155 trn(:,:,11,jpzoo) = 0.0005 * tmask(:,:,11) 156 trn(:,:,11,jpphy) = 0.0006 * tmask(:,:,11) 157 trn(:,:,11,jpno3) = 3.336 * tmask(:,:,11) 158 trn(:,:,11,jpnh4) = 0.005 * tmask(:,:,11) 159 trn(:,:,11,jpdom) = 0.004 * tmask(:,:,11) 160 161 trn(:,:,12,jpdet) = 0.002 * tmask(:,:,12) 162 trn(:,:,12,jpzoo) = 1.e-6 * tmask(:,:,12) 163 trn(:,:,12,jpphy) = 5.e-6 * tmask(:,:,12) 164 trn(:,:,12,jpno3) = 4.24 * tmask(:,:,12) 165 trn(:,:,12,jpnh4) = 0.001 * tmask(:,:,12) 166 trn(:,:,12,jpdom) = 3.e-5 * tmask(:,:,12) 194 167 195 168 DO jk=13,jpk 196 trn(:,:,jk,jpdet)=0.0 197 trn(:,:,jk,jpzoo)=0.0 198 trn(:,:,jk,jpphy)=0.0 199 trn(:,:,jk,jpnh4)=0.0 200 trn(:,:,jk,jpdom)=0.0 201 END DO 202 203 trn(:,:,13,jpno3)=5.31*tmask(:,:,13) 204 trn(:,:,14,jpno3)=6.73*tmask(:,:,14) 205 trn(:,:,15,jpno3)=8.32*tmask(:,:,15) 206 trn(:,:,16,jpno3)=10.13*tmask(:,:,16) 207 trn(:,:,17,jpno3)=11.95*tmask(:,:,17) 208 trn(:,:,18,jpno3)=13.57*tmask(:,:,18) 209 trn(:,:,19,jpno3)=15.08*tmask(:,:,19) 210 trn(:,:,20,jpno3)=16.41*tmask(:,:,20) 211 trn(:,:,21,jpno3)=17.47*tmask(:,:,21) 212 trn(:,:,22,jpno3)=18.29*tmask(:,:,22) 213 trn(:,:,23,jpno3)=18.88*tmask(:,:,23) 214 trn(:,:,24,jpno3)=19.30*tmask(:,:,24) 215 trn(:,:,25,jpno3)=19.68*tmask(:,:,25) 216 trn(:,:,26,jpno3)=19.91*tmask(:,:,26) 217 trn(:,:,27,jpno3)=19.99*tmask(:,:,27) 218 trn(:,:,28,jpno3)=20.01*tmask(:,:,28) 219 trn(:,:,29,jpno3)=20.01*tmask(:,:,29) 220 trn(:,:,30,jpno3)=20.01*tmask(:,:,30) 221 222 #elif defined key_trc_lobster1 && defined key_gyre 223 ! init NO3=f(density) by asklod AS Kremeur 2005-07 224 trn(:,:,:,jpdet)=0.1*tmask(:,:,:) 225 trn(:,:,:,jpzoo)=0.1*tmask(:,:,:) 226 trn(:,:,:,jpnh4)=0.1*tmask(:,:,:) 227 trn(:,:,:,jpphy)=0.1*tmask(:,:,:) 228 trn(:,:,:,jpdom)=1.*tmask(:,:,:) 229 DO jk=1,jpk 230 DO jj=1,jpj 231 DO ji=1,jpi 232 IF (rhd(ji,jj,jk).LE.24.5e-3) THEN 233 trn(ji,jj,jk,jpno3)=2.*tmask(ji,jj,jk) 169 trn(:,:,jk,jpdet) = 0.e0 170 trn(:,:,jk,jpzoo) = 0.e0 171 trn(:,:,jk,jpphy) = 0.e0 172 trn(:,:,jk,jpnh4) = 0.e0 173 trn(:,:,jk,jpdom) = 0.e0 174 END DO 175 176 trn(:,:,13,jpno3) = 5.31 * tmask(:,:,13) 177 trn(:,:,14,jpno3) = 6.73 * tmask(:,:,14) 178 trn(:,:,15,jpno3) = 8.32 * tmask(:,:,15) 179 trn(:,:,16,jpno3) = 10.13 * tmask(:,:,16) 180 trn(:,:,17,jpno3) = 11.95 * tmask(:,:,17) 181 trn(:,:,18,jpno3) = 13.57 * tmask(:,:,18) 182 trn(:,:,19,jpno3) = 15.08 * tmask(:,:,19) 183 trn(:,:,20,jpno3) = 16.41 * tmask(:,:,20) 184 trn(:,:,21,jpno3) = 17.47 * tmask(:,:,21) 185 trn(:,:,22,jpno3) = 18.29 * tmask(:,:,22) 186 trn(:,:,23,jpno3) = 18.88 * tmask(:,:,23) 187 trn(:,:,24,jpno3) = 19.30 * tmask(:,:,24) 188 trn(:,:,25,jpno3) = 19.68 * tmask(:,:,25) 189 trn(:,:,26,jpno3) = 19.91 * tmask(:,:,26) 190 trn(:,:,27,jpno3) = 19.99 * tmask(:,:,27) 191 trn(:,:,28,jpno3) = 20.01 * tmask(:,:,28) 192 trn(:,:,29,jpno3) = 20.01 * tmask(:,:,29) 193 trn(:,:,30,jpno3) = 20.01 * tmask(:,:,30) 194 195 # elif defined key_lobster && defined key_gyre 196 ! LOBSTER initialisation for GYRE 197 ! ---------------------- 198 ! here: init NO3=f(density) by asklod AS Kremeur 2005-07 199 trn(:,:,:,jpdet) = 0.1 * tmask(:,:,:) 200 trn(:,:,:,jpzoo) = 0.1 * tmask(:,:,:) 201 trn(:,:,:,jpnh4) = 0.1 * tmask(:,:,:) 202 trn(:,:,:,jpphy) = 0.1 * tmask(:,:,:) 203 trn(:,:,:,jpdom) = 1.0 * tmask(:,:,:) 204 DO jk = 1, jpk 205 DO jj = 1, jpj 206 DO ji = 1, jpi 207 IF( rhd(ji,jj,jk) <= 24.5e-3 ) THEN 208 trn(ji,jj,jk,jpno3) = 2. * tmask(ji,jj,jk) 234 209 ELSE 235 trn(ji,jj,jk,jpno3) =(15.55*(rhd(ji,jj,jk)*1000)-380.11)*tmask(ji,jj,jk)210 trn(ji,jj,jk,jpno3) = ( 15.55 * ( rhd(ji,jj,jk) * 1000. ) - 380.11 ) * tmask(ji,jj,jk) 236 211 ENDIF 237 212 END DO … … 239 214 END DO 240 215 216 # else 217 ! Default case 218 ! ------------ 219 DO jn = 1, jptra 220 trn(:,:,:,jn) = 0.1 * tmask(:,:,:) 221 END DO 222 223 # endif 224 225 # if defined key_dtatrc 226 ! Initialization of tracer from a file that may also be used for damping 227 CALL dta_trc( nittrc000 ) 228 DO jn = 1, jptra 229 IF( lutini(jn) ) trn(:,:,:,jn) = trdta(:,:,:,jn) * tmask(:,:,:) ! initialisation from file if required 230 END DO 231 # endif 232 233 ! before field : 234 ! ------------- 235 trb(:,:,:,:) = trn(:,:,:,:) 236 237 # if defined key_lobster 238 ! initialize the POC in sediments 239 sedpocb(:,:) = 0.e0 240 sedpocn(:,:) = 0.e0 241 sedpoca(:,:) = 0.e0 242 # endif 243 ! 244 END SUBROUTINE trc_dtr 245 241 246 #else 242 243 !! general case 244 do jn = 1, jptra245 trn(:,:,:,jn)=0.1*tmask(:,:,:) 246 enddo247 247 !!---------------------------------------------------------------------- 248 !! Dummy module : No passive tracer 249 !!---------------------------------------------------------------------- 250 CONTAINS 251 SUBROUTINE trc_dtr ! Empty routine 252 END SUBROUTINE trc_dtr 248 253 #endif 249 254 250 #if defined key_dtatrc 251 !! Initialization of tracer from a file 252 !! that may also be used for damping 253 CALL dta_trc( nittrc000 ) 254 DO jk = 1, jptra 255 IF( lutini(jk) ) THEN 256 !! initialisation from file 257 trn(:,:,:,jk) = trdta(:,:,:,jk)*tmask(:,:,:) 258 ENDIF 259 END DO 260 #endif 261 262 !! before field : 263 !! ------------- 264 trb(:,:,:,:) = trn(:,:,:,:) 265 266 #if defined key_trc_lobster1 267 !! initialize the POC in sediments 268 269 sedpocb(:,:) = 0. 270 sedpocn(:,:) = 0. 271 sedpoca(:,:) = 0. 272 #endif 273 274 END SUBROUTINE trc_dtr 275 276 #else 277 278 SUBROUTINE trc_dtr 279 !!====================== 280 !! no passive tracers 281 !!====================== 282 END SUBROUTINE trc_dtr 283 #endif 284 255 !!====================================================================== 285 256 END MODULE trcdtr
Note: See TracChangeset
for help on using the changeset viewer.