- Timestamp:
- 2019-12-05T13:18:21+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/SBC/sbcmod.F90
r11822 r12068 37 37 #endif 38 38 USE sbcice_cice ! surface boundary condition: CICE sea-ice model 39 USE sbcisf ! surface boundary condition: ice-shelf40 39 USE sbccpl ! surface boundary condition: coupled formulation 41 40 USE cpl_oasis3 ! OASIS routines for coupling … … 43 42 USE sbcrnf ! surface boundary condition: runoffs 44 43 USE sbcapr ! surface boundary condition: atmo pressure 45 USE sbcisf ! surface boundary condition: ice shelf46 44 USE sbcfwb ! surface boundary condition: freshwater budget 47 45 USE icbstp ! Icebergs … … 92 90 LOGICAL :: ll_purecpl, ll_opa, ll_not_nemo ! local logical 93 91 !! 94 NAMELIST/namsbc/ nn_fsbc , 95 & ln_usr , ln_flx , ln_blk ,&96 & ln_cpl , ln_mixcpl , nn_components,&97 & nn_ice , ln_ice_embd, 98 & ln_traqsr, ln_dm2dc ,&99 & ln_rnf , nn_fwb , ln_ssr , ln_isf , ln_apr_dyn ,&100 & ln_wave , ln_cdgw , ln_sdw , ln_tauwoc , ln_stcor ,&101 & ln_tauw , nn_lsm , nn_sdrift92 NAMELIST/namsbc/ nn_fsbc , & 93 & ln_usr , ln_flx , ln_blk , & 94 & ln_cpl , ln_mixcpl , nn_components, & 95 & nn_ice , ln_ice_embd, & 96 & ln_traqsr, ln_dm2dc , & 97 & ln_rnf , nn_fwb , ln_ssr , ln_apr_dyn, & 98 & ln_wave , ln_cdgw , ln_sdw , ln_tauwoc , ln_stcor, & 99 & ln_tauw , nn_lsm , nn_sdrift 102 100 !!---------------------------------------------------------------------- 103 101 ! … … 154 152 WRITE(numout,*) ' Patm gradient added in ocean & ice Eqs. ln_apr_dyn = ', ln_apr_dyn 155 153 WRITE(numout,*) ' runoff / runoff mouths ln_rnf = ', ln_rnf 156 WRITE(numout,*) ' iceshelf formulation ln_isf = ', ln_isf157 154 WRITE(numout,*) ' nb of iterations if land-sea-mask applied nn_lsm = ', nn_lsm 158 155 WRITE(numout,*) ' surface wave ln_wave = ', ln_wave … … 241 238 #endif 242 239 ! 243 IF( .NOT.ln_isf ) THEN !* No ice-shelf in the domain : allocate and set to zero244 IF( sbc_isf_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' )245 fwfisf (:,:) = 0._wp ; risf_tsc (:,:,:) = 0._wp246 fwfisf_b(:,:) = 0._wp ; risf_tsc_b(:,:,:) = 0._wp247 END IF248 240 IF( nn_ice == 0 ) THEN !* No sea-ice in the domain : ice fraction is always zero 249 241 IF( nn_components /= jp_iam_opa ) fr_i(:,:) = 0._wp ! except for OPA in SAS-OPA coupled case … … 335 327 IF( ln_ssr ) CALL sbc_ssr_init ! Sea-Surface Restoring initialization 336 328 ! 337 IF( ln_isf ) CALL sbc_isf_init( Kmm ) ! Compute iceshelves338 329 ! 339 330 CALL sbc_rnf_init( Kmm ) ! Runof initialization … … 404 395 rnf_b (:,: ) = rnf (:,: ) 405 396 rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 406 ENDIF407 IF( ln_isf ) THEN408 fwfisf_b (:,: ) = fwfisf (:,: )409 risf_tsc_b(:,:,:) = risf_tsc(:,:,:)410 397 ENDIF 411 398 ! … … 458 445 ENDIF 459 446 460 IF( ln_isf ) CALL sbc_isf( kt, Kmm ) ! compute iceshelves 461 462 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes 447 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes 463 448 464 449 IF( ln_ssr ) CALL sbc_ssr( kt ) ! add SST/SSS damping term … … 561 546 ! 562 547 IF(ln_ctl) THEN ! print mean trends (used for debugging) 563 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ', mask1=tmask )564 CALL prt_ctl(tab2d_1=(emp-rnf + fwfisf) , clinfo1=' emp-rnf - : ', mask1=tmask )565 CALL prt_ctl(tab2d_1=(sfx-rnf + fwfisf) , clinfo1=' sfx-rnf - : ', mask1=tmask )566 CALL prt_ctl(tab2d_1=qns 567 CALL prt_ctl(tab2d_1=qsr 568 CALL prt_ctl(tab3d_1=tmask 548 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ', mask1=tmask ) 549 CALL prt_ctl(tab2d_1=(emp-rnf) , clinfo1=' emp-rnf - : ', mask1=tmask ) 550 CALL prt_ctl(tab2d_1=(sfx-rnf) , clinfo1=' sfx-rnf - : ', mask1=tmask ) 551 CALL prt_ctl(tab2d_1=qns , clinfo1=' qns - : ', mask1=tmask ) 552 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask ) 553 CALL prt_ctl(tab3d_1=tmask , clinfo1=' tmask - : ', mask1=tmask, kdim=jpk ) 569 554 CALL prt_ctl(tab3d_1=ts(:,:,:,jp_tem,Kmm), clinfo1=' sst - : ', mask1=tmask, kdim=1 ) 570 555 CALL prt_ctl(tab3d_1=ts(:,:,:,jp_sal,Kmm), clinfo1=' sss - : ', mask1=tmask, kdim=1 )
Note: See TracChangeset
for help on using the changeset viewer.