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 5443 for branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90 – NEMO

Ignore:
Timestamp:
2015-06-19T17:18:00+02:00 (9 years ago)
Author:
davestorkey
Message:

Update 2015/dev_r5021_UKMO1_CICE_coupling branch to revision 5442 of the trunk.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r5234 r5443  
    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 
     
    8285      INTEGER ::   icpt   ! local integer 
    8386      !! 
    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 
     87      NAMELIST/namsbc/ nn_fsbc   , ln_ana    , ln_flx, ln_blk_clio, ln_blk_core, ln_mixcpl,   & 
     88         &             ln_blk_mfs, ln_apr_dyn, nn_ice, nn_ice_embd, ln_dm2dc   , ln_rnf   ,   & 
     89         &             ln_ssr    , nn_isf    , nn_fwb, ln_cdgw    , ln_wave    , ln_sdw   ,   & 
     90         &             nn_lsm    , nn_limflx , nn_components, ln_cpl 
    8791      INTEGER  ::   ios 
     92      INTEGER  ::   ierr, ierr0, ierr1, ierr2, ierr3, jpm 
     93      LOGICAL  ::   ll_purecpl 
    8894      !!---------------------------------------------------------------------- 
    8995 
     
    113119          nn_ice      =   0 
    114120      ENDIF 
    115       
     121 
    116122      IF(lwp) THEN               ! Control print 
    117123         WRITE(numout,*) '        Namelist namsbc (partly overwritten with CPP key setting)' 
     
    123129         WRITE(numout,*) '              CORE bulk  formulation                     ln_blk_core = ', ln_blk_core 
    124130         WRITE(numout,*) '              MFS  bulk  formulation                     ln_blk_mfs  = ', ln_blk_mfs 
    125          WRITE(numout,*) '              coupled    formulation (T if key_oasis3)   lk_cpl      = ', lk_cpl 
     131         WRITE(numout,*) '              ocean-atmosphere coupled formulation       ln_cpl      = ', ln_cpl 
     132         WRITE(numout,*) '              forced-coupled mixed formulation           ln_mixcpl   = ', ln_mixcpl 
     133         WRITE(numout,*) '              OASIS coupling (with atm or sas)           lk_oasis    = ', lk_oasis 
     134         WRITE(numout,*) '              components of your executable              nn_components = ', nn_components 
    126135         WRITE(numout,*) '              Multicategory heat flux formulation (LIM3) nn_limflx   = ', nn_limflx 
    127136         WRITE(numout,*) '           Misc. options of sbc : ' 
     
    150159      END SELECT 
    151160      ! 
    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       ! 
     161      IF ( nn_components /= jp_iam_nemo .AND. .NOT. lk_oasis )   & 
     162         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but key_oasis3 disabled' ) 
     163      IF ( nn_components == jp_iam_opa .AND. ln_cpl )   & 
     164         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but ln_cpl = T in OPA' ) 
     165      IF ( nn_components == jp_iam_opa .AND. ln_mixcpl )   & 
     166         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but ln_mixcpl = T in OPA' ) 
     167      IF ( ln_cpl .AND. .NOT. lk_oasis )    & 
     168         &      CALL ctl_stop( 'STOP', 'sbc_init : OASIS-coupled atmosphere model, but key_oasis3 disabled' ) 
     169      IF( ln_mixcpl .AND. .NOT. lk_oasis )    & 
     170         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) requires the cpp key key_oasis3' ) 
     171      IF( ln_mixcpl .AND. .NOT. ln_cpl )    & 
     172         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) requires ln_cpl = T' ) 
     173      IF( ln_mixcpl .AND. nn_components /= jp_iam_nemo )    & 
     174         &      CALL ctl_stop( 'the forced-coupled mixed mode (ln_mixcpl) is not yet working with sas-opa coupling via oasis' ) 
     175 
    165176      !                              ! allocate sbc arrays 
    166177      IF( sbc_oce_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_oce arrays' ) 
    167178 
    168179      !                          ! 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 
    178180      IF( nn_isf .EQ. 0 ) THEN                      ! no specific treatment in vicinity of ice shelf  
    179181         IF( sbc_isf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' ) 
    180182         fwfisf  (:,:) = 0.0_wp 
     183         fwfisf_b(:,:) = 0.0_wp 
    181184      END IF 
    182       IF( nn_ice == 0  )   fr_i(:,:) = 0.e0        ! no ice in the domain, ice fraction is always zero 
     185      IF( nn_ice == 0 .AND. nn_components /= jp_iam_opa )   fr_i(:,:) = 0.e0 ! no ice in the domain, ice fraction is always zero 
    183186 
    184187      sfx(:,:) = 0.0_wp                            ! the salt flux due to freezing/melting will be computed (i.e. will be non-zero)  
     
    190193 
    191194      !                                            ! 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 ) )   & 
     195      IF( ( nn_ice == 2 .OR. nn_ice ==3 ) .AND. .NOT.( ln_blk_clio .OR. ln_blk_core .OR. ln_cpl ) )   & 
    203196         &   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' ) 
     197      IF( nn_ice == 4 .AND. .NOT.( ln_blk_core .OR. ln_cpl ) )   & 
     198         &   CALL ctl_stop( 'CICE sea-ice model requires ln_blk_core or ln_cpl' ) 
    206199      IF( nn_ice == 4 .AND. lk_agrif )   & 
    207200         &   CALL ctl_stop( 'CICE sea-ice model not currently available with AGRIF' ) 
     
    210203      IF( ( nn_ice /= 3 ) .AND. ( nn_limflx >= 0 ) )   & 
    211204         &   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 ) ) )   & 
     205      IF( ( nn_ice == 3 ) .AND. ( ln_cpl ) .AND. ( ( nn_limflx == -1 ) .OR. ( nn_limflx == 1 ) ) )   & 
    213206         &   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 ) )   & 
     207      IF( ( nn_ice == 3 ) .AND. ( .NOT. ln_cpl ) .AND. ( nn_limflx == 2 ) )   & 
    215208         &   CALL ctl_stop( 'The chosen nn_limflx for LIM3 in forced mode cannot be 2' ) 
    216209 
    217210      IF( ln_dm2dc )   nday_qsr = -1   ! initialisation flag 
    218211 
    219       IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) )   & 
     212      IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) .AND. nn_components /= jp_iam_opa )   & 
    220213         &   CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' ) 
    221214       
    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  
    225215      IF ( ln_wave ) THEN 
    226216      !Activated wave module but neither drag nor stokes drift activated 
     
    236226         & asked coupling with drag coefficient (ln_cdgw =T) or Stokes drift (ln_sdw=T) ') 
    237227      ENDIF  
    238        
    239228      !                          ! Choice of the Surface Boudary Condition (set nsbc) 
     229      ll_purecpl = ln_cpl .AND. .NOT. ln_mixcpl 
     230      ! 
    240231      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 
     232      IF( ln_ana          ) THEN   ;   nsbc = jp_ana     ; icpt = icpt + 1   ;   ENDIF       ! analytical           formulation 
     233      IF( ln_flx          ) THEN   ;   nsbc = jp_flx     ; icpt = icpt + 1   ;   ENDIF       ! flux                 formulation 
     234      IF( ln_blk_clio     ) THEN   ;   nsbc = jp_clio    ; icpt = icpt + 1   ;   ENDIF       ! CLIO bulk            formulation 
     235      IF( ln_blk_core     ) THEN   ;   nsbc = jp_core    ; icpt = icpt + 1   ;   ENDIF       ! CORE bulk            formulation 
     236      IF( ln_blk_mfs      ) THEN   ;   nsbc = jp_mfs     ; icpt = icpt + 1   ;   ENDIF       ! MFS  bulk            formulation 
     237      IF( ll_purecpl      ) THEN   ;   nsbc = jp_purecpl ; icpt = icpt + 1   ;   ENDIF       ! Pure Coupled         formulation 
     238      IF( cp_cfg == 'gyre') THEN   ;   nsbc = jp_gyre                        ;   ENDIF       ! GYRE analytical      formulation 
     239      IF( nn_components == jp_iam_opa )   & 
     240         &                  THEN   ;   nsbc = jp_none    ; icpt = icpt + 1   ;   ENDIF       ! opa coupling via SAS module 
     241      IF( lk_esopa        )            nsbc = jp_esopa                                       ! esopa test, ALL formulations 
    249242      ! 
    250243      IF( icpt /= 1 .AND. .NOT.lk_esopa ) THEN 
     
    257250      IF(lwp) THEN 
    258251         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       ! 
     252         IF( nsbc == jp_esopa   )   WRITE(numout,*) '              ESOPA test All surface boundary conditions' 
     253         IF( nsbc == jp_gyre    )   WRITE(numout,*) '              GYRE analytical formulation' 
     254         IF( nsbc == jp_ana     )   WRITE(numout,*) '              analytical formulation' 
     255         IF( nsbc == jp_flx     )   WRITE(numout,*) '              flux formulation' 
     256         IF( nsbc == jp_clio    )   WRITE(numout,*) '              CLIO bulk formulation' 
     257         IF( nsbc == jp_core    )   WRITE(numout,*) '              CORE bulk formulation' 
     258         IF( nsbc == jp_purecpl )   WRITE(numout,*) '              pure coupled formulation' 
     259         IF( nsbc == jp_mfs     )   WRITE(numout,*) '              MFS Bulk formulation' 
     260         IF( nsbc == jp_none    )   WRITE(numout,*) '              OPA coupled to SAS via oasis' 
     261         IF( ln_mixcpl          )   WRITE(numout,*) '              + forced-coupled mixed formulation' 
     262         IF( nn_components/= jp_iam_nemo )  & 
     263            &                       WRITE(numout,*) '              + OASIS coupled SAS' 
     264      ENDIF 
     265      ! 
     266      IF( lk_oasis )   CALL sbc_cpl_init (nn_ice)   ! OASIS initialisation. must be done before: (1) first time step 
     267      !                                                     !                                            (2) the use of nn_fsbc 
     268 
     269!     nn_fsbc initialization if OPA-SAS coupling via OASIS 
     270!     sas model time step has to be declared in OASIS (mandatory) -> nn_fsbc has to be modified accordingly 
     271      IF ( nn_components /= jp_iam_nemo ) THEN 
     272 
     273         IF ( nn_components == jp_iam_opa ) nn_fsbc = cpl_freq('O_SFLX') / NINT(rdt) 
     274         IF ( nn_components == jp_iam_sas ) nn_fsbc = cpl_freq('I_SFLX') / NINT(rdt) 
     275         ! 
     276         IF(lwp)THEN 
     277            WRITE(numout,*) 
     278            WRITE(numout,*)"   OPA-SAS coupled via OASIS : nn_fsbc re-defined from OASIS namcouple ", nn_fsbc 
     279            WRITE(numout,*) 
     280         ENDIF 
     281      ENDIF 
     282 
     283      IF( MOD( nitend - nit000 + 1, nn_fsbc) /= 0 .OR.   & 
     284          MOD( nstock             , nn_fsbc) /= 0 ) THEN  
     285         WRITE(ctmp1,*) 'experiment length (', nitend - nit000 + 1, ') or nstock (', nstock,   & 
     286            &           ' is NOT a multiple of nn_fsbc (', nn_fsbc, ')' 
     287         CALL ctl_stop( ctmp1, 'Impossible to properly do model restart' ) 
     288      ENDIF 
     289      ! 
     290      IF( MOD( rday, REAL(nn_fsbc, wp) * rdt ) /= 0 )   & 
     291         &  CALL ctl_warn( 'nn_fsbc is NOT a multiple of the number of time steps in a day' ) 
     292      ! 
     293      IF( ln_dm2dc .AND. ( ( NINT(rday) / ( nn_fsbc * NINT(rdt) ) )  < 8 ) )   & 
     294         &   CALL ctl_warn( 'diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' ) 
     295 
    269296                               CALL sbc_ssm_init               ! Sea-surface mean fields initialisation 
    270297      ! 
    271298      IF( ln_ssr           )   CALL sbc_ssr_init               ! Sea-Surface Restoring initialisation 
    272299      ! 
     300                               CALL sbc_rnf_init               ! Runof initialisation 
     301      ! 
     302      IF( nn_ice == 3      )   CALL sbc_lim_init               ! LIM3 initialisation 
     303 
    273304      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  
     305       
    277306   END SUBROUTINE sbc_init 
    278307 
     
    317346                                                         ! (caution called before sbc_ssm) 
    318347      ! 
    319       CALL sbc_ssm( kt )                                 ! ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m) 
    320       !                                                  ! averaged over nf_sbc time-step 
     348      IF( nn_components /= jp_iam_sas )   CALL sbc_ssm( kt )   ! ocean sea surface variables (sst_m, sss_m, ssu_m, ssv_m) 
     349      !                                                        ! averaged over nf_sbc time-step 
    321350 
    322351      IF (ln_wave) CALL sbc_wave( kt ) 
     
    329358      CASE( jp_flx   )   ;   CALL sbc_flx     ( kt )                    ! flux formulation 
    330359      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 
     360      CASE( jp_core  )    
     361         IF( nn_components == jp_iam_sas ) & 
     362            &                CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: SAS receiving fields from OPA  
     363                             CALL sbc_blk_core( kt )                    ! bulk formulation : CORE for the ocean 
     364                                                                        ! from oce: sea surface variables (sst_m, sss_m,  ssu_m,  ssv_m) 
     365      CASE( jp_purecpl )  ;  CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! pure coupled formulation 
     366                                                                        ! 
    333367      CASE( jp_mfs   )   ;   CALL sbc_blk_mfs ( kt )                    ! bulk formulation : MFS for the ocean 
     368      CASE( jp_none  )  
     369         IF( nn_components == jp_iam_opa ) & 
     370                             CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: OPA receiving fields from SAS 
    334371      CASE( jp_esopa )                                 
    335372                             CALL sbc_ana     ( kt )                    ! ESOPA, test ALL the formulations 
     
    341378      END SELECT 
    342379 
     380      IF( ln_mixcpl )        CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing 
     381 
     382 
    343383      !                                            !==  Misc. Options  ==! 
    344384       
     
    363403      !                                                           ! (update freshwater fluxes) 
    364404!RBbug do not understand why see ticket 667 
    365       !clem-bugsal CALL lbc_lnk( emp, 'T', 1. ) 
     405!clem: it looks like it is necessary for the north fold (in certain circumstances). Don't know why. 
     406      CALL lbc_lnk( emp, 'T', 1. ) 
    366407      ! 
    367408      IF( kt == nit000 ) THEN                          !   set the forcing field at nit000 - 1    ! 
     
    404445         ! CALL iom_rstput( kt, nitrst, numrow, 'qsr_b'  , qsr  ) 
    405446         CALL iom_rstput( kt, nitrst, numrow, 'emp_b'  , emp  ) 
    406          CALL iom_rstput( kt, nitrst, numrow, 'sfx_b' , sfx ) 
     447         CALL iom_rstput( kt, nitrst, numrow, 'sfx_b'  , sfx ) 
    407448      ENDIF 
    408449 
     
    419460         CALL iom_put( "qns"   , qns        )                   ! solar heat flux 
    420461         CALL iom_put( "qsr"   ,       qsr  )                   ! solar heat flux 
    421          IF( nn_ice > 0 )   CALL iom_put( "ice_cover", fr_i )   ! ice fraction  
     462         IF( nn_ice > 0 .OR. nn_components == jp_iam_opa )   CALL iom_put( "ice_cover", fr_i )   ! ice fraction  
    422463         CALL iom_put( "taum"  , taum       )                   ! wind stress module  
    423464         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.