Changeset 11053 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/SBC/sbcmod.F90
- Timestamp:
- 2019-05-24T12:53:06+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/SBC/sbcmod.F90
r10998 r11053 341 341 IF( sbc_ice_alloc() /= 0 ) CALL ctl_stop('STOP', 'sbc_ice_alloc : unable to allocate arrays' ) 342 342 ELSEIF( nn_ice == 2 ) THEN 343 CALL ice_init( Kbb, Kmm, Kaa ) ! ICE initialization343 CALL ice_init( Kbb, Kmm, Kaa ) ! ICE initialization 344 344 ENDIF 345 345 #endif 346 IF( nn_ice == 3 ) CALL cice_sbc_init( nsbc ) ! CICE initialization347 ! 348 IF( ln_wave ) CALL sbc_wave_init ! surface wave initialisation346 IF( nn_ice == 3 ) CALL cice_sbc_init( nsbc, Kbb, Kmm ) ! CICE initialization 347 ! 348 IF( ln_wave ) CALL sbc_wave_init ! surface wave initialisation 349 349 ! 350 350 IF( lwxios ) THEN … … 442 442 CASE( 1 ) ; CALL sbc_ice_if ( kt, Kbb, Kmm ) ! Ice-cover climatology ("Ice-if" model) 443 443 #if defined key_si3 444 CASE( 2 ) ; CALL ice_stp ( kt, Kbb, nsbc )! SI3 ice model444 CASE( 2 ) ; CALL ice_stp ( kt, Kbb, Kmm, nsbc ) ! SI3 ice model 445 445 #endif 446 446 CASE( 3 ) ; CALL sbc_ice_cice ( kt, nsbc ) ! CICE ice model … … 454 454 ENDIF 455 455 456 IF( ln_isf ) CALL sbc_isf( kt, Kmm ) ! compute iceshelves457 458 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes459 460 IF( ln_ssr ) CALL sbc_ssr( kt ) ! add SST/SSS damping term461 462 IF( nn_fwb /= 0 ) CALL sbc_fwb( kt, nn_fwb, nn_fsbc ) ! control the freshwater budget456 IF( ln_isf ) CALL sbc_isf( kt, Kmm ) ! compute iceshelves 457 458 IF( ln_rnf ) CALL sbc_rnf( kt ) ! add runoffs to fresh water fluxes 459 460 IF( ln_ssr ) CALL sbc_ssr( kt ) ! add SST/SSS damping term 461 462 IF( nn_fwb /= 0 ) CALL sbc_fwb( kt, nn_fwb, nn_fsbc, Kmm ) ! control the freshwater budget 463 463 464 464 ! Special treatment of freshwater fluxes over closed seas in the model domain … … 471 471 IF ( ll_wd ) THEN ! If near WAD point limit the flux for now 472 472 zthscl = atanh(rn_wd_sbcfra) ! taper frac default is .999 473 zwdht(:,:) = ssh n(:,:) + ht_0(:,:) - rn_wdmin1 ! do this calc of water473 zwdht(:,:) = ssh(:,:,Kmm) + ht_0(:,:) - rn_wdmin1 ! do this calc of water 474 474 ! depth above wd limit once 475 475 WHERE( zwdht(:,:) <= 0.0 ) … … 557 557 ! 558 558 IF(ln_ctl) THEN ! print mean trends (used for debugging) 559 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ', mask1=tmask )560 CALL prt_ctl(tab2d_1=(emp-rnf + fwfisf) , clinfo1=' emp-rnf - : ', mask1=tmask )561 CALL prt_ctl(tab2d_1=(sfx-rnf + fwfisf) , clinfo1=' sfx-rnf - : ', mask1=tmask )562 CALL prt_ctl(tab2d_1=qns , clinfo1=' qns - : ', mask1=tmask )563 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask )564 CALL prt_ctl(tab3d_1=tmask , clinfo1=' tmask - : ', mask1=tmask, kdim=jpk )565 CALL prt_ctl(tab3d_1=ts n(:,:,:,jp_tem), clinfo1=' sst - : ', mask1=tmask, kdim=1 )566 CALL prt_ctl(tab3d_1=ts n(:,:,:,jp_sal), clinfo1=' sss - : ', mask1=tmask, kdim=1 )567 CALL prt_ctl(tab2d_1=utau , clinfo1=' utau - : ', mask1=umask, &568 & tab2d_2=vtau , clinfo2=' vtau - : ', mask2=vmask )559 CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ' , mask1=tmask ) 560 CALL prt_ctl(tab2d_1=(emp-rnf + fwfisf) , clinfo1=' emp-rnf - : ' , mask1=tmask ) 561 CALL prt_ctl(tab2d_1=(sfx-rnf + fwfisf) , clinfo1=' sfx-rnf - : ' , mask1=tmask ) 562 CALL prt_ctl(tab2d_1=qns , clinfo1=' qns - : ', mask1=tmask ) 563 CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask ) 564 CALL prt_ctl(tab3d_1=tmask , clinfo1=' tmask - : ', mask1=tmask, kdim=jpk ) 565 CALL prt_ctl(tab3d_1=ts(:,:,:,jp_tem,Kmm), clinfo1=' sst - : ', mask1=tmask, kdim=1 ) 566 CALL prt_ctl(tab3d_1=ts(:,:,:,jp_sal,Kmm), clinfo1=' sss - : ', mask1=tmask, kdim=1 ) 567 CALL prt_ctl(tab2d_1=utau , clinfo1=' utau - : ', mask1=umask, & 568 & tab2d_2=vtau , clinfo2=' vtau - : ', mask2=vmask ) 569 569 ENDIF 570 570
Note: See TracChangeset
for help on using the changeset viewer.