Changeset 8751


Ignore:
Timestamp:
2017-11-20T13:48:24+01:00 (3 years ago)
Author:
dancopsey
Message:

Added sending of first-order ice concentration in sbccpl for JULES semi-implicit
coupling of atmosphere-ice fluxes.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r8738 r8751  
    148148   INTEGER, PARAMETER ::   jps_ocyw   = 31   ! currents on grid 2 
    149149   INTEGER, PARAMETER ::   jps_wlev   = 32   ! water level  
    150    INTEGER, PARAMETER ::   jpsnd      = 32   ! total number of fields sent  
     150   INTEGER, PARAMETER ::   jps_fice1  = 33   ! first-order ice concentration (for semi-implicit coupling of atmos-ice fluxes) 
     151   INTEGER, PARAMETER ::   jpsnd      = 33   ! total number of fields sent  
    151152 
    152153   !                                  !!** namelist namsbc_cpl ** 
     
    159160   END TYPE FLD_C 
    160161   !                                   ! Send to the atmosphere   
    161    TYPE(FLD_C) ::   sn_snd_temp, sn_snd_alb, sn_snd_thick, sn_snd_crt, sn_snd_co2                         
     162   TYPE(FLD_C) ::   sn_snd_temp, sn_snd_alb, sn_snd_thick, sn_snd_crt, sn_snd_co2, sn_snd_thick1                         
    162163   !                                   ! Received from the atmosphere 
    163164   TYPE(FLD_C) ::   sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau, sn_rcv_dqnsdt, sn_rcv_qsr, sn_rcv_qns, sn_rcv_emp, sn_rcv_rnf 
     
    237238      NAMELIST/namsbc_cpl/  sn_snd_temp , sn_snd_alb  , sn_snd_thick , sn_snd_crt   , sn_snd_co2,      &  
    238239         &                  sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau   , sn_rcv_dqnsdt, sn_rcv_qsr,      &  
     240         &                  sn_snd_thick1,                                                             & 
    239241         &                  sn_snd_ifrac, sn_snd_crtw , sn_snd_wlev  , sn_rcv_hsig  , sn_rcv_phioc ,   &  
    240242         &                  sn_rcv_sdrfx, sn_rcv_sdrfy, sn_rcv_wper  , sn_rcv_wnum  , sn_rcv_wstrf ,   & 
     
    747749      ssnd(jps_hice)%clname = 'OIceTck' 
    748750      ssnd(jps_hsnw)%clname = 'OSnwTck' 
     751      ssnd(jps_fice1)%clname = 'OIceFrd' 
    749752      IF( k_ice /= 0 ) THEN 
    750753         ssnd(jps_fice)%laction = .TRUE.                  ! if ice treated in the ocean (even in climato case) 
     754         ssnd(jps_fice1)%laction = .TRUE.                 ! First-order regridded ice concentration, to be used 
     755                                                          ! producing atmos-to-ice fluxes 
    751756! Currently no namelist entry to determine sending of multi-category ice fraction so use the thickness entry for now 
    752757         IF ( TRIM( sn_snd_thick%clcat ) == 'yes' ) ssnd(jps_fice)%nct = jpl 
     758         IF ( TRIM( sn_snd_thick1%clcat ) == 'yes' ) ssnd(jps_fice1)%nct = jpl 
    753759      ENDIF 
    754760       
     
    20962102         END SELECT 
    20972103         IF( ssnd(jps_fice)%laction )   CALL cpl_snd( jps_fice, isec, ztmp3, info ) 
     2104      ENDIF 
     2105 
     2106      IF( ssnd(jps_fice1)%laction ) THEN 
     2107         SELECT CASE( sn_snd_thick1%clcat ) 
     2108         CASE( 'yes' )   ;   ztmp3(:,:,1:jpl) =  a_i(:,:,1:jpl) 
     2109         CASE( 'no'  )   ;   ztmp3(:,:,1    ) = fr_i(:,:      ) 
     2110         CASE default    ;   CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_thick1%clcat' ) 
     2111         END SELECT 
     2112         CALL cpl_snd( jps_fice1, isec, ztmp3, info ) 
    20982113      ENDIF 
    20992114       
Note: See TracChangeset for help on using the changeset viewer.