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 10009 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/diaar5.F90 – NEMO

Ignore:
Timestamp:
2018-07-29T11:23:51+02:00 (6 years ago)
Author:
gm
Message:

#1911 (ENHANCE-04): RK3 branch - step II.1 time-level dimension on ssh

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/diaar5.F90

    r9939 r10009  
    8989         ALLOCATE( zrhd(jpi,jpj,jpk) , zrhop(jpi,jpj,jpk) ) 
    9090         ALLOCATE( ztsn(jpi,jpj,jpk,jpts) ) 
    91          zarea_ssh(:,:) = area(:,:) * sshn(:,:) 
     91         zarea_ssh(:,:) = area(:,:) * ssh(:,:,Nnn) 
    9292      ENDIF 
    9393      ! 
     
    100100         CALL iom_put( 'voltot', zvol               ) 
    101101         CALL iom_put( 'sshtot', zvolssh / area_tot ) 
    102          CALL iom_put( 'sshdyn', sshn(:,:) - (zvolssh / area_tot) ) 
     102         CALL iom_put( 'sshdyn', ssh(:,:,Nnn) - (zvolssh / area_tot) ) 
    103103         ! 
    104104      ENDIF 
     
    118118               DO ji = 1, jpi 
    119119                  DO jj = 1, jpj 
    120                      zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
     120                     zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Nnn) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
    121121                  END DO 
    122122               END DO 
    123123            ELSE 
    124                zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
     124               zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Nnn) * zrhd(:,:,1) 
    125125            END IF 
    126126!!gm 
     
    147147               DO ji = 1,jpi 
    148148                  DO jj = 1,jpj 
    149                      zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
     149                     zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Nnn) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
    150150                  END DO 
    151151               END DO 
    152152            ELSE 
    153                zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
     153               zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Nnn) * zrhd(:,:,1) 
    154154            END IF 
    155155         END IF 
     
    162162         !                                         ! ocean bottom pressure 
    163163         zztmp = rho0 * grav * 1.e-4_wp               ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 
    164          zbotpres(:,:) = zztmp * ( zbotpres(:,:) + sshn(:,:) + thick0(:,:) ) 
     164         zbotpres(:,:) = zztmp * ( zbotpres(:,:) + ssh(:,:,Nnn) + thick0(:,:) ) 
    165165         CALL iom_put( 'botpres', zbotpres ) 
    166166         ! 
Note: See TracChangeset for help on using the changeset viewer.