- Timestamp:
- 2020-06-24T14:38:26+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/sbcice_cice.F90
r12614 r13151 12 12 USE oce ! ocean dynamics and tracers 13 13 USE dom_oce ! ocean space and time domain 14 !!st8 15 # if ! defined key_qco 14 16 USE domvvl 17 # else 18 USE domqco 19 # endif 20 !!st8 15 21 USE phycst, only : rcp, rho0, r1_rho0, rhos, rhoi 16 22 USE in_out_manager ! I/O manager … … 233 239 !!gm This should be put elsewhere.... (same remark for limsbc) 234 240 !!gm especially here it is assumed zstar coordinate, but it can be ztilde.... 241 !!st9 242 #if defined key_qco 243 IF( .NOT.ln_linssh ) CALL dom_qco_zgr( Kbb, Kmm, Kaa ) ! interpolation scale factor, depth and water column 244 #else 235 245 IF( .NOT.ln_linssh ) THEN 236 246 ! 237 247 DO jk = 1,jpkm1 ! adjust initial vertical scale factors 238 e3t(:,:,jk,Kmm) = e3t_0(:,:,jk)*( 1._wp + ssh(:,:,Kmm)* tmask(:,:,1)/(ht_0(:,:) + 1.0 - tmask(:,:,1)) )239 e3t(:,:,jk,Kbb) = e3t_0(:,:,jk)*( 1._wp + ssh(:,:,Kbb)* tmask(:,:,1)/(ht_0(:,:) + 1.0 - tmask(:,:,1)) )240 END DO248 e3t(:,:,jk,Kmm) = e3t_0(:,:,jk)*( 1._wp + ssh(:,:,Kmm)*r1_ht_0(:,:)*tmask(:,:,jk) ) 249 e3t(:,:,jk,Kbb) = e3t_0(:,:,jk)*( 1._wp + ssh(:,:,Kbb)*r1_ht_0(:,:)*tmask(:,:,jk) ) 250 END DO 241 251 e3t(:,:,:,Krhs) = e3t(:,:,:,Kbb) 242 252 ! Reconstruction of all vertical scale factors at now and before time-steps … … 267 277 END DO 268 278 ENDIF 279 #endif 280 !!st9 269 281 ENDIF 270 282 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.