New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11822 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES – NEMO

Ignore:
Timestamp:
2019-10-29T11:41:36+01:00 (5 years ago)
Author:
acc
Message:

Branch 2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps. Sette tested updates to branch to align with trunk changes between 10721 and 11740. Sette tests are passing but results differ from branch before these changes (except for GYRE_PISCES and VORTEX) and branch results already differed from trunk because of algorithmic fixes. Will need more checks to confirm correctness.

Location:
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES
Files:
29 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P2Z/p2zbio.F90

    r10975 r11822  
    405405      REWIND( numnatp_ref )              ! Namelist namlobphy in reference namelist : Lobster biological parameters 
    406406      READ  ( numnatp_ref, namlobphy, IOSTAT = ios, ERR = 901) 
    407 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobphy in reference namelist', lwp ) 
     407901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobphy in reference namelist' ) 
    408408      REWIND( numnatp_cfg )              ! Namelist namlobphy in configuration namelist : Lobster biological parameters 
    409409      READ  ( numnatp_cfg, namlobphy, IOSTAT = ios, ERR = 902 ) 
    410 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobphy in configuration namelist', lwp ) 
     410902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobphy in configuration namelist' ) 
    411411      IF(lwm) WRITE ( numonp, namlobphy ) 
    412412      ! 
    … …  
    422422      REWIND( numnatp_ref )              ! Namelist namlobnut in reference namelist : Lobster nutriments parameters 
    423423      READ  ( numnatp_ref, namlobnut, IOSTAT = ios, ERR = 903) 
    424 903   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobnut in reference namelist', lwp ) 
     424903   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobnut in reference namelist' ) 
    425425      REWIND( numnatp_cfg )              ! Namelist namlobnut in configuration namelist : Lobster nutriments parameters 
    426426      READ  ( numnatp_cfg, namlobnut, IOSTAT = ios, ERR = 904 ) 
    427 904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobnut in configuration namelist', lwp ) 
     427904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobnut in configuration namelist' ) 
    428428      IF(lwm) WRITE ( numonp, namlobnut ) 
    429429 
    … …  
    439439      REWIND( numnatp_ref )              ! Namelist namlobzoo in reference namelist : Lobster zooplankton parameters 
    440440      READ  ( numnatp_ref, namlobzoo, IOSTAT = ios, ERR = 905) 
    441 905   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobzoo in reference namelist', lwp ) 
     441905   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobzoo in reference namelist' ) 
    442442      REWIND( numnatp_cfg )              ! Namelist namlobzoo in configuration namelist : Lobster zooplankton parameters 
    443443      READ  ( numnatp_cfg, namlobzoo, IOSTAT = ios, ERR = 906 ) 
    444 906   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobzoo in configuration namelist', lwp ) 
     444906   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobzoo in configuration namelist' ) 
    445445      IF(lwm) WRITE ( numonp, namlobzoo ) 
    446446 
    … …  
    461461      REWIND( numnatp_ref )              ! Namelist namlobdet in reference namelist : Lobster detritus parameters 
    462462      READ  ( numnatp_ref, namlobdet, IOSTAT = ios, ERR = 907) 
    463 907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobdet in reference namelist', lwp ) 
     463907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobdet in reference namelist' ) 
    464464      REWIND( numnatp_cfg )              ! Namelist namlobdet in configuration namelist : Lobster detritus parameters 
    465465      READ  ( numnatp_cfg, namlobdet, IOSTAT = ios, ERR = 908 ) 
    466 908   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobdet in configuration namelist', lwp ) 
     466908   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobdet in configuration namelist' ) 
    467467      IF(lwm) WRITE ( numonp, namlobdet ) 
    468468 
    … …  
    476476      REWIND( numnatp_ref )              ! Namelist namlobdom in reference namelist : Lobster DOM breakdown rate 
    477477      READ  ( numnatp_ref, namlobdom, IOSTAT = ios, ERR = 909) 
    478 909   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobdom in reference namelist', lwp ) 
     478909   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlobdom in reference namelist' ) 
    479479      REWIND( numnatp_cfg )              ! Namelist namlobdom in configuration namelist : Lobster DOM breakdown rate 
    480480      READ  ( numnatp_cfg, namlobdom, IOSTAT = ios, ERR = 910 ) 
    481 910   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobdom in configuration namelist', lwp ) 
     481910   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobdom in configuration namelist' ) 
    482482      IF(lwm) WRITE ( numonp, namlobdom ) 
    483483 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P2Z/p2zopt.F90

    r10975 r11822  
    162162      REWIND( numnatp_ref )              ! Namelist namlobopt in reference namelist : Lobster options 
    163163      READ  ( numnatp_ref, namlobopt, IOSTAT = ios, ERR = 901) 
    164 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobopt in reference namelist', lwp ) 
     164901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobopt in reference namelist' ) 
    165165 
    166166      REWIND( numnatp_cfg )              ! Namelist namlobopt in configuration namelist : Lobster options 
    167167      READ  ( numnatp_cfg, namlobopt, IOSTAT = ios, ERR = 902 ) 
    168 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namlobopt in configuration namelist', lwp ) 
     168902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namlobopt in configuration namelist' ) 
    169169      IF(lwm) WRITE ( numonp, namlobopt ) 
    170170 
    … …  
    184184      REWIND( numnatp_ref )              ! Namelist namlobrat in reference namelist : Lobster ratios 
    185185      READ  ( numnatp_ref, namlobrat, IOSTAT = ios, ERR = 903) 
    186 903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobrat in reference namelist', lwp ) 
     186903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobrat in reference namelist' ) 
    187187 
    188188      REWIND( numnatp_cfg )              ! Namelist namlobrat in configuration namelist : Lobster ratios 
    189189      READ  ( numnatp_cfg, namlobrat, IOSTAT = ios, ERR = 904 ) 
    190 904   IF( ios >  0 ) CALL ctl_nam ( ios , 'namlobrat in configuration namelist', lwp ) 
     190904   IF( ios >  0 ) CALL ctl_nam ( ios , 'namlobrat in configuration namelist' ) 
    191191      IF(lwm) WRITE ( numonp, namlobrat ) 
    192192 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P2Z/p2zsed.F90

    r10975 r11822  
    135135      REWIND( numnatp_ref )              ! Namelist namlobsed in reference namelist : Lobster sediments 
    136136      READ  ( numnatp_ref, namlobsed, IOSTAT = ios, ERR = 901) 
    137 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlosed in reference namelist', lwp ) 
     137901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namlosed in reference namelist' ) 
    138138      REWIND( numnatp_cfg )              ! Namelist namlobsed in configuration namelist : Lobster sediments 
    139139      READ  ( numnatp_cfg, namlobsed, IOSTAT = ios, ERR = 902 ) 
    140 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobsed in configuration namelist', lwp ) 
     140902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namlobsed in configuration namelist' ) 
    141141      IF(lwm) WRITE ( numonp, namlobsed ) 
    142142      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zfechem.F90

    r10975 r11822  
    257257      REWIND( numnatp_ref )            ! Namelist nampisfer in reference namelist : Pisces iron chemistry 
    258258      READ  ( numnatp_ref, nampisfer, IOSTAT = ios, ERR = 901) 
    259 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisfer in reference namelist', lwp ) 
     259901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisfer in reference namelist' ) 
    260260      REWIND( numnatp_cfg )            ! Namelist nampisfer in configuration namelist : Pisces iron chemistry 
    261261      READ  ( numnatp_cfg, nampisfer, IOSTAT = ios, ERR = 902 ) 
    262 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisfer in configuration namelist', lwp ) 
     262902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisfer in configuration namelist' ) 
    263263      IF(lwm) WRITE( numonp, nampisfer ) 
    264264 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zflx.F90

    r10975 r11822  
    242242      REWIND( numnatp_ref )              ! Namelist nampisext in reference namelist : Pisces atm. conditions 
    243243      READ  ( numnatp_ref, nampisext, IOSTAT = ios, ERR = 901) 
    244 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisext in reference namelist', lwp ) 
     244901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisext in reference namelist' ) 
    245245      REWIND( numnatp_cfg )              ! Namelist nampisext in configuration namelist : Pisces atm. conditions 
    246246      READ  ( numnatp_cfg, nampisext, IOSTAT = ios, ERR = 902 ) 
    247 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisext in configuration namelist', lwp ) 
     247902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisext in configuration namelist' ) 
    248248      IF(lwm) WRITE ( numonp, nampisext ) 
    249249      ! 
    … …  
    323323         REWIND( numnatp_ref )              ! Namelist nampisatm in reference namelist : Pisces atm. sea level pressure file 
    324324         READ  ( numnatp_ref, nampisatm, IOSTAT = ios, ERR = 901) 
    325 901      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisatm in reference namelist', lwp ) 
     325901      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisatm in reference namelist' ) 
    326326         REWIND( numnatp_cfg )              ! Namelist nampisatm in configuration namelist : Pisces atm. sea level pressure file  
    327327         READ  ( numnatp_cfg, nampisatm, IOSTAT = ios, ERR = 902 ) 
    328 902      IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisatm in configuration namelist', lwp ) 
     328902      IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisatm in configuration namelist' ) 
    329329         IF(lwm) WRITE ( numonp, nampisatm ) 
    330330         ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zligand.F90

    r10975 r11822  
    128128      REWIND( numnatp_ref )              ! Namelist nampislig in reference namelist : Pisces remineralization 
    129129      READ  ( numnatp_ref, nampislig, IOSTAT = ios, ERR = 901) 
    130 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampislig in reference namelist', lwp ) 
     130901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampislig in reference namelist' ) 
    131131      REWIND( numnatp_cfg )              ! Namelist nampislig in configuration namelist : Pisces remineralization 
    132132      READ  ( numnatp_cfg, nampislig, IOSTAT = ios, ERR = 902 ) 
    133 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampislig in configuration namelist', lwp ) 
     133902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampislig in configuration namelist' ) 
    134134      IF(lwm) WRITE ( numonp, nampislig ) 
    135135      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zlim.F90

    r10975 r11822  
    255255      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 
    256256      READ  ( numnatp_ref, namp4zlim, IOSTAT = ios, ERR = 901) 
    257 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zlim in reference namelist', lwp ) 
     257901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zlim in reference namelist' ) 
    258258      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    259259      READ  ( numnatp_cfg, namp4zlim, IOSTAT = ios, ERR = 902 ) 
    260 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zlim in configuration namelist', lwp ) 
     260902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zlim in configuration namelist' ) 
    261261      IF(lwm) WRITE( numonp, namp4zlim ) 
    262262      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zlys.F90

    r10975 r11822  
    165165      REWIND( numnatp_ref )              ! Namelist nampiscal in reference namelist : Pisces CaCO3 dissolution 
    166166      READ  ( numnatp_ref, nampiscal, IOSTAT = ios, ERR = 901) 
    167 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampiscal in reference namelist', lwp ) 
     167901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampiscal in reference namelist' ) 
    168168      REWIND( numnatp_cfg )              ! Namelist nampiscal in configuration namelist : Pisces CaCO3 dissolution 
    169169      READ  ( numnatp_cfg, nampiscal, IOSTAT = ios, ERR = 902 ) 
    170 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampiscal in configuration namelist', lwp ) 
     170902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampiscal in configuration namelist' ) 
    171171      IF(lwm) WRITE( numonp, nampiscal ) 
    172172      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zmeso.F90

    r10975 r11822  
    293293      REWIND( numnatp_ref )              ! Namelist nampismes in reference namelist : Pisces mesozooplankton 
    294294      READ  ( numnatp_ref, namp4zmes, IOSTAT = ios, ERR = 901) 
    295 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zmes in reference namelist', lwp ) 
     295901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zmes in reference namelist' ) 
    296296      REWIND( numnatp_cfg )              ! Namelist nampismes in configuration namelist : Pisces mesozooplankton 
    297297      READ  ( numnatp_cfg, namp4zmes, IOSTAT = ios, ERR = 902 ) 
    298 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zmes in configuration namelist', lwp ) 
     298902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zmes in configuration namelist' ) 
    299299      IF(lwm) WRITE( numonp, namp4zmes ) 
    300300      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zmicro.F90

    r10975 r11822  
    246246      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton 
    247247      READ  ( numnatp_ref, namp4zzoo, IOSTAT = ios, ERR = 901) 
    248 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zzoo in reference namelist', lwp ) 
     248901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zzoo in reference namelist' ) 
    249249      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton 
    250250      READ  ( numnatp_cfg, namp4zzoo, IOSTAT = ios, ERR = 902 ) 
    251 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zzoo in configuration namelist', lwp ) 
     251902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zzoo in configuration namelist' ) 
    252252      IF(lwm) WRITE( numonp, namp4zzoo ) 
    253253      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zmort.F90

    r10975 r11822  
    232232      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton 
    233233      READ  ( numnatp_ref, namp4zmort, IOSTAT = ios, ERR = 901) 
    234 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zmort in reference namelist', lwp ) 
     234901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zmort in reference namelist' ) 
    235235      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton 
    236236      READ  ( numnatp_cfg, namp4zmort, IOSTAT = ios, ERR = 902 ) 
    237 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zmort in configuration namelist', lwp ) 
     237902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zmort in configuration namelist' ) 
    238238      IF(lwm) WRITE( numonp, namp4zmort ) 
    239239      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zopt.F90

    r10975 r11822  
    404404      REWIND( numnatp_ref )              ! Namelist nampisopt in reference namelist : Pisces attenuation coef. and PAR 
    405405      READ  ( numnatp_ref, nampisopt, IOSTAT = ios, ERR = 901) 
    406 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisopt in reference namelist', lwp ) 
     406901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisopt in reference namelist' ) 
    407407      REWIND( numnatp_cfg )              ! Namelist nampisopt in configuration namelist : Pisces attenuation coef. and PAR 
    408408      READ  ( numnatp_cfg, nampisopt, IOSTAT = ios, ERR = 902 ) 
    409 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisopt in configuration namelist', lwp ) 
     409902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisopt in configuration namelist' ) 
    410410      IF(lwm) WRITE ( numonp, nampisopt ) 
    411411 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zpoc.F90

    r10975 r11822  
    103103     ! particle scheme. Does not work with Kriest parameterization. 
    104104     ! ----------------------------------------------------------------------- 
     105     ztremint(:,:,:) = zremigoc(:,:,:) 
    105106     DO jk = 2, jpkm1 
    106107        DO jj = 1, jpj 
    … …  
    167168                       &   + zsizek1 ) ) * zpoc + ( prodgoc(ji,jj,jk-1) / tgfunc(ji,jj,jk-1) * ( 1.           & 
    168169                       &   - exp( -reminp(jn) * zsizek1 ) ) * exp( -reminp(jn) * zsizek ) + prodgoc(ji,jj,jk) & 
    169                        &   / tgfunc(ji,jj,jk) * ( 1. - exp( -reminp(jn) * zsizek ) ) ) * rday / rfact2 / reminp(jn)  
     170                       &   / tgfunc(ji,jj,jk) * ( 1. - exp( -reminp(jn) * zsizek ) ) ) * rday / rfact2 / reminp(jn) * alphan(jn)  
    170171                       alphat = alphat + alphag(ji,jj,jk,jn) 
    171172                       remint = remint + alphag(ji,jj,jk,jn) * reminp(jn) 
    … …  
    286287     ! layer, this spectrum is supposed to be uniform. 
    287288     ! --------------------------------------------------------------------- 
     289     ztremint(:,:,:) = zremipoc(:,:,:) 
    288290     DO jk = 1, jpkm1 
    289291        DO jj = 1, jpj 
    … …  
    498500      REWIND( numnatp_ref )              ! Namelist nampisrem in reference namelist : Pisces remineralization 
    499501      READ  ( numnatp_ref, nampispoc, IOSTAT = ios, ERR = 901) 
    500 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampispoc in reference namelist', lwp ) 
     502901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampispoc in reference namelist' ) 
    501503      REWIND( numnatp_cfg )              ! Namelist nampisrem in configuration namelist : Pisces remineralization 
    502504      READ  ( numnatp_cfg, nampispoc, IOSTAT = ios, ERR = 902 ) 
    503 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampispoc in configuration namelist', lwp ) 
     505902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampispoc in configuration namelist' ) 
    504506      IF(lwm) WRITE( numonp, nampispoc ) 
    505507 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zprod.F90

    r10975 r11822  
    226226         DO jj = 1, jpj 
    227227            DO ji = 1, jpi 
    228                zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
    229                zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
    230228               zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
    231229               zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * ( 1. - fr_i(ji,jj) ) 
    … …  
    323321     ! 
    324322     IF( ln_ligand ) THEN 
     323         zpligprod1(:,:,:) = 0._wp    ;    zpligprod2(:,:,:) = 0._wp 
    325324         DO jk = 1, jpkm1 
    326325            DO jj = 1, jpj 
    … …  
    494493      REWIND( numnatp_ref )              ! Namelist nampisprod in reference namelist : Pisces phytoplankton production 
    495494      READ  ( numnatp_ref, namp4zprod, IOSTAT = ios, ERR = 901) 
    496 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zprod in reference namelist', lwp ) 
     495901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namp4zprod in reference namelist' ) 
    497496      REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production 
    498497      READ  ( numnatp_cfg, namp4zprod, IOSTAT = ios, ERR = 902 ) 
    499 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zprod in configuration namelist', lwp ) 
     498902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namp4zprod in configuration namelist' ) 
    500499      IF(lwm) WRITE( numonp, namp4zprod ) 
    501500 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zrem.F90

    r10975 r11822  
    330330      REWIND( numnatp_ref )              ! Namelist nampisrem in reference namelist : Pisces remineralization 
    331331      READ  ( numnatp_ref, nampisrem, IOSTAT = ios, ERR = 901) 
    332 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisrem in reference namelist', lwp ) 
     332901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisrem in reference namelist' ) 
    333333      REWIND( numnatp_cfg )              ! Namelist nampisrem in configuration namelist : Pisces remineralization 
    334334      READ  ( numnatp_cfg, nampisrem, IOSTAT = ios, ERR = 902 ) 
    335 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp ) 
     335902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisrem in configuration namelist' ) 
    336336      IF(lwm) WRITE( numonp, nampisrem ) 
    337337 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zsbc.F90

    r10975 r11822  
    2727   LOGICAL , PUBLIC ::   ln_ironsed   !: boolean for Fe input from sediments 
    2828   LOGICAL , PUBLIC ::   ln_hydrofe   !: boolean for Fe input from hydrothermal vents 
    29    LOGICAL , PUBLIC ::   ln_ironice   !: boolean for Fe input from sea ice 
    3029   REAL(wp), PUBLIC ::   sedfeinput   !: Coastal release of Iron 
    3130   REAL(wp), PUBLIC ::   dustsolub    !: Solubility of the dust 
    … …  
    110109         IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_dust > 1 ) ) THEN 
    111110            CALL fld_read( kt, 1, sf_dust ) 
    112             IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN   ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) 
    113             ELSE                                               ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) ) 
    114             ENDIF 
     111            dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) * ( 1.0 - fr_i(:,:) ) 
    115112         ENDIF 
    116113      ENDIF 
    … …  
    176173      IF( ln_ndepo ) THEN 
    177174         IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_ndep > 1 ) ) THEN 
    178              zcoef = rno3 * 14E6 * ryyss 
     175             zcoef = 14. * rno3 
    179176             CALL fld_read( kt, 1, sf_ndepo ) 
    180177             nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t(:,:,1,Kmm) ) 
    181178         ENDIF 
    182179         IF( .NOT.ln_linssh ) THEN 
    183            zcoef = rno3 * 14E6 * ryyss 
    184            nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t(:,:,1,Kmm) ) 
     180            zcoef = 14. * rno3 
     181            nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t(:,:,1,Kmm) ) 
    185182         ENDIF 
    186183      ENDIF 
    … …  
    237234      REWIND( numnatp_ref )              ! Namelist nampissbc in reference namelist : Pisces external sources of nutrients 
    238235      READ  ( numnatp_ref, nampissbc, IOSTAT = ios, ERR = 901) 
    239 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampissbc in reference namelist', lwp ) 
     236901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampissbc in reference namelist' ) 
    240237      REWIND( numnatp_cfg )              ! Namelist nampissbc in configuration namelist : Pisces external sources of nutrients 
    241238      READ  ( numnatp_cfg, nampissbc, IOSTAT = ios, ERR = 902 ) 
    242 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampissbc in configuration namelist', lwp ) 
     239902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampissbc in configuration namelist' ) 
    243240      IF(lwm) WRITE ( numonp, nampissbc ) 
    244241 
    … …  
    267264      END IF 
    268265 
    269       IF( nn_ice_tr >= 0 .AND. ln_ironice ) THEN 
    270          IF(lwp) THEN 
    271             WRITE(numout,*) '   ==>>>   ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 
    272             WRITE(numout,*) '           Specify your sea ice iron concentration in nampisice instead ' 
    273             WRITE(numout,*) '           ln_ironice is forced to .FALSE. ' 
    274          ENDIF 
    275          ln_ironice = .FALSE. 
    276       ENDIF 
    277  
    278266      IF( ln_dust .OR. ln_river .OR. ln_ndepo ) THEN   ;   ll_sbc = .TRUE. 
    279267      ELSE                                             ;   ll_sbc = .FALSE. 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zsed.F90

    r10975 r11822  
    5454      INTEGER  ::  ji, jj, jk, ikt 
    5555      REAL(wp) ::  zrivalk, zrivsil, zrivno3 
    56       REAL(wp) ::  zwflux, zfminus, zfplus 
    57       REAL(wp) ::  zlim, zfact, zfactcal 
     56      REAL(wp) ::  zwflux, zlim, zfact, zfactcal 
    5857      REAL(wp) ::  zo2, zno3, zflx, zpdenit, z1pdenit, zolimit 
    5958      REAL(wp) ::  zsiloss, zcaloss, zws3, zws4, zwsc, zdep 
    … …  
    105104               zdep    = rfact2 / e3t(ji,jj,1,Kmm) 
    106105               zwflux  = fmmflx(ji,jj) / 1000._wp 
    107                zfminus = MIN( 0._wp, -zwflux ) * tr(ji,jj,1,jpfer,Kbb) * zdep 
    108                zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep 
    109                zironice(ji,jj) =  zfplus + zfminus 
     106               zironice(ji,jj) =  MAX( -0.99 * tr(ji,jj,1,jpfer,Kbb), -zwflux * icefeinput * zdep ) 
    110107            END DO 
    111108         END DO 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zsms.F90

    r11504 r11822  
    204204      REWIND( numnatp_ref )              ! Namelist nampisbio in reference namelist : Pisces variables 
    205205      READ  ( numnatp_ref, nampisbio, IOSTAT = ios, ERR = 901) 
    206 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisbio in reference namelist', lwp ) 
     206901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisbio in reference namelist' ) 
    207207      REWIND( numnatp_cfg )              ! Namelist nampisbio in configuration namelist : Pisces variables 
    208208      READ  ( numnatp_cfg, nampisbio, IOSTAT = ios, ERR = 902 ) 
    209 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisbio in configuration namelist', lwp ) 
     209902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisbio in configuration namelist' ) 
    210210      IF(lwm) WRITE( numonp, nampisbio ) 
    211211      ! 
    … …  
    235235      REWIND( numnatp_ref )              ! Namelist nampisdmp in reference namelist : Pisces damping 
    236236      READ  ( numnatp_ref, nampisdmp, IOSTAT = ios, ERR = 905) 
    237 905   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisdmp in reference namelist', lwp ) 
     237905   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampisdmp in reference namelist' ) 
    238238      REWIND( numnatp_cfg )              ! Namelist nampisdmp in configuration namelist : Pisces damping 
    239239      READ  ( numnatp_cfg, nampisdmp, IOSTAT = ios, ERR = 906 ) 
    240 906   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisdmp in configuration namelist', lwp ) 
     240906   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampisdmp in configuration namelist' ) 
    241241      IF(lwm) WRITE( numonp, nampisdmp ) 
    242242      ! 
    … …  
    250250      REWIND( numnatp_ref )              ! Namelist nampismass in reference namelist : Pisces mass conservation check 
    251251      READ  ( numnatp_ref, nampismass, IOSTAT = ios, ERR = 907) 
    252 907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampismass in reference namelist', lwp ) 
     252907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampismass in reference namelist' ) 
    253253      REWIND( numnatp_cfg )              ! Namelist nampismass in configuration namelist : Pisces mass conservation check  
    254254      READ  ( numnatp_cfg, nampismass, IOSTAT = ios, ERR = 908 ) 
    255 908   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampismass in configuration namelist', lwp ) 
     255908   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampismass in configuration namelist' ) 
    256256      IF(lwm) WRITE( numonp, nampismass ) 
    257257 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p5zlim.F90

    r10975 r11822  
    451451      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 
    452452      READ  ( numnatp_ref, namp5zlim, IOSTAT = ios, ERR = 901) 
    453 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in reference namelist', lwp ) 
     453901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in reference namelist' ) 
    454454      ! 
    455455      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    456456      READ  ( numnatp_cfg, namp5zlim, IOSTAT = ios, ERR = 902 ) 
    457 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp ) 
     457902   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist' ) 
    458458      IF(lwm) WRITE ( numonp, namp5zlim ) 
    459459      ! 
    … …  
    492492      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 
    493493      READ  ( numnatp_ref, namp5zquota, IOSTAT = ios, ERR = 903) 
    494 903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisquota in reference namelist', lwp ) 
     494903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisquota in reference namelist' ) 
    495495      ! 
    496496      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    497497      READ  ( numnatp_cfg, namp5zquota, IOSTAT = ios, ERR = 904 ) 
    498 904   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist', lwp ) 
     498904   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist' ) 
    499499      IF(lwm) WRITE ( numonp, namp5zquota ) 
    500500      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p5zmeso.F90

    r10975 r11822  
    410410      REWIND( numnatp_ref )              ! Namelist nampismes in reference namelist : Pisces mesozooplankton 
    411411      READ  ( numnatp_ref, namp5zmes, IOSTAT = ios, ERR = 901) 
    412 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in reference namelist', lwp ) 
     412901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in reference namelist' ) 
    413413      ! 
    414414      REWIND( numnatp_cfg )              ! Namelist nampismes in configuration namelist : Pisces mesozooplankton 
    415415      READ  ( numnatp_cfg, namp5zmes, IOSTAT = ios, ERR = 902 ) 
    416 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampismes in configuration namelist', lwp ) 
     416902   IF( ios >  0 ) CALL ctl_nam ( ios , 'nampismes in configuration namelist' ) 
    417417      IF(lwm) WRITE ( numonp, namp5zmes ) 
    418418      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p5zmicro.F90

    r10975 r11822  
    352352      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton 
    353353      READ  ( numnatp_ref, namp5zzoo, IOSTAT = ios, ERR = 901) 
    354 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zzoo in reference namelist', lwp ) 
     354901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zzoo in reference namelist' ) 
    355355      ! 
    356356      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton 
    357357      READ  ( numnatp_cfg, namp5zzoo, IOSTAT = ios, ERR = 902 ) 
    358 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zzoo in configuration namelist', lwp ) 
     358902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zzoo in configuration namelist' ) 
    359359      IF(lwm) WRITE ( numonp, namp5zzoo ) 
    360360      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p5zmort.F90

    r10975 r11822  
    292292      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton 
    293293      READ  ( numnatp_ref, namp5zmort, IOSTAT = ios, ERR = 901) 
    294 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zmort in reference namelist', lwp ) 
     294901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zmort in reference namelist' ) 
    295295 
    296296      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton 
    297297      READ  ( numnatp_cfg, namp5zmort, IOSTAT = ios, ERR = 902 ) 
    298 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zmort in configuration namelist', lwp ) 
     298902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zmort in configuration namelist' ) 
    299299      IF(lwm) WRITE ( numonp, namp5zmort ) 
    300300 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p5zprod.F90

    r10975 r11822  
    445445     ! 
    446446     IF( ln_ligand ) THEN 
     447         zpligprod1(:,:,:) = 0._wp    ;    zpligprod2(:,:,:) = 0._wp 
    447448         DO jk = 1, jpkm1 
    448449            DO jj = 1, jpj 
    … …  
    584585      REWIND( numnatp_ref )              ! Namelist nampisprod in reference namelist : Pisces phytoplankton production 
    585586      READ  ( numnatp_ref, namp5zprod, IOSTAT = ios, ERR = 901) 
    586 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zprod in reference namelist', lwp ) 
     587901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zprod in reference namelist' ) 
    587588 
    588589      REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production 
    589590      READ  ( numnatp_cfg, namp5zprod, IOSTAT = ios, ERR = 902 ) 
    590 902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zprod in configuration namelist', lwp ) 
     591902   IF( ios >  0 ) CALL ctl_nam ( ios , 'namp5zprod in configuration namelist' ) 
    591592      IF(lwm) WRITE ( numonp, namp5zprod ) 
    592593 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/SED/sedini.F90

    r11504 r11822  
    459459      REWIND( numnamsed_ref )              ! Namelist nam_run in reference namelist : Pisces variables 
    460460      READ  ( numnamsed_ref, nam_run, IOSTAT = ios, ERR = 901) 
    461 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in reference namelist', lwp ) 
     461901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in reference namelist' ) 
    462462 
    463463      REWIND( numnamsed_cfg )              ! Namelist nam_run in reference namelist : Pisces variables 
    464464      READ  ( numnamsed_cfg, nam_run, IOSTAT = ios, ERR = 902) 
    465 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in configuration namelist', lwp ) 
     465902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in configuration namelist' ) 
    466466 
    467467      IF (lwp) THEN 
    … …  
    475475      REWIND( numnamsed_ref )              ! Namelist nam_geom in reference namelist : Pisces variables 
    476476      READ  ( numnamsed_ref, nam_geom, IOSTAT = ios, ERR = 903) 
    477 903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in reference namelist', lwp ) 
     477903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in reference namelist' ) 
    478478 
    479479      REWIND( numnamsed_cfg )              ! Namelist nam_geom in reference namelist : Pisces variables 
    480480      READ  ( numnamsed_cfg, nam_geom, IOSTAT = ios, ERR = 904) 
    481 904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in configuration namelist', lwp ) 
     481904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in configuration namelist' ) 
    482482 
    483483      IF (lwp) THEN  
    … …  
    498498      REWIND( numnamsed_ref )              ! Namelist nam_trased in reference namelist : Pisces variables 
    499499      READ  ( numnamsed_ref, nam_trased, IOSTAT = ios, ERR = 905) 
    500 905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in reference namelist', lwp ) 
     500905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in reference namelist' ) 
    501501 
    502502      REWIND( numnamsed_cfg )              ! Namelist nam_trased in reference namelist : Pisces variables 
    503503      READ  ( numnamsed_cfg, nam_trased, IOSTAT = ios, ERR = 906) 
    504 906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in configuration namelist', lwp ) 
     504906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in configuration namelist' ) 
    505505 
    506506      DO jn = 1, jpsol 
    … …  
    531531      REWIND( numnamsed_ref )              ! Namelist nam_diased in reference namelist : Pisces variables 
    532532      READ  ( numnamsed_ref, nam_diased, IOSTAT = ios, ERR = 907) 
    533 907   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in reference namelist', lwp ) 
     533907   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in reference namelist' ) 
    534534 
    535535      REWIND( numnamsed_cfg )              ! Namelist nam_diased in reference namelist : Pisces variables 
    536536      READ  ( numnamsed_cfg, nam_diased, IOSTAT = ios, ERR = 908) 
    537 908   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in configuration namelist', lwp ) 
     537908   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in configuration namelist' ) 
    538538       
    539539      DO jn = 1, jpdia3dsed 
    … …  
    573573      REWIND( numnamsed_ref )              ! Namelist nam_inorg in reference namelist : Pisces variables 
    574574      READ  ( numnamsed_ref, nam_inorg, IOSTAT = ios, ERR = 909) 
    575 909   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in reference namelist', lwp ) 
     575909   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in reference namelist' ) 
    576576 
    577577      REWIND( numnamsed_cfg )              ! Namelist nam_inorg in reference namelist : Pisces variables 
    578578      READ  ( numnamsed_cfg, nam_inorg, IOSTAT = ios, ERR = 910) 
    579 910   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in configuration namelist', lwp ) 
     579910   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in configuration namelist' ) 
    580580 
    581581      IF (lwp) THEN 
    … …  
    599599      REWIND( numnamsed_ref )              ! Namelist nam_poc in reference namelist : Pisces variables 
    600600      READ  ( numnamsed_ref, nam_poc, IOSTAT = ios, ERR = 911) 
    601 911   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in reference namelist', lwp ) 
     601911   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in reference namelist' ) 
    602602 
    603603      REWIND( numnamsed_cfg )              ! Namelist nam_poc in reference namelist : Pisces variables 
    604604      READ  ( numnamsed_cfg, nam_poc, IOSTAT = ios, ERR = 912) 
    605 912   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in configuration namelist', lwp ) 
     605912   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in configuration namelist' ) 
    606606 
    607607      IF (lwp) THEN 
    … …  
    651651      REWIND( numnamsed_ref )              ! Namelist nam_btb in reference namelist : Pisces variables 
    652652      READ  ( numnamsed_ref, nam_btb, IOSTAT = ios, ERR = 913) 
    653 913   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in reference namelist', lwp ) 
     653913   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in reference namelist' ) 
    654654 
    655655      REWIND( numnamsed_cfg )              ! Namelist nam_btb in reference namelist : Pisces variables 
    656656      READ  ( numnamsed_cfg, nam_btb, IOSTAT = ios, ERR = 914) 
    657 914   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in configuration namelist', lwp ) 
     657914   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in configuration namelist' ) 
    658658 
    659659      IF (lwp) THEN 
    … …  
    672672      REWIND( numnamsed_ref )              ! Namelist nam_rst in reference namelist : Pisces variables 
    673673      READ  ( numnamsed_ref, nam_rst, IOSTAT = ios, ERR = 915) 
    674 915   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in reference namelist', lwp ) 
     674915   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in reference namelist' ) 
    675675 
    676676      REWIND( numnamsed_cfg )              ! Namelist nam_rst in reference namelist : Pisces variables 
    677677      READ  ( numnamsed_cfg, nam_rst, IOSTAT = ios, ERR = 916) 
    678 916   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in configuration namelist', lwp ) 
     678916   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in configuration namelist' ) 
    679679 
    680680      IF (lwp) THEN 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/SED/sedrst.F90

    r11504 r11822  
    4949            IF( ln_rst_list ) THEN 
    5050               nrst_lst = 1 
    51                nitrst = nstocklist( nrst_lst ) 
     51               nitrst = nn_stocklist( nrst_lst ) 
    5252            ELSE 
    5353               nitrst = nitend 
    5454            ENDIF 
    5555         ENDIF 
    56          IF( .NOT. ln_rst_list .AND. MOD( kt - 1, nstock ) == 0 ) THEN 
     56         IF( .NOT. ln_rst_list .AND. MOD( kt - 1, nn_stock ) == 0 ) THEN 
    5757            ! we use kt - 1 and not kt - nittrc000 to keep the same periodicity from the beginning of the experiment 
    58             nitrst = kt + nstock - 1                  ! define the next value of nitrst for restart writing 
     58            nitrst = kt + nn_stock - 1                  ! define the next value of nitrst for restart writing 
    5959            IF( nitrst > nitend )   nitrst = nitend   ! make sure we write a restart at the end of the run 
    6060         ENDIF 
    … …  
    6262         IF( kt == nittrc000 ) lrst_sed = .FALSE. 
    6363      ENDIF 
     64 
     65      IF( .NOT. ln_rst_list .AND. nn_stock == -1 )   RETURN   ! we will never do any restart 
    6466 
    6567      ! to get better performances with NetCDF format: 
    6668      ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 1) 
    6769      ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 1 
    68       IF( kt == nitrst - 2*nn_dtsed .OR. nstock == nn_dtsed .OR. ( kt == nitend - nn_dtsed .AND. .NOT. lrst_sed ) ) THEN 
     70      IF( kt == nitrst - 2*nn_dtsed .OR. nn_stock == nn_dtsed .OR. ( kt == nitend - nn_dtsed .AND. .NOT. lrst_sed ) ) THEN 
    6971         ! beware of the format used to write kt (default is i8.8, that should be large enough) 
    7072         IF( nitrst > 1.0e9 ) THEN   ;   WRITE(clkt,*       ) nitrst 
    … …  
    300302          IF( l_offline .AND. ln_rst_list ) THEN 
    301303             nrst_lst = nrst_lst + 1 
    302              nitrst = nstocklist( nrst_lst ) 
     304             nitrst = nn_stocklist( nrst_lst ) 
    303305          ENDIF 
    304306      ENDIF 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/sms_pisces.F90

    r10425 r11822  
    6767   !!* Mass conservation 
    6868   LOGICAL  ::  ln_check_mass      !: Flag to check mass conservation 
     69   LOGICAL , PUBLIC ::   ln_ironice   !: boolean for Fe input from sea ice 
    6970 
    7071   !!*  Biological fluxes for light : variables shared by pisces & lobster 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/trcice_pisces.F90

    r10213 r11822  
    119119      zpisc(jpno3,1) =  5.79e-6_wp / rno3  
    120120      zpisc(jpnh4,1) =  3.22e-7_wp / rno3 
     121      zpisc(jplgw,1) =  1.0e-9_wp 
    121122 
    122123      ! ln_p5z 
    … …  
    162163      zpisc(jpno3,2) =  3.51e-06_wp / rno3  
    163164      zpisc(jpnh4,2) =  6.15e-08_wp / rno3  
     165      zpisc(jplgw,2) =  1.0e-9_wp 
    164166 
    165167      ! ln_p5z 
    … …  
    205207      zpisc(jpno3,3) =  2.64e-5_wp / rno3   
    206208      zpisc(jpnh4,3) =  3.39e-7_wp / rno3   
     209      zpisc(jplgw,3) =  1.0e-9_wp 
    207210 
    208211      ! ln_p5z 
    … …  
    249252      zpisc(jpno3,4) = 5.36e-5_wp / rno3 
    250253      zpisc(jpnh4,4) = 7.18e-7_wp / rno3 
     254      zpisc(jplgw,4) = 1.0e-9_wp 
    251255 
    252256      ! ln_p5z 
    … …  
    266270      zpisc(jpgon,4) =  1.05e-8_wp 
    267271      zpisc(jpgop,4) =  1.05e-8_wp 
    268   
     272! 
     273!     ln_ironice and tracers in seaice are redundant. Thus, if tracers in ice 
     274!     is activated, ln_ironice should be set to false 
     275!     ------------------------------------------------------------------------ 
     276      IF( nn_ice_tr /= 0 .AND. ln_ironice ) THEN 
     277         IF(lwp) THEN 
     278            WRITE(numout,*) '   ==>>>   ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 
     279            WRITE(numout,*) '           Specify your sea ice iron concentration in nampisice instead ' 
     280            WRITE(numout,*) '           ln_ironice is forced to .FALSE. ' 
     281         ENDIF 
     282         ln_ironice = .FALSE. 
     283      ENDIF 
     284!  
    269285      DO jn = jp_pcs0, jp_pcs1 
    270286         IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1)  ! Global case 
    … …  
    310326            trc_i(:,:,jn) = trc_ice_prescr(jn) 
    311327         ENDIF 
    312         
    313328         !-- Baltic 
    314329         IF( cn_cfg == "orca" .OR. cn_cfg == "ORCA" ) THEN      
    … …  
    318333                     trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn)  
    319334               END WHERE 
    320             ELSE                                 ! prescribed tracer concentration in ice 
    321                WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    322                    54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    323                      trc_i(:,:,jn) = trc_ice_prescr(jn) 
    324                END WHERE 
    325             ENDIF ! trc_ice_ratio 
     335            ENDIF 
    326336         ENDIF 
    327337      ! 

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/trcini_pisces.F90

    r10975 r11822  
    277277 
    278278      ! Initialization of the sediment model 
    279       IF( ln_sediment) THEN 
    280          CALL sed_init 
    281       ELSE 
    282          ALLOCATE( profsed(2) ) 
    283       ENDIF 
     279      IF( ln_sediment)   CALL sed_init 
    284280 
    285281      IF(lwp) WRITE(numout,*)  

  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/trcnam_pisces.F90

    r10222 r11822  
    5757      REWIND( numnatp_ref )              ! Namelist nampisbio in reference namelist : Pisces variables 
    5858      READ  ( numnatp_ref, nampismod, IOSTAT = ios, ERR = 901) 
    59 901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampismod in reference namelist', lwp ) 
     59901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nampismod in reference namelist' ) 
    6060      REWIND( numnatp_cfg )              ! Namelist nampisbio in configuration namelist : Pisces variables 
    6161      READ  ( numnatp_cfg, nampismod, IOSTAT = ios, ERR = 902 ) 
    62 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampismod in configuration namelist', lwp ) 
     62902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nampismod in configuration namelist' ) 
    6363      IF(lwm) WRITE( numonp, nampismod ) 
    6464      ! 
Note: See TracChangeset for help on using the changeset viewer.