Changeset 4575

Timestamp:

2014-03-26T09:35:49+01:00 (10 years ago)

Author:

pabouttier

Message:

Cosmetic changes in dynldf_bilap_tam

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/dynldf_bilap_tam.F90 (modified) (10 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/dynldf_bilap_tam.F90

@@ -266 +266 @@
       !!   9.0  !  04-08  (C. Talandier) New trends organization
       !! History of the adjoint routine
+      !!   9.0  !  09-12 (F. Vigilant) adjoint of 9.0
       !!   3.4  !  12-07  (P.-A. Bouttier) 3.4 version
       !!----------------------------------------------------------------------
@@ -293 +294 @@
       zuaad = 0.0_wp
       zvaad = 0.0_wp
+
       zufad(:,:,:) = 0.0_wp
       zutad(:,:,:) = 0.0_wp
       zluad(:,:,:) = 0.0_wp
       zlvad(:,:,:) = 0.0_wp
+
       zcvad(:,:) = 0.0_wp
       zcuad(:,:) = 0.0_wp
 
       DO jk = 1, jpkm1
+
          ! Bilaplacian
          ! -----------
+
          DO jj = jpjm1, 2, -1
             DO ji = fs_jpim1, fs_2, -1   ! vector opt.
@@ -315 +320 @@
                zutad(ji  ,jj  ,jk) = zutad(ji  ,jj  ,jk) - zvaad / e2v(ji,jj)
                zutad(ji  ,jj+1,jk) = zutad(ji  ,jj+1,jk) + zvaad / e2v(ji,jj)
-               !
+
                zufad(ji  ,jj  ,jk) = zufad(ji  ,jj  ,jk) - zuaad / ze2u
                zufad(ji  ,jj-1,jk) = zufad(ji  ,jj-1,jk) + zuaad / ze2u
                zutad(ji  ,jj  ,jk) = zutad(ji  ,jj  ,jk) - zuaad / e1u(ji,jj)
                zutad(ji+1,jj  ,jk) = zutad(ji+1,jj  ,jk) + zuaad / e1u(ji,jj)
-               !
+               
                zuaad = 0.0_wp
                zvaad = 0.0_wp
             END DO
          END DO
+
          !                                             ! ===============
       END DO                                           !   End of slab
       !                                                ! ===============
+
       ! boundary conditions on the laplacian curl and div (zuf,zut)
 !!bug gm no need to do this 2 following lbc...
       CALL lbc_lnk_adj( zutad, 'T', 1.0_wp )
       CALL lbc_lnk_adj( zufad, 'F', 1.0_wp )
+
       DO jk = 1, jpkm1
+
          ! Third derivative
          ! ----------------
+
          ! Laplacian divergence
          DO jj = jpj, 2, -1
@@ -343 +353 @@
                zlvad(ji  ,jj  ,jk) = zlvad(ji  ,jj  ,jk) + zutad(ji,jj,jk) / zbt
                zlvad(ji  ,jj-1,jk) = zlvad(ji  ,jj-1,jk) - zutad(ji,jj,jk) / zbt
+
                zutad(ji,jj,jk) = 0.0_wp
             END DO
          END DO
+
          ! Laplacian Horizontal fluxes
          DO jj = jpjm1, 1, -1
@@ -353 +365 @@
             END DO
          END DO
+
          ! Laplacian curl ( * e3f if s-coordinates or z-coordinate with partial steps)
          DO jj = jpjm1, 1, -1
@@ -361 +374 @@
                zcuad(ji  ,jj  ) = zcuad(ji  ,jj  ) + zufad(ji,jj,jk)
                zcuad(ji  ,jj+1) = zcuad(ji  ,jj+1) - zufad(ji,jj,jk)
+
                zufad(ji,jj,jk) = 0.0_wp
             END DO
          END DO
+
          ! Contravariant "laplacian"
          DO jj = 1, jpj
@@ -373 +388 @@
             END DO
          END DO
+
          ! Multiply by the eddy viscosity coef. (at u- and v-points)
          zluad(:,:,jk) = zluad(:,:,jk) * fsahmu(:,:,jk)
          zlvad(:,:,jk) = zlvad(:,:,jk) * fsahmv(:,:,jk)
+
       END DO
+
       ! Boundary conditions on the laplacian  (zlu,zlv)
       CALL lbc_lnk_adj( zlvad, 'V', -1.0_wp )
       CALL lbc_lnk_adj( zluad, 'U', -1.0_wp )
+
       !                                                ! ===============
       DO jk = 1, jpkm1                                 ! Horizontal slab
@@ -385 +404 @@
          ! Laplacian
          ! ---------
+
          IF( ln_sco .OR. ln_zps ) THEN   ! s-coordinate or z-coordinate with partial steps
             DO jj = jpjm1, 2, -1
@@ -393 +413 @@
                   hdivb_ad(ji  ,jj+1,jk) = hdivb_ad(ji  ,jj+1,jk) + zlvad(ji,jj,jk) / e2v(ji,jj)
                   zlvad(ji,jj,jk) = 0.0_wp
+
                   zufad   (ji  ,jj  ,jk) = zufad   (ji  ,jj  ,jk) - zluad(ji  ,jj  ,jk) / ( e2u(ji,jj) * fse3u(ji,jj,jk) )
                   zufad   (ji  ,jj-1,jk) = zufad   (ji  ,jj-1,jk) + zluad(ji  ,jj  ,jk) / ( e2u(ji,jj) * fse3u(ji,jj,jk) )
@@ -409 +430 @@
                   hdivb_ad(ji  ,jj+1,jk) = hdivb_ad(ji  ,jj+1,jk) + zlvad(ji,jj,jk) / e2v(ji,jj)
                   zlvad(ji,jj,jk) = 0.0_wp
+                  
                   rotb_ad (ji  ,jj  ,jk) = rotb_ad (ji  ,jj  ,jk) - zluad(ji,jj,jk) / e2u(ji,jj)
                   rotb_ad (ji  ,jj-1,jk) = rotb_ad (ji  ,jj-1,jk) + zluad(ji,jj,jk) / e2u(ji,jj)