Changeset 713 for trunk/NEMO/OPA_SRC/SBC/sbcssr.F90
- Timestamp:
- 2007-10-11T11:43:16+02:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcssr.F90
r702 r713 65 65 !! 66 66 INTEGER :: ji, jj ! dummy loop indices 67 REAL(wp) :: zerp ! local scalar 67 REAL(wp) :: zerp ! local scalar for evaporation damping 68 REAL(wp) :: zqrp ! local scalar for heat flux damping 68 69 INTEGER :: ierror ! return error code 69 70 !! … … 149 150 ENDIF 150 151 ! 152 ! Initialize qrp and erp if no restoring 153 IF( nn_sst /= 1 ) qrp(:,:) = 0.e0 154 IF( nn_sss /= 1 .OR. nn_sss /= 2 ) erp(:,:) = 0.e0 151 155 ENDIF 152 156 … … 164 168 IF( nn_sst == 1 ) THEN ! Temperature restoring term 165 169 !CDIR COLLAPSE 166 qns(:,:) = qns(:,:) + dqdt * ( sst_m(:,:) - sf_sst(1)%fnow(:,:) ) 170 ! use zqrp scalar to optimize memory access (speedup the loop) 171 DO jj = 1, jpj 172 DO ji = 1, jpi 173 zqrp = dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj) ) 174 qns(ji,jj) = qns(ji,jj) + zqrp 175 qrp(ji,jj) = zqrp 176 END DO 177 END DO 167 178 ENDIF 168 179 ! 169 180 IF( nn_sss == 1 ) THEN ! Salinity damping term (salt flux, emps only) 170 181 !CDIR COLLAPSE 182 ! use zerp scalar to optimize memory access (speedup the loop) 171 183 DO jj = 1, jpj 172 184 DO ji = 1, jpi … … 175 187 & / ( sss_m(ji,jj) + 1.e-20 ) 176 188 emps(ji,jj) = emps(ji,jj) + zerp 189 erp( ji,jj) = zerp 177 190 END DO 178 191 END DO 179 192 ELSEIF( nn_sss == 2 ) THEN ! Salinity damping term (volume flux, emp and emps) 180 193 !CDIR COLLAPSE 194 ! use zerp scalar to optimize memory access (speedup the loop) 181 195 DO jj = 1, jpj 182 196 DO ji = 1, jpi … … 186 200 emp (ji,jj) = emp (ji,jj) + zerp 187 201 emps(ji,jj) = emps(ji,jj) + zerp 202 erp (ji,jj) = zerp 188 203 END DO 189 204 END DO
Note: See TracChangeset
for help on using the changeset viewer.