- Timestamp:
- 2019-12-10T15:16:33+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg_ts.F90
r12145 r12150 31 31 USE dom_oce ! ocean space and time domain 32 32 USE sbc_oce ! surface boundary condition: ocean 33 USE isf_oce ! ice shelf variable (fwfisf) 33 34 USE zdf_oce ! vertical physics: variables 34 35 USE zdfdrg ! vertical physics: top/bottom drag coef. 35 USE sbcisf ! ice shelf variable (fwfisf)36 36 USE sbcapr ! surface boundary condition: atmospheric pressure 37 37 USE dynadv , ONLY: ln_dynadv_vec … … 337 337 ! ! --------------------------------------------------- ! 338 338 IF (ln_bt_fw) THEN ! FORWARD integration: use kt+1/2 fluxes (NOW+1/2) 339 zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf(:,:))339 zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) 340 340 ELSE ! CENTRED integration: use kt-1/2 + kt+1/2 fluxes (NOW) 341 341 zztmp = r1_rau0 * r1_2 342 zssh_frc(:,:) = zztmp * ( emp(:,:) + emp_b(:,:) - rnf(:,:) - rnf_b(:,:) + fwfisf(:,:) + fwfisf_b(:,:) ) 342 zssh_frc(:,:) = zztmp * ( emp(:,:) + emp_b(:,:) & 343 & - rnf(:,:) - rnf_b(:,:) & 344 & + fwfisf_cav(:,:) + fwfisf_cav_b(:,:) & 345 & + fwfisf_par(:,:) + fwfisf_par_b(:,:) ) 343 346 ENDIF 344 347 ! != Add Stokes drift divergence =! (if exist) … … 346 349 zssh_frc(:,:) = zssh_frc(:,:) + div_sd(:,:) 347 350 ENDIF 351 ! 352 ! ! ice sheet coupling 353 IF ( ln_isf .AND. ln_isfcpl ) THEN 354 ! 355 ! ice sheet coupling 356 IF( ln_rstart .AND. kt == nit000 ) THEN 357 zssh_frc(:,:) = zssh_frc(:,:) + risfcpl_ssh(:,:) 358 END IF 359 ! 360 ! conservation option 361 IF( ln_isfcpl_cons ) THEN 362 zssh_frc(:,:) = zssh_frc(:,:) + risfcpl_cons_ssh(:,:) 363 END IF 364 ! 365 END IF 348 366 ! 349 367 #if defined key_asminc
Note: See TracChangeset
for help on using the changeset viewer.