Changeset 14715

Timestamp:

2021-04-15T13:41:16+02:00 (3 years ago)

Author:

clem

Message:

solve ticket #2641

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/OCE/TRA/traqsr.F90 (modified) (1 diff)

NEMO/releases/r4.0/r4.0-HEAD/src/OCE/TRA/traqsr.F90

@@ -281 +281 @@
          zetot(:,:,nksr+1:jpk) = 0._wp     ! below ~400m set to zero
          DO jk = nksr, 1, -1
-            zetot(:,:,jk) = zetot(:,:,jk+1) + qsr_hc(:,:,jk) * rau0_rcp
-         END DO         
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1
+                  zetot(ji,jj,jk) = zetot(ji,jj,jk+1) + qsr_hc(ji,jj,jk) * rau0_rcp
+               ENDDO
+            ENDDO
+         END DO
+         CALL lbc_lnk( 'traqsr', zetot, 'T', 1._wp )         
          CALL iom_put( 'qsr3d', zetot )   ! 3D distribution of shortwave Radiation
          DEALLOCATE( zetot )