Changeset 7617 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
- Timestamp:
- 2017-01-27T18:17:01+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
r7180 r7617 69 69 INTEGER :: ji, jj, jk, jic, jn 70 70 REAL(wp) :: zdep, zlam1a, zlam1b, zlamfac 71 REAL(wp) :: zkeq, zfeequi, zfesatur, zfecoll, fe3sol, fe3sol171 REAL(wp) :: zkeq, zfeequi, zfesatur, zfecoll, zfecollc, fe3sol, fe3sol1 72 72 REAL(wp) :: zdenom1, zscave, zaggdfea, zaggdfeb, zcoag 73 73 REAL(wp) :: ztrc, zdust … … 78 78 REAL(wp), POINTER, DIMENSION(:,:,:) :: zFeL1, zFeL2, zTL2, zFe2, zFeP 79 79 REAL(wp), POINTER, DIMENSION(:,: ) :: zstrn, zstrn2 80 REAL(wp), POINTER, DIMENSION(:,:,:) :: zcoll3d, zscav3d 80 81 REAL(wp) :: zzFeL1, zzFeL2, zzFe2, zzFeP, zzFe3, zzstrn2 81 82 REAL(wp) :: zrum, zcodel, zargu, zlight … … 84 85 REAL(wp) :: zxs, zfunc, zp, zq, zd, zr, zphi, zfff, zp3, zq2 85 86 REAL(wp) :: ztfe, zoxy, zhplus 86 REAL(wp) :: zstep 87 REAL(wp) :: zstep, zrfact2 87 88 #if defined key_ligand 88 89 REAL(wp) :: zaggliga, zaggligb 89 90 REAL(wp) :: dissol, zligco 91 REAL(wp), POINTER, DIMENSION(:,:,:) :: zlcoll3d 90 92 #endif 91 93 CHARACTER (len=25) :: charout … … 96 98 ! Allocate temporary workspace 97 99 CALL wrk_alloc( jpi, jpj, jpk, zFe3, zFeL1, zTL1, ztotlig, precip ) 100 CALL wrk_alloc( jpi, jpj, jpk, zcoll3d, zscav3d ) 101 #if defined key_ligand 102 CALL wrk_alloc( jpi, jpj, jpk, zlcoll3d ) 103 #endif 98 104 zFe3 (:,:,:) = 0. 99 105 zFeL1(:,:,:) = 0. … … 327 333 ! It requires certainly some more work as it is very poorly constrained. 328 334 ! ---------------------------------------------------------------- 335 zfecollc = trb(ji,jj,jk,jpdoc) / ( 40.E-6 + trb(ji,jj,jk,jpdoc) ) * zfecoll 329 336 zlam1a = ( 0.369 * 0.3 * trb(ji,jj,jk,jpdoc) + 102.4 * trb(ji,jj,jk,jppoc) ) * xdiss(ji,jj,jk) & 330 337 & + ( 114. * 0.3 * trb(ji,jj,jk,jpdoc) ) 331 zaggdfea = zlam1a * zstep * zfecoll 338 zaggdfea = zlam1a * zstep * zfecollc 332 339 ! 333 340 #if defined key_kriest … … 335 342 #else 336 343 zlam1b = 3.53E3 * trb(ji,jj,jk,jpgoc) * xdiss(ji,jj,jk) 337 zaggdfeb = zlam1b * zstep * zfecoll 344 zaggdfeb = zlam1b * zstep * zfecollc 338 345 #endif 339 346 ! precipitation of Fe3+, creation of nanoparticles … … 346 353 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zscave * zdenom2 + zaggdfeb 347 354 #endif 355 zscav3d(ji,jj,jk) = zscave 356 zcoll3d(ji,jj,jk) = zaggdfea + zaggdfeb 348 357 #if defined key_ligand 349 358 zaggliga = zlam1a * zstep * zligco … … 355 364 tra(ji,jj,jk,jpfep) = tra(ji,jj,jk,jpfep) + precip(ji,jj,jk) 356 365 tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) - zaggliga - zaggligb 366 zlcoll3d(ji,jj,jk) = zaggliga + zaggligb 357 367 #endif 358 368 END DO … … 379 389 IF( ln_diatrc .AND. lk_iomput ) THEN 380 390 IF( jnt == nrdttrc ) THEN 391 zrfact2 = 1.e3 * rfact2r ! conversion from mol/L/timestep into mol/m3/s 381 392 CALL iom_put("Fe3" , zFe3 (:,:,:) * tmask(:,:,:) ) ! Fe3+ 382 393 CALL iom_put("FeL1" , zFeL1 (:,:,:) * tmask(:,:,:) ) ! FeL1 383 394 CALL iom_put("TL1" , zTL1 (:,:,:) * tmask(:,:,:) ) ! TL1 384 395 CALL iom_put("Totlig" , ztotlig(:,:,:) * tmask(:,:,:) ) ! TL 385 CALL iom_put("Biron" , biron (:,:,:) * 1e9 * tmask(:,:,:) ) ! biron 396 CALL iom_put("Biron" , biron (:,:,:) * 1e9 * tmask(:,:,:) ) ! biron 397 CALL iom_put("FESCAV" , zscav3d(:,:,:) * 1e9 * tmask(:,:,:) * zrfact2 ) 398 CALL iom_put("FECOLL" , zcoll3d(:,:,:) * 1e9 * tmask(:,:,:) * zrfact2 ) 399 #if defined key_ligand 400 CALL iom_put("LGWCOLL", zlcoll3d(:,:,:) * 1e9 * tmask(:,:,:) * zrfact2 ) 401 #endif 386 402 IF( ln_fechem ) THEN 387 403 CALL iom_put("Fe2" , zFe2 (:,:,:) * tmask(:,:,:) ) ! Fe2+ … … 400 416 ! 401 417 CALL wrk_dealloc( jpi, jpj, jpk, zFe3, zFeL1, zTL1, ztotlig, precip ) 418 CALL wrk_dealloc( jpi, jpj, jpk, zscav3d, zcoll3d ) 419 #if defined key_ligand 420 CALL wrk_dealloc( jpi, jpj, jpk, zlcoll3d ) 421 #endif 422 402 423 IF( ln_fechem ) THEN 403 424 CALL wrk_dealloc( jpi, jpj, zstrn, zstrn2 )
Note: See TracChangeset
for help on using the changeset viewer.