Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90
- Timestamp:
- 2010-12-27T18:33:53+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90
- Property svn:keywords set to Id
r1836 r2528 21 21 USE iom 22 22 USE trcdta 23 USE lib_mpp 24 USE lib_fortran 23 25 24 26 IMPLICIT NONE … … 118 120 IF(lwp) WRITE(numout,*) 119 121 120 IF( cp_cfg == "orca" .AND. .NOT. lk_ trc_c1d ) THEN ! ORCA condiguration (not 1D) !122 IF( cp_cfg == "orca" .AND. .NOT. lk_c1d ) THEN ! ORCA condiguration (not 1D) ! 121 123 ! ! --------------------------- ! 122 124 ! set total alkalinity, phosphate, nitrate & silicate 123 125 124 zalksum = 0.e0125 zpo4sum = 0.e0126 zno3sum = 0.e0127 zsilsum = 0.e0128 DO jk = 1, jpk129 DO jj = 1, jpj130 DO ji = 1, jpi131 zvol = cvol(ji,jj,jk)132 # if defined key_off_degrad133 zvol = zvol * facvol(ji,jj,jk)134 # endif135 zalksum = zalksum + trn(ji,jj,jk,jptal) * zvol136 zpo4sum = zpo4sum + trn(ji,jj,jk,jppo4) * zvol137 zno3sum = zno3sum + trn(ji,jj,jk,jpno3) * zvol138 zsilsum = zsilsum + trn(ji,jj,jk,jpsil) * zvol139 END DO140 END DO141 END DO142 IF( lk_mpp ) CALL mpp_sum( zalksum ) ! sum over the global domain143 IF( lk_mpp ) CALL mpp_sum( zpo4sum ) ! sum over the global domain144 IF( lk_mpp ) CALL mpp_sum( zno3sum ) ! sum over the global domain145 IF( lk_mpp ) CALL mpp_sum( zsilsum ) ! sum over the global domain146 126 zarea = 1. / areatot * 1.e6 147 zalksum = zalksum * zarea 148 zpo4sum = zpo4sum * zarea / 122. 149 zno3sum = zno3sum * zarea / 7.6 150 zsilsum = zsilsum * zarea 127 # if defined key_degrad 128 zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:) * facvol(:,:,:) ) * zarea 129 zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:) * facvol(:,:,:) ) * zarea / 122. 130 zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:) * facvol(:,:,:) ) * zarea / 7.6 131 zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:) * facvol(:,:,:) ) * zarea 132 # else 133 zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:) ) * zarea 134 zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:) ) * zarea / 122. 135 zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:) ) * zarea / 7.6 136 zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:) ) * zarea 137 # endif 151 138 152 139 IF(lwp) WRITE(numout,*) ' TALK mean : ', zalksum … … 263 250 #if defined key_dtatrc 264 251 ! Restore close seas values to initial data 265 CALL trc_dta( nit trc000 )252 CALL trc_dta( nit000 ) 266 253 DO jn = 1, jptra 267 254 IF( lutini(jn) ) THEN
Note: See TracChangeset
for help on using the changeset viewer.