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 15516 – NEMO

Changeset 15516


Ignore:
Timestamp:
2021-11-16T15:40:29+01:00 (6 months ago)
Author:
vancop
Message:

fixes in diagnostics

Location:
NEMO/branches/2020/SI3-03_VP_rheology/src/ICE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_evp.F90

    r15499 r15516  
    134134      REAL(wp) ::   zintb, zintn                                        ! dummy argument 
    135135      REAL(wp) ::   zfac_x, zfac_y 
    136       REAL(wp) ::   zshear, zdum1, zdum2 
     136      REAL(wp) ::   zdum1, zdum2 
    137137      ! 
    138138      REAL(wp), DIMENSION(jpi,jpj) ::   zp_delt                         ! P/delta at T points 
    139       REAL(wp), DIMENSION(jpi,jpj) ::   zdeltastar_t                    ! Delta* at T-points 
    140       REAL(wp), DIMENSION(jpi,jpj) ::   zten_i                          ! Tension 
     139      REAL(wp), DIMENSION(jpi,jpj) ::   zdeltat                         ! Delta at T-points 
     140      REAL(wp), DIMENSION(jpi,jpj) ::   ztens, zshear                   ! Tension and shear 
    141141      REAL(wp), DIMENSION(jpi,jpj) ::   zbeta                           ! beta coef from Kimmritz 2017 
    142142      ! 
     
    410410               zdelta = SQRT( zdiv2 + ( zdt2 + zds2 ) * z1_ecc2 )   
    411411                
    412                zdeltastar_t(ji,jj) = zdelta + rn_creepl                     ! store delta* at previous iterate (for ellipse computation) 
     412               zdeltat(ji,jj) = zdelta  ! store delta at previous iterate (for ellipse computation) 
    413413 
    414414               ! P/delta at T points 
     
    758758            zdt2 = zdt * zdt 
    759759             
    760             zten_i(ji,jj) = zdt 
     760            ztens(ji,jj) = zdt 
    761761             
    762762            ! shear**2 at T points (doc eq. A16) 
     
    764764               &   + zds(ji,jj-1) * zds(ji,jj-1) * e1e2f(ji,jj-1) + zds(ji-1,jj-1) * zds(ji-1,jj-1) * e1e2f(ji-1,jj-1)  & 
    765765               &   ) * 0.25_wp * r1_e1e2t(ji,jj) 
     766 
     767            zshear(ji,jj) = SQRT( zds2 ) 
    766768             
    767             ! shear at T points 
     769            ! maximum shear rate at T points (includes tension, output only) 
    768770            pshear_i(ji,jj) = SQRT( zdt2 + zds2 ) 
    769771 
     
    827829               zfac             =   strength(ji,jj) / ( pdelta_i(ji,jj) + rn_creepl )  
    828830               zsig1            =   zfac * ( pdivu_i(ji,jj) - pdelta_i(ji,jj) ) 
    829                zsig2            =   zfac * z1_ecc2 * zten_i(ji,jj) 
    830                zsig12           =   zfac * z1_ecc2 * pshear_i(ji,jj) 
     831               zsig2            =   zfac * z1_ecc2 * ztens(ji,jj) 
     832               zsig12           =   zfac * z1_ecc2 * zshear(ji,jj) * 0.5_wp 
    831833                
    832834               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008) 
     
    863865               !                        following Lemieux & Dupont (2020) 
    864866               zfac             =   zp_delt(ji,jj) 
    865                zsig1            =   zfac * ( pdivu_i(ji,jj) - zdeltastar_t(ji,jj) ) 
    866                zsig2            =   zfac * z1_ecc2 * zten_i(ji,jj) 
    867                zsig12           =   zfac * z1_ecc2 * pshear_i(ji,jj) 
     867               zsig1            =   zfac * ( pdivu_i(ji,jj) - zdeltat(ji,jj) ) 
     868               zsig2            =   zfac * z1_ecc2 * ztens(ji,jj) 
     869               zsig12           =   zfac * z1_ecc2 * zshear(ji,jj) * 0.5_wp 
    868870                
    869871               ! Stress invariants (sigma_I, sigma_II, Coon 1974, Feltham 2008), T-point 
     
    10041006 
    10051007      ! time 
    1006       it = ( kt - 1 ) * kitermax + kiter 
     1008      it = ( kt - nit000 ) * kitermax + kiter 
    10071009       
    10081010      ! convergence 
     
    10241026         istatus = NF90_PUT_VAR( ncvgid, nvarid, (/zresm/), (/it/), (/1/) ) 
    10251027         ! close file 
    1026          IF( kt == nitend )   istatus = NF90_CLOSE(ncvgid) 
     1028         IF( kt == nitend - nn_fsbc + 1 .AND. kiter == kitermax )    istatus = NF90_CLOSE( ncvgid ) 
    10271029      ENDIF 
    10281030       
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90

    r15515 r15516  
    11281128 
    11291129               CALL rhg_cvg_vp( kt, i_out, i_inn, i_inn_tot, nn_nout_vp, nn_ninn_vp, nn_nvp,        & 
    1130                       &         u_ice, v_ice, zu_b, zv_b, zu_b_outer, zv_b_outer,                   & 
     1130                      &         u_ice, v_ice, zu_b, zv_b, zu_c, zv_c_outer                 ,        & 
    11311131                      &         zmt, za_iU, za_iV, zuerr_max, zverr_max, zglob_area,                & 
    11321132                      &         zrhsu, zAU, zBU, zCU, zDU, zEU, zFU,                                & 
Note: See TracChangeset for help on using the changeset viewer.