- Timestamp:
- 2014-07-23T18:27:21+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO2_ice_shelves/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90
r4666 r4726 188 188 189 189 ! note that mbkt is set to 1 over land ==> use surface tmask 190 zprt(:,:) = lmask(:,:) * REAL( mbkt(:,:) , wp )190 zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) , wp ) 191 191 CALL iom_rstput( 0, 0, inum4, 'mbathy', zprt, ktype = jp_i2 ) ! ! nb of ocean T-points 192 zprt(:,:) = lmask(:,:) * REAL( mikt(:,:) , wp )192 zprt(:,:) = ssmask(:,:) * REAL( mikt(:,:) , wp ) 193 193 CALL iom_rstput( 0, 0, inum4, 'misf', zprt, ktype = jp_i2 ) ! ! nb of ocean T-points 194 zprt(:,:) = lmask(:,:) * REAL( icedep(:,:) , wp )194 zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp ) 195 195 CALL iom_rstput( 0, 0, inum4, 'isfdraft', zprt, ktype = jp_r4 ) ! ! nb of ocean T-points 196 196 … … 227 227 DO jj = 1,jpj 228 228 DO ji = 1,jpi 229 e3tp(ji,jj) = e3t_0(ji,jj,mbkt(ji,jj)) * lmask(ji,jj)230 e3wp(ji,jj) = e3w_0(ji,jj,mbkt(ji,jj)) * lmask(ji,jj)229 e3tp(ji,jj) = e3t_0(ji,jj,mbkt(ji,jj)) * ssmask(ji,jj) 230 e3wp(ji,jj) = e3w_0(ji,jj,mbkt(ji,jj)) * ssmask(ji,jj) 231 231 END DO 232 232 END DO … … 252 252 DO jj = 1,jpj 253 253 DO ji = 1,jpi 254 zprt(ji,jj) = gdept_0(ji,jj,mbkt(ji,jj) ) * lmask(ji,jj)255 zprw(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1) * lmask(ji,jj)254 zprt(ji,jj) = gdept_0(ji,jj,mbkt(ji,jj) ) * ssmask(ji,jj) 255 zprw(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1) * ssmask(ji,jj) 256 256 END DO 257 257 END DO
Note: See TracChangeset
for help on using the changeset viewer.