- Timestamp:
- 2019-12-05T10:57:46+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfcavgam.F90
r11987 r12062 65 65 ! gamma is constant (specified in namelist) 66 66 ! nothing to do 67 CASE (' ad15', 'hj99')67 CASE ('vel', 'vel_stab') 68 68 ! compute velocity in tbl 69 69 CALL isf_tbl(un(:,:,:) ,zutbl(:,:),'U', miku, rhisf_tbl_cav) … … 89 89 pgt(:,:) = rn_gammat0 90 90 pgs(:,:) = rn_gammas0 91 CASE ( ' ad15' ) ! gamma is proportional to u*92 CALL gammats_ AD15 ( zutbl, zvtbl, rCd0_top, rn_vtide**2,pgt, pgs )93 CASE ( ' hj99' ) ! gamma depends of stability of boundary layer and u*94 CALL gammats_ HJ99(pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs )91 CASE ( 'vel' ) ! gamma is proportional to u* 92 CALL gammats_vel ( zutbl, zvtbl, rCd0_top, rn_vtide**2, pgt, pgs ) 93 CASE ( 'vel_stab' ) ! gamma depends of stability of boundary layer and u* 94 CALL gammats_vel_stab (pttbl, pstbl, zutbl, zvtbl, rCd0_top, r_ke0_top, pqoce, pqfwf, pgt, pgs ) 95 95 CASE DEFAULT 96 96 CALL ctl_stop('STOP','method to compute gamma (cn_gammablk) is unknown (should not see this)') … … 115 115 !!----------------------------------------------------------------------------------------------------- 116 116 ! 117 SUBROUTINE gammats_ AD15(putbl, pvtbl, pCd, pke2, & ! <<== in117 SUBROUTINE gammats_vel( putbl, pvtbl, pCd, pke2, & ! <<== in 118 118 & pgt, pgs ) ! ==>> out gammats [m/s] 119 119 !!---------------------------------------------------------------------- … … 145 145 CALL iom_put('isfustar',zustar(:,:)) 146 146 ! 147 END SUBROUTINE gammats_ AD15148 149 SUBROUTINE gammats_ HJ99( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, & ! <<== in150 & pgt , pgs ) ! ==>> out gammats [m/s]147 END SUBROUTINE gammats_vel 148 149 SUBROUTINE gammats_vel_stab( pttbl, pstbl, putbl, pvtbl, pCd, pke2, pqoce, pqfwf, & ! <<== in 150 & pgt , pgs ) ! ==>> out gammats [m/s] 151 151 !!---------------------------------------------------------------------- 152 152 !! ** Purpose : compute the coefficient echange coefficient … … 252 252 END DO 253 253 254 END SUBROUTINE gammats_ HJ99254 END SUBROUTINE gammats_vel_stab 255 255 256 256 END MODULE isfcavgam
Note: See TracChangeset
for help on using the changeset viewer.