Changeset 713 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2007-10-11T11:43:16+02:00 (16 years ago)
- Location:
- trunk/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90
r702 r713 18 18 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qns !: sea heat flux: non solar [W/m2] 19 19 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qsr !: sea heat flux: solar [W/m2] 20 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qrp !: heat flux damping [w/m2] 20 21 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: emp !: freshwater budget: volume flux [Kg/m2/s] 21 22 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: emps !: freshwater budget: concentration/dillution [Kg/m2/s] 22 23 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: erp !: evaporation damping [kg/m2/s] 24 23 25 !!---------------------------------------------------------------------- 24 26 !! Sea Surface Mean fields -
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.