- Timestamp:
- 2019-11-26T16:06:55+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo_TEST_MERGE/src/OCE/SBC/sbcmod.F90
r11536 r11970 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 … … 91 89 LOGICAL :: ll_purecpl, ll_opa, ll_not_nemo ! local logical 92 90 !! 93 NAMELIST/namsbc/ nn_fsbc , 94 & ln_usr , ln_flx , ln_blk ,&95 & ln_cpl , ln_mixcpl , nn_components,&96 & nn_ice , ln_ice_embd, 97 & ln_traqsr, ln_dm2dc ,&98 & ln_rnf , nn_fwb , ln_ssr , ln_isf , ln_apr_dyn ,&99 & ln_wave , ln_cdgw , ln_sdw , ln_tauwoc , ln_stcor ,&100 & ln_tauw , nn_lsm , nn_sdrift91 NAMELIST/namsbc/ nn_fsbc , & 92 & ln_usr , ln_flx , ln_blk , & 93 & ln_cpl , ln_mixcpl , nn_components, & 94 & nn_ice , ln_ice_embd, & 95 & ln_traqsr, ln_dm2dc , & 96 & ln_rnf , nn_fwb , ln_ssr , ln_apr_dyn, & 97 & ln_wave , ln_cdgw , ln_sdw , ln_tauwoc , ln_stcor, & 98 & ln_tauw , nn_lsm , nn_sdrift 101 99 !!---------------------------------------------------------------------- 102 100 ! … … 153 151 WRITE(numout,*) ' Patm gradient added in ocean & ice Eqs. ln_apr_dyn = ', ln_apr_dyn 154 152 WRITE(numout,*) ' runoff / runoff mouths ln_rnf = ', ln_rnf 155 WRITE(numout,*) ' iceshelf formulation ln_isf = ', ln_isf156 153 WRITE(numout,*) ' nb of iterations if land-sea-mask applied nn_lsm = ', nn_lsm 157 154 WRITE(numout,*) ' surface wave ln_wave = ', ln_wave … … 239 236 #endif 240 237 ! 241 IF( .NOT.ln_isf ) THEN !* No ice-shelf in the domain : allocate and set to zero242 IF( sbc_isf_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' )243 fwfisf (:,:) = 0._wp ; risf_tsc (:,:,:) = 0._wp244 fwfisf_b(:,:) = 0._wp ; risf_tsc_b(:,:,:) = 0._wp245 END IF246 238 IF( nn_ice == 0 ) THEN !* No sea-ice in the domain : ice fraction is always zero 247 239 IF( nn_components /= jp_iam_opa ) fr_i(:,:) = 0._wp ! except for OPA in SAS-OPA coupled case … … 333 325 IF( ln_ssr ) CALL sbc_ssr_init ! Sea-Surface Restoring initialization 334 326 ! 335 IF( ln_isf ) CALL sbc_isf_init ! Compute iceshelves336 !337 327 CALL sbc_rnf_init ! Runof initialization 338 328 ! … … 401 391 rnf_b (:,: ) = rnf (:,: ) 402 392 rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 403 ENDIF404 IF( ln_isf ) THEN405 fwfisf_b (:,: ) = fwfisf (:,: )406 risf_tsc_b(:,:,:) = risf_tsc(:,:,:)407 393 ENDIF 408 394 ! … … 454 440 IF( .NOT. ln_passive_mode ) CALL lbc_lnk( 'sbcmod', emp, 'T', 1. ) ! ensure restartability with icebergs 455 441 ENDIF 456 457 IF( ln_isf ) CALL sbc_isf( kt ) ! compute iceshelves458 442 459 443 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes … … 558 542 ! 559 543 IF(ln_ctl) THEN ! print mean trends (used for debugging) 560 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i- : ', mask1=tmask )561 CALL prt_ctl(tab2d_1=(emp-rnf + fwfisf), clinfo1=' emp-rnf- : ', mask1=tmask )562 CALL prt_ctl(tab2d_1=(sfx-rnf + fwfisf), clinfo1=' sfx-rnf- : ', mask1=tmask )544 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ', mask1=tmask ) 545 CALL prt_ctl(tab2d_1=(emp-rnf) , clinfo1=' emp-rnf - : ', mask1=tmask ) 546 CALL prt_ctl(tab2d_1=(sfx-rnf) , clinfo1=' sfx-rnf - : ', mask1=tmask ) 563 547 CALL prt_ctl(tab2d_1=qns , clinfo1=' qns - : ', mask1=tmask ) 564 548 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask )
Note: See TracChangeset
for help on using the changeset viewer.