Ignore:
Timestamp:
2013-08-14T10:54:56+02:00 (8 years ago)
Author:
rfurner
Message:

Unnecessary lbc_lnk calls removed from domzgr and precision corrected in some lines, see ticket

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r3998 r4006  
    12531253         DO jj = 1, jpj 
    12541254            DO ji = 1, jpi 
    1255                ztaper = EXP( -(gphit(ji,jj)/8._wp)**2 ) 
     1255               ztaper = EXP( -(gphit(ji,jj)/8._wp)**2._wp ) 
    12561256               hbatt(ji,jj) = rn_sbot_max * ztaper + hbatt(ji,jj) * ( 1._wp - ztaper ) 
    12571257            END DO 
     
    13681368      fsde3w(:,:,:) = gdep3w(:,:,:) 
    13691369      ! 
    1370       where (e3t   (:,:,:).eq.0.0)  e3t(:,:,:) = 1.0 
    1371       where (e3u   (:,:,:).eq.0.0)  e3u(:,:,:) = 1.0 
    1372       where (e3v   (:,:,:).eq.0.0)  e3v(:,:,:) = 1.0 
    1373       where (e3f   (:,:,:).eq.0.0)  e3f(:,:,:) = 1.0 
    1374       where (e3w   (:,:,:).eq.0.0)  e3w(:,:,:) = 1.0 
    1375       where (e3uw  (:,:,:).eq.0.0)  e3uw(:,:,:) = 1.0 
    1376       where (e3vw  (:,:,:).eq.0.0)  e3vw(:,:,:) = 1.0 
     1370      where (e3t   (:,:,:).eq.0.0)  e3t(:,:,:) = 1._wp 
     1371      where (e3u   (:,:,:).eq.0.0)  e3u(:,:,:) = 1._wp 
     1372      where (e3v   (:,:,:).eq.0.0)  e3v(:,:,:) = 1._wp 
     1373      where (e3f   (:,:,:).eq.0.0)  e3f(:,:,:) = 1._wp 
     1374      where (e3w   (:,:,:).eq.0.0)  e3w(:,:,:) = 1._wp 
     1375      where (e3uw  (:,:,:).eq.0.0)  e3uw(:,:,:) = 1._wp 
     1376      where (e3vw  (:,:,:).eq.0.0)  e3vw(:,:,:) = 1._wp 
    13771377 
    13781378#if defined key_agrif 
     
    17521752      ENDDO 
    17531753      ! 
    1754       CALL lbc_lnk(e3t ,'T',1.) ; CALL lbc_lnk(e3u ,'T',1.) 
    1755       CALL lbc_lnk(e3v ,'T',1.) ; CALL lbc_lnk(e3f ,'T',1.) 
    1756       CALL lbc_lnk(e3w ,'T',1.) 
    1757       CALL lbc_lnk(e3uw,'T',1.) ; CALL lbc_lnk(e3vw,'T',1.) 
    1758       ! 
    17591754      !                                               ! ============= 
    17601755 
     
    18531848      !!---------------------------------------------------------------------- 
    18541849      ! 
    1855       pf =   (   TANH( rn_theta * ( -(pk-0.5_wp) / REAL(jpkm1) + rn_thetb )  )   & 
     1850      pf =   (   TANH( rn_theta * ( -(pk-0.5_wp) / REAL(jpkm1,wp) + rn_thetb )  )   & 
    18561851         &     - TANH( rn_thetb * rn_theta                                )  )   & 
    18571852         & * (   COSH( rn_theta                           )                      & 
     
    18791874      ! 
    18801875      IF ( rn_theta == 0 ) then      ! uniform sigma 
    1881          pf1 = - ( pk1 - 0.5_wp ) / REAL( jpkm1 ) 
     1876         pf1 = - ( pk1 - 0.5_wp ) / REAL( jpkm1,wp ) 
    18821877      ELSE                        ! stretched sigma 
    1883          pf1 =   ( 1._wp - pbb ) * ( SINH( rn_theta*(-(pk1-0.5_wp)/REAL(jpkm1)) ) ) / SINH( rn_theta )              & 
    1884             &  + pbb * (  (TANH( rn_theta*( (-(pk1-0.5_wp)/REAL(jpkm1)) + 0.5_wp) ) - TANH( 0.5_wp * rn_theta )  )  & 
     1878         pf1 =   ( 1._wp - pbb ) * ( SINH( rn_theta*(-(pk1-0.5_wp)/REAL(jpkm1,wp)) ) ) / SINH( rn_theta )              & 
     1879            &  + pbb * (  (TANH( rn_theta*( (-(pk1-0.5_wp)/REAL(jpkm1,wp)) + 0.5_wp) ) - TANH( 0.5_wp * rn_theta )  )  & 
    18851880            &        / ( 2._wp * TANH( 0.5_wp * rn_theta ) )  ) 
    18861881      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.