Ignore:
Timestamp:
2015-07-15T17:46:12+02:00 (5 years ago)
Author:
andrewryan
Message:

merged in latest version of trunk alongside changes to SAO_SRC to be compatible with latest OBS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r5034 r5600  
    1313   !!            3.4  ! 2011-11  (C. Harris) CICE added as an option 
    1414   !!            3.5  ! 2012-11  (A. Coward, G. Madec) Rethink of heat, mass and salt surface fluxes 
     15   !!            3.6  ! 2014-11  (P. Mathiot, C. Harris) add ice shelves melting                     
    1516   !!---------------------------------------------------------------------- 
    1617 
     
    2324   USE phycst           ! physical constants 
    2425   USE sbc_oce          ! Surface boundary condition: ocean fields 
     26   USE trc_oce          ! shared ocean-passive tracers variables 
    2527   USE sbc_ice          ! Surface boundary condition: ice fields 
    2628   USE sbcdcy           ! surface boundary condition: diurnal cycle 
     
    3739   USE sbcice_cice      ! surface boundary condition: CICE    sea-ice model 
    3840   USE sbccpl           ! surface boundary condition: coupled florulation 
     41   USE cpl_oasis3       ! OASIS routines for coupling 
    3942   USE sbcssr           ! surface boundary condition: sea surface restoring 
    4043   USE sbcrnf           ! surface boundary condition: runoffs 
     
    5053   USE timing           ! Timing 
    5154   USE sbcwave          ! Wave module 
     55   USE bdy_par          ! Require lk_bdy 
    5256 
    5357   IMPLICIT NONE 
     
    8286      INTEGER ::   icpt   ! local integer 
    8387      !! 
    84       NAMELIST/namsbc/ nn_fsbc   , ln_ana    , ln_flx,  ln_blk_clio, ln_blk_core,           & 
    85          &             ln_blk_mfs, ln_apr_dyn, nn_ice,  nn_ice_embd, ln_dm2dc   , ln_rnf,   & 
    86          &             ln_ssr    ,  nn_isf , nn_fwb    , ln_cdgw , ln_wave , ln_sdw, nn_lsm, nn_limflx 
     88      NAMELIST/namsbc/ nn_fsbc   , ln_ana    , ln_flx, ln_blk_clio, ln_blk_core, ln_mixcpl,   & 
     89         &             ln_blk_mfs, ln_apr_dyn, nn_ice, nn_ice_embd, ln_dm2dc   , ln_rnf   ,   & 
     90         &             ln_ssr    , nn_isf    , nn_fwb, ln_cdgw    , ln_wave    , ln_sdw   ,   & 
     91         &             nn_lsm    , nn_limflx , nn_components, ln_cpl 
    8792      INTEGER  ::   ios 
     93      INTEGER  ::   ierr, ierr0, ierr1, ierr2, ierr3, jpm 
     94      LOGICAL  ::   ll_purecpl 
    8895      !!---------------------------------------------------------------------- 
    8996 
     
    113120          nn_ice      =   0 
    114121      ENDIF 
    115       
     122 
    116123      IF(lwp) THEN               ! Control print 
    117124         WRITE(numout,*) '        Namelist namsbc (partly overwritten with CPP key setting)' 
     
    123130         WRITE(numout,*) '              CORE bulk  formulation                     ln_blk_core = ', ln_blk_core 
    124131         WRITE(numout,*) '              MFS  bulk  formulation                     ln_blk_mfs  = ', ln_blk_mfs 
    125          WRITE(numout,*) '              coupled    formulation (T if key_oasis3)   lk_cpl      = ', lk_cpl 
     132         WRITE(numout,*) '              ocean-atmosphere coupled formulation       ln_cpl      = ', ln_cpl 
     133         WRITE(numout,*) '              forced-coupled mixed formulation           ln_mixcpl   = ', ln_mixcpl 
     134         WRITE(numout,*) '              OASIS coupling (with atm or sas)           lk_oasis    = ', lk_oasis 
     135         WRITE(numout,*) '              components of your executable              nn_components = ', nn_components 
    126136         WRITE(numout,*) '              Multicategory heat flux formulation (LIM3) nn_limflx   = ', nn_limflx 
    127137         WRITE(numout,*) '           Misc. options of sbc : ' 
     
    150160      END SELECT 
    151161      ! 
    152 #if defined key_top && ! defined key_offline 
    153       ltrcdm2dc = (ln_dm2dc .AND. ln_blk_core .AND. nn_ice==2) 
    154       IF( ltrcdm2dc )THEN 
    155          IF(lwp)THEN 
    156             WRITE(numout,*)"analytical diurnal cycle, core bulk formulation and LIM2 use: " 
    157             WRITE(numout,*)"Diurnal cycle on physics but not in passive tracers" 
    158          ENDIF 
    159       ENDIF 
    160 #else  
    161       ltrcdm2dc =  .FALSE. 
    162 #endif 
    163  
    164       ! 
     162      IF ( nn_components /= jp_iam_nemo .AND. .NOT. lk_oasis )   & 
     163         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but key_oasis3 disabled' ) 
     164      IF ( nn_components == jp_iam_opa .AND. ln_cpl )   & 
     165         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but ln_cpl = T in OPA' ) 
     166      IF ( nn_components == jp_iam_opa .AND. ln_mixcpl )   & 
     167         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but ln_mixcpl = T in OPA' ) 
     168      IF ( ln_cpl .AND. .NOT. lk_oasis )    & 
     169         &      CALL ctl_stop( 'STOP', 'sbc_init : OASIS-coupled atmosphere model, but key_oasis3 disabled' ) 
     170      IF( ln_mixcpl .AND. .NOT. lk_oasis )    & 
     171         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) requires the cpp key key_oasis3' ) 
     172      IF( ln_mixcpl .AND. .NOT. ln_cpl )    & 
     173         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) requires ln_cpl = T' ) 
     174      IF( ln_mixcpl .AND. nn_components /= jp_iam_nemo )    & 
     175         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) is not yet working with sas-opa coupling via oasis' ) 
     176 
    165177      !                              ! allocate sbc arrays 
    166178      IF( sbc_oce_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_oce arrays' ) 
    167179 
    168180      !                          ! Checks: 
    169       IF( .NOT. ln_rnf ) THEN                      ! no specific treatment in vicinity of river mouths  
    170          ln_rnf_mouth  = .false.                       
    171          IF( sbc_rnf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_rnf arrays' ) 
    172          nkrnf         = 0 
    173          rnf     (:,:) = 0.0_wp 
    174          rnf_b   (:,:) = 0.0_wp 
    175          rnfmsk  (:,:) = 0.0_wp 
    176          rnfmsk_z(:)   = 0.0_wp 
    177       ENDIF 
    178       IF( nn_isf .EQ. 0 ) THEN                      ! no specific treatment in vicinity of ice shelf  
     181      IF( nn_isf .EQ. 0 ) THEN                      ! variable initialisation if no ice shelf  
    179182         IF( sbc_isf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' ) 
    180183         fwfisf  (:,:) = 0.0_wp 
     184         fwfisf_b(:,:) = 0.0_wp 
     185         rdivisf       = 0.0_wp 
    181186      END IF 
    182       IF( nn_ice == 0  )   fr_i(:,:) = 0.e0        ! no ice in the domain, ice fraction is always zero 
     187      IF( nn_ice == 0 .AND. nn_components /= jp_iam_opa )   fr_i(:,:) = 0.e0 ! no ice in the domain, ice fraction is always zero 
    183188 
    184189      sfx(:,:) = 0.0_wp                            ! the salt flux due to freezing/melting will be computed (i.e. will be non-zero)  
     
    190195 
    191196      !                                            ! restartability    
    192       IF( MOD( nitend - nit000 + 1, nn_fsbc) /= 0 .OR.   & 
    193           MOD( nstock             , nn_fsbc) /= 0 ) THEN  
    194          WRITE(ctmp1,*) 'experiment length (', nitend - nit000 + 1, ') or nstock (', nstock,   & 
    195             &           ' is NOT a multiple of nn_fsbc (', nn_fsbc, ')' 
    196          CALL ctl_stop( ctmp1, 'Impossible to properly do model restart' ) 
    197       ENDIF 
    198       ! 
    199       IF( MOD( rday, REAL(nn_fsbc, wp) * rdt ) /= 0 )   & 
    200          &  CALL ctl_warn( 'nn_fsbc is NOT a multiple of the number of time steps in a day' ) 
    201       ! 
    202       IF( ( nn_ice == 2 .OR. nn_ice ==3 ) .AND. .NOT.( ln_blk_clio .OR. ln_blk_core .OR. lk_cpl ) )   & 
     197      IF( ( nn_ice == 2 .OR. nn_ice ==3 ) .AND. .NOT.( ln_blk_clio .OR. ln_blk_core .OR. ln_cpl ) )   & 
    203198         &   CALL ctl_stop( 'LIM sea-ice model requires a bulk formulation or coupled configuration' ) 
    204       IF( nn_ice == 4 .AND. .NOT.( ln_blk_core .OR. lk_cpl ) )   & 
    205          &   CALL ctl_stop( 'CICE sea-ice model requires ln_blk_core or lk_cpl' ) 
     199      IF( nn_ice == 4 .AND. .NOT.( ln_blk_core .OR. ln_cpl ) )   & 
     200         &   CALL ctl_stop( 'CICE sea-ice model requires ln_blk_core or ln_cpl' ) 
    206201      IF( nn_ice == 4 .AND. lk_agrif )   & 
    207202         &   CALL ctl_stop( 'CICE sea-ice model not currently available with AGRIF' ) 
     
    210205      IF( ( nn_ice /= 3 ) .AND. ( nn_limflx >= 0 ) )   & 
    211206         &   WRITE(numout,*) 'The nn_limflx>=0 option has no effect if sea ice model is not LIM3' 
    212       IF( ( nn_ice == 3 ) .AND. ( lk_cpl ) .AND. ( ( nn_limflx == -1 ) .OR. ( nn_limflx == 1 ) ) )   & 
     207      IF( ( nn_ice == 3 ) .AND. ( ln_cpl ) .AND. ( ( nn_limflx == -1 ) .OR. ( nn_limflx == 1 ) ) )   & 
    213208         &   CALL ctl_stop( 'The chosen nn_limflx for LIM3 in coupled mode must be 0 or 2' ) 
    214       IF( ( nn_ice == 3 ) .AND. ( .NOT. lk_cpl ) .AND. ( nn_limflx == 2 ) )   & 
     209      IF( ( nn_ice == 3 ) .AND. ( .NOT. ln_cpl ) .AND. ( nn_limflx == 2 ) )   & 
    215210         &   CALL ctl_stop( 'The chosen nn_limflx for LIM3 in forced mode cannot be 2' ) 
    216211 
    217212      IF( ln_dm2dc )   nday_qsr = -1   ! initialisation flag 
    218213 
    219       IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) )   & 
     214      IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) .AND. nn_components /= jp_iam_opa )   & 
    220215         &   CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' ) 
    221216       
    222       IF( ln_dm2dc .AND. ( ( NINT(rday) / ( nn_fsbc * NINT(rdt) ) )  < 8 ) )   & 
    223          &   CALL ctl_warn( 'diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' ) 
    224  
    225217      IF ( ln_wave ) THEN 
    226218      !Activated wave module but neither drag nor stokes drift activated 
     
    236228         & asked coupling with drag coefficient (ln_cdgw =T) or Stokes drift (ln_sdw=T) ') 
    237229      ENDIF  
    238        
    239230      !                          ! Choice of the Surface Boudary Condition (set nsbc) 
     231      ll_purecpl = ln_cpl .AND. .NOT. ln_mixcpl 
     232      ! 
    240233      icpt = 0 
    241       IF( ln_ana          ) THEN   ;   nsbc = jp_ana    ; icpt = icpt + 1   ;   ENDIF       ! analytical      formulation 
    242       IF( ln_flx          ) THEN   ;   nsbc = jp_flx    ; icpt = icpt + 1   ;   ENDIF       ! flux            formulation 
    243       IF( ln_blk_clio     ) THEN   ;   nsbc = jp_clio   ; icpt = icpt + 1   ;   ENDIF       ! CLIO bulk       formulation 
    244       IF( ln_blk_core     ) THEN   ;   nsbc = jp_core   ; icpt = icpt + 1   ;   ENDIF       ! CORE bulk       formulation 
    245       IF( ln_blk_mfs      ) THEN   ;   nsbc = jp_mfs    ; icpt = icpt + 1   ;   ENDIF       ! MFS  bulk       formulation 
    246       IF( lk_cpl          ) THEN   ;   nsbc = jp_cpl    ; icpt = icpt + 1   ;   ENDIF       ! Coupled         formulation 
    247       IF( cp_cfg == 'gyre') THEN   ;   nsbc = jp_gyre                       ;   ENDIF       ! GYRE analytical formulation 
    248       IF( lk_esopa        )            nsbc = jp_esopa                                      ! esopa test, ALL formulations 
     234      IF( ln_ana          ) THEN   ;   nsbc = jp_ana     ; icpt = icpt + 1   ;   ENDIF       ! analytical           formulation 
     235      IF( ln_flx          ) THEN   ;   nsbc = jp_flx     ; icpt = icpt + 1   ;   ENDIF       ! flux                 formulation 
     236      IF( ln_blk_clio     ) THEN   ;   nsbc = jp_clio    ; icpt = icpt + 1   ;   ENDIF       ! CLIO bulk            formulation 
     237      IF( ln_blk_core     ) THEN   ;   nsbc = jp_core    ; icpt = icpt + 1   ;   ENDIF       ! CORE bulk            formulation 
     238      IF( ln_blk_mfs      ) THEN   ;   nsbc = jp_mfs     ; icpt = icpt + 1   ;   ENDIF       ! MFS  bulk            formulation 
     239      IF( ll_purecpl      ) THEN   ;   nsbc = jp_purecpl ; icpt = icpt + 1   ;   ENDIF       ! Pure Coupled         formulation 
     240      IF( cp_cfg == 'gyre') THEN   ;   nsbc = jp_gyre                        ;   ENDIF       ! GYRE analytical      formulation 
     241      IF( nn_components == jp_iam_opa )   & 
     242         &                  THEN   ;   nsbc = jp_none    ; icpt = icpt + 1   ;   ENDIF       ! opa coupling via SAS module 
     243      IF( lk_esopa        )            nsbc = jp_esopa                                       ! esopa test, ALL formulations 
    249244      ! 
    250245      IF( icpt /= 1 .AND. .NOT.lk_esopa ) THEN 
     
    257252      IF(lwp) THEN 
    258253         WRITE(numout,*) 
    259          IF( nsbc == jp_esopa )   WRITE(numout,*) '              ESOPA test All surface boundary conditions' 
    260          IF( nsbc == jp_gyre  )   WRITE(numout,*) '              GYRE analytical formulation' 
    261          IF( nsbc == jp_ana   )   WRITE(numout,*) '              analytical formulation' 
    262          IF( nsbc == jp_flx   )   WRITE(numout,*) '              flux formulation' 
    263          IF( nsbc == jp_clio  )   WRITE(numout,*) '              CLIO bulk formulation' 
    264          IF( nsbc == jp_core  )   WRITE(numout,*) '              CORE bulk formulation' 
    265          IF( nsbc == jp_cpl   )   WRITE(numout,*) '              coupled formulation' 
    266          IF( nsbc == jp_mfs   )   WRITE(numout,*) '              MFS Bulk formulation' 
    267       ENDIF 
    268       ! 
     254         IF( nsbc == jp_esopa   )   WRITE(numout,*) '              ESOPA test All surface boundary conditions' 
     255         IF( nsbc == jp_gyre    )   WRITE(numout,*) '              GYRE analytical formulation' 
     256         IF( nsbc == jp_ana     )   WRITE(numout,*) '              analytical formulation' 
     257         IF( nsbc == jp_flx     )   WRITE(numout,*) '              flux formulation' 
     258         IF( nsbc == jp_clio    )   WRITE(numout,*) '              CLIO bulk formulation' 
     259         IF( nsbc == jp_core    )   WRITE(numout,*) '              CORE bulk formulation' 
     260         IF( nsbc == jp_purecpl )   WRITE(numout,*) '              pure coupled formulation' 
     261         IF( nsbc == jp_mfs     )   WRITE(numout,*) '              MFS Bulk formulation' 
     262         IF( nsbc == jp_none    )   WRITE(numout,*) '              OPA coupled to SAS via oasis' 
     263         IF( ln_mixcpl          )   WRITE(numout,*) '              + forced-coupled mixed formulation' 
     264         IF( nn_components/= jp_iam_nemo )  & 
     265            &                       WRITE(numout,*) '              + OASIS coupled SAS' 
     266      ENDIF 
     267      ! 
     268      IF( lk_oasis )   CALL sbc_cpl_init (nn_ice)   ! OASIS initialisation. must be done before: (1) first time step 
     269      !                                                     !                                            (2) the use of nn_fsbc 
     270 
     271!     nn_fsbc initialization if OPA-SAS coupling via OASIS 
     272!     sas model time step has to be declared in OASIS (mandatory) -> nn_fsbc has to be modified accordingly 
     273      IF ( nn_components /= jp_iam_nemo ) THEN 
     274 
     275         IF ( nn_components == jp_iam_opa ) nn_fsbc = cpl_freq('O_SFLX') / NINT(rdt) 
     276         IF ( nn_components == jp_iam_sas ) nn_fsbc = cpl_freq('I_SFLX') / NINT(rdt) 
     277         ! 
     278         IF(lwp)THEN 
     279            WRITE(numout,*) 
     280            WRITE(numout,*)"   OPA-SAS coupled via OASIS : nn_fsbc re-defined from OASIS namcouple ", nn_fsbc 
     281            WRITE(numout,*) 
     282         ENDIF 
     283      ENDIF 
     284 
     285      IF( MOD( nitend - nit000 + 1, nn_fsbc) /= 0 .OR.   & 
     286          MOD( nstock             , nn_fsbc) /= 0 ) THEN  
     287         WRITE(ctmp1,*) 'experiment length (', nitend - nit000 + 1, ') or nstock (', nstock,   & 
     288            &           ' is NOT a multiple of nn_fsbc (', nn_fsbc, ')' 
     289         CALL ctl_stop( ctmp1, 'Impossible to properly do model restart' ) 
     290      ENDIF 
     291      ! 
     292      IF( MOD( rday, REAL(nn_fsbc, wp) * rdt ) /= 0 )   & 
     293         &  CALL ctl_warn( 'nn_fsbc is NOT a multiple of the number of time steps in a day' ) 
     294      ! 
     295      IF( ln_dm2dc .AND. ( ( NINT(rday) / ( nn_fsbc * NINT(rdt) ) )  < 8 ) )   & 
     296         &   CALL ctl_warn( 'diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' ) 
     297 
    269298                               CALL sbc_ssm_init               ! Sea-surface mean fields initialisation 
    270299      ! 
    271300      IF( ln_ssr           )   CALL sbc_ssr_init               ! Sea-Surface Restoring initialisation 
    272301      ! 
     302                               CALL sbc_rnf_init               ! Runof initialisation 
     303      ! 
     304      IF( nn_ice == 3      )   CALL sbc_lim_init               ! LIM3 initialisation 
     305 
    273306      IF( nn_ice == 4      )   CALL cice_sbc_init( nsbc )      ! CICE initialisation 
    274       ! 
    275       IF( nsbc   == jp_cpl )   CALL sbc_cpl_init (nn_ice)      ! OASIS initialisation. must be done before first time step 
    276  
     307       
    277308   END SUBROUTINE sbc_init 
    278309 
     
    314345      !                                            ! ---------------------------------------- ! 
    315346      ! 
    316       IF( ln_apr_dyn ) CALL sbc_apr( kt )                ! atmospheric pressure provided at kt+0.5*nn_fsbc 
     347      IF ( .NOT. lk_bdy ) then 
     348         IF( ln_apr_dyn ) CALL sbc_apr( kt )                ! atmospheric pressure provided at kt+0.5*nn_fsbc 
     349      ENDIF 
    317350                                                         ! (caution called before sbc_ssm) 
    318351      ! 
    319       CALL sbc_ssm( kt )                                 ! ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m) 
    320       !                                                  ! averaged over nf_sbc time-step 
     352      IF( nn_components /= jp_iam_sas )   CALL sbc_ssm( kt )   ! ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m) 
     353      !                                                        ! averaged over nf_sbc time-step 
    321354 
    322355      IF (ln_wave) CALL sbc_wave( kt ) 
     
    329362      CASE( jp_flx   )   ;   CALL sbc_flx     ( kt )                    ! flux formulation 
    330363      CASE( jp_clio  )   ;   CALL sbc_blk_clio( kt )                    ! bulk formulation : CLIO for the ocean 
    331       CASE( jp_core  )   ;   CALL sbc_blk_core( kt )                    ! bulk formulation : CORE for the ocean 
    332       CASE( jp_cpl   )   ;   CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! coupled formulation 
     364      CASE( jp_core  )    
     365         IF( nn_components == jp_iam_sas ) & 
     366            &                CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: SAS receiving fields from OPA  
     367                             CALL sbc_blk_core( kt )                    ! bulk formulation : CORE for the ocean 
     368                                                                        ! from oce: sea surface variables (sst_m, sss_m,  ssu_m,  ssv_m) 
     369      CASE( jp_purecpl )  ;  CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! pure coupled formulation 
     370                                                                        ! 
    333371      CASE( jp_mfs   )   ;   CALL sbc_blk_mfs ( kt )                    ! bulk formulation : MFS for the ocean 
     372      CASE( jp_none  )  
     373         IF( nn_components == jp_iam_opa ) & 
     374                             CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: OPA receiving fields from SAS 
    334375      CASE( jp_esopa )                                 
    335376                             CALL sbc_ana     ( kt )                    ! ESOPA, test ALL the formulations 
     
    341382      END SELECT 
    342383 
     384      IF( ln_mixcpl )        CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing 
     385 
     386 
    343387      !                                            !==  Misc. Options  ==! 
    344388       
     
    363407      !                                                           ! (update freshwater fluxes) 
    364408!RBbug do not understand why see ticket 667 
    365       !clem-bugsal CALL lbc_lnk( emp, 'T', 1. ) 
     409!clem: it looks like it is necessary for the north fold (in certain circumstances). Don't know why. 
     410      CALL lbc_lnk( emp, 'T', 1. ) 
    366411      ! 
    367412      IF( kt == nit000 ) THEN                          !   set the forcing field at nit000 - 1    ! 
     
    404449         ! CALL iom_rstput( kt, nitrst, numrow, 'qsr_b'  , qsr  ) 
    405450         CALL iom_rstput( kt, nitrst, numrow, 'emp_b'  , emp  ) 
    406          CALL iom_rstput( kt, nitrst, numrow, 'sfx_b' , sfx ) 
     451         CALL iom_rstput( kt, nitrst, numrow, 'sfx_b'  , sfx ) 
    407452      ENDIF 
    408453 
     
    419464         CALL iom_put( "qns"   , qns        )                   ! solar heat flux 
    420465         CALL iom_put( "qsr"   ,       qsr  )                   ! solar heat flux 
    421          IF( nn_ice > 0 )   CALL iom_put( "ice_cover", fr_i )   ! ice fraction  
     466         IF( nn_ice > 0 .OR. nn_components == jp_iam_opa )   CALL iom_put( "ice_cover", fr_i )   ! ice fraction  
    422467         CALL iom_put( "taum"  , taum       )                   ! wind stress module  
    423468         CALL iom_put( "wspd"  , wndm       )                   ! wind speed  module over free ocean or leads in presence of sea-ice 
Note: See TracChangeset for help on using the changeset viewer.