Changeset 1678 for trunk/NEMO/TOP_SRC/PISCES
- Timestamp:
- 2009-10-26T12:38:26+01:00 (15 years ago)
- Location:
- trunk/NEMO/TOP_SRC/PISCES
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90
r1457 r1678 64 64 REAL(wp) :: zcoef1, zcoef2 65 65 #endif 66 #if defined key_trc_dia3d && defined key_kriest && defined key_iomput67 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d68 #endif69 66 CHARACTER (len=25) :: charout 70 67 … … 130 127 #endif 131 128 132 133 # if defined key_trc_dia3d && defined key_kriest134 #if ! defined key_iomput135 trc3d(:,:,:,jp_pcs0_3d + 10) = tra(:,:,:,jpcal) * xnegtr(:,:,:) * 1.e3 * rfact2r * tmask(:,:,:)136 #else137 zw3d(:,:,:) = tra(:,:,:,jpcal) * xnegtr(:,:,:) * 1.e3 * rfact2r * tmask(:,:,:)138 IF( jnt == nrdttrc ) CALL iom_put( "PBSi", zw3d )139 # endif140 #endif141 129 ! 142 130 IF(ln_ctl) THEN ! print mean trends (used for debugging) -
trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90
r1534 r1678 74 74 REAL(wp), DIMENSION(jpi,jpj) :: zkgco2, zkgo2, zh2co3 75 75 #if defined key_trc_diaadd && defined key_iomput 76 REAL(wp), DIMENSION(jpi,jpj) :: zcflx, zoflx, zkg, zd elc76 REAL(wp), DIMENSION(jpi,jpj) :: zcflx, zoflx, zkg, zdpco2, zdpo2 77 77 #endif 78 78 CHARACTER (len=25) :: charout … … 171 171 172 172 ! Compute O2 flux 173 zfld16 = atcox * chemc(ji,jj,2) * tmask(ji,jj,1) * zkgo2(ji,jj)173 zfld16 = atcox * chemc(ji,jj,2) * tmask(ji,jj,1) * zkgo2(ji,jj) 174 174 zflu16 = trn(ji,jj,1,jpoxy) * tmask(ji,jj,1) * zkgo2(ji,jj) 175 175 tra(ji,jj,1,jpoxy) = tra(ji,jj,1,jpoxy) + ( zfld16 - zflu16 ) / fse3t(ji,jj,1) … … 181 181 trc2d(ji,jj,jp_pcs0_2d + 1) = ( zfld16 - zflu16 ) * 1000. * tmask(ji,jj,1) 182 182 trc2d(ji,jj,jp_pcs0_2d + 2) = zkgco2(ji,jj) * tmask(ji,jj,1) 183 trc2d(ji,jj,jp_pcs0_2d + 3) = atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) * tmask(ji,jj,1) 183 trc2d(ji,jj,jp_pcs0_2d + 3) = ( atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) ) & 184 & * tmask(ji,jj,1) 184 185 # else 185 186 zcflx(ji,jj) = ( zfld - zflu ) * 1000. * tmask(ji,jj,1) 186 187 zoflx(ji,jj) = ( zfld16 - zflu16 ) * 1000. * tmask(ji,jj,1) 187 188 zkg (ji,jj) = zkgco2(ji,jj) * tmask(ji,jj,1) 188 zdelc(ji,jj) = atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) * tmask(ji,jj,1) 189 zdpco2(ji,jj) = ( atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) ) & 190 & * tmask(ji,jj,1) 191 zdpo2 (ji,jj) = ( atcox - trn(ji,jj,1,jpoxy) / ( chemc(ji,jj,2) + rtrn ) ) & 192 & * tmask(ji,jj,1) 189 193 # endif 190 194 #endif … … 224 228 225 229 # if defined key_trc_diaadd && defined key_iomput 226 CALL iom_put( "Cflx", zcflx ) 227 CALL iom_put( "Oflx", zoflx ) 228 CALL iom_put( "Kg" , zkg ) 229 CALL iom_put( "DelC", zdelc ) 230 CALL iom_put( "Cflx" , zcflx ) 231 CALL iom_put( "Oflx" , zoflx ) 232 CALL iom_put( "Kg" , zkg ) 233 CALL iom_put( "Dpco2", zdpco2 ) 234 CALL iom_put( "Dpo2" , zdpo2 ) 230 235 #endif 231 236 -
trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90
r1457 r1678 63 63 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zco3 64 64 #if defined key_trc_dia3d && defined key_iomput 65 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 65 REAL(wp) :: zrfact2 66 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zcaldiss, zw3d 66 67 #endif 67 68 CHARACTER (len=25) :: charout … … 72 73 zco3(:,:,:) = 0. 73 74 75 # if defined key_trc_dia3d && defined key_iomput 76 zcaldiss(:,:,:) = 0. 77 # endif 74 78 ! ------------------------------------------- 75 79 ! COMPUTE [CO3--] and [H+] CONCENTRATIONS … … 146 150 zremco3 = zdispot / rmoss 147 151 zco3(ji,jj,jk) = zco3(ji,jj,jk) + zremco3 * rfact 148 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2.*zremco3 149 tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zremco3 150 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zremco3 151 152 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2. * zremco3 153 tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zremco3 154 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zremco3 155 156 # if defined key_trc_dia3d && defined key_iomput 157 zcaldiss(ji,jj,jk) = zremco3 ! calcite dissolution 158 # endif 152 159 END DO 153 160 END DO … … 160 167 trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon * tmask(:,:,:) 161 168 # else 162 zw3d(:,:,:) = hi (:,:,:) * tmask(:,:,:) 169 zrfact2 = 1.e3 * rfact2r 170 zw3d(:,:,:) = hi (:,:,:) * tmask(:,:,:) 163 171 CALL iom_put( "PH", zw3d ) 164 zw3d(:,:,:) = zco3(:,:,:) * tmask(:,:,:)172 zw3d(:,:,:) = zco3(:,:,:) * tmask(:,:,:) 165 173 CALL iom_put( "CO3", zw3d ) 166 zw3d(:,:,:) = aksp(:,:,:) / calcon * tmask(:,:,:)174 zw3d(:,:,:) = aksp(:,:,:) / calcon * tmask(:,:,:) 167 175 CALL iom_put( "CO3sat", zw3d ) 176 zw3d(:,:,:) = zcaldiss(:,:,:) * zrfact2 * tmask(:,:,:) 177 CALL iom_put( "Dcal", zw3d ) 168 178 # endif 169 179 # endif -
trunk/NEMO/TOP_SRC/PISCES/p4zmeso.F90
r1152 r1678 21 21 USE p4zint 22 22 USE p4zsink 23 USE iom 23 24 24 25 IMPLICIT NONE … … 73 74 REAL(wp),DIMENSION(jpi,jpj,jpk) :: zgrazfff,zgrazffe 74 75 CHARACTER (len=25) :: charout 76 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput 77 REAL(wp) :: zrfact2 78 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 79 #endif 80 75 81 !!--------------------------------------------------------------------- 76 82 77 83 78 84 IF( ( kt * jnt ) == nittrc000 ) CALL p4z_meso_init ! Initialization (first time-step only) 85 86 zrespz2 (:,:,:) = 0. 87 ztortz2 (:,:,:) = 0. 88 zgrazd (:,:,:) = 0. 89 zgrazz (:,:,:) = 0. 90 zgrazpof(:,:,:) = 0. 91 zgrazn (:,:,:) = 0. 92 zgrazpoc(:,:,:) = 0. 93 zgraznf (:,:,:) = 0. 94 zgrazf (:,:,:) = 0. 95 zgrazfff(:,:,:) = 0. 96 zgrazffe(:,:,:) = 0. 79 97 80 98 zstep = rfact2 / rjjss ! Time step duration for biology … … 151 169 ! ---------------------------------- 152 170 # if ! defined key_kriest 171 # if ! defined key_off_degrad 153 172 zgrazffe(ji,jj,jk) = grazflux * zstep * wsbio4(ji,jj,jk) & 154 # if defined key_off_degrad 155 & * facvol(ji,jj,jk) & 173 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jpgoc) * trn(ji,jj,jk,jpmes) 174 # else 175 zgrazffe(ji,jj,jk) = grazflux * zstep * wsbio4(ji,jj,jk) * facvol(ji,jj,jk) & 176 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jpgoc) * trn(ji,jj,jk,jpmes) 156 177 # endif 157 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jpgoc) * trn(ji,jj,jk,jpmes)158 159 178 zgrazfff(ji,jj,jk) = zgrazffe(ji,jj,jk) & 160 & * trn(ji,jj,jk,jpbfe) / (trn(ji,jj,jk,jpgoc) + rtrn)179 & * trn(ji,jj,jk,jpbfe) / (trn(ji,jj,jk,jpgoc) + rtrn) 161 180 # else 162 ! KRIEST3181 !!--------------------------- KRIEST3 ------------------------------------------- 163 182 !! zgrazffe(ji,jj,jk) = 0.5 * 1.3e-2 / 5.5e-7 * 0.3 * zstep * wsbio3(ji,jj,jk) & 164 183 !! & * tgfunc(ji,jj,jk) * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpmes) & … … 167 186 # endif 168 187 !! & / (trn(ji,jj,jk,jppoc) * 1.e7 + 0.1) 169 188 !!--------------------------- KRIEST3 ------------------------------------------- 189 190 # if ! defined key_off_degrad 170 191 zgrazffe(ji,jj,jk) = grazflux * zstep * wsbio3(ji,jj,jk) & 171 # if defined key_off_degrad 172 & * facvol(ji,jj,jk) & 192 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpmes) 193 # else 194 zgrazffe(ji,jj,jk) = grazflux * zstep * wsbio3(ji,jj,jk) * facvol(ji,jj,jk) & 195 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpmes) 173 196 # endif 174 & * tgfunc2(ji,jj,jk) * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpmes)175 197 176 198 zgrazfff(ji,jj,jk) = zgrazffe(ji,jj,jk) & 177 & * trn(ji,jj,jk,jpsfe) / (trn(ji,jj,jk,jppoc) + rtrn)199 & * trn(ji,jj,jk,jpsfe) / (trn(ji,jj,jk,jppoc) + rtrn) 178 200 # endif 179 201 END DO … … 253 275 tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) - zgrazf(ji,jj,jk) 254 276 255 zprcaca = xfracal(ji,jj,jk) * part * unass2 * zgrazn(ji,jj,jk) 256 277 zprcaca = xfracal(ji,jj,jk) * unass2 * zgrazn(ji,jj,jk) 278 #if defined key_trc_dia3d 279 prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo) 280 #endif 281 zprcaca = part * zprcaca 257 282 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprcaca 258 283 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca … … 284 309 END DO 285 310 ! 311 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput 312 zrfact2 = 1.e3 * rfact2r 313 zw3d(:,:,:) = ( zgrazd(:,:,:) + zgrazz(:,:,:) + zgrazn(:,:,:) & 314 & + zgrazpoc(:,:,:) + zgrazffe(:,:,:) ) * zrfact2 * tmask(:,:,:) 315 IF( jnt == nrdttrc ) CALL iom_put( "Graz2" , zw3d ) 316 317 zw3d(:,:,:) = prodcal(:,:,:) * zrfact2 * tmask(:,:,:) 318 IF( jnt == nrdttrc ) CALL iom_put( "Pcal" , zw3d ) 319 #endif 320 286 321 IF(ln_ctl) THEN ! print mean trends (used for debugging) 287 322 WRITE(charout, FMT="('meso')") -
trunk/NEMO/TOP_SRC/PISCES/p4zmicro.F90
r1152 r1678 21 21 USE p4zint 22 22 USE p4zsink 23 USE iom 23 24 24 25 IMPLICIT NONE … … 69 70 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazmf, zgrazsf, zgrazpf 70 71 CHARACTER (len=25) :: charout 72 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput 73 REAL(wp) :: zrfact2 74 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 75 #endif 71 76 72 77 !!--------------------------------------------------------------------- … … 74 79 IF( ( kt * jnt ) == nittrc000 ) CALL p4z_micro_init ! Initialization (first time-step only) 75 80 76 77 zstep = rfact2 / rjjss ! Time step duration for biology 81 zrespz (:,:,:) = 0. 82 ztortz (:,:,:) = 0. 83 zgrazp (:,:,:) = 0. 84 zgrazm (:,:,:) = 0. 85 zgrazsd(:,:,:) = 0. 86 zgrazmf(:,:,:) = 0. 87 zgrazsf(:,:,:) = 0. 88 zgrazpf(:,:,:) = 0. 89 90 91 zstep = rfact2 / rjjss ! Time step duration for biology 78 92 79 93 DO jk = 1, jpkm1 … … 202 216 & + unass * ( zgrazpf(ji,jj,jk) + zgrazsf (ji,jj,jk)) & 203 217 & - (1.-unass) * zgrazmf(ji,jj,jk) 204 zprcaca = xfracal(ji,jj,jk) * part * unass * zgrazp(ji,jj,jk) 218 zprcaca = xfracal(ji,jj,jk) * unass * zgrazp(ji,jj,jk) 219 #if defined key_trc_dia3d 220 prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo) 221 #endif 222 zprcaca = part * zprcaca 205 223 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprcaca 206 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca224 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca 207 225 tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca 208 226 #if defined key_kriest … … 213 231 END DO 214 232 ! 233 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput 234 zrfact2 = 1.e3 * rfact2r 235 zw3d(:,:,:) = ( zgrazp(:,:,:) + zgrazm(:,:,:) + zgrazsd(:,:,:) ) * zrfact2 * tmask(:,:,:) 236 IF( jnt == nrdttrc ) CALL iom_put( "Graz" , zw3d ) 237 #endif 238 215 239 IF(ln_ctl) THEN ! print mean trends (used for debugging) 216 240 WRITE(charout, FMT="('micro')") -
trunk/NEMO/TOP_SRC/PISCES/p4zmort.F90
r1152 r1678 89 89 90 90 91 #if defined key_trc_dia3d 92 prodcal(:,:,:) = 0. !: Initialisation of calcite production variable 93 #endif 94 91 95 DO jk = 1, jpkm1 92 96 DO jj = 1, jpj … … 126 130 tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zmortp * zfactfe 127 131 zprcaca = xfracal(ji,jj,jk) * zmortp 132 #if defined key_trc_dia3d 133 prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo) 134 #endif 128 135 zfracal = 0.5 * xfracal(ji,jj,jk) 129 136 tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprcaca -
trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90
r1457 r1678 240 240 zw2d(:,: ) = heup(:,: ) * tmask(:,:,1) 241 241 zw3d(:,:,:) = etot(:,:,:) * tmask(:,:,:) 242 IF( jnt == nrdttrc ) CALL iom_put( " heup", zw2d )242 IF( jnt == nrdttrc ) CALL iom_put( "Heup", zw2d ) 243 243 IF( jnt == nrdttrc ) CALL iom_put( "PAR" , zw3d ) 244 244 # endif -
trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90
r1511 r1678 89 89 REAL(wp), DIMENSION(jpi,jpj) :: zsidep 90 90 REAL(wp), DIMENSION(jpi,jpj,jpk) :: znitrpot, zirondep 91 #if defined key_trc_diaadd && defined key_iomput 92 REAL(wp), DIMENSION(jpi,jpj) :: zw2d 91 #if defined key_diaadd || defined key_trc_dia3d 92 REAL(wp) :: zrfact2 93 # if defined key_iomput 94 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 95 REAL(wp), DIMENSION(jpi,jpj) :: zw2d 96 # endif 93 97 #endif 94 98 CHARACTER (len=25) :: charout … … 267 271 IF( lk_mpp ) CALL mpp_sum( zdenitot ) ! sum over the global domain 268 272 269 ! Potential nitrogen fi cation dependant on temperature and iron273 ! Potential nitrogen fixation dependant on temperature and iron 270 274 ! ------------------------------------------------------------- 271 275 … … 318 322 END DO 319 323 320 #if defined key_trc_diaadd 324 #if defined key_trc_diaadd || defined key_trc_dia3d 325 zrfact2 = 1.e+3 * rfact2r 321 326 # if ! defined key_iomput 322 trc2d(:,:,jp_pcs0_2d + 11) = zirondep(:,:,1) * 1.e+3 * rfact2r* fse3t(:,:,1) * tmask(:,:,1)323 trc2d(:,:,jp_pcs0_2d + 12) = znitrpot(:,:,1) * 1.e-7 * 1.e+3 / rfact2* fse3t(:,:,1) * tmask(:,:,1)327 trc2d(:,:,jp_pcs0_2d + 11) = zirondep(:,:,1) * zrfact2 * fse3t(:,:,1) * tmask(:,:,1) 328 trc2d(:,:,jp_pcs0_2d + 12) = znitrpot(:,:,1) * 1.e-7 * zrfact2 * fse3t(:,:,1) * tmask(:,:,1) 324 329 # else 325 330 ! write diagnostics 326 zw2d(:,:) = zirondep(:,:,1) * 1.e+3 * rfact2r * fse3t(:,:,1) * tmask(:,:,1) 327 IF( jnt == nrdttrc ) CALL iom_put( "Fedep", zw2d ) 328 zw2d(:,:) = znitrpot(:,:,1) * 1.e-7 * 1.e+3 / rfact2 * fse3t(:,:,1) * tmask(:,:,1) 329 IF( jnt == nrdttrc ) CALL iom_put( "Nfix", zw2d ) 331 zw2d(:,:) = ( zirondep(:,:,1) + ironsed(:,:,1) * rfact2 ) & 332 & * zrfact2 * fse3t(:,:,1) * tmask(:,:,1) 333 IF( jnt == nrdttrc ) CALL iom_put( "Irondep", zw2d ) 334 zw3d(:,:,:) = znitrpot(:,:,:) * 1.e-7 * zrfact2 * fse3t(:,:,:) * tmask(:,:,:) 335 IF( jnt == nrdttrc ) CALL iom_put( "Nfix", zw3d ) 330 336 # endif 331 337 -
trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90
r1457 r1678 101 101 INTEGER :: iksed1 102 102 #if defined key_iomput 103 REAL(wp), DIMENSION(jpi,jpj) :: zw2d104 103 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 105 104 #endif … … 302 301 trc3d(:,:,:,jp_pcs0_3d + 17) = wsbio4 (:,:,:) * tmask(:,:,:) 303 302 #else 304 zw 2d(:,: ) = sinking (:,:,iksed1) * zrfact2 * tmask(:,:,1)305 IF( jnt == nrdttrc ) CALL iom_put( "PMO" , zw2d )306 zw 2d(:,: ) = sinking2(:,:,iksed1) * zrfact2 * tmask(:,:,1)307 IF( jnt == nrdttrc ) CALL iom_put( "PMO2", zw 2d )308 zw 2d(:,: ) = sinkfer (:,:,iksed1) * zrfact2 * tmask(:,:,1)309 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe1", zw 2d )310 zw 2d(:,: ) = sinksil (:,:,iksed1) * zrfact2 * tmask(:,:,1)311 IF( jnt == nrdttrc ) CALL iom_put( "ExpSi", zw 2d )312 zw 2d(:,: ) = sinkcal (:,:,iksed1) * zrfact2 * tmask(:,:,1)313 IF( jnt == nrdttrc ) CALL iom_put( "ExpCaCO3", zw 2d )314 zw3d(:,:,:) = sinking (:,:,:) 303 zw3d(:,:,:) = sinking (:,:,:) * zrfact2 * tmask(:,:,:) 304 IF( jnt == nrdttrc ) CALL iom_put( "PMO" , zw3d ) 305 zw3d(:,:,:) = sinking2(:,:,:) * zrfact2 * tmask(:,:,:) 306 IF( jnt == nrdttrc ) CALL iom_put( "PMO2", zw3d ) 307 zw3d(:,:,:) = sinkfer (:,:,:) * zrfact2 * tmask(:,:,:) 308 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe1", zw3d ) 309 zw3d(:,:,:) = sinksil (:,:,:) * zrfact2 * tmask(:,:,:) 310 IF( jnt == nrdttrc ) CALL iom_put( "ExpSi", zw3d ) 311 zw3d(:,:,:) = sinkcal (:,:,:) * zrfact2 * tmask(:,:,:) 312 IF( jnt == nrdttrc ) CALL iom_put( "ExpCaCO3", zw3d ) 313 zw3d(:,:,:) = sinking (:,:,:) * zrfact2 * tmask(:,:,:) 315 314 IF( jnt == nrdttrc ) CALL iom_put( "POCFlx", zw3d ) 316 zw3d(:,:,:) = sinking2(:,:,:) 315 zw3d(:,:,:) = sinking2(:,:,:) * zrfact2 * tmask(:,:,:) 317 316 IF( jnt == nrdttrc ) CALL iom_put( "GOCFlx", zw3d ) 318 zw3d(:,:,:) = sinksil (:,:,:) 317 zw3d(:,:,:) = sinksil (:,:,:) * zrfact2 * tmask(:,:,:) 319 318 IF( jnt == nrdttrc ) CALL iom_put( "SiFlx", zw3d ) 320 zw3d(:,:,:) = sinkcal (:,:,:) 319 zw3d(:,:,:) = sinkcal (:,:,:) * zrfact2 * tmask(:,:,:) 321 320 IF( jnt == nrdttrc ) CALL iom_put( "CaCO3Flx", zw3d ) 322 zw3d(:,:,:) = znum3d (:,:,:) 321 zw3d(:,:,:) = znum3d (:,:,:) * tmask(:,:,:) 323 322 IF( jnt == nrdttrc ) CALL iom_put( "xnum", zw3d ) 324 zw3d(:,:,:) = wsbio3 (:,:,:) 323 zw3d(:,:,:) = wsbio3 (:,:,:) * tmask(:,:,:) 325 324 IF( jnt == nrdttrc ) CALL iom_put( "W1", zw3d ) 326 zw3d(:,:,:) = wsbio4 (:,:,:) 325 zw3d(:,:,:) = wsbio4 (:,:,:) * tmask(:,:,:) 327 326 IF( jnt == nrdttrc ) CALL iom_put( "W2", zw3d ) 328 327 # endif … … 493 492 #endif 494 493 #if defined key_iomput 495 REAL(wp), DIMENSION(jpi,jpj ) :: zw2d494 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 496 495 #endif 497 496 CHARACTER (len=25) :: charout … … 623 622 trc2d(:,:,jp_pcs0_2d + 9) = sinkcal (:,:,iksed1) * zrfact2 * tmask(:,:,1) 624 623 # else 625 ! write diagnostics 626 zw2d(:,:) = sinking (:,:,iksed1) * zrfact2 * tmask(:,:,1) 627 IF( jnt == nrdttrc ) CALL iom_put( "PMO", zw2d ) 628 zw2d(:,:) = sinking2(:,:,iksed1) * zrfact2 * tmask(:,:,1) 629 IF( jnt == nrdttrc ) CALL iom_put( "PMO2", zw2d ) 630 zw2d(:,:) = sinkfer (:,:,iksed1) * zrfact2 * tmask(:,:,1) 631 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe1", zw2d ) 632 zw2d(:,:) = sinkfer2(:,:,iksed1) * zrfact2 * tmask(:,:,1) 633 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe2", zw2d ) 634 zw2d(:,:) = sinksil (:,:,iksed1) * zrfact2 * tmask(:,:,1) 635 IF( jnt == nrdttrc ) CALL iom_put( "ExpSi", zw2d ) 636 zw2d(:,:) = sinkcal (:,:,iksed1) * zrfact2 * tmask(:,:,1) 637 IF( jnt == nrdttrc ) CALL iom_put( "ExpCaCO3", zw2d ) 624 zw3d(:,:,:) = sinking (:,:,:) * zrfact2 * tmask(:,:,:) 625 IF( jnt == nrdttrc ) CALL iom_put( "ExpPOC" , zw3d ) 626 zw3d(:,:,:) = sinking2(:,:,:) * zrfact2 * tmask(:,:,:) 627 IF( jnt == nrdttrc ) CALL iom_put( "ExpGOC", zw3d ) 628 zw3d(:,:,:) = sinkfer (:,:,:) * zrfact2 * tmask(:,:,:) 629 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe1", zw3d ) 630 zw3d(:,:,:) = sinkfer2(:,:,:) * zrfact2 * tmask(:,:,:) 631 IF( jnt == nrdttrc ) CALL iom_put( "ExpFe2", zw3d ) 632 zw3d(:,:,:) = sinksil (:,:,:) * zrfact2 * tmask(:,:,:) 633 IF( jnt == nrdttrc ) CALL iom_put( "ExpSi", zw3d ) 634 zw3d(:,:,:) = sinkcal (:,:,:) * zrfact2 * tmask(:,:,:) 635 IF( jnt == nrdttrc ) CALL iom_put( "Expcal", zw3d ) 638 636 # endif 639 637 #endif -
trunk/NEMO/TOP_SRC/PISCES/sms_pisces.F90
r1445 r1678 60 60 REAL(wp), DIMENSION(jpi,jpj,jpk) :: xlimbac !: ?? 61 61 REAL(wp), DIMENSION(jpi,jpj,jpk) :: xdiss !: ?? 62 #if defined key_trc_dia3d 63 REAL(wp), DIMENSION(jpi,jpj,jpk) :: prodcal !: Calcite production 64 #endif 62 65 63 66 !!* Variable for chemistry of the CO2 cycle
Note: See TracChangeset
for help on using the changeset viewer.