- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90
r5385 r6808 59 59 REAL(wp) :: xconv = 0.01_wp / 3600._wp !: coefficients for conversion 60 60 61 !!* Substitution62 # include "top_substitute.h90"63 61 !!---------------------------------------------------------------------- 64 62 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 122 120 123 121 DO jm = 1, 10 124 !CDIR NOVERRCHK125 122 DO jj = 1, jpj 126 !CDIR NOVERRCHK127 123 DO ji = 1, jpi 128 124 … … 155 151 ! ------------------------------------------- 156 152 157 !CDIR NOVERRCHK158 153 DO jj = 1, jpj 159 !CDIR NOVERRCHK160 154 DO ji = 1, jpi 161 155 ztc = MIN( 35., tsn(ji,jj,1,jp_tem) ) … … 186 180 oce_co2(ji,jj) = ( zfld - zflu ) * rfact2 * e1e2t(ji,jj) * tmask(ji,jj,1) * 1000. 187 181 ! compute the trend 188 tra(ji,jj,1,jpdic) = tra(ji,jj,1,jpdic) + ( zfld - zflu ) * rfact2 / fse3t(ji,jj,1)182 tra(ji,jj,1,jpdic) = tra(ji,jj,1,jpdic) + ( zfld - zflu ) * rfact2 / e3t_n(ji,jj,1) 189 183 190 184 ! Compute O2 flux … … 192 186 zflu16 = trb(ji,jj,1,jpoxy) * tmask(ji,jj,1) * zkgo2(ji,jj) 193 187 zoflx(ji,jj) = zfld16 - zflu16 194 tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + zoflx(ji,jj) * rfact2 / fse3t(ji,jj,1)188 tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + zoflx(ji,jj) * rfact2 / e3t_n(ji,jj,1) 195 189 END DO 196 190 END DO … … 400 394 401 395 !!====================================================================== 402 END MODULE 396 END MODULE p4zflx
Note: See TracChangeset
for help on using the changeset viewer.