New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4582 – NEMO

Changeset 4582


Ignore:
Timestamp:
2014-03-26T10:30:54+01:00 (10 years ago)
Author:
pabouttier
Message:

Correct a wrong adjoint computation in ssh_wzv_adj, see Ticket #1276

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/sshwzv_tam.F90

    r3611 r4582  
    269269      !                                                ! integrate from the bottom the hor. divergence 
    270270      DO jk = 1, jpkm1 
    271          hdivn_ad(:,:,jk  ) = hdivn_ad(:,:,jk  ) - fse3t_n(:,:,jk) * wn_ad(:,:,jk) 
    272          wn_ad(   :,:,jk+1) = wn_ad(:,:,jk+1) + wn_ad(:,:,jk) 
     271         hdivn_ad(:,:,jk  ) = hdivn_ad(:,:,jk  ) - wn_ad(:,:,jk) * fse3t_n(:,:,jk) 
     272         wn_ad(   :,:,jk+1) = wn_ad(   :,:,jk+1) + wn_ad(:,:,jk) 
    273273         wn_ad(   :,:,jk  ) = 0.0_wp 
    274274      END DO 
     
    290290         ENDIF 
    291291         !                                                ! Sea surface elevation time stepping 
    292          sshb_ad(:,:) = sshb_ad(:,:) + ssha_ad(:,:)* tmask(:,:,1) 
    293          emp_ad( :,:) = emp_ad(:,:)  - z2dt * z1_rau0 * ssha_ad(:,:) * tmask(:,:,1) 
    294          emp_b_ad( :,:) = emp_b_ad(:,:) - z2dt * z1_rau0 * ssha_ad(:,:) * tmask(:,:,1) 
    295          zhdivad(:,:) = zhdivad(:,:) - z2dt * tmask(:,:,1) * ssha_ad(:,:) 
     292         ssha_ad( :,:) = ssha_ad( :,:) * tmask(:,:,1) 
     293         sshb_ad( :,:) = sshb_ad( :,:) + ssha_ad(:,:) 
     294         ssha_ad( :,:) = ssha_ad( :,:) * z2dt 
     295         zhdivad( :,:) = zhdivad( :,:) - ssha_ad(:,:) 
     296         ssha_ad( :,:) = ssha_ad( :,:) * z1_rau0  
     297         emp_ad(  :,:) = emp_ad(  :,:) - ssha_ad(:,:) 
     298         emp_b_ad(:,:) = emp_b_ad(:,:) - ssha_ad(:,:) 
    296299         ssha_ad(:,:) = 0.0_wp 
    297300 
     
    307310      ! 
    308311      ! 
    309       wn_ad(:,:,jpk) = 0._wp 
     312      IF( kt == nit000 ) wn_ad(:,:,jpk) = 0._wp 
    310313      ! 
    311314      ! 
Note: See TracChangeset for help on using the changeset viewer.