Changeset 14325


Ignore:
Timestamp:
2021-01-21T09:24:32+01:00 (5 months 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.