Changeset 763 for branches/dev_001_GM/NEMO/TOP_SRC/SMS/trcfreons.F90
- Timestamp:
- 2007-12-13T14:52:50+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/TOP_SRC/SMS/trcfreons.F90
r719 r763 1 1 MODULE trcfreons 2 !!============================================================== 3 !! *** MODULE trcfreons *** 4 !! Passive tracer : CFC main model 5 !!============================================================== 2 !!====================================================================== 3 !! *** MODULE trcfreons *** 4 !! TOP : CFC main model 5 !!====================================================================== 6 !! History : - ! 1999-10 (JC. Dutay) original code 7 !! 1.0 ! 2004-03 (C. Ethe) free form + modularity 8 !!---------------------------------------------------------------------- 6 9 #if defined key_cfc 7 !!-------------------------------------------------------------- 8 !! 'key_cfc' CFC model 9 !!-------------------------------------------------------------- 10 !! * Modules used 10 !!---------------------------------------------------------------------- 11 !! 'key_cfc' CFC chemical model 12 !!---------------------------------------------------------------------- 13 !! trc_freons : compute and add CFC suface forcing to CFC trends 14 !! trc_freons_cst : sets constants for CFC surface forcing computation 15 !!---------------------------------------------------------------------- 11 16 USE daymod 12 17 USE sms … … 14 19 USE trc 15 20 16 17 21 IMPLICIT NONE 18 22 PRIVATE 19 23 20 !! * Routine accessibility 21 PUBLIC trc_freons 22 23 !! * Module variables 24 REAL(wp), DIMENSION(jptra) :: & ! coefficient for solubility of CFC11 in mol/l/atm 25 soa1, soa2, soa3, soa4, & 26 sob1, sob2, sob3 27 28 REAL(wp), DIMENSION(jptra) :: & ! coefficients for schmidt number in degre Celcius 29 sca1, sca2, sca3, sca4 30 31 REAL(wp) :: & ! coefficients for conversion 32 xconv1 = 1.0 , & ! conversion from to 33 xconv2 = 0.01/3600., & ! conversion from cm/h to m/s: 34 xconv3 = 1.0e+3 , & ! conversion from mol/l/atm to mol/m3/atm 35 xconv4 = 1.0e-12 ! conversion from mol/m3/atm to mol/m3/pptv 24 PUBLIC trc_freons ! called in ??? 25 26 REAL(wp), DIMENSION(jptra) :: soa1, soa2, soa3, soa4 ! coefficient for solubility of CFC [mol/l/atm] 27 REAL(wp), DIMENSION(jptra) :: sob1, sob2, sob3 ! " " 28 REAL(wp), DIMENSION(jptra) :: sca1, sca2, sca3, sca4 ! coefficients for schmidt number in degre Celcius 29 30 ! ! coefficients for conversion 31 REAL(wp) :: xconv1 = 1.0 ! conversion from to 32 REAL(wp) :: xconv2 = 0.01/3600. ! conversion from cm/h to m/s: 33 REAL(wp) :: xconv3 = 1.0e+3 ! conversion from mol/l/atm to mol/m3/atm 34 REAL(wp) :: xconv4 = 1.0e-12 ! conversion from mol/m3/atm to mol/m3/pptv 36 35 37 36 !! * Substitutions 38 37 # include "passivetrc_substitute.h90" 39 40 !!---------------------------------------------------------------------- 41 !! TOP 1.0 , LOCEAN-IPSL (2005) 42 !! $Header$ 43 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 38 !!---------------------------------------------------------------------- 39 !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005) 40 !! $Id:$ 41 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 44 42 !!---------------------------------------------------------------------- 45 43 … … 62 60 !! - the input function is in pico-mol/m3/s and the 63 61 !! freons concentration in pico-mol/m3 64 !!65 !! History :66 !! 8.1 ! 99-10 (JC. Dutay) original code67 !! 9.0 ! 04-03 (C. Ethe) free form + modularity68 62 !!---------------------------------------------------------------------- 69 !! * Arguments70 63 INTEGER, INTENT( in ) :: kt ! ocean time-step index 71 72 !! * Local declarations 73 INTEGER :: & 74 ji, jj, jn, jm 75 76 INTEGER :: & 77 iyear_beg, iyear_end, & 78 imonth, im1, im2 79 80 REAL(wp) :: & 81 ztap, zdtap, & 82 zt1, zt2, zt3, zv2 83 84 REAL(wp), DIMENSION(jphem,jptra) :: & 85 zpatm ! atmospheric function 86 87 REAL(wp) :: & 88 zsol, & ! solubility 89 zsch ! schmidt number 90 64 !! 65 INTEGER :: ji, jj, jn, jm 66 INTEGER :: iyear_beg, iyear_end 67 INTEGER :: imonth, im1, im2 68 69 REAL(wp) :: ztap, zdtap 70 REAL(wp) :: zt1, zt2, zt3, zv2 71 REAL(wp) :: zsol ! solubility 72 REAL(wp) :: zsch ! schmidt number 91 73 92 REAL(wp), DIMENSION(jpi,jpj,jptra) :: & 93 zca_cfc, & ! concentration 94 zak_cfc ! transfert coefficients 95 74 REAL(wp), DIMENSION(jphem,jptra) :: zpatm ! atmospheric function 75 REAL(wp), DIMENSION(jpi,jpj,jptra) :: zca_cfc ! concentration 76 REAL(wp), DIMENSION(jpi,jpj,jptra) :: zak_cfc ! transfert coefficients 96 77 !!---------------------------------------------------------------------- 97 98 78 99 79 IF( kt == nittrc000 ) CALL trc_freons_cst … … 117 97 118 98 119 120 121 99 ! Temporal and spatial interpolation at time k 122 100 ! -------------------------------------------------- … … 125 103 zpatm(jm,jn) = ( p_cfc(iyear_beg, jm, jn) * FLOAT (im1) & 126 104 & + p_cfc(iyear_end, jm, jn) * FLOAT (im2) ) / 12. 127 ENDDO 128 END DO 129 130 DO jn = 1, jptra 131 DO jj = 1, jpj 132 DO ji = 1, jpi 133 pp_cfc(ji,jj,jn) = xphem(ji,jj) * zpatm(1,jn) & 134 & + ( 1.- xphem(ji,jj) ) * zpatm(2,jn) 135 END DO 136 END DO 137 ENDDO 105 END DO 106 END DO 107 108 DO jn = 1, jptra 109 pp_cfc(:,:,jn) = xphem(:,:) * zpatm(1,jn) & 110 & + ( 1.- xphem(:,:) ) * zpatm(2,jn) 111 END DO 138 112 139 113 … … 160 134 END DO 161 135 END DO 162 END DO136 END DO 163 137 164 138 … … 181 155 END DO 182 156 END DO 183 END DO157 END DO 184 158 185 159 … … 199 173 END DO 200 174 END DO 201 END DO175 END DO 202 176 203 177 … … 212 186 END DO 213 187 END DO 214 END DO188 END DO 215 189 216 190 ! -------------------------------------------- … … 223 197 END DO 224 198 END DO 225 ENDDO 226 227 199 END DO 200 ! 228 201 END SUBROUTINE trc_freons 202 229 203 230 204 SUBROUTINE trc_freons_cst … … 232 206 !! *** trc_freons_cst *** 233 207 !! 234 !! Purpose : sets constants for CFC model 235 !! --------- 236 !! 237 !! 238 !! History : 239 !! 8.2 ! 04-06 (JC. Dutay) original code 240 !! 9.0 ! 05-10 (C. Ethe) Modularity 208 !! ** Purpose : sets constants for CFC model 241 209 !!--------------------------------------------------------------------- 242 !! TOP 1.0 , LOCEAN-IPSL (2005) 243 !! $Header$ 244 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 245 !!---------------------------------------------------------------- 246 !! Local declarations 247 INTEGER :: jn 210 INTEGER :: jn 211 !!--------------------------------------------------------------------- 248 212 249 213 DO jn = 1, jptra … … 290 254 WRITE(numout,*) 'coefficient for schmidt of tracer',ctrcnm(jn) 291 255 WRITE(numout,*) sca1(jn), sca2(jn),sca3(jn), sca4(jn) 292 END DO293 256 END DO 257 ! 294 258 END SUBROUTINE trc_freons_cst 259 295 260 #else 296 261 !!---------------------------------------------------------------------- 297 !! D efault option Dummy module262 !! Dummy module No CFC model 298 263 !!---------------------------------------------------------------------- 299 264 CONTAINS
Note: See TracChangeset
for help on using the changeset viewer.
