Changeset 15075 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zbc.F90

Timestamp:

2021-07-02T17:15:57+02:00 (3 years ago)

Author:

aumont

Message:

major update of the sediment module

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zbc.F90 (modified) (4 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zbc.F90

@@ -24 +24 @@
    LOGICAL , PUBLIC ::   ln_ironsed   !: boolean for Fe input from sediments
    LOGICAL , PUBLIC ::   ln_hydrofe   !: boolean for Fe input from hydrothermal vents
+   LOGICAL , PUBLIC ::   ln_dust_inp  !: boolean for Fe input from hydrothermal vents
    REAL(wp), PUBLIC ::   sedfeinput   !: Coastal release of Iron
    REAL(wp), PUBLIC ::   icefeinput   !: Iron concentration in sea ice
@@ -97 +98 @@
 
          zwdust = 0.03 / ( wdust / rday ) / ( 250. * rday )
-         DO_3D( 1, 1, 1, 1, 2, jpkm1 )
-            zdustdep = dust(ji,jj) * zwdust * rfact * EXP( -gdept(ji,jj,jk,Kmm) /( 250. * wdust ) )
-            !
-            tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zdustdep * mfrac / mMass_Fe 
-            tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) + zdustdep * 1.e-3 / mMass_P
-            tr(ji,jj,jk,jpsil,Krhs) = tr(ji,jj,jk,jpsil,Krhs) + zdustdep * 0.269 / mMass_Si
-         END_3D
+
+         ! Atmospheric input of Iron dissolves in the water column
+         IF ( ln_trc_sbc(jpfer) ) THEN
+            DO_3D( 1, 1, 1, 1, 2, jpkm1 )
+               zdustdep = dust(ji,jj) * zwdust * rfact * EXP( -gdept(ji,jj,jk,Kmm) /( 250. * wdust ) )
+               !
+               tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zdustdep * mfrac / mMass_Fe 
+            END_3D
+
+            IF( lk_iomput ) THEN
+                ! surface downward dust depo of iron
+                jl = n_trc_indsbc(jpfer)
+                CALL iom_put( "Irondep", ( rf_trsfac(jl) * sf_trcsbc(jl)%fnow(:,:,1) / rn_sbc_time ) * 1.e+3 * tmask(:,:,1) )
+
+            ENDIF
+
+         ENDIF
+
+         ! Atmospheric input of PO4 dissolves in the water column
+         IF ( ln_trc_sbc(jppo4) ) THEN
+            DO_3D( 1, 1, 1, 1, 2, jpkm1 )
+               zdustdep = dust(ji,jj) * zwdust * rfact * EXP( -gdept(ji,jj,jk,Kmm) /( 250. * wdust ) )
+               !
+               tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) + zdustdep * 1.e-3 / mMass_P
+            END_3D
+         ENDIF
+
+         ! Atmospheric input of Si dissolves in the water column
+         IF ( ln_trc_sbc(jpsil) ) THEN
+            DO_3D( 1, 1, 1, 1, 2, jpkm1 )
+               zdustdep = dust(ji,jj) * zwdust * rfact * EXP( -gdept(ji,jj,jk,Kmm) /( 250. * wdust ) )
+               !
+               tr(ji,jj,jk,jpsil,Krhs) = tr(ji,jj,jk,jpsil,Krhs) + zdustdep * 0.269 / mMass_Si
+            END_3D
+         ENDIF
+
          !
          IF( lk_iomput ) THEN
-             ! surface downward dust depo of iron
-             jl = n_trc_indsbc(jpfer)
-             CALL iom_put( "Irondep", ( rf_trsfac(jl) * sf_trcsbc(jl)%fnow(:,:,1) / rn_sbc_time ) * 1.e+3 * tmask(:,:,1) )
-
              ! dust concentration at surface
              CALL iom_put( "pdust"  , dust(:,:) / ( wdust / rday ) * tmask(:,:,1) ) ! dust concentration at surface
@@ -225 +251 @@
       !!
       NAMELIST/nampisbc/cn_dir, sn_dust, sn_ironsed, sn_hydrofe, &
-        &                ln_ironsed, ln_ironice, ln_hydrofe,    &
-        &                sedfeinput, distcoast, icefeinput, wdust, mfrac,  &
+        &                ln_ironsed, ln_ironice, ln_hydrofe, ln_dust_inp,   &
+        &                sedfeinput, distcoast, icefeinput, wdust, mfrac,   &
         &                hratio, lgw_rath
       !!----------------------------------------------------------------------
@@ -267 +293 @@
       END IF
 
-      ll_bc    = ( ln_trcbc .AND. lltrcbc )  .OR. ln_hydrofe .OR. ln_ironsed .OR. ln_ironice
-      ll_dust  =  ln_trc_sbc(jpfer)   
+      ll_bc    = ( ln_trcbc .AND. lltrcbc )  .OR. ln_hydrofe .OR. ln_ironsed .OR. ln_ironice .OR. ln_dust_inp
+      ll_dust  =  ln_dust_inp
       ll_ndepo =  ln_trc_sbc(jpno3) .OR. ln_trc_sbc(jpnh4)   
       ll_river =  ln_trc_cbc(jpno3)