- Timestamp:
- 2012-04-13T15:48:38+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_3352_UKMO8_CICE/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r3294 r3356 678 678 679 679 680 FUNCTION tfreez( psal ) RESULT( ptf )680 FUNCTION tfreez( psal, ppress ) RESULT( ptf ) 681 681 !!---------------------------------------------------------------------- 682 682 !! *** ROUTINE eos_init *** … … 691 691 !!---------------------------------------------------------------------- 692 692 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 693 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ), OPTIONAL :: ppress ! pressure [dBar] 693 694 ! Leave result array automatic rather than making explicitly allocated 694 695 REAL(wp), DIMENSION(jpi,jpj) :: ptf ! freezing temperature [Celcius] 696 REAL(wp), DIMENSION(jpi,jpj) :: zpress ! pressure [dBar] 695 697 !!---------------------------------------------------------------------- 696 698 ! 697 699 ptf(:,:) = ( - 0.0575_wp + 1.710523e-3_wp * SQRT( psal(:,:) ) & 698 700 & - 2.154996e-4_wp * psal(:,:) ) * psal(:,:) 701 702 IF( PRESENT(ppress) ) THEN 703 ptf(:,:) = ptf(:,:) - 7.53e-4_wp * ppress(:,:) 704 ENDIF 705 699 706 ! 700 707 END FUNCTION tfreez
Note: See TracChangeset
for help on using the changeset viewer.