Changeset 7646 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r6505 r7646 204 204 !! 205 205 !! ln_teos10 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 206 !! Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Cel cius, sa=35.5 g/kg206 !! Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Celsius, sa=35.5 g/kg 207 207 !! 208 208 !! ln_eos80 : polynomial EOS-80 equation of state is used for rho(t,s,z). 209 !! Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Cel cius, sp=35.5 psu209 !! Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Celsius, sp=35.5 psu 210 210 !! 211 211 !! ln_seos : simplified equation of state … … 221 221 !! TEOS-10 Manual, 2010 222 222 !!---------------------------------------------------------------------- 223 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]223 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 224 224 ! ! 2 : salinity [psu] 225 225 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: prd ! in situ density [-] … … 316 316 !! 317 317 !!---------------------------------------------------------------------- 318 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]318 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 319 319 ! ! 2 : salinity [psu] 320 320 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: prd ! in situ density [-] … … 481 481 !! 482 482 !!---------------------------------------------------------------------- 483 REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]483 REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 484 484 ! ! 2 : salinity [psu] 485 485 REAL(wp), DIMENSION(jpi,jpj) , INTENT(in ) :: pdep ! depth [m] … … 907 907 !! 908 908 !!---------------------------------------------------------------------- 909 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! pot. temperature and salinity [Cel cius,psu]910 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pab ! thermal/haline expansion coef. [Cel cius-1,psu-1]909 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! pot. temperature and salinity [Celsius,psu] 910 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pab ! thermal/haline expansion coef. [Celsius-1,psu-1] 911 911 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: pn2 ! Brunt-Vaisala frequency squared [1/s^2] 912 912 ! … … 944 944 !! *** ROUTINE eos_pt_from_ct *** 945 945 !! 946 !! ** Purpose : Compute pot.temp. from cons. temp. [Cel cius]946 !! ** Purpose : Compute pot.temp. from cons. temp. [Celsius] 947 947 !! 948 948 !! ** Method : rational approximation (5/3th order) of TEOS-10 algorithm … … 952 952 !! Rational approximation to TEOS10 algorithm (rms error on WOA13 values: 4.0e-5 degC) 953 953 !!---------------------------------------------------------------------- 954 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ctmp ! Cons. Temp [Celcius]955 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu]954 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ctmp ! Cons. Temp [Celsius] 955 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 956 956 ! Leave result array automatic rather than making explicitly allocated 957 REAL(wp), DIMENSION(jpi,jpj) :: ptmp ! potential temperature [Cel cius]957 REAL(wp), DIMENSION(jpi,jpj) :: ptmp ! potential temperature [Celsius] 958 958 ! 959 959 INTEGER :: ji, jj ! dummy loop indices … … 1003 1003 !! *** ROUTINE eos_fzp *** 1004 1004 !! 1005 !! ** Purpose : Compute the freezing point temperature [Cel cius]1006 !! 1007 !! ** Method : UNESCO freezing point (ptf) in Cel cius is given by1005 !! ** Purpose : Compute the freezing point temperature [Celsius] 1006 !! 1007 !! ** Method : UNESCO freezing point (ptf) in Celsius is given by 1008 1008 !! ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 1009 1009 !! checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m … … 1013 1013 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 1014 1014 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ), OPTIONAL :: pdep ! depth [m] 1015 REAL(wp), DIMENSION(jpi,jpj), INTENT(out ) :: ptf ! freezing temperature [Cel cius]1015 REAL(wp), DIMENSION(jpi,jpj), INTENT(out ) :: ptf ! freezing temperature [Celsius] 1016 1016 ! 1017 1017 INTEGER :: ji, jj ! dummy loop indices … … 1056 1056 !! *** ROUTINE eos_fzp *** 1057 1057 !! 1058 !! ** Purpose : Compute the freezing point temperature [Cel cius]1059 !! 1060 !! ** Method : UNESCO freezing point (ptf) in Cel cius is given by1058 !! ** Purpose : Compute the freezing point temperature [Celsius] 1059 !! 1060 !! ** Method : UNESCO freezing point (ptf) in Celsius is given by 1061 1061 !! ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 1062 1062 !! checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m … … 1066 1066 REAL(wp), INTENT(in ) :: psal ! salinity [psu] 1067 1067 REAL(wp), INTENT(in ), OPTIONAL :: pdep ! depth [m] 1068 REAL(wp), INTENT(out) :: ptf ! freezing temperature [Cel cius]1068 REAL(wp), INTENT(out) :: ptf ! freezing temperature [Celsius] 1069 1069 ! 1070 1070 REAL(wp) :: zs ! local scalars
Note: See TracChangeset
for help on using the changeset viewer.