Changeset 4627 for trunk/NEMOGCM/NEMO/OPA_SRC
- Timestamp:
- 2014-05-07T15:58:44+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r4624 r4627 60 60 PUBLIC sbc_ice_cice ! routine called by sbc 61 61 62 INTEGER , PARAMETER :: ji_off = INT ( (jpiglo - nx_global) / 2 )63 INTEGER , PARAMETER :: jj_off = INT ( (jpjglo - ny_global) / 2 )62 INTEGER :: ji_off 63 INTEGER :: jj_off 64 64 65 65 INTEGER , PARAMETER :: jpfld = 13 ! maximum number of files to read … … 158 158 IF(lwp) WRITE(numout,*)'cice_sbc_init' 159 159 160 ji_off = INT ( (jpiglo - nx_global) / 2 ) 161 jj_off = INT ( (jpjglo - ny_global) / 2 ) 162 160 163 ! Initialize CICE 161 164 CALL CICE_Initialize … … 220 223 sshb(:,:) = sshb(:,:) - snwice_mass(:,:) * r1_rau0 221 224 ! 222 ! Note: Changed the initial values of sshb and sshn=> need to recompute ssh[u,v,f]_[b,n]223 ! which were previously set in domvvl224 IF ( lk_vvl ) THEN ! Is this necessary? embd 2 should be restricted to vvl only???225 DO jj = 1, jpjm1226 DO ji = 1, jpim1 ! caution: use of Vector Opt. not possible227 zcoefu = 0.5 * umask(ji,jj,1) / ( e1u(ji,jj) * e2u(ji,jj) )228 zcoefv = 0.5 * vmask(ji,jj,1) / ( e1v(ji,jj) * e2v(ji,jj) )229 zcoeff = 0.25 * umask(ji,jj,1) * umask(ji,jj+1,1)230 sshu_b(ji,jj) = zcoefu * ( e1t(ji ,jj) * e2t(ji ,jj) * sshb(ji ,jj) &231 & + e1t(ji+1,jj) * e2t(ji+1,jj) * sshb(ji+1,jj) )232 sshv_b(ji,jj) = zcoefv * ( e1t(ji,jj ) * e2t(ji,jj ) * sshb(ji,jj ) &233 & + e1t(ji,jj+1) * e2t(ji,jj+1) * sshb(ji,jj+1) )234 sshu_n(ji,jj) = zcoefu * ( e1t(ji ,jj) * e2t(ji ,jj) * sshn(ji ,jj) &235 & + e1t(ji+1,jj) * e2t(ji+1,jj) * sshn(ji+1,jj) )236 sshv_n(ji,jj) = zcoefv * ( e1t(ji,jj ) * e2t(ji,jj ) * sshn(ji,jj ) &237 & + e1t(ji,jj+1) * e2t(ji,jj+1) * sshn(ji,jj+1) )238 END DO239 END DO240 CALL lbc_lnk( sshu_b, 'U', 1. ) ; CALL lbc_lnk( sshu_n, 'U', 1. )241 CALL lbc_lnk( sshv_b, 'V', 1. ) ; CALL lbc_lnk( sshv_n, 'V', 1. )242 DO jj = 1, jpjm1243 DO ji = 1, jpim1 ! NO Vector Opt.244 sshf_n(ji,jj) = 0.5 * umask(ji,jj,1) * umask(ji,jj+1,1) &245 & / ( e1f(ji,jj ) * e2f(ji,jj ) ) &246 & * ( e1u(ji,jj ) * e2u(ji,jj ) * sshu_n(ji,jj ) &247 & + e1u(ji,jj+1) * e2u(ji,jj+1) * sshu_n(ji,jj+1) )248 END DO249 END DO250 CALL lbc_lnk( sshf_n, 'F', 1. )251 ENDIF252 225 ENDIF 253 226
Note: See TracChangeset
for help on using the changeset viewer.