Changeset 7753 for trunk/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_sbc.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_sbc.F90
r7698 r7753 109 109 ztrp= - 40.e0 ! retroaction term on heat fluxes (W/m2/K) 110 110 zconv = 3.16e-5 ! convertion factor: 1 m/yr => 3.16e-5 mm/s 111 !$OMP PARALLEL DO schedule(static) private(jj, ji, t_star)112 111 DO jj = 1, jpj 113 112 DO ji = 1, jpi … … 138 137 139 138 ! freshwater (mass flux) and update of qns with heat content of emp 140 !$OMP PARALLEL DO schedule(static) private(jj, ji) 141 DO jj = 1, jpj 142 DO ji = 1, jpi 143 emp (ji,jj) = emp(ji,jj) - zsumemp * tmask(ji,jj,1) ! freshwater flux (=0 in domain average) 144 sfx (ji,jj) = 0.0_wp ! no salt flux 145 qns (ji,jj) = qns(ji,jj) - emp(ji,jj) * sst_m(ji,jj) * rcp ! evap and precip are at SST 146 END DO 147 END DO 139 emp (:,:) = emp(:,:) - zsumemp * tmask(:,:,1) ! freshwater flux (=0 in domain average) 140 sfx (:,:) = 0.0_wp ! no salt flux 141 qns (:,:) = qns(:,:) - emp(:,:) * sst_m(:,:) * rcp ! evap and precip are at SST 148 142 149 143 … … 172 166 ztau_sais = 0.015 173 167 ztaun = ztau - ztau_sais * COS( (ztime - ztimemax) / (ztimemin - ztimemax) * rpi ) 174 !$OMP PARALLEL175 !$OMP DO schedule(static) private(jj, ji)176 168 DO jj = 1, jpj 177 169 DO ji = 1, jpi … … 185 177 ! module of wind stress and wind speed at T-point 186 178 zcoef = 1. / ( zrhoa * zcdrag ) 187 !$OMP DO schedule(static) private(jj, ji, ztx, zty, zmod)188 179 DO jj = 2, jpjm1 189 180 DO ji = fs_2, fs_jpim1 ! vect. opt. … … 195 186 END DO 196 187 END DO 197 !$OMP END PARALLEL198 188 CALL lbc_lnk( taum(:,:), 'T', 1. ) ; CALL lbc_lnk( wndm(:,:), 'T', 1. ) 199 189
Note: See TracChangeset
for help on using the changeset viewer.