Changeset 10276 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2018-11-02T18:15:27+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r10181 r10276 251 251 CALL lbc_lnk ( fr_iv , 'V', 1. ) 252 252 253 !!*****JUST ADDED BUT NEEDS CONFIRMING!!**** 253 254 254 ! Snow and ice thickness 255 CALL cice2nemo(vice,thick_i,'T', 1. ) !!Do this for snow depth too? set thick_s. Duplicated below for embedded ice? 256 257 ! vice is volume per unit area of grid cell = thickness 255 ! vice, vsno is volume per unit area of grid cell = thickness 256 CALL cice2nemo(vice,thick_i,'T', 1. ) 257 CALL cice2nemo(vsno,thick_s,'T', 1. ) 258 258 259 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN ! Confirm what this if loop is for259 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 260 260 DO jl = 1,ncat 261 261 CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. ) … … 263 263 ENDDO 264 264 ENDIF 265 265 266 266 267 ! T point to U point 267 268 ! T point to V point 269 270 ! Sea ice thickness 268 271 thick_iu(:,:)=0.0 269 272 thick_iv(:,:)=0.0 … … 278 281 CALL lbc_lnk ( thick_iv , 'V', 1. ) 279 282 280 !!********************************************************** 283 ! Snow depth 284 thick_su(:,:)=0.0 285 thick_sv(:,:)=0.0 286 DO jj=1,jpjm1 287 DO ji=1,jpim1 288 thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1) 289 thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1) 290 ENDDO 291 ENDDO 292 293 CALL lbc_lnk ( thick_su , 'U', 1. ) 294 CALL lbc_lnk ( thick_sv , 'V', 1. ) 281 295 282 296
Note: See TracChangeset
for help on using the changeset viewer.