Changeset 1073
- Timestamp:
- 2008-06-05T14:15:34+02:00 (16 years ago)
- Location:
- trunk/NEMO/TOP_SRC/PISCES
- Files:
-
- 2 added
- 17 edited
-
p4zbio.F90 (modified) (2 diffs)
-
p4zche.F90 (modified) (1 diff)
-
p4zflx.F90 (modified) (2 diffs)
-
p4zint.F90 (modified) (1 diff)
-
p4zlim.F90 (modified) (1 diff)
-
p4zlys.F90 (modified) (2 diffs)
-
p4zmeso.F90 (modified) (2 diffs)
-
p4zmicro.F90 (modified) (2 diffs)
-
p4zmort.F90 (modified) (2 diffs)
-
p4zopt.F90 (modified) (1 diff)
-
p4zprod.F90 (modified) (10 diffs)
-
p4zrem.F90 (modified) (2 diffs)
-
p4zsed.F90 (modified) (2 diffs)
-
p4zsink.F90 (modified) (2 diffs)
-
sms_pisces.F90 (added)
-
trcctl_pisces.F90 (added)
-
trcini_pisces.F90 (modified) (1 diff)
-
trclsm_pisces.F90 (modified) (1 diff)
-
trcsms_pisces.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90
r935 r1073 16 16 USE oce_trc ! 17 17 USE trp_trc ! 18 USE sms 18 USE sms_pisces ! 19 19 USE p4zsink ! 20 20 USE p4zopt ! … … 30 30 PRIVATE 31 31 32 PUBLIC p4z_bio ! called in p4zprg.F9032 PUBLIC p4z_bio 33 33 34 34 !! * Shared module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zche.F90
r935 r1073 19 19 USE oce_trc ! 20 20 USE trp_trc ! 21 USE sms 21 USE sms_pisces ! 22 22 23 23 IMPLICIT NONE 24 24 PRIVATE 25 25 26 PUBLIC p4z_che ! called in p4zprg.F9026 PUBLIC p4z_che 27 27 28 28 !! * Shared module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90
r1035 r1073 20 20 USE oce_trc ! 21 21 USE trp_trc 22 USE sms 22 USE sms_pisces 23 23 USE prtctl_trc 24 24 USE p4zche … … 27 27 PRIVATE 28 28 29 PUBLIC p4z_flx ! called in p4zprg.F9029 PUBLIC p4z_flx 30 30 31 31 REAL(wp) :: & ! pre-industrial atmospheric [co2] (ppm) -
trunk/NEMO/TOP_SRC/PISCES/p4zint.F90
r935 r1073 15 15 USE oce_trc ! 16 16 USE trp_trc 17 USE sms 17 USE sms_pisces 18 18 19 19 IMPLICIT NONE 20 20 PRIVATE 21 21 22 PUBLIC p4z_int ! called in p4zprg.F9022 PUBLIC p4z_int 23 23 24 24 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & !: -
trunk/NEMO/TOP_SRC/PISCES/p4zlim.F90
r935 r1073 17 17 USE oce_trc ! 18 18 USE trp_trc ! 19 USE sms 19 USE sms_pisces ! 20 20 21 21 IMPLICIT NONE 22 22 PRIVATE 23 23 24 PUBLIC p4z_lim ! called in p4zprg.F9024 PUBLIC p4z_lim 25 25 26 26 !! * Shared module variables 27 REAL(wp), PUBLIC :: conc0 = 2.e-6_wp , & !: 28 & conc1 = 10.e-6_wp , & !: 29 & conc2 = 2.e-11_wp , & !: 30 & conc2m = 8.E-11_wp , & !: 31 & conc3 = 1.e-10_wp , & !: 32 & conc3m = 4.e-10_wp , & !: 33 & concnnh4 = 1.e-7_wp , & !: 34 & concdnh4 = 5.e-7_wp , & !: 35 & xksi1 = 2.E-6_wp , & !: 36 & xksi2 = 3.33E-6_wp , & !: 37 & xkdoc = 417.E-6_wp , & !: 38 & caco3r = 0.3_wp !: 27 REAL(wp), PUBLIC :: & 28 conc0 = 2.e-6_wp , & !: 29 conc1 = 10.e-6_wp , & !: 30 conc2 = 2.e-11_wp , & !: 31 conc2m = 8.E-11_wp , & !: 32 conc3 = 1.e-10_wp , & !: 33 conc3m = 4.e-10_wp , & !: 34 concnnh4 = 1.e-7_wp , & !: 35 concdnh4 = 5.e-7_wp , & !: 36 xksi1 = 2.E-6_wp , & !: 37 xksi2 = 3.33E-6_wp , & !: 38 xkdoc = 417.E-6_wp , & !: 39 caco3r = 0.3_wp !: 39 40 40 41 -
trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90
r935 r1073 20 20 USE oce_trc ! 21 21 USE trp_trc 22 USE sms 22 USE sms_pisces 23 23 USE prtctl_trc 24 24 … … 29 29 30 30 !! * Shared module variables 31 REAL(wp), PUBLIC :: kdca = 0.327e3_wp , & !: 32 & nca = 1.0_wp !: 31 REAL(wp), PUBLIC :: & 32 kdca = 0.327e3_wp , & !: 33 nca = 1.0_wp !: 33 34 34 35 !! * Module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zmeso.F90
r935 r1073 17 17 USE oce_trc ! 18 18 USE trp_trc ! 19 USE sms 19 USE sms_pisces ! 20 20 USE prtctl_trc 21 21 USE p4zint … … 29 29 !! * Shared module variables 30 30 REAL(wp), PUBLIC :: & 31 xprefc = 1.0_wp, & !:32 xprefp = 0.2_wp, & !:33 xprefz = 1.0_wp, & !:31 xprefc = 1.0_wp , & !: 32 xprefp = 0.2_wp , & !: 33 xprefz = 1.0_wp , & !: 34 34 xprefpoc = 0.0_wp , & !: 35 resrat2 = 0.005_wp, & !:36 mzrat2 = 0.03_wp, & !:35 resrat2 = 0.005_wp , & !: 36 mzrat2 = 0.03_wp , & !: 37 37 grazrat2 = 0.7_wp , & !: 38 xkgraz2 = 20E-6_wp, & !:39 unass2 = 0.3_wp , &!:40 sigma2 = 0.6_wp , &!:41 epsher2 = 0.33_wp , &!:38 xkgraz2 = 20E-6_wp , & !: 39 unass2 = 0.3_wp , & !: 40 sigma2 = 0.6_wp , & !: 41 epsher2 = 0.33_wp , & !: 42 42 grazflux = 5.E3_wp 43 43 -
trunk/NEMO/TOP_SRC/PISCES/p4zmicro.F90
r935 r1073 17 17 USE oce_trc ! 18 18 USE trp_trc ! 19 USE sms 19 USE sms_pisces ! 20 20 USE prtctl_trc 21 21 USE p4zint … … 32 32 xpref2p = 0.5_wp , & !: 33 33 xpref2d = 0.5_wp , & !: 34 resrat = 0.03_wp , & !:35 mzrat = 0.0_wp, & !:36 grazrat = 4.0_wp , & !:37 xkgraz = 20E-6_wp , & !:38 unass = 0.3_wp , & !:39 sigma1 = 0.6_wp , & !:40 epsher = 0.33_wp34 resrat = 0.03_wp , & !: 35 mzrat = 0.0_wp , & !: 36 grazrat = 4.0_wp , & !: 37 xkgraz = 20E-6_wp , & !: 38 unass = 0.3_wp , & !: 39 sigma1 = 0.6_wp , & !: 40 epsher = 0.33_wp 41 41 42 42 -
trunk/NEMO/TOP_SRC/PISCES/p4zmort.F90
r935 r1073 17 17 USE oce_trc ! 18 18 USE trp_trc ! 19 USE sms 19 USE sms_pisces ! 20 20 USE p4zsink 21 21 USE prtctl_trc … … 24 24 PRIVATE 25 25 26 PUBLIC p4z_mort ! called in p4zbio.F9026 PUBLIC p4z_mort 27 27 28 28 29 29 !! * Shared module variables 30 REAL(wp), PUBLIC :: wchl = 0.001_wp , & !: 31 & wchld = 0.02_wp , & !: 32 & mprat = 0.01_wp , & !: 33 & mprat2 = 0.01_wp , & !: 34 & mpratm = 0.01_wp !: 30 REAL(wp), PUBLIC :: & 31 wchl = 0.001_wp , & !: 32 wchld = 0.02_wp , & !: 33 mprat = 0.01_wp , & !: 34 mprat2 = 0.01_wp , & !: 35 mpratm = 0.01_wp !: 35 36 36 37 !! * Module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90
r935 r1073 16 16 USE oce_trc ! 17 17 USE trp_trc 18 USE sms 18 USE sms_pisces 19 19 20 20 IMPLICIT NONE 21 21 PRIVATE 22 22 23 PUBLIC p4z_opt ! called in p4zprg.F9023 PUBLIC p4z_opt 24 24 25 25 !! * Shared module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90
r935 r1073 16 16 USE oce_trc ! 17 17 USE trp_trc ! 18 USE sms ! 19 USE p4zday ! 18 USE sms_pisces ! 20 19 USE prtctl_trc 21 20 USE p4zint … … 29 28 30 29 !! * Shared module variables 31 REAL(wp), PUBLIC :: pislope = 3.0_wp , & !: 32 & pislope2 = 3.0_wp , & !: 33 & excret = 10.e-5_wp , & !: 34 & excret2 = 0.05_wp , & !: 35 & chlcnm = 0.033_wp , & !: 36 & chlcdm = 0.05_wp , & !: 37 & fecnm = 10.E-6_wp , & !: 38 & fecdm = 15.E-6_wp , & !: 39 & grosip = 0.151_wp 30 REAL(wp), PUBLIC :: & 31 pislope = 3.0_wp , & !: 32 pislope2 = 3.0_wp , & !: 33 excret = 10.e-5_wp , & !: 34 excret2 = 0.05_wp , & !: 35 chlcnm = 0.033_wp , & !: 36 chlcdm = 0.05_wp , & !: 37 fecnm = 10.E-6_wp , & !: 38 fecdm = 15.E-6_wp , & !: 39 grosip = 0.151_wp 40 40 41 41 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & … … 43 43 44 44 REAL(wp) :: & 45 tpp = 0. !: Total primary production 45 texcret , & !: 1 - excret 46 texcret2 , & !: 1 - excret2 47 rpis180 , & !: rpi / 180 48 tpp = 0. !: Total primary production 46 49 47 50 !!* Substitution … … 69 72 REAL(wp) :: zprdiachl, zprbiochl, zsilim, ztn, zadap, zadap2 70 73 REAL(wp) :: zlim, zsilfac2, zsiborn, zprod, zetot2, zmax, zproreg, zproreg2 71 REAL(wp) :: zmxltst, zmxlday, zlim1 , zexcret, zexcret274 REAL(wp) :: zmxltst, zmxlday, zlim1 72 75 REAL(wp) :: zpislopen , zpislope2n 76 REAL(wp) :: zrum, zcodel, zargu 73 77 REAL(wp), DIMENSION(jpi,jpj) :: zmixnano , zmixdiat, zstrn 74 78 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zpislopead , zpislopead2 … … 98 102 nspyr = INT( raass / rdt ) 99 103 100 zexcret = 1. - excret101 zexcret2 = 1. - excret2102 104 103 105 ! Computation of the optimal production … … 110 112 # endif 111 113 112 CALL p4z_day ! Computation of the day length 114 ! compute the day length depending on latitude and the day 115 !-------------------------------------------------------- 116 IF(lwp) write(numout,*) 117 IF(lwp) write(numout,*) 'p4zday : - Julian day ', nday_year 118 IF(lwp) write(numout,*) '~~~~~~' 119 120 IF( nleapy == 1 .AND. MOD( nyear, 4 ) == 0 ) THEN 121 zrum = FLOAT( nday_year - 80 ) / 366. 122 ELSE 123 zrum = FLOAT( nday_year - 80 ) / 365. 124 ENDIF 125 zcodel = ASIN( SIN( zrum * rpi * 2. ) * SIN( rpis180 * 23.5 ) ) 126 127 ! day length in hours 128 zstrn(:,:) = 0. 129 DO jj = 1, jpj 130 DO ji = 1, jpi 131 zargu = TAN( zcodel ) * TAN( gphit(ji,jj) * rpis180 ) 132 zargu = MAX( -1., MIN( 1., zargu ) ) 133 zstrn(ji,jj) = MAX( 0.0, 24. - 2. * ACOS( zargu ) / rpis180 / 15. ) 134 END DO 135 END DO 136 113 137 114 138 !CDIR NOVERRCHK … … 210 234 ! Computation of the fractionnal day length 211 235 ! ----------------------------------------- 212 zstrn(:,:) = strn(:,:) 213 214 DO jj = 1, jpj 215 DO ji = 1, jpi 216 236 ! zstrn(:,:) = strn(:,:) 237 238 WHERE( zstrn(:,:) < 1.e0 ) zstrn(:,:) = 24. 239 zstrn(:,:) = 24. / zstrn(:,:) 240 ! DO jj = 1, jpj 241 ! DO ji = 1, jpi 242 ! 217 243 ! Computation of the maximum light intensity 218 244 ! ------------------------------------------ 219 IF( zstrn(ji,jj) < 1.e0 ) zstrn(ji,jj) = 24.220 zstrn(ji,jj) = 24./zstrn(ji,jj)221 END DO222 END DO245 ! IF( zstrn(ji,jj) < 1.e0 ) zstrn(ji,jj) = 24. 246 ! zstrn(ji,jj) = 24. / zstrn(ji,jj) 247 ! END DO 248 ! END DO 223 249 224 250 !CDIR NOVERRCHK … … 304 330 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) - zpronew(ji,jj,jk) - zpronewd(ji,jj,jk) 305 331 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) - zproreg - zproreg2 306 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zprorca(ji,jj,jk) * zexcret307 tra(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) + zprochln(ji,jj,jk) * zexcret308 tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) + zprofen(ji,jj,jk) * zexcret309 tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) + zprorcad(ji,jj,jk) * zexcret2310 tra(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) + zprochld(ji,jj,jk) * zexcret2311 tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) + zprofed(ji,jj,jk) * zexcret2312 tra(ji,jj,jk,jpbsi) = tra(ji,jj,jk,jpbsi) + zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) * zexcret2332 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zprorca(ji,jj,jk) * texcret 333 tra(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) + zprochln(ji,jj,jk) * texcret 334 tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) + zprofen(ji,jj,jk) * texcret 335 tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) + zprorcad(ji,jj,jk) * texcret2 336 tra(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) + zprochld(ji,jj,jk) * texcret2 337 tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) + zprofed(ji,jj,jk) * texcret2 338 tra(ji,jj,jk,jpbsi) = tra(ji,jj,jk,jpbsi) + zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) * texcret2 313 339 tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + & 314 340 & excret2 * zprorcad(ji,jj,jk) + excret * zprorca(ji,jj,jk) … … 316 342 & + ( o2ut + o2nit ) * ( zpronew(ji,jj,jk) + zpronewd(ji,jj,jk) ) 317 343 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) & 318 & - zexcret * zprofen(ji,jj,jk) - zexcret2 * zprofed(ji,jj,jk)344 & - texcret * zprofen(ji,jj,jk) - texcret2 * zprofed(ji,jj,jk) 319 345 tra(ji,jj,jk,jpsil) = tra(ji,jj,jk,jpsil) & 320 & - zexcret2 * zprorcad(ji,jj,jk) * zysopt(ji,jj,jk)346 & - texcret2 * zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) 321 347 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprorca(ji,jj,jk) - zprorcad(ji,jj,jk) 322 348 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) & … … 405 431 ENDIF 406 432 433 rpis180 = rpi / 180. 434 texcret = 1.0 - excret 435 texcret2 = 1.0 - excret2 436 407 437 END SUBROUTINE p4z_prod_init 408 438 -
trunk/NEMO/TOP_SRC/PISCES/p4zrem.F90
r935 r1073 17 17 USE oce_trc ! 18 18 USE trp_trc ! 19 USE sms 19 USE sms_pisces ! 20 20 USE prtctl_trc 21 21 USE p4zint … … 31 31 32 32 !! * Shared module variables 33 REAL(wp), PUBLIC :: xremik = 0.3_wp , & !: 34 & xremip = 0.025_wp , & !: 35 & nitrif = 0.05_wp , & !: 36 & xsirem = 0.015_wp , & !: 37 & xlam1 = 0.005_wp , & !: 38 & oxymin = 1.e-6_wp !: 33 REAL(wp), PUBLIC :: & 34 xremik = 0.3_wp , & !: 35 xremip = 0.025_wp , & !: 36 nitrif = 0.05_wp , & !: 37 xsirem = 0.015_wp , & !: 38 xlam1 = 0.005_wp , & !: 39 oxymin = 1.e-6_wp !: 40 39 41 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & !: 40 42 & denitr !: denitrification array -
trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90
r935 r1073 18 18 USE oce_trc ! 19 19 USE trp_trc 20 USE sms 20 USE sms_pisces 21 21 USE lib_mpp 22 22 USE prtctl_trc … … 34 34 PRIVATE 35 35 36 PUBLIC p4z_sed ! called in p4zprg.F9036 PUBLIC p4z_sed 37 37 38 38 !! * Shared module variables 39 LOGICAL, PUBLIC :: bdustfer = .FALSE. , & !: 40 & briver = .FALSE. , & !: 41 & bndepo = .FALSE. , & !: 42 & bsedinput = .FALSE. !: 43 44 REAL(wp), PUBLIC :: sedfeinput = 1.E-9_wp , & !: 45 & dustsolub = 0.014_wp !: 39 LOGICAL, PUBLIC :: & 40 bdustfer = .FALSE. , & !: 41 briver = .FALSE. , & !: 42 bndepo = .FALSE. , & !: 43 bsedinput = .FALSE. !: 44 45 REAL(wp), PUBLIC :: & 46 sedfeinput = 1.E-9_wp , & !: 47 dustsolub = 0.014_wp !: 46 48 47 49 !! * Module variables -
trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90
r935 r1073 13 13 USE oce_trc ! 14 14 USE trp_trc 15 USE sms 16 USE p4zsink2 ! 15 USE sms_pisces 17 16 USE prtctl_trc 18 17 … … 615 614 #endif 616 615 616 SUBROUTINE p4z_sink2( pwsink, psinkflx, jp_tra ) 617 !!--------------------------------------------------------------------- 618 !! *** ROUTINE p4z_sink2 *** 619 !! 620 !! ** Purpose : Compute the sedimentation terms for the various sinking 621 !! particles. The scheme used to compute the trends is based 622 !! on MUSCL. 623 !! 624 !! ** Method : - this ROUTINE compute not exactly the advection but the 625 !! transport term, i.e. div(u*tra). 626 !!--------------------------------------------------------------------- 627 INTEGER , INTENT(in ) :: jp_tra ! tracer index index 628 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj,jpk) :: pwsink ! sinking speed 629 REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk) :: psinkflx ! sinking fluxe 630 !! 631 INTEGER :: ji, jj, jk, jn 632 REAL(wp) :: zigma,zew,zstep,zign, zflx 633 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztraz, zakz 634 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zwsink2 635 !!--------------------------------------------------------------------- 636 637 zstep = rfact2 / 2. 638 639 ztraz(:,:,:) = 0.e0 640 zakz (:,:,:) = 0.e0 641 642 DO jk = 1, jpkm1 643 # if defined key_off_degrad 644 zwsink2(:,:,jk+1) = -pwsink(:,:,jk) / rjjss * tmask(:,:,jk+1) * facvol(:,:,jk) 645 # else 646 zwsink2(:,:,jk+1) = -pwsink(:,:,jk) / rjjss * tmask(:,:,jk+1) 647 # endif 648 END DO 649 zwsink2(:,:,1) = 0.e0 650 651 652 ! Vertical advective flux 653 DO jn = 1, 2 654 ! first guess of the slopes interior values 655 DO jk = 2, jpkm1 656 ztraz(:,:,jk) = ( trn(:,:,jk-1,jp_tra) - trn(:,:,jk,jp_tra) ) * tmask(:,:,jk) 657 END DO 658 ztraz(:,:,1 ) = 0.0 659 ztraz(:,:,jpk) = 0.0 660 661 ! slopes 662 DO jk = 2, jpkm1 663 DO jj = 1,jpj 664 DO ji = 1, jpi 665 zign = 0.25 + SIGN( 0.25, ztraz(ji,jj,jk) * ztraz(ji,jj,jk+1) ) 666 zakz(ji,jj,jk) = ( ztraz(ji,jj,jk) + ztraz(ji,jj,jk+1) ) * zign 667 END DO 668 END DO 669 END DO 670 671 ! Slopes limitation 672 DO jk = 2, jpkm1 673 DO jj = 1, jpj 674 DO ji = 1, jpi 675 zakz(ji,jj,jk) = SIGN( 1., zakz(ji,jj,jk) ) * & 676 & MIN( ABS( zakz(ji,jj,jk) ), 2. * ABS(ztraz(ji,jj,jk+1)), 2. * ABS(ztraz(ji,jj,jk) ) ) 677 END DO 678 END DO 679 END DO 680 681 ! vertical advective flux 682 DO jk = 1, jpkm1 683 DO jj = 1, jpj 684 DO ji = 1, jpi 685 zigma = zwsink2(ji,jj,jk+1) * zstep / fse3w(ji,jj,jk+1) 686 zew = zwsink2(ji,jj,jk+1) 687 psinkflx(ji,jj,jk+1) = -zew * ( trn(ji,jj,jk,jp_tra) - 0.5 * ( 1 + zigma ) * zakz(ji,jj,jk) ) * zstep 688 END DO 689 END DO 690 END DO 691 ! 692 ! Boundary conditions 693 psinkflx(:,:,1 ) = 0.e0 694 psinkflx(:,:,jpk) = 0.e0 695 696 DO jk=1,jpkm1 697 DO jj = 1,jpj 698 DO ji = 1, jpi 699 zflx = ( psinkflx(ji,jj,jk) - psinkflx(ji,jj,jk+1) ) / fse3t(ji,jj,jk) 700 trn(ji,jj,jk,jp_tra) = trn(ji,jj,jk,jp_tra) + zflx 701 END DO 702 END DO 703 END DO 704 705 ENDDO 706 707 DO jk=1,jpkm1 708 DO jj = 1,jpj 709 DO ji = 1, jpi 710 zflx = ( psinkflx(ji,jj,jk) - psinkflx(ji,jj,jk+1) ) / fse3t(ji,jj,jk) 711 trb(ji,jj,jk,jp_tra) = trb(ji,jj,jk,jp_tra) + 2. * zflx 712 END DO 713 END DO 714 END DO 715 716 trn(:,:,:,jp_tra) = trb(:,:,:,jp_tra) 717 psinkflx(:,:,:) = 2. * psinkflx(:,:,:) 718 719 ! 720 END SUBROUTINE p4z_sink2 721 617 722 #else 618 723 !!====================================================================== -
trunk/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
r1007 r1073 17 17 !!---------------------------------------------------------------------- 18 18 USE par_trc ! TOP parameters 19 USE sms 19 USE sms_pisces ! Source Minus Sink variables 20 20 USE trc 21 21 USE oce_trc ! ocean variables -
trunk/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90
r1007 r1073 18 18 USE par_trc ! TOP parameters 19 19 USE trc ! TOP variables 20 USE sms 20 USE sms_pisces ! sms trends 21 21 22 22 -
trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90
r935 r1073 15 15 USE oce_trc ! 16 16 USE trp_trc 17 USE sms 17 USE sms_pisces 18 18 USE lbclnk 19 19 USE lib_mpp … … 53 53 !! 54 54 INTEGER :: jnt, jn 55 INTEGER :: iyy, imm, idd56 57 58 55 !!--------------------------------------------------------------------- 59 56 60 61 iyy = ndastp/10000 62 imm = (ndastp - iyy*10000)/100 63 idd = (ndastp - iyy*10000 - imm*100) 64 65 IF( ndayflxtr /= idd) THEN ! New days 57 IF( ndayflxtr /= nday ) THEN ! New days 66 58 ! 67 ndayflxtr = idd59 ndayflxtr = nday 68 60 69 61 CALL p4z_che ! computation of chemical constants
Note: See TracChangeset
for help on using the changeset viewer.
