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 11774 for NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPATAM_SRC/TAM/trj_tam.F90 – NEMO

Ignore:
Timestamp:
2019-10-23T16:30:17+02:00 (4 years ago)
Author:
smueller
Message:

Correction of the state-variable advancement in adjoint mode (application of the patch attached to ticket #1443)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPATAM_SRC/TAM/trj_tam.F90

    r3662 r11774  
    578578         &                   '  = ', zwtr1, zwtr2 
    579579 
    580       IF ( kstp /= nit000-1 ) THEN 
     580      IF ( ( kstp /= nit000-1 ).AND.( kdir == 1 ) ) THEN 
    581581         tsb(:,:,:,:) = tsn(:,:,:,:) 
    582582         ub(:,:,:) = un(:,:,:) 
     
    591591      tsn(:,:,:,jp_tem)     = zwtr1 * tnr1    (:,:,:) + zwtr2 * tnr2    (:,:,:) 
    592592      tsn(:,:,:,jp_sal)     = zwtr1 * snr1    (:,:,:) + zwtr2 * snr2    (:,:,:) 
     593      IF ( kdir == -1 ) THEN 
     594         zwtr1  = ( stpr2 - zstp + 1  ) * zden 
     595         zwtr2  = ( zstp - 1 - stpr1 ) * zden 
     596         ub(:,:,:)     = zwtr1 * unr1    (:,:,:) + zwtr2 * unr2    (:,:,:) 
     597         vb(:,:,:)     = zwtr1 * vnr1    (:,:,:) + zwtr2 * vnr2    (:,:,:) 
     598         tsb(:,:,:,jp_tem)     = zwtr1 * tnr1    (:,:,:) + zwtr2 * tnr2    (:,:,:) 
     599         tsb(:,:,:,jp_sal)     = zwtr1 * snr1    (:,:,:) + zwtr2 * snr2    (:,:,:) 
     600         IF(lwp)WRITE(numout,*) ' before lin. interp. coeff.', & 
     601           &                   '  = ', zwtr1, zwtr2 
     602 
     603         zwtr1  = ( stpr2 - zstp      ) * zden 
     604         zwtr2  = ( zstp  - stpr1     ) * zden       
     605      END IF 
    593606      IF ( kstp == nit000-1 ) THEN 
    594607         tsb(:,:,:,:) = tsn(:,:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.