Changeset 10975 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zligand.F90
- Timestamp:
- 2019-05-13T18:34:33+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zligand.F90
r10416 r10975 33 33 CONTAINS 34 34 35 SUBROUTINE p4z_ligand( kt, knt )35 SUBROUTINE p4z_ligand( kt, knt, Kbb, Krhs ) 36 36 !!--------------------------------------------------------------------- 37 37 !! *** ROUTINE p4z_ligand *** … … 39 39 !! ** Purpose : Compute remineralization/scavenging of organic ligands 40 40 !!--------------------------------------------------------------------- 41 INTEGER, INTENT(in) :: kt, knt ! ocean time step 41 INTEGER, INTENT(in) :: kt, knt ! ocean time step 42 INTEGER, INTENT(in) :: Kbb, Krhs ! time level indices 42 43 ! 43 44 INTEGER :: ji, jj, jk … … 62 63 ! This is based on the idea that as LGW is lower 63 64 ! there is a larger fraction of refractory OM 64 zlgwr = max( rlgs , rlgw * exp( -2 * (tr b(ji,jj,jk,jplgw)*1e9) ) ) ! years65 zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * blim(ji,jj,jk) * tr b(ji,jj,jk,jplgw)65 zlgwr = max( rlgs , rlgw * exp( -2 * (tr(ji,jj,jk,jplgw,Kbb)*1e9) ) ) ! years 66 zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * blim(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb) 66 67 ! photochem loss of weak ligand 67 zlgwpr = prlgw * xstep * etot(ji,jj,jk) * tr b(ji,jj,jk,jplgw) * (1. - fr_i(ji,jj))68 tr a(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zlgwp - zlgwr - zlgwpr68 zlgwpr = prlgw * xstep * etot(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb) * (1. - fr_i(ji,jj)) 69 tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs) + zlgwp - zlgwr - zlgwpr 69 70 zligrem(ji,jj,jk) = zlgwr 70 71 zligpr(ji,jj,jk) = zlgwpr … … 97 98 WRITE(charout, FMT="('ligand1')") 98 99 CALL prt_ctl_trc_info(charout) 99 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)100 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 100 101 ENDIF 101 102 !
Note: See TracChangeset
for help on using the changeset viewer.