Changeset 7190


Ignore:
Timestamp:
2016-11-04T12:19:19+01:00 (4 years ago)
Author:
aumont
Message:

various important bug fixes

Location:
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zagg.F90

    r6841 r7190  
    128128               tra(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) - zaggdop - zaggdop2 - zaggdop3 
    129129               ! 
    130                conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zagg + zaggdoc + zaggdoc3 
     130               conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zaggpoc + zaggdoc + zaggdoc3 
    131131               prodgoc(ji,jj,jk) = prodgoc(ji,jj,jk) + zaggpoc + zaggdoc2 
    132132               ! 
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zbio.F90

    r6966 r7190  
    2828   USE p5zpoc          !  Remineralisation of organic particles 
    2929   USE p5zagg          !  Aggregation of particles 
     30   USE p4zfechem       !  Iron chemistry 
    3031   USE p4zligand       !  Remineralization of ligands 
    3132   USE prtctl_trc      !  print control for debugging 
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zpoc.F90

    r7180 r7190  
    289289                   END DO 
    290290                   ! Mean remineralization rate in the mixed layer 
    291                    zremipoc(ji,jj,jk) = MIN(xremip, MAX(0., remint )) 
     291                   zremipoc(ji,jj,jk) = MIN(xremipc, MAX(0., remint )) 
    292292                ENDIF 
    293293              ENDIF 
     
    370370                  END DO 
    371371                  ! Mean remineralization rate in the water column 
    372                   zremipoc(ji,jj,jk) = MIN(xremip, MAX(0., remint )) 
     372                  zremipoc(ji,jj,jk) = MIN(xremipc, MAX(0., remint )) 
    373373                ENDIF 
    374374              ENDIF 
     
    493493         ! 
    494494         alphan(1) = gamain(reminup, rshape, ifault) 
    495          reminp(1) = gamain(reminup, rshape+1.0, ifault) * xremip / alphan(1) 
     495         reminp(1) = gamain(reminup, rshape+1.0, ifault) * xremipc / alphan(1) 
    496496         DO jn = 2, jcpoc-1 
    497497            reminup = 1./ 400. * exp(float(jn) * remindelta) 
     
    499499            alphan(jn) = gamain(reminup, rshape, ifault) - gamain(remindown, rshape, ifault) 
    500500            reminp(jn) = gamain(reminup, rshape+1.0, ifault) - gamain(remindown, rshape+1.0, ifault) 
    501             reminp(jn) = reminp(jn) * xremip / alphan(jn) 
     501            reminp(jn) = reminp(jn) * xremipc / alphan(jn) 
    502502         END DO 
    503503         remindown = 1. / 400. * exp(float(jcpoc-1) * remindelta) 
    504504         alphan(jcpoc) = 1.0 - gamain(remindown, rshape, ifault) 
    505          reminp(jcpoc) = reminp(jcpoc) * xremip / alphan(jcpoc) 
     505         reminp(jcpoc) = 1.0 - gamain(remindown, rshape+1.0, ifault) 
     506         reminp(jcpoc) = reminp(jcpoc) * xremipc / alphan(jcpoc) 
    506507      ELSE 
    507508         alphan(jcpoc) = 1. 
    508          reminp(jcpoc) = xremip 
     509         reminp(jcpoc) = xremipc 
    509510      ENDIF 
    510511 
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zprod.F90

    r7184 r7190  
    9090      CHARACTER (len=25) :: charout 
    9191      REAL(wp), POINTER, DIMENSION(:,:  ) :: zmixnano, zmixpico, zmixdiat, zstrn 
    92       REAL(wp), POINTER, DIMENSION(:,:,:) :: zpislopeadn, zpislopeadp, zpislopead 
     92      REAL(wp), POINTER, DIMENSION(:,:,:) :: zpislopeadn, zpislopeadp, zpislopeadd 
    9393      REAL(wp), POINTER, DIMENSION(:,:,:) :: zprbio, zprpic, zprdia, zysopt 
    9494      REAL(wp), POINTER, DIMENSION(:,:,:) :: zprchln, zprchlp, zprchld 
     
    109109      CALL wrk_alloc( jpi, jpj,      zmixnano, zmixpico, zmixdiat, zstrn ) 
    110110      CALL wrk_alloc( jpi, jpj, jpk, zmxl_fac, zmxl_chl ) 
    111       CALL wrk_alloc( jpi, jpj, jpk, zpislopeadn, zpislopeadp, zpislopead, zysopt )  
     111      CALL wrk_alloc( jpi, jpj, jpk, zpislopeadn, zpislopeadp, zpislopeadd, zysopt )  
    112112      CALL wrk_alloc( jpi, jpj, jpk, zprdia, zprpic, zprbio, zprorcan, zprorcap, zprorcad ) 
    113113      CALL wrk_alloc( jpi, jpj, jpk, zprofed, zprofep, zprofen ) 
     
    548548      ! 
    549549      NAMELIST/nampisprod/ pislopen, pislopep, pisloped, xadap, bresp, excretn, excretp, excretd,  & 
    550          &                 thetannm, thetanpm, thetandm, chlcmin, grosip            & 
     550         &                 thetannm, thetanpm, thetandm, chlcmin, grosip 
    551551 
    552552      INTEGER :: ios                 ! Local integer output status for namelist read 
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zrem.F90

    r7180 r7190  
    2424   USE p4zche          !  chemical model 
    2525   USE p5zlim          ! Phytoplankton limitation factors 
     26   USE p5zsink         !  Sinking of particles 
    2627   USE p5zprod         !  Production by phytoplankton 
    2728   USE prtctl_trc      !  print control for debugging 
     
    267268               tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) - zosil 
    268269               tra(ji,jj,jk,jpsil) = tra(ji,jj,jk,jpsil) + zosil 
    269                ENDIF 
    270270            END DO 
    271271         END DO 
  • branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zsms.F90

    r6966 r7190  
    133133         ! 
    134134         CALL p5z_bio (kt, knt)    ! Biology 
    135          CALL p4z_lys( kt, jnt )   ! Compute CaCO3 saturation 
     135         CALL p4z_lys( kt, knt )   ! Compute CaCO3 saturation 
    136136         CALL p4z_flx( kt, knt )   ! Compute surface fluxes 
    137137         CALL p5z_sed (kt, knt)    ! Sedimentation 
     
    217217#if defined key_ligand 
    218218      NAMELIST/nampisbio/ nrdttrc, wsbio, xkmort, ferat3, no3rat3, po4rat3, wsbio2,    & 
    219       &                   wsbio2max, wsbio2scale, niter1max, niter2max, wfep 
     219      &                   wsbio2max, wsbio2scale, niter1max, niter2max, wfep, ldocp, ldocz, lthet 
    220220#else 
    221221      NAMELIST/nampisbio/ nrdttrc, wsbio, xkmort, ferat3, no3rat3, po4rat3, wsbio2,    & 
Note: See TracChangeset for help on using the changeset viewer.