- Timestamp:
- 2017-12-19T09:26:25+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zlim.F90
r7646 r9124 154 154 !!--------------------------------------------------------------------- 155 155 ! 156 IF( nn_timing == 1 )CALL timing_start('p5z_lim')156 IF( ln_timing ) CALL timing_start('p5z_lim') 157 157 ! 158 158 zratchl = 6.0 … … 448 448 ENDIF 449 449 ! 450 IF( nn_timing == 1) CALL timing_stop('p5z_lim')450 IF( ln_timing ) CALL timing_stop('p5z_lim') 451 451 ! 452 452 END SUBROUTINE p5z_lim 453 453 454 454 455 SUBROUTINE p5z_lim_init 455 456 456 !!---------------------------------------------------------------------- 457 457 !! *** ROUTINE p5z_lim_init *** … … 465 465 !! 466 466 !!---------------------------------------------------------------------- 467 467 INTEGER :: ios ! Local integer output status for namelist read 468 !! 468 469 NAMELIST/namp5zlim/ concnno3, concpno3, concdno3, concnnh4, concpnh4, concdnh4, & 469 470 & concnfer, concpfer, concdfer, concbfe, concnpo4, concppo4, & … … 471 472 & xsizephy, xsizern, xsizerp, xsizerd, xksi1, xksi2, xkdoc, & 472 473 & caco3r, oxymin 473 474 ! 474 475 NAMELIST/namp5zquota/ qnnmin, qnnmax, qpnmin, qpnmax, qnpmin, qnpmax, qppmin, & 475 476 & qppmax, qndmin, qndmax, qpdmin, qpdmax, qfnmax, qfpmax, qfdmax, & 476 477 & qfnopt, qfpopt, qfdopt 477 478 479 INTEGER :: ios ! Local integer output status for namelist read 480 478 !!---------------------------------------------------------------------- 479 ! 481 480 REWIND( numnatp_ref ) ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 482 481 READ ( numnatp_ref, namp5zlim, IOSTAT = ios, ERR = 901) 483 482 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in reference namelist', lwp ) 484 483 ! 485 484 REWIND( numnatp_cfg ) ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters 486 485 READ ( numnatp_cfg, namp5zlim, IOSTAT = ios, ERR = 902 ) 487 902 IF( ios /=0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp )486 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp ) 488 487 IF(lwm) WRITE ( numonp, namp5zlim ) 489 488 ! 490 489 IF(lwp) THEN ! control print 491 490 WRITE(numout,*) ' ' … … 523 522 READ ( numnatp_ref, namp5zquota, IOSTAT = ios, ERR = 903) 524 523 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisquota in reference namelist', lwp ) 525 524 ! 526 525 REWIND( numnatp_cfg ) ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters 527 526 READ ( numnatp_cfg, namp5zquota, IOSTAT = ios, ERR = 904 ) 528 904 IF( ios /=0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist', lwp )527 904 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist', lwp ) 529 528 IF(lwm) WRITE ( numonp, namp5zquota ) 530 529 ! 531 530 IF(lwp) THEN ! control print 532 531 WRITE(numout,*) ' ' … … 552 551 WRITE(numout,*) ' Maximal Fe quota for diatoms qfdmax = ', qfdmax 553 552 ENDIF 554 553 ! 555 554 zpsino3 = 2.3 * rno3 556 555 zpsinh4 = 1.8 * rno3 … … 560 559 ! 561 560 END SUBROUTINE p5z_lim_init 561 562 562 563 563 INTEGER FUNCTION p5z_lim_alloc() … … 569 569 !!---------------------------------------------------------------------- 570 570 ierr(:) = 0 571 571 ! 572 572 !* Biological arrays for phytoplankton growth 573 573 ALLOCATE( xpicono3(jpi,jpj,jpk), xpiconh4(jpi,jpj,jpk), & … … 585 585 & fvpuptk (jpi,jpj,jpk), xlimpic (jpi,jpj,jpk), & 586 586 & xlimsi (jpi,jpj,jpk), STAT=ierr(1) ) 587 587 ! 588 588 !* Minimum/maximum quotas of phytoplankton 589 589 ALLOCATE( xqnnmin (jpi,jpj,jpk), xqnnmax(jpi,jpj,jpk), & … … 593 593 & xqndmin (jpi,jpj,jpk), xqndmax(jpi,jpj,jpk), & 594 594 & xqpdmin (jpi,jpj,jpk), xqpdmax(jpi,jpj,jpk), STAT=ierr(2) ) 595 !595 ! 596 596 p5z_lim_alloc = MAXVAL( ierr ) 597 597 !
Note: See TracChangeset
for help on using the changeset viewer.