Changeset 13299 for NEMO/branches
- Timestamp:
- 2020-07-13T14:34:12+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/ticket2487
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/ticket2487/cfgs/SHARED/namelist_ref
r12288 r13299 213 213 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 214 214 ! ! =F levitating ice (no pressure, mass and salt exchanges) 215 ln_ice_sladj = .false. ! sea-level adjustment due to initial ice and snow mass 215 216 ! Misc. options of sbc : 216 217 ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/ticket2487/src/ICE/iceistate.F90
r13294 r13299 18 18 USE oce ! dynamics and tracers variables 19 19 USE dom_oce ! ocean domain 20 USE sbc_oce , ONLY : sst_m, sss_m, ln_ice_embd 20 USE sbc_oce , ONLY : sst_m, sss_m, ln_ice_embd , ln_ice_sladj 21 21 USE sbc_ice , ONLY : tn_ice, snwice_mass, snwice_mass_b 22 22 USE eosbn2 ! equation of state … … 379 379 at_i(:,:) = SUM( a_i, dim=3 ) 380 380 ! 381 !----------------- -----------------------------382 ! 3) Snow-ice mass (case ice is fully embedded)383 !----------------- -----------------------------381 !----------------- 382 ! 3) Snow-ice mass 383 !----------------- 384 384 snwice_mass (:,:) = tmask(:,:,1) * SUM( rhos * v_s(:,:,:) + rhoi * v_i(:,:,:), dim=3 ) ! snow+ice mass 385 385 snwice_mass_b(:,:) = snwice_mass(:,:) … … 428 428 END DO 429 429 ENDIF 430 ELSE 430 ELSEIF( ln_ice_sladj ) THEN ! adjustment of sea level due to initial snow+ice mass 431 431 z1_area = 1.0_wp / glob_sum( 'iceistate', e1e2t(:,:) ) 432 432 zsshadj = glob_sum( 'iceistate', e1e2t(:,:) * snwice_mass(:,:) ) * r1_rau0 * z1_area … … 438 438 rsshadj = zsshadj 439 439 #endif 440 WRITE(ctmp1,'(A34,F10.6,A32)') 'iceistate: mean ssh adjusted by ', -1.0_wp * zsshadj, ' m to compensate for the initial'441 CALL ctl_warn( ctmp1, ' ice+snow mass' )440 IF(lwp) WRITE(numout,'(A35,F10.6,A21)') 'iceistate: sea level adjusted by ', -1.0_wp * zsshadj, ' m to compensate for' 441 IF(lwp) WRITE(numout,*) ' the initial snow+ice mass' 442 442 sshn(:,:) = sshn(:,:) - zsshadj 443 443 sshb(:,:) = sshb(:,:) - zsshadj -
NEMO/branches/2020/ticket2487/src/OCE/SBC/sbc_oce.F90
r12132 r13299 50 50 ! !: =F levitating ice (no presure effect) with mass and salt exchanges 51 51 ! !: =T embedded sea-ice (pressure effect + mass and salt exchanges) 52 LOGICAL , PUBLIC :: ln_ice_sladj !: flag to enable initial sea-level adjustment due to initial sea-ice and snow mass 53 ! (this option has no effect if ln_ice_embd equals .TRUE.) 52 54 INTEGER , PUBLIC :: nn_components !: flag for sbc module (including sea-ice) coupling mode (see component definition below) 53 55 INTEGER , PUBLIC :: nn_fwb !: FreshWater Budget: -
NEMO/branches/2020/ticket2487/src/OCE/SBC/sbcice_cice.F90
r13294 r13299 271 271 END DO 272 272 ENDIF 273 ELSE 273 ELSEIF( ln_ice_sladj ) THEN ! adjustment of the sea level due to initial snow+ice mass 274 274 z1_area = 1.0_wp / glob_sum( 'sbcice_cice', e1e2t(:,:) ) 275 275 zsshadj = glob_sum( 'sbcice_cice', e1e2t(:,:) * snwice_mass(:,:) ) * r1_rau0 * z1_area … … 281 281 rsshadj = zsshadj 282 282 #endif 283 WRITE(ctmp1,'(A36,F10.6,A24)') 'sbcice_cice: mean ssh adjusted by ', -1.0_wp * zsshadj, ' m to compensate for the'284 CALL ctl_warn( ctmp1, ' initial ice+snow mass' )283 IF(lwp) WRITE(numout,'(A35,F10.6,A21)') 'iceistate: sea level adjusted by ', -1.0_wp * zsshadj, ' m to compensate for' 284 IF(lwp) WRITE(numout,*) ' the initial snow+ice mass' 285 285 sshn(:,:) = sshn(:,:) - zsshadj 286 286 sshb(:,:) = sshb(:,:) - zsshadj -
NEMO/branches/2020/ticket2487/src/OCE/SBC/sbcmod.F90
r12276 r13299 94 94 & ln_usr , ln_flx , ln_blk , & 95 95 & ln_cpl , ln_mixcpl, nn_components, & 96 & nn_ice , ln_ice_embd, 96 & nn_ice , ln_ice_embd, ln_ice_sladj, & 97 97 & ln_traqsr, ln_dm2dc , & 98 98 & ln_rnf , nn_fwb , ln_ssr , ln_isf , ln_apr_dyn , & … … 146 146 WRITE(numout,*) ' ice management in the sbc (=0/1/2/3) nn_ice = ', nn_ice 147 147 WRITE(numout,*) ' ice embedded into ocean ln_ice_embd = ', ln_ice_embd 148 WRITE(numout,*) ' sea-level adj. due to initial ice/snow ln_ice_sladj = ', ln_ice_sladj 148 149 WRITE(numout,*) ' Misc. options of sbc : ' 149 150 WRITE(numout,*) ' Light penetration in temperature Eq. ln_traqsr = ', ln_traqsr
Note: See TracChangeset
for help on using the changeset viewer.