- Timestamp:
- 2020-07-08T13:34:28+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_2494_2375/src/OCE/ICB/icbdyn.F90
r10570 r13265 14 14 USE dom_oce ! NEMO ocean domain 15 15 USE phycst ! NEMO physical constants 16 USE in_out_manager ! IO parameters 16 17 ! 17 18 USE icb_oce ! define iceberg arrays … … 155 156 pt%xi = zxi_n 156 157 pt%yj = zyj_n 157 158 ! update actual position159 pt%lon = icb_utl_bilin_x(glamt, pt%xi, pt%yj )160 pt%lat = icb_utl_bilin(gphit, pt%xi, pt%yj, 'T' )161 158 162 159 berg => berg%next ! switch to the next berg … … 270 267 ! Interpolate gridded fields to berg 271 268 nknberg = berg%number(1) 272 CALL icb_utl_interp( pxi, pe1, zuo, zui, zua, zssh_x, & 273 & pyj, pe2, zvo, zvi, zva, zssh_y, zsst, zcn, zhi, zff ) 269 CALL icb_utl_interp( pxi, pyj, pe1=pe1, pe2=pe2, & ! scale factor 270 & puo=zuo, pui=zui, pua=zua, & ! oce/ice/atm velocities 271 & pvo=zvo, pvi=zvi, pva=zva, & ! oce/ice/atm velocities 272 & pssh_i=zssh_x, pssh_j=zssh_y, & ! ssh gradient 273 & phi=zhi, pff=zff ) ! ice thickness and coriolis 274 275 IF (lwp) THEN 276 WRITE(numout,*) 'icbdyn ', pxi, pyj 277 WRITE(numout,*) pe1, zuo, zui, zua, zssh_x 278 WRITE(numout,*) pe2, zvo, zvi, zva, zssh_y 279 WRITE(numout,*) zhi, zff 280 WRITE(numout,*) '' 281 END IF 282 274 283 275 284 zM = berg%current_point%mass
Note: See TracChangeset
for help on using the changeset viewer.