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 14325 – NEMO

Changeset 14325


Ignore:
Timestamp:
2021-01-21T09:24:32+01:00 (4 years ago)
Author:
aumont
Message:

bug in the CacO3 dissolution scheme

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedinorg.F90

    r14323 r14325  
    163163            zco3eq(ji)     = MAX( rtrn, zco3eq(ji) )  
    164164            zundsat(ji,jk) = ( zco3eq(ji) - co3por(ji,jk) ) / zco3eq(ji) 
     165            saturco3(ji,jk) = zundsat(ji,jk) 
    165166         ENDDO 
    166167      ENDDO 
    … …  
    172173               zbeta = 1.0 - reac_cal * dtsed2 * zundsat(ji,jk) + reac_cal * dtsed2 * zsolid1 
    173174               zgamma = -reac_cal * dtsed2 * zundsat(ji,jk) 
    174                zundsat(ji,jk) = ( -zbeta + SQRT( zbeta**2 - 4.0 * zgamma ) ) / ( 2.0 * reac_sil * znusil * dtsed2 ) 
     175               zundsat(ji,jk) = ( -zbeta + SQRT( zbeta**2 - 4.0 * zgamma ) ) / ( 2.0 * reac_cal * dtsed2 ) 
    175176               saturco3(ji,jk) = zundsat(ji,jk) 
    176177               zreasat = reac_cal * dtsed2 * zundsat(ji,jk) * zsolid1 / ( 1. + reac_cal * dtsed2 * zundsat(ji,jk) ) 
Note: See TracChangeset for help on using the changeset viewer.