Changeset 605 for trunk/NEMO/OPA_SRC/LDF/ldfeiv.F90
- Timestamp:
- 2007-02-20T17:12:23+01:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/LDF/ldfeiv.F90
r461 r605 124 124 ! Compute aeiw by multiplying Ro^2 and T^-1 125 125 aeiw(ji,jj) = zross(ji,jj) * zross(ji,jj) * SQRT( zah(ji,jj) / zhw(ji,jj) ) * tmask(ji,jj,1) 126 ! Take the minimum between aeiw and 1000m^2/s for depth levels 127 ! lower than 20 (21 in w- point) 128 IF( mbathy(ji,jj) <= 21. ) aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000. ) 126 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R02 127 ! Take the minimum between aeiw and aeiv0 for depth levels 128 ! lower than 20 (21 in w- point) 129 IF( mbathy(ji,jj) <= 21. ) aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000. ) 130 ENDIF 129 131 END DO 130 132 … … 135 137 END DO 136 138 137 ! ORCA R05: Take the minimum between aeiw and 1000m2/s139 ! ORCA R05: Take the minimum between aeiw and aeiv0 138 140 IF( cp_cfg == "orca" .AND. jp_cfg == 05 ) THEN ! ORCA R05 139 141 DO ji = 2, jpim1 140 aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000.)142 aeiw(ji,jj) = MIN( aeiw(ji,jj), aeiv0 ) 141 143 END DO 142 144 ENDIF … … 176 178 zaht = ( 1. - MIN( 1., ABS( ff(ji,jj) / zf20 ) ) ) * ( aht0 - zaht_min ) & 177 179 & + aht0 * upsrnfh(ji,jj) ! enhanced near river mouths 178 ahtu(ji,jj) = MAX( zaht_min, aeiu(ji,jj) ) + zaht179 ahtv(ji,jj) = MAX( zaht_min, aeiv(ji,jj) ) + zaht180 ahtw(ji,jj) = MAX( zaht_min, aeiw(ji,jj) ) + zaht180 ahtu(ji,jj) = MAX( MAX( zaht_min, aeiu(ji,jj) ) + zaht, aht0 ) 181 ahtv(ji,jj) = MAX( MAX( zaht_min, aeiv(ji,jj) ) + zaht, aht0 ) 182 ahtw(ji,jj) = MAX( MAX( zaht_min, aeiw(ji,jj) ) + zaht, aht0 ) 181 183 END DO 182 184 END DO … … 298 300 ! Compute aeiw by multiplying Ro^2 and T^-1 299 301 aeiw(ji,jj) = zross(ji,jj) * zross(ji,jj) * SQRT( zah(ji,jj) / zhw(ji,jj) ) * tmask(ji,jj,1) 300 ! Take the minimum between aeiw and 1000m^2/s for depth levels 301 ! lower than 20 (21 in w- point) 302 IF( mbathy(ji,jj) <= 21. ) aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000. ) 302 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R02 303 ! Take the minimum between aeiw and aeiv0 for depth levels 304 ! lower than 20 (21 in w- point) 305 IF( mbathy(ji,jj) <= 21. ) aeiw(ji,jj) = MIN( aeiw(ji,jj), 1000. ) 306 ENDIF 303 307 END DO 304 308 END DO 305 309 306 310 ! Decrease the coefficient in the tropics (20N-20S) 307 311 zf20 = 2. * omega * sin( rad * 20. ) 308 312 DO jj = 2, jpjm1 309 313 DO ji = fs_2, fs_jpim1 ! vector opt. … … 312 316 END DO 313 317 314 ! ORCA R05: Take the minimum between aeiw and 1000m2/s318 ! ORCA R05: Take the minimum between aeiw and aeiv0 315 319 IF( cp_cfg == "orca" .AND. jp_cfg == 05 ) THEN 316 320 DO jj = 2, jpjm1
Note: See TracChangeset
for help on using the changeset viewer.