Changeset 10975 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/CFC
- Timestamp:
- 2019-05-13T18:34:33+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/CFC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/CFC/trcini_cfc.F90
r10068 r10975 31 31 CONTAINS 32 32 33 SUBROUTINE trc_ini_cfc 33 SUBROUTINE trc_ini_cfc( Kmm ) 34 34 !!---------------------------------------------------------------------- 35 35 !! *** trc_ini_cfc *** … … 39 39 !! ** Method : - Read the namcfc namelist and check the parameter values 40 40 !!---------------------------------------------------------------------- 41 INTEGER, INTENT(in) :: Kmm ! time level indices 41 42 INTEGER :: ji, jj, jn, jl, jm, js, io, ierr 42 INTEGER :: iskip = 6 ! number of 1st descriptor lines43 INTEGER :: iskip = 6 ! number of 1st descriptor lines 43 44 REAL(wp) :: zyy, zyd 44 45 CHARACTER(len = 20) :: cltra … … 90 91 DO jl = 1, jp_cfc 91 92 jn = jp_cfc0 + jl - 1 92 tr n(:,:,:,jn) = 0._wp93 tr(:,:,:,jn,Kmm) = 0._wp 93 94 END DO 94 95 ENDIF -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/CFC/trcsms_cfc.F90
r10966 r10975 54 54 CONTAINS 55 55 56 SUBROUTINE trc_sms_cfc( kt, K mm)56 SUBROUTINE trc_sms_cfc( kt, Kbb, Kmm, Krhs ) 57 57 !!---------------------------------------------------------------------- 58 58 !! *** ROUTINE trc_sms_cfc *** … … 70 70 !! CFC concentration in pico-mol/m3 71 71 !!---------------------------------------------------------------------- 72 INTEGER, INTENT(in) :: kt ! ocean time-step index73 INTEGER, INTENT(in) :: K mm! ocean time level72 INTEGER, INTENT(in) :: kt ! ocean time-step index 73 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! ocean time level 74 74 ! 75 75 INTEGER :: ji, jj, jn, jl, jm … … 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 = ( ts n(ji,jj,1,jp_tem) + 273.16 ) * 0.01131 ztap = ( ts(ji,jj,1,jp_tem,Kmm) + 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 + ts n(ji,jj,1,jp_sal) * zdtap )134 & + soa(4,jl) * ztap * ztap + ts(ji,jj,1,jp_sal,Kmm) * zdtap ) 135 135 ELSE 136 136 zsol = 0.e0 … … 143 143 ! Computation of speed transfert 144 144 ! Schmidt number revised in Wanninkhof (2014) 145 zt1 = ts n(ji,jj,1,jp_tem)145 zt1 = ts(ji,jj,1,jp_tem,Kmm) 146 146 zt2 = zt1 * zt1 147 147 zt3 = zt1 * zt2 … … 155 155 156 156 ! Input function : speed *( conc. at equil - concen at surface ) 157 ! tr nin pico-mol/l idem qtr; ak in en m/a158 qtr_cfc(ji,jj,jl) = -zak_cfc * ( tr b(ji,jj,1,jn) - zca_cfc ) &157 ! tr(:,:,:,:,Kmm) in pico-mol/l idem qtr; ak in en m/a 158 qtr_cfc(ji,jj,jl) = -zak_cfc * ( tr(ji,jj,1,jn,Kbb) - zca_cfc ) & 159 159 & * tmask(ji,jj,1) * ( 1. - fr_i(ji,jj) ) 160 160 ! Add the surface flux to the trend 161 tr a(ji,jj,1,jn) = tra(ji,jj,1,jn) + qtr_cfc(ji,jj,jl) / e3t_n(ji,jj,1)161 tr(ji,jj,1,jn,Krhs) = tr(ji,jj,1,jn,Krhs) + qtr_cfc(ji,jj,jl) / e3t(ji,jj,1,Kmm) 162 162 163 163 ! cumulation of surface flux at each time step … … 192 192 IF( l_trdtrc ) THEN 193 193 DO jn = jp_cfc0, jp_cfc1 194 CALL trd_trc( tr a(:,:,:,jn), jn, jptra_sms, kt, Kmm ) ! save trends194 CALL trd_trc( tr(:,:,:,jn,Krhs), jn, jptra_sms, kt, Kmm ) ! save trends 195 195 END DO 196 196 END IF -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/CFC/trcwri_cfc.F90
r10069 r10975 20 20 CONTAINS 21 21 22 SUBROUTINE trc_wri_cfc 22 SUBROUTINE trc_wri_cfc( Kmm ) 23 23 !!--------------------------------------------------------------------- 24 24 !! *** ROUTINE trc_wri_trc *** … … 26 26 !! ** Purpose : output passive tracers fields 27 27 !!--------------------------------------------------------------------- 28 INTEGER, INTENT(in) :: Kmm ! time level indices 28 29 CHARACTER (len=20) :: cltra 29 30 INTEGER :: jn … … 34 35 DO jn = jp_cfc0, jp_cfc1 35 36 cltra = TRIM( ctrcnm(jn) ) ! short title for tracer 36 CALL iom_put( cltra, tr n(:,:,:,jn) )37 CALL iom_put( cltra, tr(:,:,:,jn,Kmm) ) 37 38 END DO 38 39 !
Note: See TracChangeset
for help on using the changeset viewer.