Ignore:
Timestamp:
2018-11-22T12:39:19+01:00 (2 years ago)
Author:
cetlod
Message:

Bugfix in sediment model, see ticket #2171

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90

    r10222 r10356  
    678678         &         - 0.07711*zsal + 0.0041249*zsal15 
    679679 
     680         ! CONVERT FROM DIFFERENT PH SCALES 
     681         total2free  = 1.0/(1.0 + zst/zcks) 
     682         free2SWS    = 1. + zst/zcks + zft/(zckf*total2free) 
     683         total2SWS   = total2free * free2SWS 
     684         SWS2total   = 1.0 / total2SWS 
     685 
     686 
    680687         ! K1, K2 OF CARBONIC ACID, KB OF BORIC ACID, KW (H2O) (LIT.?) 
    681688         zak1    = 10**(zck1) * total2SWS 
     
    743750         aksis(ji) = zaksi * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 ) 
    744751 
    745          ! CONVERT FROM DIFFERENT PH SCALES 
    746          total2free  = 1.0/(1.0 + zst/aks3s(ji)) 
    747          free2SWS    = 1. + zst/aks3s(ji) + zft/akf3s(ji) 
    748          total2SWS   = total2free * free2SWS 
    749          SWS2total   = 1.0 / total2SWS 
    750  
    751752         ! Convert to total scale 
    752753         ak1s(ji)  = ak1s(ji)  * SWS2total 
Note: See TracChangeset for help on using the changeset viewer.