Changeset 4609 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z

Timestamp:

2014-03-31T15:06:20+02:00 (10 years ago)

Author:

cetlod

Message:

bugfix on PISCES/P2Z, see ticket #1298

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

@@ -80 +80 @@
 
       ! Allocate temporary workspace
-                       CALL wrk_alloc( jpi, jpj, jpk, zwork, ztra )
-      IF( ln_diatrc )  CALL wrk_alloc( jpi, jpj, zw2d )
+      CALL wrk_alloc( jpi, jpj, jpk, zwork, ztra )
       IF( l_trdtrc ) THEN
          CALL wrk_alloc( jpi, jpj, jpk, ztrbio )
@@ -110 +109 @@
 
       IF( ln_diatrc ) THEN 
-         DO jk = 1, jpkm1
+         CALL wrk_alloc( jpi, jpj, zw2d )
+         zw2d(:,:) =  ztra(:,:,1) * fse3t(:,:,1) * 86400.
+         DO jk = 2, jpkm1
             zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400.
          END DO
@@ -121 +122 @@
       ENDIF
       !
-      IF( ln_diabio )  trbio(:,:,:,jp_pcs0_trd + 7) = ztra(:,:,:)
+      IF( ln_diabio .AND. .NOT. lk_iomput )  trbio(:,:,:,jp_pcs0_trd + 7) = ztra(:,:,:)
       CALL wrk_dealloc( jpi, jpj, jpk, zwork, ztra )
       !