Changeset 11962 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p6.F90
- Timestamp:
- 2019-11-25T23:31:07+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p6.F90
r11851 r11962 43 43 PRIVATE 44 44 45 PUBLIC :: COARE3P6_INIT, TURB_COARE3P645 PUBLIC :: SBCBLK_ALGO_COARE3P6_INIT, TURB_COARE3P6 46 46 47 47 !! COARE own values for given constants: … … 55 55 56 56 57 SUBROUTINE coare3p6_init(l_use_cs, l_use_wl)57 SUBROUTINE sbcblk_algo_coare3p6_init(l_use_cs, l_use_wl) 58 58 !!--------------------------------------------------------------------- 59 !! *** FUNCTION coare3p6_init ***59 !! *** FUNCTION sbcblk_algo_coare3p6_init *** 60 60 !! 61 61 !! INPUT : … … 71 71 ierr = 0 72 72 ALLOCATE ( Tau_ac(jpi,jpj) , Qnt_ac(jpi,jpj), dT_wl(jpi,jpj), Hz_wl(jpi,jpj), STAT=ierr ) 73 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P6_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' )73 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P6_INIT => allocation of Tau_ac, Qnt_ac, dT_wl & Hz_wl failed!' ) 74 74 Tau_ac(:,:) = 0._wp 75 75 Qnt_ac(:,:) = 0._wp … … 80 80 ierr = 0 81 81 ALLOCATE ( dT_cs(jpi,jpj), STAT=ierr ) 82 IF( ierr > 0 ) CALL ctl_stop( ' COARE3P6_INIT => allocation of dT_cs failed!' )82 IF( ierr > 0 ) CALL ctl_stop( ' SBCBLK_ALGO_COARE3P6_INIT => allocation of dT_cs failed!' ) 83 83 dT_cs(:,:) = -0.25_wp ! First guess of skin correction 84 84 END IF 85 END SUBROUTINE coare3p6_init85 END SUBROUTINE sbcblk_algo_coare3p6_init 86 86 87 87 … … 185 185 REAL(wp), DIMENSION(jpi,jpj) :: ztmp0, ztmp1, ztmp2 186 186 ! 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: & 188 & zeta_t, & ! stability parameter at height zt 189 & zsst, & ! to back up the initial bulk SST 190 & pdTc, & ! SST increment "dT" for cool-skin correction [K] 191 & pdTw, & ! SST increment "dT" for warm layer correction [K] 192 & zHwl ! depth of warm-layer [m] 187 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zeta_t ! stability parameter at height zt 188 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zsst ! to back up the initial bulk SST 189 ! 193 190 CHARACTER(len=40), PARAMETER :: crtnm = 'turb_coare3p6@sbcblk_algo_coare3p6' 194 191 !!---------------------------------------------------------------------------------- 195 IF ( kt == nit000 ) CALL COARE3P6_INIT(l_use_cs, l_use_wl)192 IF ( kt == nit000 ) CALL SBCBLK_ALGO_COARE3P6_INIT(l_use_cs, l_use_wl) 196 193 197 194 l_zt_equal_zu = .FALSE.
Note: See TracChangeset
for help on using the changeset viewer.