- Timestamp:
- 2017-12-26T17:32:56+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r9168 r9169 1272 1272 CASE( np_teos10 ) !== polynomial TEOS-10 ==! 1273 1273 IF(lwp) WRITE(numout,*) 1274 IF(lwp) WRITE(numout,*) ' 1274 IF(lwp) WRITE(numout,*) ' ==>>> use of TEOS-10 equation of state (cons. temp. and abs. salinity)' 1275 1275 ! 1276 1276 l_useCT = .TRUE. ! model temperature is Conservative temperature … … 1464 1464 ! 1465 1465 IF(lwp) WRITE(numout,*) 1466 IF(lwp) WRITE(numout,*) ' 1466 IF(lwp) WRITE(numout,*) ' ==>>> use of EOS-80 equation of state (pot. temp. and pract. salinity)' 1467 1467 ! 1468 1468 l_useCT = .FALSE. ! model temperature is Potential temperature … … 1655 1655 IF(lwp) THEN 1656 1656 WRITE(numout,*) 1657 WRITE(numout,*) ' use of simplified eos: rhd(dT=T-10,dS=S-35,Z) = ' 1658 WRITE(numout,*) ' [-a0*(1+lambda1/2*dT+mu1*Z)*dT + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS]/rau0' 1659 WRITE(numout,*) 1660 WRITE(numout,*) ' thermal exp. coef. rn_a0 = ', rn_a0 1661 WRITE(numout,*) ' saline cont. coef. rn_b0 = ', rn_b0 1662 WRITE(numout,*) ' cabbeling coef. rn_lambda1 = ', rn_lambda1 1663 WRITE(numout,*) ' cabbeling coef. rn_lambda2 = ', rn_lambda2 1664 WRITE(numout,*) ' thermobar. coef. rn_mu1 = ', rn_mu1 1665 WRITE(numout,*) ' thermobar. coef. rn_mu2 = ', rn_mu2 1666 WRITE(numout,*) ' 2nd cabbel. coef. rn_nu = ', rn_nu 1667 WRITE(numout,*) ' Caution: rn_beta0=0 incompatible with ddm parameterization ' 1657 WRITE(numout,*) ' ==>>> use of simplified eos: ' 1658 WRITE(numout,*) ' rhd(dT=T-10,dS=S-35,Z) = [-a0*(1+lambda1/2*dT+mu1*Z)*dT ' 1659 WRITE(numout,*) ' + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS] / rau0' 1660 WRITE(numout,*) ' with the following coefficients :' 1661 WRITE(numout,*) ' thermal exp. coef. rn_a0 = ', rn_a0 1662 WRITE(numout,*) ' saline cont. coef. rn_b0 = ', rn_b0 1663 WRITE(numout,*) ' cabbeling coef. rn_lambda1 = ', rn_lambda1 1664 WRITE(numout,*) ' cabbeling coef. rn_lambda2 = ', rn_lambda2 1665 WRITE(numout,*) ' thermobar. coef. rn_mu1 = ', rn_mu1 1666 WRITE(numout,*) ' thermobar. coef. rn_mu2 = ', rn_mu2 1667 WRITE(numout,*) ' 2nd cabbel. coef. rn_nu = ', rn_nu 1668 WRITE(numout,*) ' Caution: rn_beta0=0 incompatible with ddm parameterization ' 1668 1669 ENDIF 1669 1670 l_useCT = .TRUE. ! Use conservative temperature … … 1682 1683 IF(lwp) THEN 1683 1684 IF( l_useCT ) THEN 1684 WRITE(numout,*) ' model uses Conservative Temperature' 1685 WRITE(numout,*) ' Important: model must be initialized with CT and SA fields' 1685 WRITE(numout,*) 1686 WRITE(numout,*) ' ==>>> model uses Conservative Temperature' 1687 WRITE(numout,*) ' Important: model must be initialized with CT and SA fields' 1686 1688 ELSE 1687 WRITE(numout,*) ' model does not use Conservative Temperature' 1689 WRITE(numout,*) 1690 WRITE(numout,*) ' ==>>> model does not use Conservative Temperature' 1688 1691 ENDIF 1689 1692 ENDIF 1690 1693 ! 1691 1694 IF(lwp) WRITE(numout,*) 1692 IF(lwp) WRITE(numout,*) ' volumic mass of reference rau0 = ', rau0 , ' kg/m^3' 1693 IF(lwp) WRITE(numout,*) ' 1. / rau0 r1_rau0 = ', r1_rau0, ' m^3/kg' 1694 IF(lwp) WRITE(numout,*) ' ocean specific heat rcp = ', rcp , ' J/Kelvin' 1695 IF(lwp) WRITE(numout,*) ' rau0 * rcp rau0_rcp = ', rau0_rcp 1696 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) r1_rau0_rcp = ', r1_rau0_rcp 1695 IF(lwp) WRITE(numout,*) ' Associated physical constant' 1696 IF(lwp) WRITE(numout,*) ' volumic mass of reference rau0 = ', rau0 , ' kg/m^3' 1697 IF(lwp) WRITE(numout,*) ' 1. / rau0 r1_rau0 = ', r1_rau0, ' m^3/kg' 1698 IF(lwp) WRITE(numout,*) ' ocean specific heat rcp = ', rcp , ' J/Kelvin' 1699 IF(lwp) WRITE(numout,*) ' rau0 * rcp rau0_rcp = ', rau0_rcp 1700 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) r1_rau0_rcp = ', r1_rau0_rcp 1697 1701 ! 1698 1702 END SUBROUTINE eos_init
Note: See TracChangeset
for help on using the changeset viewer.