Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/CFC
- Timestamp:
- 2010-12-27T18:33:53+01:00 (14 years ago)
- Location:
- trunk/NEMOGCM/NEMO/TOP_SRC/CFC
- Files:
-
- 2 deleted
- 4 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/CFC/par_cfc.F90
r2047 r2528 6 6 !! History : 2.0 ! 2007-12 (C. Ethe, G. Madec) revised architecture 7 7 !!---------------------------------------------------------------------- 8 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)8 !! NEMO/TOP 3.3 , NEMO Consortium (2010) 9 9 !! $Id$ 10 !! Software governed by the CeCILL licence ( modipsl/doc/NEMO_CeCILL.txt)10 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 11 11 !!---------------------------------------------------------------------- 12 12 USE par_lobster, ONLY : jp_lobster !: number of tracers in LOBSTER -
trunk/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r2047 r2528 30 30 31 31 !!---------------------------------------------------------------------- 32 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)32 !! NEMO/TOP 3.3 , NEMO Consortium (2010) 33 33 !! $Id$ 34 !! Software governed by the CeCILL licence ( modipsl/doc/NEMO_CeCILL.txt)34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 35 !!---------------------------------------------------------------------- 36 36 -
trunk/NEMOGCM/NEMO/TOP_SRC/CFC/trcrst_cfc.F90
- Property svn:keywords set to Id
-
trunk/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r2423 r2528 18 18 USE par_trc ! TOP parameters 19 19 USE trc ! TOP variables 20 USE trdm ld_trc_oce21 USE trdm ld_trc20 USE trdmod_oce 21 USE trdmod_trc 22 22 USE iom 23 23 … … 34 34 INTEGER , PUBLIC :: npyear ! Number of years read in CFC1112 file 35 35 36 REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, 2) :: p_cfc ! partial hemispheric pressure for CFC37 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: xphem ! spatial interpolation factor for patm36 REAL(wp), PUBLIC, DIMENSION(jpyear,jphem, 2 ) :: p_cfc ! partial hemispheric pressure for CFC 37 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: xphem ! spatial interpolation factor for patm 38 38 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp_cfc) :: qtr_cfc ! flux at surface 39 39 REAL(wp), PUBLIC, DIMENSION(jpi,jpj ,jp_cfc) :: qint_cfc ! cumulative flux … … 52 52 # include "top_substitute.h90" 53 53 !!---------------------------------------------------------------------- 54 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)54 !! NEMO/TOP 3.3 , NEMO Consortium (2010) 55 55 !! $Id$ 56 !! Software governed by the CeCILL licence ( modipsl/doc/NEMO_CeCILL.txt)56 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 57 57 !!---------------------------------------------------------------------- 58 58 … … 93 93 !!---------------------------------------------------------------------- 94 94 95 IF( kt == nit trc000 ) CALL trc_cfc_cst95 IF( kt == nit000 ) CALL trc_cfc_cst 96 96 97 97 ! Temporal interpolation … … 129 129 ! coefficient for solubility for CFC-11/12 in mol/l/atm 130 130 IF( tmask(ji,jj,1) .GE. 0.5 ) THEN 131 ztap = ( t n(ji,jj,1) + 273.16 ) * 0.01131 ztap = ( tsn(ji,jj,1,jp_tem) + 273.16 ) * 0.01 132 132 zdtap = sob(1,jl) + ztap * ( sob(2,jl) + ztap * sob(3,jl) ) 133 133 zsol = EXP( soa(1,jl) + soa(2,jl) / ztap + soa(3,jl) * LOG( ztap ) & 134 & + soa(4,jl) * ztap * ztap + sn(ji,jj,1) * zdtap )134 & + soa(4,jl) * ztap * ztap + tsn(ji,jj,1,jp_sal) * zdtap ) 135 135 ELSE 136 136 zsol = 0.e0 … … 143 143 ! Computation of speed transfert 144 144 ! Schmidt number 145 zt1 = t n(ji,jj,1)145 zt1 = tsn(ji,jj,1,jp_tem) 146 146 zt2 = zt1 * zt1 147 147 zt3 = zt1 * zt2 … … 156 156 ! trn in pico-mol/l idem qtr; ak in en m/s 157 157 qtr_cfc(ji,jj,jl) = -zak_cfc * ( trb(ji,jj,1,jn) - zca_cfc ) & 158 #if defined key_ off_degrad158 #if defined key_degrad 159 159 & * facvol(ji,jj,1) & 160 160 #endif … … 173 173 ! !----------------! 174 174 175 #if defined key_ trc_diaadd175 #if defined key_diatrc 176 176 ! Save diagnostics , just for CFC11 177 177 # if ! defined key_iomput … … 187 187 DO jn = jp_cfc0, jp_cfc1 188 188 ztrcfc(:,:,:) = tra(:,:,:,jn) 189 CALL trd_mod_trc( ztrcfc, jn, jptr c_trd_sms, kt ) ! save trends189 CALL trd_mod_trc( ztrcfc, jn, jptra_trd_sms, kt ) ! save trends 190 190 END DO 191 191 END IF
Note: See TracChangeset
for help on using the changeset viewer.