Ignore:
Timestamp:
2014-02-06T17:47:57+01:00 (7 years ago)
Author:
jchanut
Message:

Missing Asselin correction with Agrif, #1203

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/NST_SRC/agrif_top_update.F90

    r3680 r4491  
    7878            ENDDO 
    7979         ELSE 
     80            IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN 
     81            ! Add asselin part 
     82               DO jn = n1, n2 
     83                  DO jk = k1, k2 
     84                     DO jj = j1, j2 
     85                        DO ji = i1, i2 
     86                           IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 
     87                              trb(ji,jj,jk,jn) = trb(ji,jj,jk,jn) &  
     88                                 & + atfp * ( tabres(ji,jj,jk,jn) & 
     89                                               - trn(ji,jj,jk,jn) ) * tmask(ji,jj,jk) 
     90                           ENDIF 
     91                        ENDDO 
     92                     ENDDO 
     93                  ENDDO 
     94               ENDDO 
     95            ENDIF 
     96 
    8097            DO jn = n1, n2 
    8198               DO jk = k1, k2 
Note: See TracChangeset for help on using the changeset viewer.