- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
TOP/PISCES/P4Z/p4zpoc.F90 (modified) (16 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zpoc.F90
r11536 r11949 44 44 CONTAINS 45 45 46 SUBROUTINE p4z_poc( kt, knt )46 SUBROUTINE p4z_poc( kt, knt, Kbb, Kmm, Krhs ) 47 47 !!--------------------------------------------------------------------- 48 48 !! *** ROUTINE p4z_poc *** … … 52 52 !! ** Method : - ??? 53 53 !!--------------------------------------------------------------------- 54 INTEGER, INTENT(in) :: kt, knt ! ocean time step and ??? 54 INTEGER, INTENT(in) :: kt, knt ! ocean time step and ??? 55 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 55 56 ! 56 57 INTEGER :: ji, jj, jk, jn … … 113 114 ! ------------------------------------------------------------ 114 115 ! 115 IF( gdept _n(ji,jj,jk) > zdep ) THEN116 IF( gdept(ji,jj,jk,Kmm) > zdep ) THEN 116 117 alphat = 0. 117 118 remint = 0. 118 119 ! 119 zsizek1 = e3t _n(ji,jj,jk-1) / 2. / (wsbio4(ji,jj,jk-1) + rtrn) * tgfunc(ji,jj,jk-1)120 zsizek = e3t _n(ji,jj,jk) / 2. / (wsbio4(ji,jj,jk) + rtrn) * tgfunc(ji,jj,jk)120 zsizek1 = e3t(ji,jj,jk-1,Kmm) / 2. / (wsbio4(ji,jj,jk-1) + rtrn) * tgfunc(ji,jj,jk-1) 121 zsizek = e3t(ji,jj,jk,Kmm) / 2. / (wsbio4(ji,jj,jk) + rtrn) * tgfunc(ji,jj,jk) 121 122 ! 122 IF ( gdept _n(ji,jj,jk-1) <= zdep ) THEN123 IF ( gdept(ji,jj,jk-1,Kmm) <= zdep ) THEN 123 124 ! 124 125 ! The first level just below the mixed layer needs a … … 131 132 ! POC concentration is computed using the lagrangian 132 133 ! framework. It is only used for the lability param 133 zpoc = tr b(ji,jj,jk-1,jpgoc) + consgoc(ji,jj,jk) * rday / rfact2 &134 & * e3t _n(ji,jj,jk) / 2. / (wsbio4(ji,jj,jk) + rtrn)134 zpoc = tr(ji,jj,jk-1,jpgoc,Kbb) + consgoc(ji,jj,jk) * rday / rfact2 & 135 & * e3t(ji,jj,jk,Kmm) / 2. / (wsbio4(ji,jj,jk) + rtrn) 135 136 zpoc = MAX(0., zpoc) 136 137 ! … … 158 159 ! --------------------------------------------------- 159 160 ! 160 zpoc = tr b(ji,jj,jk-1,jpgoc) + consgoc(ji,jj,jk-1) * rday / rfact2 &161 & * e3t _n(ji,jj,jk-1) / 2. / (wsbio4(ji,jj,jk-1) + rtrn) + consgoc(ji,jj,jk) &162 & * rday / rfact2 * e3t _n(ji,jj,jk) / 2. / (wsbio4(ji,jj,jk) + rtrn)161 zpoc = tr(ji,jj,jk-1,jpgoc,Kbb) + consgoc(ji,jj,jk-1) * rday / rfact2 & 162 & * e3t(ji,jj,jk-1,Kmm) / 2. / (wsbio4(ji,jj,jk-1) + rtrn) + consgoc(ji,jj,jk) & 163 & * rday / rfact2 * e3t(ji,jj,jk,Kmm) / 2. / (wsbio4(ji,jj,jk) + rtrn) 163 164 zpoc = max(0., zpoc) 164 165 ! … … 198 199 ! -------------------------------------------------------- 199 200 zremig = zremigoc(ji,jj,jk) * xstep * tgfunc(ji,jj,jk) 200 zorem2 = zremig * tr b(ji,jj,jk,jpgoc)201 zorem2 = zremig * tr(ji,jj,jk,jpgoc,Kbb) 201 202 orem(ji,jj,jk) = zorem2 202 zorem3(ji,jj,jk) = zremig * solgoc * tr b(ji,jj,jk,jpgoc)203 zofer2 = zremig * tr b(ji,jj,jk,jpbfe)204 zofer3 = zremig * solgoc * tr b(ji,jj,jk,jpbfe)203 zorem3(ji,jj,jk) = zremig * solgoc * tr(ji,jj,jk,jpgoc,Kbb) 204 zofer2 = zremig * tr(ji,jj,jk,jpbfe,Kbb) 205 zofer3 = zremig * solgoc * tr(ji,jj,jk,jpbfe,Kbb) 205 206 206 207 ! ------------------------------------- 207 tr a(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zorem3(ji,jj,jk)208 tr a(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) - zorem2 - zorem3(ji,jj,jk)209 tr a(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zofer3210 tr a(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zofer2 - zofer3211 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zorem2212 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer2208 tr(ji,jj,jk,jppoc,Krhs) = tr(ji,jj,jk,jppoc,Krhs) + zorem3(ji,jj,jk) 209 tr(ji,jj,jk,jpgoc,Krhs) = tr(ji,jj,jk,jpgoc,Krhs) - zorem2 - zorem3(ji,jj,jk) 210 tr(ji,jj,jk,jpsfe,Krhs) = tr(ji,jj,jk,jpsfe,Krhs) + zofer3 211 tr(ji,jj,jk,jpbfe,Krhs) = tr(ji,jj,jk,jpbfe,Krhs) - zofer2 - zofer3 212 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zorem2 213 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zofer2 213 214 zfolimi(ji,jj,jk) = zofer2 214 215 END DO … … 222 223 ! -------------------------------------------------------- 223 224 zremig = zremigoc(ji,jj,jk) * xstep * tgfunc(ji,jj,jk) 224 zopoc2 = zremig * tr b(ji,jj,jk,jpgoc)225 zopoc2 = zremig * tr(ji,jj,jk,jpgoc,Kbb) 225 226 orem(ji,jj,jk) = zopoc2 226 zorem3(ji,jj,jk) = zremig * solgoc * tr b(ji,jj,jk,jpgoc)227 zopon2 = xremipn / xremipc * zremig * tr b(ji,jj,jk,jpgon)228 zopop2 = xremipp / xremipc * zremig * tr b(ji,jj,jk,jpgop)229 zofer2 = xremipn / xremipc * zremig * tr b(ji,jj,jk,jpbfe)227 zorem3(ji,jj,jk) = zremig * solgoc * tr(ji,jj,jk,jpgoc,Kbb) 228 zopon2 = xremipn / xremipc * zremig * tr(ji,jj,jk,jpgon,Kbb) 229 zopop2 = xremipp / xremipc * zremig * tr(ji,jj,jk,jpgop,Kbb) 230 zofer2 = xremipn / xremipc * zremig * tr(ji,jj,jk,jpbfe,Kbb) 230 231 231 232 ! ------------------------------------- 232 tr a(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zorem3(ji,jj,jk)233 tr a(ji,jj,jk,jppon) = tra(ji,jj,jk,jppon) + solgoc * zopon2234 tr a(ji,jj,jk,jppop) = tra(ji,jj,jk,jppop) + solgoc * zopop2235 tr a(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + solgoc * zofer2236 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zopoc2237 tr a(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + zopon2238 tr a(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zopop2239 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer2240 tr a(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) - zopoc2 - zorem3(ji,jj,jk)241 tr a(ji,jj,jk,jpgon) = tra(ji,jj,jk,jpgon) - zopon2 * (1. + solgoc)242 tr a(ji,jj,jk,jpgop) = tra(ji,jj,jk,jpgop) - zopop2 * (1. + solgoc)243 tr a(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zofer2 * (1. + solgoc)233 tr(ji,jj,jk,jppoc,Krhs) = tr(ji,jj,jk,jppoc,Krhs) + zorem3(ji,jj,jk) 234 tr(ji,jj,jk,jppon,Krhs) = tr(ji,jj,jk,jppon,Krhs) + solgoc * zopon2 235 tr(ji,jj,jk,jppop,Krhs) = tr(ji,jj,jk,jppop,Krhs) + solgoc * zopop2 236 tr(ji,jj,jk,jpsfe,Krhs) = tr(ji,jj,jk,jpsfe,Krhs) + solgoc * zofer2 237 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zopoc2 238 tr(ji,jj,jk,jpdon,Krhs) = tr(ji,jj,jk,jpdon,Krhs) + zopon2 239 tr(ji,jj,jk,jpdop,Krhs) = tr(ji,jj,jk,jpdop,Krhs) + zopop2 240 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zofer2 241 tr(ji,jj,jk,jpgoc,Krhs) = tr(ji,jj,jk,jpgoc,Krhs) - zopoc2 - zorem3(ji,jj,jk) 242 tr(ji,jj,jk,jpgon,Krhs) = tr(ji,jj,jk,jpgon,Krhs) - zopon2 * (1. + solgoc) 243 tr(ji,jj,jk,jpgop,Krhs) = tr(ji,jj,jk,jpgop,Krhs) - zopop2 * (1. + solgoc) 244 tr(ji,jj,jk,jpbfe,Krhs) = tr(ji,jj,jk,jpbfe,Krhs) - zofer2 * (1. + solgoc) 244 245 zfolimi(ji,jj,jk) = zofer2 245 246 END DO … … 251 252 WRITE(charout, FMT="('poc1')") 252 253 CALL prt_ctl_trc_info(charout) 253 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)254 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 254 255 ENDIF 255 256 … … 272 273 DO ji = 1, jpi 273 274 zdep = hmld(ji,jj) 274 IF (tmask(ji,jj,jk) == 1. .AND. gdept _n(ji,jj,jk) <= zdep ) THEN275 totprod(ji,jj) = totprod(ji,jj) + prodpoc(ji,jj,jk) * e3t _n(ji,jj,jk) * rday/ rfact2275 IF (tmask(ji,jj,jk) == 1. .AND. gdept(ji,jj,jk,Kmm) <= zdep ) THEN 276 totprod(ji,jj) = totprod(ji,jj) + prodpoc(ji,jj,jk) * e3t(ji,jj,jk,Kmm) * rday/ rfact2 276 277 ! The temperature effect is included here 277 totthick(ji,jj) = totthick(ji,jj) + e3t _n(ji,jj,jk)* tgfunc(ji,jj,jk)278 totcons(ji,jj) = totcons(ji,jj) - conspoc(ji,jj,jk) * e3t _n(ji,jj,jk) * rday/ rfact2 &279 & / ( tr b(ji,jj,jk,jppoc) + rtrn )278 totthick(ji,jj) = totthick(ji,jj) + e3t(ji,jj,jk,Kmm)* tgfunc(ji,jj,jk) 279 totcons(ji,jj) = totcons(ji,jj) - conspoc(ji,jj,jk) * e3t(ji,jj,jk,Kmm) * rday/ rfact2 & 280 & / ( tr(ji,jj,jk,jppoc,Kbb) + rtrn ) 280 281 ENDIF 281 282 END DO … … 294 295 alphat = 0.0 295 296 remint = 0.0 296 IF( gdept _n(ji,jj,jk) <= zdep ) THEN297 IF( gdept(ji,jj,jk,Kmm) <= zdep ) THEN 297 298 DO jn = 1, jcpoc 298 299 ! For each lability class, the system is supposed to be … … 331 332 IF (tmask(ji,jj,jk) == 1.) THEN 332 333 zdep = hmld(ji,jj) 333 IF( gdept _n(ji,jj,jk) > zdep ) THEN334 IF( gdept(ji,jj,jk,Kmm) > zdep ) THEN 334 335 alphat = 0. 335 336 remint = 0. 336 337 ! 337 338 ! the scale factors are corrected with temperature 338 zsizek1 = e3t _n(ji,jj,jk-1) / 2. / (wsbio3(ji,jj,jk-1) + rtrn) * tgfunc(ji,jj,jk-1)339 zsizek = e3t _n(ji,jj,jk) / 2. / (wsbio3(ji,jj,jk) + rtrn) * tgfunc(ji,jj,jk)339 zsizek1 = e3t(ji,jj,jk-1,Kmm) / 2. / (wsbio3(ji,jj,jk-1) + rtrn) * tgfunc(ji,jj,jk-1) 340 zsizek = e3t(ji,jj,jk,Kmm) / 2. / (wsbio3(ji,jj,jk) + rtrn) * tgfunc(ji,jj,jk) 340 341 ! 341 342 ! Special treatment of the level just below the MXL … … 343 344 ! --------------------------------------------------- 344 345 ! 345 IF ( gdept _n(ji,jj,jk-1) <= zdep ) THEN346 IF ( gdept(ji,jj,jk-1,Kmm) <= zdep ) THEN 346 347 ! 347 348 ! Computation of the POC concentration using the 348 349 ! lagrangian algorithm 349 zpoc = tr b(ji,jj,jk-1,jppoc) + conspoc(ji,jj,jk) * rday / rfact2 &350 & * e3t _n(ji,jj,jk) / 2. / (wsbio3(ji,jj,jk) + rtrn)350 zpoc = tr(ji,jj,jk-1,jppoc,Kbb) + conspoc(ji,jj,jk) * rday / rfact2 & 351 & * e3t(ji,jj,jk,Kmm) / 2. / (wsbio3(ji,jj,jk) + rtrn) 351 352 zpoc = max(0., zpoc) 352 353 ! … … 368 369 ! -------------------------------------------------------- 369 370 ! 370 zpoc = tr b(ji,jj,jk-1,jppoc) + conspoc(ji,jj,jk-1) * rday / rfact2 &371 & * e3t _n(ji,jj,jk-1) / 2. / (wsbio3(ji,jj,jk-1) + rtrn) + conspoc(ji,jj,jk) &372 & * rday / rfact2 * e3t _n(ji,jj,jk) / 2. / (wsbio3(ji,jj,jk) + rtrn)371 zpoc = tr(ji,jj,jk-1,jppoc,Kbb) + conspoc(ji,jj,jk-1) * rday / rfact2 & 372 & * e3t(ji,jj,jk-1,Kmm) / 2. / (wsbio3(ji,jj,jk-1) + rtrn) + conspoc(ji,jj,jk) & 373 & * rday / rfact2 * e3t(ji,jj,jk,Kmm) / 2. / (wsbio3(ji,jj,jk) + rtrn) 373 374 zpoc = max(0., zpoc) 374 375 ! … … 411 412 ! -------------------------------------------------------- 412 413 zremip = zremipoc(ji,jj,jk) * xstep * tgfunc(ji,jj,jk) 413 zorem = zremip * tr b(ji,jj,jk,jppoc)414 zofer = zremip * tr b(ji,jj,jk,jpsfe)415 416 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zorem414 zorem = zremip * tr(ji,jj,jk,jppoc,Kbb) 415 zofer = zremip * tr(ji,jj,jk,jpsfe,Kbb) 416 417 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zorem 417 418 orem(ji,jj,jk) = orem(ji,jj,jk) + zorem 418 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer419 tr a(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zorem420 tr a(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zofer419 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zofer 420 tr(ji,jj,jk,jppoc,Krhs) = tr(ji,jj,jk,jppoc,Krhs) - zorem 421 tr(ji,jj,jk,jpsfe,Krhs) = tr(ji,jj,jk,jpsfe,Krhs) - zofer 421 422 zfolimi(ji,jj,jk) = zfolimi(ji,jj,jk) + zofer 422 423 ENDIF … … 431 432 ! -------------------------------------------------------- 432 433 zremip = zremipoc(ji,jj,jk) * xstep * tgfunc(ji,jj,jk) 433 zopoc = zremip * tr b(ji,jj,jk,jppoc)434 zopoc = zremip * tr(ji,jj,jk,jppoc,Kbb) 434 435 orem(ji,jj,jk) = orem(ji,jj,jk) + zopoc 435 zopon = xremipn / xremipc * zremip * tr b(ji,jj,jk,jppon)436 zopop = xremipp / xremipc * zremip * tr b(ji,jj,jk,jppop)437 zofer = xremipn / xremipc * zremip * tr b(ji,jj,jk,jpsfe)438 439 tr a(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zopoc440 tr a(ji,jj,jk,jppon) = tra(ji,jj,jk,jppon) - zopon441 tr a(ji,jj,jk,jppop) = tra(ji,jj,jk,jppop) - zopop442 tr a(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zofer443 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zopoc444 tr a(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + zopon445 tr a(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zopop446 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer436 zopon = xremipn / xremipc * zremip * tr(ji,jj,jk,jppon,Kbb) 437 zopop = xremipp / xremipc * zremip * tr(ji,jj,jk,jppop,Kbb) 438 zofer = xremipn / xremipc * zremip * tr(ji,jj,jk,jpsfe,Kbb) 439 440 tr(ji,jj,jk,jppoc,Krhs) = tr(ji,jj,jk,jppoc,Krhs) - zopoc 441 tr(ji,jj,jk,jppon,Krhs) = tr(ji,jj,jk,jppon,Krhs) - zopon 442 tr(ji,jj,jk,jppop,Krhs) = tr(ji,jj,jk,jppop,Krhs) - zopop 443 tr(ji,jj,jk,jpsfe,Krhs) = tr(ji,jj,jk,jpsfe,Krhs) - zofer 444 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zopoc 445 tr(ji,jj,jk,jpdon,Krhs) = tr(ji,jj,jk,jpdon,Krhs) + zopon 446 tr(ji,jj,jk,jpdop,Krhs) = tr(ji,jj,jk,jpdop,Krhs) + zopop 447 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zofer 447 448 zfolimi(ji,jj,jk) = zfolimi(ji,jj,jk) + zofer 448 449 END DO … … 463 464 WRITE(charout, FMT="('poc2')") 464 465 CALL prt_ctl_trc_info(charout) 465 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)466 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 466 467 ENDIF 467 468 !
Note: See TracChangeset
for help on using the changeset viewer.
