- Timestamp:
- 2016-08-29T17:55:22+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r6140 r6900 195 195 !! 196 196 !! nn_eos = -1 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 197 !! Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Cel cius, sa=35.5 g/kg197 !! Check value: rho = 1028.21993233072 kg/m^3 for z=3000 dbar, ct=3 Celsius, sa=35.5 g/kg 198 198 !! 199 199 !! nn_eos = 0 : polynomial EOS-80 equation of state is used for rho(t,s,z). 200 !! Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Cel cius, sp=35.5 psu200 !! Check value: rho = 1028.35011066567 kg/m^3 for z=3000 dbar, pt=3 Celsius, sp=35.5 psu 201 201 !! 202 202 !! nn_eos = 1 : simplified equation of state … … 212 212 !! TEOS-10 Manual, 2010 213 213 !!---------------------------------------------------------------------- 214 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]214 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 215 215 ! ! 2 : salinity [psu] 216 216 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: prd ! in situ density [-] … … 307 307 !! 308 308 !!---------------------------------------------------------------------- 309 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]309 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 310 310 ! ! 2 : salinity [psu] 311 311 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: prd ! in situ density [-] … … 472 472 !! 473 473 !!---------------------------------------------------------------------- 474 REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Cel cius]474 REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in ) :: pts ! 1 : potential temperature [Celsius] 475 475 ! ! 2 : salinity [psu] 476 476 REAL(wp), DIMENSION(jpi,jpj) , INTENT(in ) :: pdep ! depth [m] … … 897 897 !! 898 898 !!---------------------------------------------------------------------- 899 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! pot. temperature and salinity [Cel cius,psu]900 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pab ! thermal/haline expansion coef. [Cel cius-1,psu-1]899 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pts ! pot. temperature and salinity [Celsius,psu] 900 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in ) :: pab ! thermal/haline expansion coef. [Celsius-1,psu-1] 901 901 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT( out) :: pn2 ! Brunt-Vaisala frequency squared [1/s^2] 902 902 ! … … 934 934 !! *** ROUTINE eos_pt_from_ct *** 935 935 !! 936 !! ** Purpose : Compute pot.temp. from cons. temp. [Cel cius]936 !! ** Purpose : Compute pot.temp. from cons. temp. [Celsius] 937 937 !! 938 938 !! ** Method : rational approximation (5/3th order) of TEOS-10 algorithm … … 942 942 !! Rational approximation to TEOS10 algorithm (rms error on WOA13 values: 4.0e-5 degC) 943 943 !!---------------------------------------------------------------------- 944 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ctmp ! Cons. Temp [Celcius] 944 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ctmp ! Cons. Temp 945 [Celsius] 945 946 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 946 947 ! Leave result array automatic rather than making explicitly allocated 947 REAL(wp), DIMENSION(jpi,jpj) :: ptmp ! potential temperature [Cel cius]948 REAL(wp), DIMENSION(jpi,jpj) :: ptmp ! potential temperature [Celsius] 948 949 ! 949 950 INTEGER :: ji, jj ! dummy loop indices … … 993 994 !! *** ROUTINE eos_fzp *** 994 995 !! 995 !! ** Purpose : Compute the freezing point temperature [Cel cius]996 !! 997 !! ** Method : UNESCO freezing point (ptf) in Cel cius is given by996 !! ** Purpose : Compute the freezing point temperature [Celsius] 997 !! 998 !! ** Method : UNESCO freezing point (ptf) in Celsius is given by 998 999 !! ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 999 1000 !! checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m … … 1003 1004 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 1004 1005 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ), OPTIONAL :: pdep ! depth [m] 1005 REAL(wp), DIMENSION(jpi,jpj), INTENT(out ) :: ptf ! freezing temperature [Cel cius]1006 REAL(wp), DIMENSION(jpi,jpj), INTENT(out ) :: ptf ! freezing temperature [Celsius] 1006 1007 ! 1007 1008 INTEGER :: ji, jj ! dummy loop indices … … 1044 1045 !! *** ROUTINE eos_fzp *** 1045 1046 !! 1046 !! ** Purpose : Compute the freezing point temperature [Cel cius]1047 !! 1048 !! ** Method : UNESCO freezing point (ptf) in Cel cius is given by1047 !! ** Purpose : Compute the freezing point temperature [Celsius] 1048 !! 1049 !! ** Method : UNESCO freezing point (ptf) in Celsius is given by 1049 1050 !! ptf(t,z) = (-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s - 7.53e-4*z 1050 1051 !! checkvalue: tf=-2.588567 Celsius for s=40psu, z=500m … … 1054 1055 REAL(wp), INTENT(in ) :: psal ! salinity [psu] 1055 1056 REAL(wp), INTENT(in ), OPTIONAL :: pdep ! depth [m] 1056 REAL(wp), INTENT(out) :: ptf ! freezing temperature [Cel cius]1057 REAL(wp), INTENT(out) :: ptf ! freezing temperature [Celsius] 1057 1058 ! 1058 1059 REAL(wp) :: zs ! local scalars
Note: See TracChangeset
for help on using the changeset viewer.