Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/friction_loc.f90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/friction_loc.f90
r5136 r5158 67 67 IF (pole_sud) jje = jj_end 68 68 69 doj = jjb, jje70 doi = 1, iip169 DO j = jjb, jje 70 DO i = 1, iip1 71 71 u2(i, j) = ucov(i, j, 1) * ucov(i, j, 1) * unscu2(i, j) 72 72 enddo … … 78 78 IF (pole_sud) jje = jj_end - 1 79 79 80 doj = jjb, jje81 doi = 1, iip180 DO j = jjb, jje 81 DO i = 1, iip1 82 82 v2(i, j) = vcov(i, j, 1) * vcov(i, j, 1) * unscv2(i, j) 83 83 enddo … … 90 90 IF (pole_sud) jje = jj_end - 1 91 91 92 doj = jjb, jje93 doi = 2, iip192 DO j = jjb, jje 93 DO i = 2, iip1 94 94 modv(i, j) = sqrt(0.5 * (u2(i - 1, j) + u2(i, j) + v2(i, j - 1) + v2(i, j))) 95 95 enddo … … 103 103 vpoln = 0. 104 104 105 doi = 2, iip1105 DO i = 2, iip1 106 106 zco = cos(rlonv(i)) * (rlonu(i) - rlonu(i - 1)) 107 107 zsi = sin(rlonv(i)) * (rlonu(i) - rlonu(i - 1)) … … 111 111 enddo 112 112 vpn = sqrt(upoln * upoln + vpoln * vpoln) / pi 113 doi = 1, iip1113 DO i = 1, iip1 114 114 ! modv(i,1)=vpn 115 115 modv(i, 1) = modv(i, 2) … … 121 121 upols = 0. 122 122 vpols = 0. 123 doi = 2, iip1123 DO i = 2, iip1 124 124 zco = cos(rlonv(i)) * (rlonu(i) - rlonu(i - 1)) 125 125 zsi = sin(rlonv(i)) * (rlonu(i) - rlonu(i - 1)) … … 129 129 enddo 130 130 vps = sqrt(upols * upols + vpols * vpols) / pi 131 doi = 1, iip1131 DO i = 1, iip1 132 132 ! modv(i,jjp1)=vps 133 133 modv(i, jjp1) = modv(i, jjm) … … 143 143 IF (pole_sud) jje = jj_end - 1 144 144 145 doj = jjb, jje146 doi = 1, iim145 DO j = jjb, jje 146 DO i = 1, iim 147 147 ucov(i, j, 1) = ucov(i, j, 1) & 148 148 - cfric * pdt * 0.5 * (modv(i + 1, j) + modv(i, j)) * ucov(i, j, 1) … … 155 155 IF (pole_sud) jje = jj_end - 1 156 156 157 doj = jjb, jje158 doi = 1, iip1157 DO j = jjb, jje 158 DO i = 1, iip1 159 159 vcov(i, j, 1) = vcov(i, j, 1) & 160 160 - cfric * pdt * 0.5 * (modv(i, j + 1) + modv(i, j)) * vcov(i, j, 1) … … 173 173 174 174 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 175 dol = 1, llm175 DO l = 1, llm 176 176 ucov(1:iip1, jjb:jje, l) = ucov(1:iip1, jjb:jje, l) * & 177 177 (1. - pdt * kfrict(l)) … … 185 185 186 186 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 187 dol = 1, llm187 DO l = 1, llm 188 188 vcov(1:iip1, jjb:jje, l) = vcov(1:iip1, jjb:jje, l) * & 189 189 (1. - pdt * kfrict(l))
Note: See TracChangeset
for help on using the changeset viewer.