Changeset 2113 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA
- Timestamp:
- 2010-09-22T21:20:46+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/trasbc.F90
r2104 r2113 152 152 DO jj = 1, jpj 153 153 DO ji = 1, jpi 154 rnf_dep(ji,jj) = 0.e0155 DO jk = 1, rnf_mod_dep(ji,jj) ! recalculates rnf_dep to be the depth156 rnf_dep(ji,jj) = rnf_dep(ji,jj) + fse3t(ji,jj,jk) ! in metres to the bottom of the relevant grid box157 END DO158 154 zdep = 1. / rnf_dep(ji,jj) 159 155 zse3t= 1. / fse3t(ji,jj,1) 160 IF( rnf_t mp(ji,jj) == -999 ) rnf_tmp(ji,jj) = tsn(ji,jj,1,jp_tem) ! if not specified set runoff temp to be sst156 IF( rnf_tem(ji,jj) == -999 ) rnf_tem(ji,jj) = tsn(ji,jj,1,jp_tem) ! if not specified set runoff temp to be sst 161 157 162 158 IF( rnf(ji,jj) > 0.e0 ) THEN … … 181 177 182 178 ! direct flux 183 zta = rnf_t mp(ji,jj) * rnf(ji,jj) * zsrau * zdep179 zta = rnf_tem(ji,jj) * rnf(ji,jj) * zsrau * zdep 184 180 zsa = rnf_sal(ji,jj) * rnf(ji,jj) * zsrau * zdep 185 181 … … 190 186 ELSE 191 187 DO jk = 1, rnf_mod_dep(ji,jj) 192 zta = ( rnf_t mp(ji,jj) - tsn(ji,jj,jk,jp_tem) ) * rnf(ji,jj) * zsrau * zdep188 zta = ( rnf_tem(ji,jj) - tsn(ji,jj,jk,jp_tem) ) * rnf(ji,jj) * zsrau * zdep 193 189 zsa = ( rnf_sal(ji,jj) - tsn(ji,jj,jk,jp_sal) ) * rnf(ji,jj) * zsrau * zdep 194 190 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) + zta ! add the trend to the general tracer trend
Note: See TracChangeset
for help on using the changeset viewer.