Changeset 15450 for NEMO/trunk/src/TOP/PISCES/SED/sedwri.F90

Timestamp:

2021-10-27T16:32:08+02:00 (3 years ago)

Author:

cetlod

Message:

Some updates to make the PISCES/SED module usable. Totally orthogonal with no effect on other parts of the code

File:

• NEMO/trunk/src/TOP/PISCES/SED/sedwri.F90 (modified) (6 diffs)

NEMO/trunk/src/TOP/PISCES/SED/sedwri.F90

@@ -4 +4 @@
    !!         Sediment diagnostics :  write sediment output files
    !!======================================================================
-   USE par_sed
    USE sed
    USE sedarr
    USE lib_mpp         ! distribued memory computing library
    USE iom
+   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 
    IMPLICIT NONE
@@ -38 +38 @@
       CHARACTER(len = 20)  ::  cltra 
       REAL(wp)  :: zrate
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zdta, zflx
+      REAL(wp), DIMENSION(jpoce, jpksed)     :: zdta
+      REAL(wp), DIMENSION(jpoce, jptrased+1) :: zflx
+      REAL(wp), DIMENSION(jpi, jpj, jpksed, jptrased)   :: trcsedi
+      REAL(wp), DIMENSION(jpi, jpj, jpksed, jpdia3dsed) :: flxsedi3d
+      REAL(wp), DIMENSION(jpi, jpj, jpdia2dsed) :: flxsedi2d
 
       !!-------------------------------------------------------------------
@@ -54 +58 @@
       IF (lwp) WRITE(numsed,*) ' '
       
-      ALLOCATE( zdta(jpoce,jpksed) )    ;   ALLOCATE( zflx(jpoce,jpwatp1) )
-
       ! Initialize variables
       ! --------------------
@@ -88 +90 @@
       CALL unpack_arr( jpoce, flxsedi3d(1:jpi,1:jpj,1:jpksed,2)  , iarroce(1:jpoce), &
          &                   co3por(1:jpoce,1:jpksed)  )
+
+      CALL unpack_arr( jpoce, flxsedi3d(1:jpi,1:jpj,1:jpksed,3)  , iarroce(1:jpoce), &
+         &                   saturco3(1:jpoce,1:jpksed)  )
+
       
 !      flxsedi3d = 0.
@@ -95 +101 @@
          DO ji = 1, jpoce
             zflx(ji,jw) = ( pwcp(ji,1,jw) - pwcp_dta(ji,jw) ) &
-               &         * 1.e3 / 1.e2 * dzkbot(ji) / rDt_trc
+               &         * 1.e3 * ( 1.e-2 * dzkbot(ji) ) / 1.E4 / rDt_trc
+         ENDDO
+      ENDDO
+
+      ! Calculation of fluxes g/cm2/s
+      DO js = 1, jpsol
+         zrate =  1.0 / rDt_trc
+         DO ji = 1, jpoce
+            zflx(ji,jpwat+js) = zflx(ji,jpwat+js) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate
          ENDDO
       ENDDO
@@ -101 +115 @@
       ! Calculation of accumulation rate per dt
       DO js = 1, jpsol
-         zrate =  1.0 / ( denssol * por1(jpksed) ) / rDt_trc
+         zrate =  1.0 / rDt_trc
          DO ji = 1, jpoce
-            zflx(ji,jpwatp1) = zflx(ji,jpwatp1) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate
+            zflx(ji,jptrased+1) = zflx(ji,jptrased+1) + ( tosed(ji,js) - fromsed(ji,js) ) * zrate
          ENDDO
       ENDDO
 
-      DO jn = 1, jpdia2dsed - 1 
+      DO jn = 1, jpdia2dsed - 2 
          CALL unpack_arr( jpoce, flxsedi2d(1:jpi,1:jpj,jn), iarroce(1:jpoce), zflx(1:jpoce,jn)  )
       END DO
+
       zflx(:,1) = dzdep(:) / dtsed
-      CALL unpack_arr( jpoce, flxsedi2d(1:jpi,1:jpj,jpdia2dsed), iarroce(1:jpoce), zflx(1:jpoce,1) )
+      CALL unpack_arr( jpoce, flxsedi2d(1:jpi,1:jpj,jpdia2dsed-1), iarroce(1:jpoce), zflx(1:jpoce,1) )
 
-       ! Start writing data
-       ! ---------------------
-       DO jn = 1, jptrased
-          cltra = sedtrcd(jn) ! short title for 3D diagnostic
-          CALL iom_put( cltra, trcsedi(:,:,:,jn) )
-       END DO
+      CALL unpack_arr( jpoce, flxsedi2d(1:jpi,1:jpj,jpdia2dsed), iarroce(1:jpoce), rstepros(1:jpoce) )
+      !
+!      CALL lbc_lnk( 'sedwri', trcsedi(:,:,:,:), 'T', 1._wp )
+!      CALL lbc_lnk( 'sedwri', flxsedi3d(:,:,:,:), 'T', 1._wp )
+!      CALL lbc_lnk( 'sedwri', flxsedi2d(:,:,:), 'T', 1._wp )
 
-       DO jn = 1, jpdia3dsed
-          cltra = seddia3d(jn) ! short title for 3D diagnostic
-          CALL iom_put( cltra, flxsedi3d(:,:,:,jn) )
-       END DO
+      ! Start writing data
+      ! ---------------------
+      DO jn = 1, jptrased
+         cltra = sedtrcd(jn) ! short title for 3D diagnostic
+         CALL iom_put( cltra, trcsedi(:,:,:,jn) )
+      END DO
 
-       DO jn = 1, jpdia2dsed
-          cltra = seddia2d(jn) ! short title for 2D diagnostic
-          CALL iom_put( cltra, flxsedi2d(:,:,jn) )
-       END DO
+      DO jn = 1, jpdia3dsed
+         cltra = seddia3d(jn) ! short title for 3D diagnostic
+         CALL iom_put( cltra, flxsedi3d(:,:,:,jn) )
+      END DO
 
-
-      DEALLOCATE( zdta )    ;   DEALLOCATE( zflx )
+      DO jn = 1, jpdia2dsed
+         cltra = seddia2d(jn) ! short title for 2D diagnostic
+         CALL iom_put( cltra, flxsedi2d(:,:,jn) )
+      END DO
 
       IF( ln_timing )  CALL timing_stop('sed_wri')