- Timestamp:
- 2018-03-30T18:08:29+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90
r8533 r9450 71 71 REAL(wp) :: zgrazp, zgrazm, zgrazsd 72 72 REAL(wp) :: zgrazmf, zgrazsf, zgrazpf 73 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d 73 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d, zfezoo 74 REAL(wp), POINTER, DIMENSION(:,:,:) :: zzligprod 74 75 CHARACTER (len=25) :: charout 75 76 !!--------------------------------------------------------------------- … … 77 78 IF( nn_timing == 1 ) CALL timing_start('p4z_micro') 78 79 ! 79 CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 80 CALL wrk_alloc( jpi, jpj, jpk, zgrazing, zfezoo ) 81 IF (ln_ligand) THEN 82 CALL wrk_alloc( jpi, jpj, jpk, zzligprod ) 83 zzligprod(:,:,:) = 0._wp 84 ENDIF 80 85 ! 81 86 DO jk = 1, jpkm1 … … 139 144 tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zgrarem - zgrarsig 140 145 ! 141 IF( ln_ligand ) tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + (zgrarem - zgrarsig) * ldocz 146 IF( ln_ligand ) THEN 147 tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + (zgrarem - zgrarsig) * ldocz 148 zzligprod(ji,jj,jk) = (zgrarem - zgrarsig) * ldocz 149 ENDIF 142 150 ! 143 151 tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * zgrarsig 144 152 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgrafer 153 zfezoo(ji,jj,jk) = zgrafer 145 154 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zgrapoc 146 155 prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + zgrapoc … … 184 193 CALL iom_put( "GRAZ1", zw3d ) 185 194 ENDIF 195 IF( iom_use( "FEZOO" ) ) THEN 196 zw3d(:,:,:) = zfezoo(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ! 197 CALL iom_put( "FEZOO", zw3d ) 198 ENDIF 199 IF( iom_use( "LPRODZ" ) .AND. ln_ligand ) THEN 200 zw3d(:,:,:) = zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 201 CALL iom_put( "LPRODZ" , zw3d ) 202 ENDIF 186 203 CALL wrk_dealloc( jpi, jpj, jpk, zw3d ) 187 204 ENDIF … … 194 211 ENDIF 195 212 ! 196 CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 213 CALL wrk_dealloc( jpi, jpj, jpk, zgrazing, zfezoo ) 214 IF (ln_ligand) CALL wrk_dealloc( jpi, jpj, jpk, zzligprod ) 197 215 ! 198 216 IF( nn_timing == 1 ) CALL timing_stop('p4z_micro')
Note: See TracChangeset
for help on using the changeset viewer.