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

Changeset 2773


Ignore:
Timestamp:
2011-06-01T12:49:52+02:00 (13 years ago)
Author:
cetlod
Message:

PISCES trunk, bug in p4zrem.F90 : wrk_3d_2 is used twice

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90

    r2715 r2773  
    6161      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released 
    6262      USE wrk_nemo, ONLY:   ztempbac => wrk_2d_1 
    63       USE wrk_nemo, ONLY:   zdepbac  => wrk_3d_2 , zfesatur => wrk_3d_2 , zolimi => wrk_3d_4 
     63      USE wrk_nemo, ONLY:   zdepbac  => wrk_3d_2 , zolimi => wrk_3d_3 
    6464      ! 
    6565      INTEGER, INTENT(in) ::   kt ! ocean time step 
     
    6767      INTEGER  ::   ji, jj, jk 
    6868      REAL(wp) ::   zremip, zremik , zlam1b 
    69       REAL(wp) ::   zkeq  , zfeequi, zsiremin 
     69      REAL(wp) ::   zkeq  , zfeequi, zsiremin, zfesatur 
    7070      REAL(wp) ::   zsatur, zsatur2, znusil 
    7171      REAL(wp) ::   zbactfer, zorem, zorem2, zofer 
     
    7878      !!--------------------------------------------------------------------- 
    7979 
    80       IF(  wrk_in_use(2, 1)  .OR.  wrk_in_use(3, 2,3,4)  ) THEN 
     80      IF(  wrk_in_use(2, 1)  .OR.  wrk_in_use(3, 2,3)  ) THEN 
    8181         CALL ctl_stop('p4z_rem: requested workspace arrays unavailable')   ;   RETURN 
    8282      ENDIF 
     
    8484       ! Initialisation of temprary arrys 
    8585       zdepbac (:,:,:) = 0._wp 
    86        zfesatur(:,:,:) = 0._wp 
    8786       zolimi  (:,:,:) = 0._wp 
    8887       ztempbac(:,:)   = 0._wp 
     
    276275               zsatur  = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 
    277276               zsatur  = MAX( rtrn, zsatur ) 
    278                zsatur2 = zsatur * ( 1. + tn(ji,jj,jk) / 400.)**4 
    279                znusil  = 0.225  * ( 1. + tn(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9 
     277               zsatur2 = zsatur * ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**4 
     278               znusil  = 0.225  * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2**9 
    280279               zsiremin = xsirem * zstep * znusil 
    281280               zosil = zsiremin * trn(ji,jj,jk,jpdsi) 
     
    294293       ENDIF 
    295294 
    296       zfesatur(:,:,:) = 0.6e-9 
     295      zfesatur = 0.6e-9 
    297296!CDIR NOVERRCHK 
    298297      DO jk = 1, jpkm1 
     
    323322               !     ------------------------------------------------------------ 
    324323               zkeq = fekeq(ji,jj,jk) 
    325                zfeequi = ( -( 1. + zfesatur(ji,jj,jk) * zkeq - zkeq * trn(ji,jj,jk,jpfer) )               & 
    326                   &        + SQRT( ( 1. + zfesatur(ji,jj,jk) * zkeq - zkeq * trn(ji,jj,jk,jpfer) )**2       & 
     324               zfeequi = ( -( 1. + zfesatur * zkeq - zkeq * trn(ji,jj,jk,jpfer) )               & 
     325                  &        + SQRT( ( 1. + zfesatur * zkeq - zkeq * trn(ji,jj,jk,jpfer) )**2       & 
    327326                  &               + 4. * trn(ji,jj,jk,jpfer) * zkeq) ) / ( 2. * zkeq ) 
    328327 
     
    395394      ! 
    396395      IF(  wrk_not_released(2, 1)     .OR.   & 
    397            wrk_not_released(3, 2,3,4)  )   CALL ctl_stop('p4z_rem: failed to release workspace arrays') 
     396           wrk_not_released(3, 2,3)  )   CALL ctl_stop('p4z_rem: failed to release workspace arrays') 
    398397      ! 
    399398   END SUBROUTINE p4z_rem 
Note: See TracChangeset for help on using the changeset viewer.