Changeset 6140 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
r5836 r6140 30 30 PUBLIC p4z_fechem_init ! called in trcsms_pisces.F90 31 31 32 !! * Shared module variables33 LOGICAL :: ln_fechem !: boolean for complex iron chemistryfollowing Tagliabue and voelker34 LOGICAL :: ln_ligvar !: boolean for variable ligand concentration following Tagliabue and voelker35 REAL(wp), PUBLIC :: xlam1 !: scavenging rate of Iron36 REAL(wp), PUBLIC :: xlamdust !: scavenging rate of Iron by dust37 REAL(wp), PUBLIC :: ligand !: ligand concentration in the ocean 38 32 LOGICAL :: ln_fechem !: boolean for complex iron chemistry following Tagliabue and voelker 33 LOGICAL :: ln_ligvar !: boolean for variable ligand concentration following Tagliabue and voelker 34 REAL(wp), PUBLIC :: xlam1 !: scavenging rate of Iron 35 REAL(wp), PUBLIC :: xlamdust !: scavenging rate of Iron by dust 36 REAL(wp), PUBLIC :: ligand !: ligand concentration in the ocean 37 38 !!gm Not DOCTOR norm !!! 39 39 REAL(wp) :: kl1, kl2, kb1, kb2, ks, kpr, spd, con, kth 40 40 41 !! * Substitutions42 # include "domzgr_substitute.h90"43 41 !!---------------------------------------------------------------------- 44 42 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 61 59 !! and one particulate form (ln_fechem) 62 60 !!--------------------------------------------------------------------- 63 ! 64 INTEGER, INTENT(in) :: kt, knt ! ocean time step 61 INTEGER, INTENT(in) :: kt, knt ! ocean time step 65 62 ! 66 63 INTEGER :: ji, jj, jk, jic 64 CHARACTER (len=25) :: charout 67 65 REAL(wp) :: zdep, zlam1a, zlam1b, zlamfac 68 66 REAL(wp) :: zkeq, zfeequi, zfesatur, zfecoll … … 79 77 REAL(wp) :: ztfe, zoxy 80 78 REAL(wp) :: zstep 81 CHARACTER (len=25) :: charout82 79 !!--------------------------------------------------------------------- 83 80 ! 84 81 IF( nn_timing == 1 ) CALL timing_start('p4z_fechem') 85 82 ! 86 ! Allocate temporary workspace 87 CALL wrk_alloc( jpi, jpj, jpk, zFe3, zFeL1, zTL1, ztotlig ) 83 CALL wrk_alloc( jpi,jpj,jpk, zFe3, zFeL1, zTL1, ztotlig ) 88 84 zFe3 (:,:,:) = 0. 89 85 zFeL1(:,:,:) = 0. 90 86 zTL1 (:,:,:) = 0. 91 87 IF( ln_fechem ) THEN 92 CALL wrk_alloc( jpi, jpj, jpk,zFe2, zFeL2, zTL2, zFeP )88 CALL wrk_alloc( jpi,jpj,jpk, zFe2, zFeL2, zTL2, zFeP ) 93 89 zFe2 (:,:,:) = 0. 94 90 zFeL2(:,:,:) = 0. … … 253 249 zlamfac = MAX( 0.e0, ( gphit(ji,jj) + 55.) / 30. ) 254 250 zlamfac = MIN( 1. , zlamfac ) 255 zdep = MIN( 1., 1000. / fsdept(ji,jj,jk) ) 251 !!gm very small BUG : it is unlikely but possible that gdept_n = 0 ..... 252 zdep = MIN( 1., 1000. / gdept_n(ji,jj,jk) ) 256 253 zlam1b = xlam1 * MAX( 0.e0, ( trb(ji,jj,jk,jpfer) * 1.e9 - ztotlig(ji,jj,jk) ) ) 257 254 zcoag = zfeequi * zlam1b * zstep + 1E-4 * ( 1. - zlamfac ) * zdep * zstep * trb(ji,jj,jk,jpfer)
Note: See TracChangeset
for help on using the changeset viewer.