New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12068 for NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/SBC/sbcmod.F90 – NEMO

Ignore:
Timestamp:
2019-12-05T13:18:21+01:00 (4 years ago)
Author:
davestorkey
Message:

2019/UKMO_MERGE_2019 : Merging in changes from ENHANCE-02_ISF_nemo.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/SBC/sbcmod.F90

    r11822 r12068  
    3737#endif 
    3838   USE sbcice_cice    ! surface boundary condition: CICE sea-ice model 
    39    USE sbcisf         ! surface boundary condition: ice-shelf 
    4039   USE sbccpl         ! surface boundary condition: coupled formulation 
    4140   USE cpl_oasis3     ! OASIS routines for coupling 
     
    4342   USE sbcrnf         ! surface boundary condition: runoffs 
    4443   USE sbcapr         ! surface boundary condition: atmo pressure  
    45    USE sbcisf         ! surface boundary condition: ice shelf 
    4644   USE sbcfwb         ! surface boundary condition: freshwater budget 
    4745   USE icbstp         ! Icebergs 
     
    9290      LOGICAL ::   ll_purecpl, ll_opa, ll_not_nemo   ! local logical 
    9391      !! 
    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_sdrift 
     92      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 
    102100      !!---------------------------------------------------------------------- 
    103101      ! 
     
    154152         WRITE(numout,*) '         Patm gradient added in ocean & ice Eqs.    ln_apr_dyn    = ', ln_apr_dyn 
    155153         WRITE(numout,*) '         runoff / runoff mouths                     ln_rnf        = ', ln_rnf 
    156          WRITE(numout,*) '         iceshelf formulation                       ln_isf        = ', ln_isf 
    157154         WRITE(numout,*) '         nb of iterations if land-sea-mask applied  nn_lsm        = ', nn_lsm 
    158155         WRITE(numout,*) '         surface wave                               ln_wave       = ', ln_wave 
     
    241238#endif 
    242239      ! 
    243       IF( .NOT.ln_isf ) THEN        !* No ice-shelf in the domain : allocate and set to zero 
    244          IF( sbc_isf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' ) 
    245          fwfisf  (:,:)   = 0._wp   ;   risf_tsc  (:,:,:) = 0._wp 
    246          fwfisf_b(:,:)   = 0._wp   ;   risf_tsc_b(:,:,:) = 0._wp 
    247       END IF 
    248240      IF( nn_ice == 0 ) THEN        !* No sea-ice in the domain : ice fraction is always zero 
    249241         IF( nn_components /= jp_iam_opa )   fr_i(:,:) = 0._wp    ! except for OPA in SAS-OPA coupled case 
     
    335327      IF( ln_ssr      )   CALL sbc_ssr_init              ! Sea-Surface Restoring initialization 
    336328      ! 
    337       IF( ln_isf      )   CALL sbc_isf_init( Kmm )       ! Compute iceshelves 
    338329      ! 
    339330                          CALL sbc_rnf_init( Kmm )       ! Runof initialization 
     
    404395            rnf_b    (:,:  ) = rnf    (:,:  ) 
    405396            rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 
    406          ENDIF 
    407          IF( ln_isf )  THEN 
    408             fwfisf_b  (:,:  ) = fwfisf  (:,:  )                
    409             risf_tsc_b(:,:,:) = risf_tsc(:,:,:)               
    410397         ENDIF 
    411398        ! 
     
    458445      ENDIF 
    459446 
    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 
    463448 
    464449      IF( ln_ssr         )   CALL sbc_ssr( kt )                        ! add SST/SSS damping term 
     
    561546      ! 
    562547      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                 , clinfo1=' qns      - : ', mask1=tmask ) 
    567          CALL prt_ctl(tab2d_1=qsr                 , clinfo1=' qsr      - : ', mask1=tmask ) 
    568          CALL prt_ctl(tab3d_1=tmask               , clinfo1=' tmask    - : ', mask1=tmask, kdim=jpk ) 
     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 ) 
    569554         CALL prt_ctl(tab3d_1=ts(:,:,:,jp_tem,Kmm), clinfo1=' sst      - : ', mask1=tmask, kdim=1   ) 
    570555         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.