- Timestamp:
- 2018-11-22T12:39:19+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90
r10222 r10356 678 678 & - 0.07711*zsal + 0.0041249*zsal15 679 679 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 680 687 ! K1, K2 OF CARBONIC ACID, KB OF BORIC ACID, KW (H2O) (LIT.?) 681 688 zak1 = 10**(zck1) * total2SWS … … 743 750 aksis(ji) = zaksi * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 ) 744 751 745 ! CONVERT FROM DIFFERENT PH SCALES746 total2free = 1.0/(1.0 + zst/aks3s(ji))747 free2SWS = 1. + zst/aks3s(ji) + zft/akf3s(ji)748 total2SWS = total2free * free2SWS749 SWS2total = 1.0 / total2SWS750 751 752 ! Convert to total scale 752 753 ak1s(ji) = ak1s(ji) * SWS2total
Note: See TracChangeset
for help on using the changeset viewer.