- Timestamp:
- 2018-07-29T11:23:51+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/diaar5.F90
r9939 r10009 89 89 ALLOCATE( zrhd(jpi,jpj,jpk) , zrhop(jpi,jpj,jpk) ) 90 90 ALLOCATE( ztsn(jpi,jpj,jpk,jpts) ) 91 zarea_ssh(:,:) = area(:,:) * ssh n(:,:)91 zarea_ssh(:,:) = area(:,:) * ssh(:,:,Nnn) 92 92 ENDIF 93 93 ! … … 100 100 CALL iom_put( 'voltot', zvol ) 101 101 CALL iom_put( 'sshtot', zvolssh / area_tot ) 102 CALL iom_put( 'sshdyn', ssh n(:,:) - (zvolssh / area_tot) )102 CALL iom_put( 'sshdyn', ssh(:,:,Nnn) - (zvolssh / area_tot) ) 103 103 ! 104 104 ENDIF … … 118 118 DO ji = 1, jpi 119 119 DO jj = 1, jpj 120 zbotpres(ji,jj) = zbotpres(ji,jj) + ssh n(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) 121 121 END DO 122 122 END DO 123 123 ELSE 124 zbotpres(:,:) = zbotpres(:,:) + ssh n(:,:) * zrhd(:,:,1)124 zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Nnn) * zrhd(:,:,1) 125 125 END IF 126 126 !!gm … … 147 147 DO ji = 1,jpi 148 148 DO jj = 1,jpj 149 zbotpres(ji,jj) = zbotpres(ji,jj) + ssh n(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) 150 150 END DO 151 151 END DO 152 152 ELSE 153 zbotpres(:,:) = zbotpres(:,:) + ssh n(:,:) * zrhd(:,:,1)153 zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Nnn) * zrhd(:,:,1) 154 154 END IF 155 155 END IF … … 162 162 ! ! ocean bottom pressure 163 163 zztmp = rho0 * grav * 1.e-4_wp ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 164 zbotpres(:,:) = zztmp * ( zbotpres(:,:) + ssh n(:,:) + thick0(:,:) )164 zbotpres(:,:) = zztmp * ( zbotpres(:,:) + ssh(:,:,Nnn) + thick0(:,:) ) 165 165 CALL iom_put( 'botpres', zbotpres ) 166 166 !
Note: See TracChangeset
for help on using the changeset viewer.