Changeset 10332 for NEMO/trunk/src/OCE/ICB
- Timestamp:
- 2018-11-19T13:02:15+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/ICB/icbutl.F90
r10068 r10332 24 24 #if defined key_si3 25 25 USE ice, ONLY: u_ice, v_ice, hm_i ! SI3 variables 26 USE icevar ! ice_var_sshdyn 27 USE sbc_ice, ONLY: snwice_mass, snwice_mass_b 26 28 #endif 27 29 … … 60 62 !! ** Method : - blah blah 61 63 !!---------------------------------------------------------------------- 62 64 #if defined key_si3 65 REAL(wp), DIMENSION(jpi,jpj) :: zssh_lead_m ! ocean surface (ssh_m) if ice is not embedded 66 ! ! ocean surface in leads if ice is embedded 67 #endif 63 68 ! copy nemo forcing arrays into iceberg versions with extra halo 64 69 ! only necessary for variables not on T points … … 84 89 ui_e(:,:) = 0._wp ; ui_e(1:jpi, 1:jpj) = u_ice(:,:) 85 90 vi_e(:,:) = 0._wp ; vi_e(1:jpi, 1:jpj) = v_ice(:,:) 91 ! 92 ! compute ssh slope using ssh_lead if embedded 93 zssh_lead_m(:,:) = ice_var_sshdyn(ssh_m, snwice_mass, snwice_mass_b) 94 ssh_e(:,:) = 0._wp ; ssh_e(1:jpi, 1:jpj) = zssh_lead_m(:,:) * tmask(:,:,1) 86 95 ! 87 96 CALL lbc_lnk_icb( hicth, 'T', +1._wp, 1, 1 ) 88 97 CALL lbc_lnk_icb( ui_e , 'U', -1._wp, 1, 1 ) 89 98 CALL lbc_lnk_icb( vi_e , 'V', -1._wp, 1, 1 ) 99 #else 100 ssh_e(:,:) = 0._wp ; ssh_e(1:jpi, 1:jpj) = ssh_m(:,:) * tmask(:,:,1) 90 101 #endif 91 102 92 103 !! special for ssh which is used to calculate slope 93 104 !! so fudge some numbers all the way around the boundary 94 95 ssh_e(:,:) = 0._wp ; ssh_e(1:jpi, 1:jpj) = ssh_m(:,:) * tmask(:,:,1)96 105 ssh_e(0 , :) = ssh_e(1 , :) 97 106 ssh_e(jpi+1, :) = ssh_e(jpi, :)
Note: See TracChangeset
for help on using the changeset viewer.