Changeset 5845 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90
- Timestamp:
- 2015-10-31T08:40:45+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90
r5836 r5845 34 34 REAL(wp), PUBLIC :: xhr ! coeff for martin''s remineralisation profile 35 35 36 !! * Substitutions37 # include "domzgr_substitute.h90"38 36 !!---------------------------------------------------------------------- 39 37 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 102 100 DO jj = 1, jpj 103 101 DO ji = 1, jpi 104 ztra(ji,jj,jk) = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / fse3t(ji,jj,jk)102 ztra(ji,jj,jk) = - ( zwork(ji,jj,jk) - zwork(ji,jj,jk+1) ) / e3t_n(ji,jj,jk) 105 103 tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + ztra(ji,jj,jk) 106 104 END DO … … 111 109 IF( iom_use( "TDETSED" ) ) THEN 112 110 CALL wrk_alloc( jpi, jpj, zw2d ) 113 zw2d(:,:) = ztra(:,:,1) * fse3t(:,:,1) * 86400.111 zw2d(:,:) = ztra(:,:,1) * e3t_n(:,:,1) * 86400. 114 112 DO jk = 2, jpkm1 115 zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400.113 zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * e3t_n(:,:,jk) * 86400. 116 114 END DO 117 115 CALL iom_put( "TDETSED", zw2d ) … … 121 119 IF( ln_diatrc ) THEN 122 120 CALL wrk_alloc( jpi, jpj, zw2d ) 123 zw2d(:,:) = ztra(:,:,1) * fse3t(:,:,1) * 86400.121 zw2d(:,:) = ztra(:,:,1) * e3t_n(:,:,1) * 86400. 124 122 DO jk = 2, jpkm1 125 zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400.123 zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * e3t_n(:,:,jk) * 86400. 126 124 END DO 127 125 trc2d(:,:,jp_pcs0_2d + 7) = zw2d(:,:)
Note: See TracChangeset
for help on using the changeset viewer.