Possible typos in p4zsed.F90
Hi guys, not sure if you experienced the same problem but I noticed weird values in some biogeochemical parameters and I think this is related to a few errors in the new version. See below the original (commented out) and my modifications (below). I would appreciate any comments regarding the validity of this fix (or maybe I just got lucky... as the weird values are gone):
DO ji = IRANGE
zdep = rfact2 / fse3t(ji,jj,KSURF)
! irondep(ji,jj,1) = dustsolub & ! & * mfrac * dust(ji,jj) / 55.85 * zdep & ! & + 3.e-10 / ryyss * zdep !! ! irondep(ji,jj,1) = (dustsolub & & * mfrac * dust(ji,jj) / (55.85 * rmtss) & & + 3.e-10 / ryyss) * zdep ! !
! sidep (ji,jj) = 8.8 * 0.075 & ! & * mfrac * dust(ji,jj) / 28.1 * zdep sidep (ji,jj) = 8.8 * 0.075 & & * mfrac * dust(ji,jj) / 28.1 /rmtss * zdep
!
! po4dep (ji,jj) = 0.1 * 0.021 & ! & * mfrac * dust(ji,jj) / 31. / po4r * zdep po4dep (ji,jj) = 0.1 * 0.021 & & * mfrac * dust(ji,jj) / 31. /rmtss / po4r * zdep
!
tra(ji,jj,1,jpsil) = tra(ji,jj,1,jpsil) + sidep (ji,jj)
tra(ji,jj,1,jppo4) = tra(ji,jj,1,jppo4) + po4dep(ji,jj)
tra(ji,jj,1,jpfer) = tra(ji,jj,1,jpfer) + irondep(ji,jj,1)
END DO
END DO