Changeset 14072 for NEMO/trunk/src/OCE/DIU/diu_coolskin.F90
- Timestamp:
- 2020-12-04T08:48:38+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DIU/diu_coolskin.F90
r13722 r14072 20 20 USE lib_mpp 21 21 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 22 22 23 23 IMPLICIT NONE 24 24 PRIVATE … … 34 34 REAL(wp), PRIVATE, PARAMETER :: pp_rhoa = 1.20421_wp ! density of air (at 20C) 35 35 REAL(wp), PRIVATE, PARAMETER :: pp_cda = 1.45e-3_wp ! assumed air-sea drag coefficient for calculating wind speed 36 36 37 37 ! Key variables 38 38 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: x_csdsst ! Cool skin delta SST … … 46 46 !! $Id$ 47 47 !! Software governed by the CeCILL license (see ./LICENSE) 48 !!---------------------------------------------------------------------- 49 CONTAINS 50 48 !!---------------------------------------------------------------------- 49 CONTAINS 50 51 51 SUBROUTINE diurnal_sst_coolskin_init 52 52 !!---------------------------------------------------------------------- … … 55 55 !! ** Purpose : initialise the cool skin model 56 56 !! 57 !! ** Method : 57 !! ** Method : 58 58 !! 59 59 !! ** Reference : 60 !! 60 !! 61 61 !!---------------------------------------------------------------------- 62 62 ALLOCATE( x_csdsst(jpi,jpj), x_csthick(jpi,jpj) ) … … 73 73 !! ** Purpose : Time-step the Artale cool skin model 74 74 !! 75 !! ** Method : 75 !! ** Method : 76 76 !! 77 !! ** Reference : 77 !! ** Reference : 78 78 !!---------------------------------------------------------------------- 79 79 ! Dummy variables … … 82 82 REAL(wp), INTENT(IN), DIMENSION(jpi,jpj) :: psrho ! Water density (kg/m^3) 83 83 REAL(wp), INTENT(IN) :: pDt ! Time-step 84 84 85 85 ! Local variables 86 REAL(wp), DIMENSION(jpi,jpj) :: z_fv ! Friction velocity 86 REAL(wp), DIMENSION(jpi,jpj) :: z_fv ! Friction velocity 87 87 REAL(wp), DIMENSION(jpi,jpj) :: z_gamma ! Dimensionless function of wind speed 88 88 REAL(wp), DIMENSION(jpi,jpj) :: z_lamda ! Sauders (dimensionless) proportionality constant … … 91 91 REAL(wp) :: z_zty ! Temporary v wind stress 92 92 REAL(wp) :: z_zmod ! Temporary total wind stress 93 93 94 94 INTEGER :: ji,jj 95 95 !!---------------------------------------------------------------------- … … 105 105 ELSE 106 106 z_fv(ji,jj) = 0. 107 z_wspd(ji,jj) = 0. 107 z_wspd(ji,jj) = 0. 108 108 ENDIF 109 109 !
Note: See TracChangeset
for help on using the changeset viewer.