Changeset 1255 for trunk/NEMO/TOP_SRC/PISCES
- Timestamp:
- 2009-01-13T11:20:17+01:00 (15 years ago)
- Location:
- trunk/NEMO/TOP_SRC/PISCES
- Files:
-
- 5 edited
-
p4zprod.F90 (modified) (2 diffs)
-
p4zrem.F90 (modified) (2 diffs)
-
p4zsed.F90 (modified) (3 diffs)
-
trclsm_pisces.F90 (modified) (6 diffs)
-
trcsms_pisces.F90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90
r1180 r1255 73 73 REAL(wp) :: zmxltst, zmxlday, zlim1 74 74 REAL(wp) :: zpislopen , zpislope2n 75 REAL(wp) :: zrum, zcodel, zargu 75 REAL(wp) :: zrum, zcodel, zargu, zvol 76 76 REAL(wp), DIMENSION(jpi,jpj) :: zmixnano , zmixdiat, zstrn 77 77 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zpislopead , zpislopead2 … … 355 355 DO jj = 1, jpj 356 356 DO ji = 1, jpi 357 tpp = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) ) &357 zvol = cvol(ji,jj,jk) 358 358 #if defined key_off_degrad 359 & * facvol(ji,jj,jk) &359 zvol = zvol * facvol(ji,jj,jk) 360 360 #endif 361 & * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) * tmask_i(ji,jj)361 tpp = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) ) * zvol 362 362 END DO 363 363 END DO 364 364 END DO 365 366 IF( lk_mpp ) CALL mpp_sum( tpp ) 365 367 366 368 IF( MOD( kt, nspyr ) == 0 ) THEN -
trunk/NEMO/TOP_SRC/PISCES/p4zrem.F90
r1152 r1255 65 65 REAL(wp) :: zkeq , zfeequi, zsiremin 66 66 REAL(wp) :: zsatur, zsatur2, znusil 67 REAL(wp) :: zbactfer, zorem, zorem2, zofer, zofer2 68 REAL(wp) :: zosil, zdenom, zdenom1, zdenom2, zscave, zaggdfe 67 REAL(wp) :: zbactfer, zorem, zorem2, zofer 68 REAL(wp) :: zosil, zdenom1, zscave, zaggdfe 69 #if ! defined key_kriest 70 REAL(wp) :: zofer2, zdenom, zdenom2 71 #endif 69 72 REAL(wp) :: zlamfac, zstep, zonitr 70 73 REAL(wp), DIMENSION(jpi,jpj) :: ztempbac … … 302 305 303 306 #if defined key_kriest 304 zdenom1 = trn(ji,jj,jk,jppoc) / &307 zdenom1 = trn(ji,jj,jk,jppoc) / & 305 308 & ( trn(ji,jj,jk,jppoc) + trn(ji,jj,jk,jpdsi) + trn(ji,jj,jk,jpcal) + rtrn ) 306 309 #else -
trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90
r1180 r1255 269 269 DO jj = 1,jpj 270 270 DO ji = 1,jpi 271 zdenitot = zdenitot + denitr(ji,jj,jk) * rdenit * e1t(ji,jj) * e2t(ji,jj) & 272 & *fse3t(ji,jj,jk) * tmask(ji,jj,jk) * tmask_i(ji,jj) * xnegtr(ji,jj,jk) 271 zdenitot = zdenitot + denitr(ji,jj,jk) * rdenit * cvol(ji,jj,jk) * xnegtr(ji,jj,jk) 273 272 END DO 274 273 END DO … … 302 301 DO jj = 1, jpj 303 302 DO ji = 1, jpi 304 znitrpottot = znitrpottot + znitrpot(ji,jj,jk) * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) & 305 & * tmask(ji,jj,jk) * tmask_i(ji,jj) 303 znitrpottot = znitrpottot + znitrpot(ji,jj,jk) * cvol(ji,jj,jk) 306 304 END DO 307 305 END DO … … 596 594 DO jj = 2 , jpjm1 597 595 DO ji = 2, jpim1 598 zcoef = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,1) * tmask(ji,jj,1) * tmask_i(ji,jj) * raass596 zcoef = cvol(ji,jj,1) * raass 599 597 rivpo4input = rivpo4input + rivinp(ji,jj) * zcoef 600 598 rivalkinput = rivalkinput + cotdep(ji,jj) * zcoef -
trunk/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90
r1146 r1255 46 46 !!---------------------------------------------------------------------- 47 47 CHARACTER (len=32) :: clname 48 INTEGER :: jn49 48 !! 50 49 #if defined key_trc_diaadd 50 INTEGER :: jl, jn 51 51 ! definition of additional diagnostic as a structure 52 52 TYPE DIAG … … 127 127 nwritedia = 10 ! default values 128 128 129 DO jn = jp_pcs0_2d, jp_pcs1_2d 129 DO jl = 1, jp_pisces_2d 130 jn = jp_pcs0_2d + jl - 1 130 131 WRITE(ctrc2d(jn),'("2D_",I1)') jn ! short name 131 132 WRITE(ctrc2l(jn),'("2D DIAGNOSTIC NUMBER ",I2)') jn ! long name … … 133 134 END DO 134 135 ! ! 3D output arrays 135 DO jn = jp_pcs0_3d, jp_pcs1_3d 136 DO jl = 1, jp_pisces_3d 137 jn = jp_pcs0_3d + jl - 1 136 138 WRITE(ctrc3d(jn),'("3D_",I1)') jn ! short name 137 139 WRITE(ctrc3l(jn),'("3D DIAGNOSTIC NUMBER ",I2)') jn ! long name … … 142 144 READ ( numnat, nampisdia ) 143 145 144 DO jn = jp_pcs0_2d, jp_pcs1_2d 145 ctrc2d(jn) = pisdia2d(jn)%snamedia 146 ctrc2l(jn) = pisdia2d(jn)%lnamedia 147 ctrc2u(jn) = pisdia2d(jn)%unitdia 146 DO jl = 1, jp_pisces_2d 147 jn = jp_pcs0_2d + jl - 1 148 ctrc2d(jn) = pisdia2d(jl)%snamedia 149 ctrc2l(jn) = pisdia2d(jl)%lnamedia 150 ctrc2u(jn) = pisdia2d(jl)%unitdia 148 151 END DO 149 152 150 DO jn = jp_pcs0_3d, jp_pcs1_3d 151 ctrc3d(jn) = pisdia3d(jn)%snamedia 152 ctrc3l(jn) = pisdia3d(jn)%lnamedia 153 ctrc3u(jn) = pisdia3d(jn)%unitdia 153 DO jl = 1, jp_pisces_3d 154 jn = jp_pcs0_3d + jl - 1 155 ctrc3d(jn) = pisdia3d(jl)%snamedia 156 ctrc3l(jn) = pisdia3d(jl)%lnamedia 157 ctrc3u(jn) = pisdia3d(jl)%unitdia 154 158 END DO 155 159 … … 158 162 WRITE(numout,*) ' Namelist : natadd' 159 163 WRITE(numout,*) ' frequency of outputs for additional arrays nwritedia = ', nwritedia 160 DO jn = jp_pcs0_3d, jp_pcs1_3d 164 DO jl = 1, jp_pisces_3d 165 jn = jp_pcs0_3d + jl - 1 161 166 WRITE(numout,*) ' 3d output field No : ',jn 162 167 WRITE(numout,*) ' short name : ', TRIM(ctrc3d(jn)) … … 166 171 END DO 167 172 168 DO jn = jp_pcs0_2d, jp_pcs1_2d 173 DO jl = 1, jp_pisces_2d 174 jn = jp_pcs0_2d + jl - 1 169 175 WRITE(numout,*) ' 2d output field No : ',jn 170 176 WRITE(numout,*) ' short name : ', TRIM(ctrc2d(jn)) -
trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90
r1185 r1255 25 25 USE p4zlys ! 26 26 USE p4zflx ! 27 28 USE trdmld_trc_oce 29 USE trdmld_trc 27 30 28 31 USE sedmodel … … 55 58 !! 56 59 INTEGER :: jnt, jn 60 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztrpis ! used for pisces sms trends 57 61 !!--------------------------------------------------------------------- 58 62 … … 62 66 63 67 CALL p4z_che ! computation of chemical constants 64 65 68 CALL p4z_int ! computation of various rates for biogeochemistry 66 69 ! … … 71 74 ! 72 75 CALL p4z_bio (kt, jnt) ! Compute soft tissue production (POC) 73 74 75 76 CALL p4z_sed (kt, jnt) ! compute soft tissue remineralisation 76 77 77 ! 78 78 trb(:,:,:,:) = trn(:,:,:,:) … … 81 81 82 82 CALL p4z_lys( kt ) ! Compute CaCO3 saturation 83 84 83 CALL p4z_flx( kt ) ! Compute surface fluxes 85 84 … … 90 89 CALL lbc_lnk( tra(:,:,:,jn), 'T', 1. ) 91 90 END DO 91 92 IF( l_trdtrc ) THEN 93 DO jn = jp_pcs0, jp_pcs1 94 ztrpis(:,:,:) = tra(:,:,:,jn) 95 CALL trd_mod_trc( ztrpis, jn, jptrc_trd_sms, kt ) ! save trends 96 END DO 97 END IF 92 98 93 99 #if defined key_sed
Note: See TracChangeset
for help on using the changeset viewer.
