- Timestamp:
- 2013-11-26T17:47:17+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn.F90
r4292 r4317 97 97 98 98 IF (lk_vvl) THEN 99 phur 1(:,:) = 0.100 phvr 1(:,:) = 0.99 phur(:,:) = 0. 100 phvr(:,:) = 0. 101 101 DO jk = 1, jpkm1 102 phur 1(:,:) = phur1(:,:) + fse3u_a(:,:,jk) * umask(:,:,jk)103 phvr 1(:,:) = phvr1(:,:) + fse3v_a(:,:,jk) * vmask(:,:,jk)102 phur(:,:) = phur(:,:) + fse3u_a(:,:,jk) * umask(:,:,jk) 103 phvr(:,:) = phvr(:,:) + fse3v_a(:,:,jk) * vmask(:,:,jk) 104 104 pua2d(:,:) = pua2d(:,:) + fse3u_a(:,:,jk) * umask(:,:,jk) * ua(:,:,jk) 105 105 pva2d(:,:) = pva2d(:,:) + fse3v_a(:,:,jk) * vmask(:,:,jk) * va(:,:,jk) 106 106 END DO 107 phur 1(:,:) = umask(:,:,1) / ( phur1(:,:) + 1. - umask(:,:,1) )108 phvr 1(:,:) = vmask(:,:,1) / ( phvr1(:,:) + 1. - vmask(:,:,1) )109 pua2d(:,:) = pua2d(:,:) * phur 1(:,:)110 pva2d(:,:) = pva2d(:,:) * phvr 1(:,:)107 phur(:,:) = umask(:,:,1) / ( phur(:,:) + 1. - umask(:,:,1) ) 108 phvr(:,:) = vmask(:,:,1) / ( phvr(:,:) + 1. - vmask(:,:,1) ) 109 pua2d(:,:) = pua2d(:,:) * phur(:,:) 110 pva2d(:,:) = pva2d(:,:) * phvr(:,:) 111 111 ELSE 112 112 DO jk = 1, jpkm1 … … 119 119 120 120 DO jk = 1 , jpkm1 121 ua(:,:,jk) = ua(:,:,jk) - pua2d(:,:)122 va(:,:,jk) = va(:,:,jk) - pva2d(:,:)121 ua(:,:,jk) = (ua(:,:,jk) - pua2d(:,:)) * umask(:,:,jk) 122 va(:,:,jk) = (va(:,:,jk) - pva2d(:,:)) * vmask(:,:,jk) 123 123 END DO 124 124 … … 152 152 153 153 DO jk = 1 , jpkm1 154 ub(:,:,jk) = ub(:,:,jk) - pub2d(:,:)155 vb(:,:,jk) = vb(:,:,jk) - pvb2d(:,:)154 ub(:,:,jk) = (ub(:,:,jk) - pub2d(:,:)) * umask(:,:,jk) 155 vb(:,:,jk) = (vb(:,:,jk) - pvb2d(:,:)) * vmask(:,:,jk) 156 156 END DO 157 157 END IF … … 183 183 184 184 CALL wrk_dealloc(jpi,jpj,pua2d,pva2d) 185 IF ( ll_orlanski ) CALL wrk_dealloc(jpi,jpj,pub2d,pvb2d )185 IF ( ll_orlanski ) CALL wrk_dealloc(jpi,jpj,pub2d,pvb2d,phur1,phvr1) 186 186 187 187 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn')
Note: See TracChangeset
for help on using the changeset viewer.