Changeset 13200 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90

Timestamp:

2020-07-02T08:58:10+02:00 (4 years ago)

Author:

aumont

Message:

minor bug fixes for calcite production and diagnostics

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsed.F90

@@ -159 +159 @@
          DO jk = 2, jpkm1
             zirondep(:,:,jk) = dust(:,:) * mfrac * zwdust * rfact2 * EXP( -gdept_n(:,:,jk) / (250. * wdust) )
-            zpdep   (:,:,jk) = zirondep(:,:,jk) * 0.023
+!            zpdep   (:,:,jk) = zirondep(:,:,jk) * 0.023
+            zpdep   (:,:,jk) = zirondep(:,:,jk) * 0.38 / po4r
          END DO
          ! Solubilization of particles in the water column (Si, P, Fe)
@@ -173 +174 @@
                 &  CALL iom_put( "Irondep", zirondep(:,:,1) * 1.e+3 * rfact2r * e3t_n(:,:,1) * tmask(:,:,1) ) ! surface downward dust depo of iron
                 IF( iom_use( "pdust" ) )   &
-                &  CALL iom_put( "pdust"  , dust(:,:) / ( wdust * rday )  * tmask(:,:,1) ) ! dust concentration at surface
+                &  CALL iom_put( "pdust"  , dust(:,:) / ( wdust / rday )  * tmask(:,:,1) ) ! dust concentration at surface
             ENDIF
          ENDIF