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 9049 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rhg_evp.F90 – NEMO

Ignore:
Timestamp:
2017-12-14T14:55:17+01:00 (5 years ago)
Author:
clem
Message:

group one communcation in diags

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rhg_evp.F90

    r9019 r9049  
    112112      INTEGER ::   jter         ! local integers 
    113113      ! 
    114       REAL(wp) ::   zrhoco                                                   ! rau0 * rn_cio 
    115       REAL(wp) ::   zdtevp, z1_dtevp                                         ! time step for subcycling 
    116       REAL(wp) ::   ecc2, z1_ecc2                                            ! square of yield ellipse eccenticity 
    117       REAL(wp) ::   zalph1, z1_alph1, zalph2, z1_alph2                       ! alpha coef from Bouillon 2009 or Kimmritz 2017 
    118       REAL(wp) ::   zm1, zm2, zm3, zmassU, zmassV                            ! ice/snow mass 
    119       REAL(wp) ::   zdelta, zp_delf, zds2, zdt, zdt2, zdiv, zdiv2            ! temporary scalars 
    120       REAL(wp) ::   zTauO, zTauB, zTauE, zvel                                ! temporary scalars 
    121       ! 
    122       REAL(wp) ::   zresm                                                    ! Maximal error on ice velocity 
    123       REAL(wp) ::   zintb, zintn                                             ! dummy argument 
     114      REAL(wp) ::   zrhoco                                              ! rau0 * rn_cio 
     115      REAL(wp) ::   zdtevp, z1_dtevp                                    ! time step for subcycling 
     116      REAL(wp) ::   ecc2, z1_ecc2                                       ! square of yield ellipse eccenticity 
     117      REAL(wp) ::   zalph1, z1_alph1, zalph2, z1_alph2                  ! alpha coef from Bouillon 2009 or Kimmritz 2017 
     118      REAL(wp) ::   zm1, zm2, zm3, zmassU, zmassV                       ! ice/snow mass 
     119      REAL(wp) ::   zdelta, zp_delf, zds2, zdt, zdt2, zdiv, zdiv2       ! temporary scalars 
     120      REAL(wp) ::   zTauO, zTauB, zTauE, zvel                           ! temporary scalars 
     121      ! 
     122      REAL(wp) ::   zresm                                               ! Maximal error on ice velocity 
     123      REAL(wp) ::   zintb, zintn                                        ! dummy argument 
    124124      REAL(wp) ::   zfac_x, zfac_y 
    125125      REAL(wp) ::   zshear, zdum1, zdum2 
    … …  
    142142      REAL(wp), DIMENSION(jpi,jpj) ::   zu_ice, zv_ice, zresr           ! check convergence 
    143143      REAL(wp), DIMENSION(jpi,jpj) ::   zpice                           ! array used for the calculation of ice surface slope: 
    144       !                                                                      !   ocean surface (ssh_m) if ice is not embedded 
    145       !                                                                       !   ice top surface if ice is embedded    
     144      !                                                                 !    ocean surface (ssh_m) if ice is not embedded 
     145      !                                                                 !    ice top surface if ice is embedded    
    146146      REAL(wp), DIMENSION(jpi,jpj) ::   zCorx, zCory                    ! Coriolis stress array 
    147147      REAL(wp), DIMENSION(jpi,jpj) ::   ztaux_oi, ztauy_oi              ! Ocean-to-ice stress array 
    … …  
    810810         END DO 
    811811          
    812          CALL lbc_lnk_multi(   zdiag_sig1   , 'T',  1., zdiag_sig2   , 'T',  1.,   & 
    813             &                  zdiag_dssh_dx, 'U', -1., zdiag_dssh_dy, 'V', -1.,   & 
    814             &                  zdiag_corstrx, 'U', -1., zdiag_corstry, 'V', -1.,   &  
    815             &                  zdiag_intstrx, 'U', -1., zdiag_intstry, 'V', -1.    ) 
    816           
    817          CALL lbc_lnk_multi(   zdiag_utau_oi, 'U', -1., zdiag_vtau_oi, 'V', -1.    ) 
    818           
    819          CALL lbc_lnk_multi(   zdiag_xmtrp_ice, 'U', -1., zdiag_xmtrp_snw, 'U', -1.,   & 
    820             &                  zdiag_xatrp    , 'U', -1., zdiag_ymtrp_ice, 'V', -1.,   & 
    821             &                  zdiag_ymtrp_snw, 'V', -1., zdiag_yatrp    , 'V', -1.    ) 
     812         CALL lbc_lnk_multi( zdiag_sig1   , 'T',  1., zdiag_sig2   , 'T',  1.,   & 
     813            &                zdiag_dssh_dx, 'U', -1., zdiag_dssh_dy, 'V', -1.,   & 
     814            &                zdiag_corstrx, 'U', -1., zdiag_corstry, 'V', -1.,   &  
     815            &                zdiag_intstrx, 'U', -1., zdiag_intstry, 'V', -1.    ) 
     816                   
     817         CALL lbc_lnk_multi( zdiag_utau_oi  , 'U', -1., zdiag_vtau_oi  , 'V', -1.,   & 
     818            &                zdiag_xmtrp_ice, 'U', -1., zdiag_xmtrp_snw, 'U', -1.,   & 
     819            &                zdiag_xatrp    , 'U', -1., zdiag_ymtrp_ice, 'V', -1.,   & 
     820            &                zdiag_ymtrp_snw, 'V', -1., zdiag_yatrp    , 'V', -1.    ) 
    822821          
    823822         IF( iom_use('normstr' ) )   CALL iom_put( 'normstr'  ,  zdiag_sig1(:,:)      )   ! Normal stress 
Note: See TracChangeset for help on using the changeset viewer.