- Timestamp:
- 2015-12-04T17:05:58+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90
r5866 r6004 29 29 USE sbc_ice ! surface boundary condition: ice 30 30 USE sbc_oce ! surface boundary condition: ocean 31 USE sbccpl 31 USE sbccpl ! surface boundary condition: coupled interface 32 32 USE oce , ONLY : sshn, sshb, snwice_mass, snwice_mass_b, snwice_fmass 33 33 USE albedo ! albedo parameters 34 ! 34 35 USE lbclnk ! ocean lateral boundary condition - MPP exchanges 35 36 USE lib_mpp ! MPP library … … 43 44 PRIVATE 44 45 45 PUBLIC lim_sbc_init_2 46 PUBLIC lim_sbc_flx_2 47 PUBLIC lim_sbc_tau_2 46 PUBLIC lim_sbc_init_2 ! called by ice_init_2 47 PUBLIC lim_sbc_flx_2 ! called by sbc_ice_lim_2 48 PUBLIC lim_sbc_tau_2 ! called by sbc_ice_lim_2 48 49 49 50 REAL(wp) :: r1_rdtice ! = 1. / rdt_ice … … 52 53 REAL(wp) :: rone = 1._wp ! - - 53 54 ! 54 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: soce_0, sice_0 ! constant SSS and ice salinity used in levitating sea-ice case55 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: soce_0, sice_0 ! fix SSS and ice salinity used in levitating case 0 55 56 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: utau_oce, vtau_oce ! air-ocean surface i- & j-stress [N/m2] 56 57 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: tmod_io ! modulus of the ice-ocean relative velocity [m/s] … … 101 102 !!--------------------------------------------------------------------- 102 103 INTEGER, INTENT(in) :: kt ! number of iteration 103 ! !104 ! 104 105 INTEGER :: ji, jj ! dummy loop indices 105 106 INTEGER :: ii0, ii1, ij0, ij1 ! local integers … … 113 114 REAL(wp), POINTER, DIMENSION(:,:,:) :: zalb, zalbp ! 2D/3D workspace 114 115 !!--------------------------------------------------------------------- 115 116 ! 116 117 CALL wrk_alloc( jpi, jpj, zqnsoce ) 117 118 CALL wrk_alloc( jpi, jpj, 1, zalb, zalbp ) 118 119 SELECT CASE( nn_ice_embd ) 120 CASE( 0 ) ; zswitch = 1 ! (0) standard levitating sea-ice : salt exchange only121 CASE( 1, 2 ) ; zswitch = 0! (1) levitating sea-ice: salt and volume exchange but no pressure effect122 123 END SELECT !119 ! 120 SELECT CASE( nn_ice_embd ) ! levitating or embedded sea-ice option 121 CASE( 0 ) ; zswitch = 1 ! (0) old levitating sea-ice : salt exchange only 122 CASE( 1, 2 ) ; zswitch = 0 ! (1) levitating sea-ice: salt and volume exchange but no pressure effect 123 ! ! (2) embedded sea-ice : salt and volume fluxes and pressure 124 END SELECT 124 125 125 126 !------------------------------------------! … … 302 303 INTEGER , INTENT(in) :: kt ! ocean time-step index 303 304 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pu_oce, pv_oce ! surface ocean currents 304 ! !305 ! 305 306 INTEGER :: ji, jj ! dummy loop indices 306 307 REAL(wp) :: zfrldu, zat_u, zu_i, zutau_ice, zu_t, zmodt ! local scalar … … 434 435 !! ** input : Namelist namicedia 435 436 !!------------------------------------------------------------------- 436 !437 INTEGER :: jk ! local integer437 INTEGER :: jk ! local integer 438 !!------------------------------------------------------------------- 438 439 ! 439 440 IF(lwp) WRITE(numout,*) … … 471 472 !!gm I really don't like this staff here... Find a way to put that elsewhere or differently 472 473 !!gm 473 IF( .NOT. 474 IF( .NOT.ln_linssh ) THEN 474 475 475 476 do jk = 1,jpkm1 ! adjust initial vertical scale factors
Note: See TracChangeset
for help on using the changeset viewer.