Changeset 1004 for trunk/NEMO/TOP_SRC/CFC
- Timestamp:
- 2008-05-30T11:47:14+02:00 (16 years ago)
- Location:
- trunk/NEMO/TOP_SRC/CFC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/CFC/par_cfc.F90
r933 r1004 24 24 LOGICAL, PUBLIC, PARAMETER :: lk_cfc = .TRUE. !: CFC flag 25 25 INTEGER, PUBLIC, PARAMETER :: jp_cfc = 2 !: number of passive tracers 26 INTEGER, PUBLIC, PARAMETER :: jp_cfc_2d = 1!: additional 2d output arrays ('key_trc_diaadd')27 INTEGER, PUBLIC, PARAMETER :: jp_cfc_3d = 1!: additional 3d output arrays ('key_trc_diaadd')26 INTEGER, PUBLIC, PARAMETER :: jp_cfc_2d = 0 !: additional 2d output arrays ('key_trc_diaadd') 27 INTEGER, PUBLIC, PARAMETER :: jp_cfc_3d = 0 !: additional 3d output arrays ('key_trc_diaadd') 28 28 INTEGER, PUBLIC, PARAMETER :: jp_cfc_trd = 0 !: number of sms trends for CFC 29 29 -
trunk/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r933 r1004 52 52 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~' 53 53 54 54 55 ! Initialization of boundaries conditions 55 56 ! --------------------------------------- … … 73 74 ENDIF 74 75 DO jn = jp_cfc0, jp_cfc1 75 qint(:,:,jn) = 0.e0 76 trn(:,:,:,jn) = 0.e0 77 qint(:,: ,jn) = 0.e0 76 78 END DO 77 79 ENDIF … … 134 136 END DO 135 137 ! 138 139 IF(lwp) WRITE(numout,*) 'Initialization of CFC tracers done' 140 IF(lwp) WRITE(numout,*) ' ' 141 136 142 END SUBROUTINE trc_ini_cfc 137 143 -
trunk/NEMO/TOP_SRC/CFC/trclsm_cfc.F90
r933 r1004 42 42 !! ** input : Namelist namcfc 43 43 !!---------------------------------------------------------------------- 44 CHARACTER (len=32) :: clname = 'namelist. trc.sms'44 CHARACTER (len=32) :: clname = 'namelist.sms.cfc' 45 45 INTEGER :: numnat 46 46 !! -
trunk/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r933 r1004 36 36 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp_cfc0:jp_cfc1) :: qint ! flux function 37 37 38 REAL(wp), DIMENSION(jp_cfc ) :: soa1, soa2, soa3, soa4 ! coefficient for solubility of CFC [mol/l/atm]39 REAL(wp), DIMENSION(jp_cfc ) :: sob1, sob2, sob3 ! " "40 REAL(wp), DIMENSION(jp_cfc ) :: sca1, sca2, sca3, sca4 ! coefficients for schmidt number in degre Celcius38 REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) :: soa1, soa2, soa3, soa4 ! coefficient for solubility of CFC [mol/l/atm] 39 REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) :: sob1, sob2, sob3 ! " " 40 REAL(wp), DIMENSION(jp_cfc0:jp_cfc1) :: sca1, sca2, sca3, sca4 ! coefficients for schmidt number in degre Celcius 41 41 42 42 ! ! coefficients for conversion … … 45 45 REAL(wp) :: xconv3 = 1.0e+3 ! conversion from mol/l/atm to mol/m3/atm 46 46 REAL(wp) :: xconv4 = 1.0e-12 ! conversion from mol/m3/atm to mol/m3/pptv 47 48 REAL(wp) :: xrhoa = 1.22 ! Air density kg/m349 REAL(wp) :: xcd = 1.5e-3 ! drag coefficient50 47 51 48 !! * Substitutions … … 88 85 REAL(wp) :: zca_cfc ! concentration at equilibrium 89 86 REAL(wp) :: zak_cfc ! transfert coefficients 90 91 REAL(wp) :: ztx, zty, ztau92 REAL(wp), DIMENSION(jpi,jpj) :: zws ! wind speed93 87 94 88 REAL(wp), DIMENSION(jphem,jp_cfc) :: zpatm ! atmospheric function … … 110 104 ENDIF 111 105 iyear_end = iyear_beg + 1 112 113 114 ! Estimation of wind speed as a function of wind stress115 !CDIR NOVERRCHK116 DO jj = 1, jpj117 !CDIR NOVERRCHK118 DO ji = 1, jpi119 ztx = utau(ji,jj) * umask(ji,jj,1)120 zty = vtau(ji,jj) * vmask(ji,jj,1)121 ztau = SQRT( ztx * ztx + zty * zty )122 zws(ji,jj) = SQRT ( ztau / ( xrhoa * xcd ) )123 ENDDO124 ENDDO125 106 126 107 ! !------------! … … 164 145 165 146 ! speed transfert : formulae of wanninkhof 1992 166 zv2 = zws(ji,jj) * zws(ji,jj)147 zv2 = wndm(ji,jj) * wndm(ji,jj) 167 148 zsch = zsch / 660. 168 149 zak_cfc = ( 0.39 * xconv2 * zv2 / SQRT(zsch) ) * tmask(ji,jj,1) … … 189 170 END SUBROUTINE trc_sms_cfc 190 171 191 192 172 SUBROUTINE trc_cfc_cst 193 173 !!--------------------------------------------------------------------- … … 196 176 !! ** Purpose : sets constants for CFC model 197 177 !!--------------------------------------------------------------------- 198 INTEGER :: jn 199 !!--------------------------------------------------------------------- 200 201 DO jn = jp_cfc0, jp_cfc1 202 IF ( jn == jp11 ) THEN 203 ! coefficient for solubility of CFC11 in mol/l/atm 204 soa1(jn) = -229.9261 205 soa2(jn) = 319.6552 206 soa3(jn) = 119.4471 207 soa4(jn) = -1.39165 208 sob1(jn) = -0.142382 209 sob2(jn) = 0.091459 210 sob3(jn) = -0.0157274 211 212 ! coefficients for schmidt number in degre Celcius 213 sca1(jn) = 3501.8 214 sca2(jn) = -210.31 215 sca3(jn) = 6.1851 216 sca4(jn) = -0.07513 217 218 ELSE IF( jn == jp12 ) THEN 219 220 ! coefficient for solubility of CFC12 in mol/l/atm 221 soa1(jn) = -218.0971 222 soa2(jn) = 298.9702 223 soa3(jn) = 113.8049 224 soa4(jn) = -1.39165 225 sob1(jn) = -0.143566 226 sob2(jn) = 0.091015 227 sob3(jn) = -0.0153924 228 229 ! coefficients for schmidt number in degre Celcius 230 sca1(jn) = 3845.4 231 sca2(jn) = -228.95 232 sca3(jn) = 6.1908 233 sca4(jn) = -0.067430 234 ENDIF 235 236 WRITE(numout,*) 'coefficient for solubility of tracer',ctrcnm(jn) 237 WRITE(numout,*) soa1(jn), soa2(jn),soa3(jn), soa4(jn), & 238 & sob1(jn), sob2(jn),sob3(jn) 239 WRITE(numout,*) 240 WRITE(numout,*) 'coefficient for schmidt of tracer',ctrcnm(jn) 241 WRITE(numout,*) sca1(jn), sca2(jn),sca3(jn), sca4(jn) 242 END DO 243 ! 178 179 ! coefficient for solubility of CFC11/CFC12 in mol/l/atm 180 181 soa1(jp11) = -229.9261 182 soa2(jp11) = 319.6552 183 soa3(jp11) = 119.4471 184 soa4(jp11) = -1.39165 185 sob1(jp11) = -0.142382 186 sob2(jp11) = 0.091459 187 sob3(jp11) = -0.0157274 188 189 soa1(jp12) = -218.0971 190 soa2(jp12) = 298.9702 191 soa3(jp12) = 113.8049 192 soa4(jp12) = -1.39165 193 sob1(jp12) = -0.143566 194 sob2(jp12) = 0.091015 195 sob3(jp12) = -0.0153924 196 197 198 ! coefficients for schmidt number in degre Celcius 199 sca1(jp11) = 3501.8 200 sca2(jp11) = -210.31 201 sca3(jp11) = 6.1851 202 sca4(jp11) = -0.07513 203 204 sca1(jp12) = 3845.4 205 sca2(jp12) = -228.95 206 sca3(jp12) = 6.1908 207 sca4(jp12) = -0.067430 208 244 209 END SUBROUTINE trc_cfc_cst 245 210
Note: See TracChangeset
for help on using the changeset viewer.