Changeset 7359
- Timestamp:
- 2016-11-28T18:44:31+01:00 (8 years ago)
- Location:
- branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r7353 r7359 412 412 srcv(jpr_itx1)%clgrid = 'U' ! ocean stress used after its transformation 413 413 srcv(jpr_ity1)%clgrid = 'V' ! i.e. it is always at U- & V-points for i- & j-comp. resp. 414 ENDIF 414 415 ENDIF 415 416 -
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90
r7353 r7359 96 96 zk = ABS(zsp0)/MAX(ABS(5.97_wp*ztransp),0.0000001_wp) 97 97 ! Depth attenuation 98 zfac = EXP(-2.0_wp*zk* fsdept(ji,jj,jk))/(1.0_wp+8.0_wp*zk*fsdept(ji,jj,jk))98 zfac = EXP(-2.0_wp*zk*gdept_n(ji,jj,jk))/(1.0_wp+8.0_wp*zk*gdept_n(ji,jj,jk)) 99 99 ! 100 100 usd3dt(ji,jj,jk) = zfac * zusd2dt(ji,jj) * tmask(ji,jj,jk) … … 139 139 ! 140 140 DO jk = jpkm1, 1, -1 ! integrate from the bottom the e3t * hor. divergence 141 wsd3d(:,:,jk) = wsd3d(:,:,jk+1) - fse3t_n(:,:,jk) * ze3hdiv(:,:,jk)141 wsd3d(:,:,jk) = wsd3d(:,:,jk+1) - e3t_n(:,:,jk) * ze3hdiv(:,:,jk) 142 142 END DO 143 143 #if defined key_bdy -
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90
r7351 r7359 108 108 IF( ln_wave .AND. ln_sdw ) THEN 109 109 DO jk = 1, jpkm1 110 zun(:,:,jk) = e2u(:,:) * fse3u(:,:,jk) * &110 zun(:,:,jk) = e2u(:,:) * e3u_n(:,:,jk) * & 111 111 & ( un(:,:,jk) + usd3d(:,:,jk) ) ! eulerian transport + Stokes Drift 112 zvn(:,:,jk) = e1v(:,:) * fse3v(:,:,jk) * &112 zvn(:,:,jk) = e1v(:,:) * e3v_n(:,:,jk) * & 113 113 & ( vn(:,:,jk) + vsd3d(:,:,jk) ) 114 114 zwn(:,:,jk) = e1e2t(:,:) * & … … 117 117 ELSE 118 118 DO jk = 1, jpkm1 119 zun(:,:,jk) = e2u (:,:) * fse3u(:,:,jk) * un(:,:,jk) ! eulerian transport only120 zvn(:,:,jk) = e1v (:,:) * fse3v(:,:,jk) * vn(:,:,jk)119 zun(:,:,jk) = e2u (:,:) * e3u_n(:,:,jk) * un(:,:,jk) ! eulerian transport only 120 zvn(:,:,jk) = e1v (:,:) * e3v_n(:,:,jk) * vn(:,:,jk) 121 121 zwn(:,:,jk) = e1e2t(:,:) * wn(:,:,jk) 122 122 END DO -
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfqiao.F90
r7353 r7359 70 70 DO jj = 1, jpjm1 71 71 DO ji = 1, fs_jpim1 72 qbv(ji,jj,jk) = 1.0 * 0.353553 * swh(ji,jj) * tsd2d(ji,jj) * &73 & EXP(3.0 * wnum(ji,jj) * &74 & (-MIN( fsdepw(ji ,jj ,jk), fsdepw(ji+1,jj ,jk), &75 & fsdepw(ji ,jj+1,jk), fsdepw(ji+1,jj+1,jk)))) &72 qbv(ji,jj,jk) = 1.0 * 0.353553 * swh(ji,jj) * tsd2d(ji,jj) * & 73 & EXP(3.0 * wnum(ji,jj) * & 74 & (-MIN( gdepw_n(ji ,jj ,jk), gdepw_n(ji+1,jj ,jk), & 75 & gdepw_n(ji ,jj+1,jk), gdepw_n(ji+1,jj+1,jk)))) & 76 76 & * wmask(ji,jj,jk) 77 77 END DO -
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90
r7351 r7359 205 205 DO jj = 2, jpjm1 206 206 DO ji = fs_2, fs_jpim1 207 IF( fsdept(ji,jj,jk) < ekm_dep(ji,jj) ) THEN207 IF( gdepw_n(ji,jj,jk) < ekm_dep(ji,jj) ) THEN 208 208 avmv(ji,jj,jk) = MAX( avmv(ji,jj,jk), rn_wvmix ) 209 209 avmu(ji,jj,jk) = MAX( avmu(ji,jj,jk), rn_wvmix )
Note: See TracChangeset
for help on using the changeset viewer.