Changeset 10181 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2018-10-09T11:29:47+02: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
r9987 r10181 57 57 flatn_f,fsurfn_f,fcondtopn_f, & 58 58 #ifdef key_asminc 59 daice_da, fresh_da,fsalt_da,&59 daice_da,dsit_da,fresh_da,fsalt_da, & 60 60 #endif 61 61 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl, & … … 250 250 CALL lbc_lnk ( fr_iu , 'U', 1. ) 251 251 CALL lbc_lnk ( fr_iv , 'V', 1. ) 252 253 !!*****JUST ADDED BUT NEEDS CONFIRMING!!**** 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 258 259 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN ! Confirm what this if loop is for 260 DO jl = 1,ncat 261 CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. ) 262 CALL cice2nemo(vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1. ) 263 ENDDO 264 ENDIF 265 266 ! T point to U point 267 ! T point to V point 268 thick_iu(:,:)=0.0 269 thick_iv(:,:)=0.0 270 DO jj=1,jpjm1 271 DO ji=1,jpim1 272 thick_iu(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji+1,jj))*umask(ji,jj,1) 273 thick_iv(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji,jj+1))*vmask(ji,jj,1) 274 ENDDO 275 ENDDO 276 277 CALL lbc_lnk ( thick_iu , 'U', 1. ) 278 CALL lbc_lnk ( thick_iv , 'V', 1. ) 279 280 !!********************************************************** 281 252 282 253 283 ! ! embedded sea ice … … 310 340 nfresh_da(:,:) = 0.0 311 341 nfsalt_da(:,:) = 0.0 312 ndaice_da(:,:) = 0.0 342 ndaice_da(:,:) = 0.0 343 ndsit_da(:,:) = 0.0 313 344 #endif 314 345 ! … … 469 500 ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1) 470 501 Call nemo2cice(ztmp,daice_da,'T', 1. ) 502 !Ice thickness change (from assimilation) 503 ztmp(:,:)=ndsit_da(:,:)*tmask(:,:,1) 504 Call nemo2cice(ztmp,dsit_da,'T', 1. ) 471 505 #endif 472 506
Note: See TracChangeset
for help on using the changeset viewer.