- Timestamp:
- 2019-12-05T18:41:39+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/ISF/isfcavgam.F90
r12068 r12077 10 10 !! isfcav_gammats : compute exchange coeficient gamma 11 11 !!---------------------------------------------------------------------- 12 USE isf 12 USE isf_oce 13 13 USE isfutils, ONLY: debug 14 14 USE isftbl , ONLY: isf_tbl … … 46 46 !! 47 47 !! ** Method : select the gamma formulation 48 !! 3 method available (cst, AD15 and HJ99)48 !! 3 method available (cst, vel and vel_stab) 49 49 !!--------------------------------------------------------------------- 50 50 !!-------------------------- OUT ------------------------------------- … … 66 66 ! gamma is constant (specified in namelist) 67 67 ! nothing to do 68 CASE (' ad15', 'hj99')68 CASE ('vel', 'vel_stab') 69 69 ! compute velocity in tbl 70 70 CALL isf_tbl(Kmm, uu(:,:,:,Kmm) ,zutbl(:,:),'U', miku, rhisf_tbl_cav) … … 90 90 pgt(:,:) = rn_gammat0 91 91 pgs(:,:) = rn_gammas0 92 CASE ( ' ad15' ) ! gamma is proportional to u*93 CALL gammats_ AD15( zutbl, zvtbl, rCd0_top, r_ke0_top, pgt, pgs )94 CASE ( ' hj99' ) ! gamma depends of stability of boundary layer and u*95 CALL gammats_ HJ99(Kmm, pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs )92 CASE ( 'vel' ) ! gamma is proportional to u* 93 CALL gammats_vel ( zutbl, zvtbl, rCd0_top, r_ke0_top, pgt, pgs ) 94 CASE ( 'vel_stab' ) ! gamma depends of stability of boundary layer and u* 95 CALL gammats_vel_stab (Kmm, pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 96 96 CASE DEFAULT 97 97 CALL ctl_stop('STOP','method to compute gamma (cn_gammablk) is unknown (should not see this)') … … 116 116 !!----------------------------------------------------------------------------------------------------- 117 117 ! 118 SUBROUTINE gammats_ AD15(putbl, pvtbl, pCd, pke2, & ! <<== in118 SUBROUTINE gammats_vel( putbl, pvtbl, pCd, pke2, & ! <<== in 119 119 & pgt, pgs ) ! ==>> out gammats [m/s] 120 120 !!---------------------------------------------------------------------- … … 123 123 !! ** Method : gamma is velocity dependent ( gt= gt0 * Ustar ) 124 124 !! 125 !! ** Reference : Jenkins et al., 2010, JPO, p2298-2312 126 !! Asay-Davis et al., Geosci. Model Dev., 9, 2471-2497, 2016 125 !! ** Reference : Asay-Davis et al., Geosci. Model Dev., 9, 2471-2497, 2016 127 126 !!--------------------------------------------------------------------- 128 127 !!-------------------------- OUT ------------------------------------- … … 146 145 CALL iom_put('isfustar',zustar(:,:)) 147 146 ! 148 END SUBROUTINE gammats_ AD15149 150 SUBROUTINE gammats_ HJ99( Kmm, pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, & ! <<== in151 & pgt , pgs ) ! ==>> out gammats [m/s]147 END SUBROUTINE gammats_vel 148 149 SUBROUTINE gammats_vel_stab( Kmm, pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, & ! <<== in 150 & pgt , pgs ) ! ==>> out gammats [m/s] 152 151 !!---------------------------------------------------------------------- 153 152 !! ** Purpose : compute the coefficient echange coefficient … … 186 185 ! 187 186 ! compute ustar 188 zustar(:,:) = SQRT( pCd * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + r_ke0_top) )187 zustar(:,:) = SQRT( pCd * ( putbl(:,:) * putbl(:,:) + pvtbl(:,:) * pvtbl(:,:) + pke2 ) ) 189 188 ! 190 189 ! output ustar … … 254 253 END DO 255 254 256 END SUBROUTINE gammats_ HJ99255 END SUBROUTINE gammats_vel_stab 257 256 258 257 END MODULE isfcavgam
Note: See TracChangeset
for help on using the changeset viewer.