Ignore:
Timestamp:
2019-08-20T14:14:56+02:00 (14 months ago)
Author:
acc
Message:

Branch: dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap. Minor bugfix in step.F90 to enable AGRIF SETTE tests to run. Also merged prettification changes to NST routines from the dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps branch. Neither AGRIF_DEMO nor VORTEX restart perfectly (drifting after 8 and 121 timesteps, respectively).

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_update.F90

    r10989 r11463  
    8484               DO jj=j1,j2 
    8585                  DO ji=i1,i2 
    86                      tabres(ji,jj,jk,jn) = (tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) & 
    87                                            * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
     86                     tabres(ji,jj,jk,jn) = (    tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) & 
     87                       &                   * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
    8888                  END DO 
    8989               END DO 
     
    9393            DO jj=j1,j2 
    9494               DO ji=i1,i2 
    95                   tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) & 
    96                                            + (tmask(ji,jj,jk)-1)*999._wp 
     95                  tabres(ji,jj,jk,n2) =   tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) & 
     96                    &                  + (tmask(ji,jj,jk)-1)*999._wp 
    9797               END DO 
    9898            END DO 
     
    138138                     DO ji=i1,i2 
    139139                        IF( tabres_child(ji,jj,jk,jn) .NE. 0. ) THEN 
    140                            tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb) &  
    141                                  & + atfp * ( tabres_child(ji,jj,jk,jn) & 
    142                                  &          - tr(ji,jj,jk,jn,Kmm) ) * tmask(ji,jj,jk) 
     140                           tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb)                     &  
     141                             &                  + atfp * ( tabres_child(ji,jj,jk,jn)     & 
     142                             &                                     - tr(ji,jj,jk,jn,Kmm) & 
     143                             &                           ) *      tmask(ji,jj,jk) 
    143144                        ENDIF 
    144145                     ENDDO 
     
    183184                  DO ji=i1,i2 
    184185!> jc tmp 
    185                      tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm)  * e3t(ji,jj,jk,Kmm) / e3t_0(ji,jj,jk) 
    186 !                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) 
     186                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) / e3t_0(ji,jj,jk) 
     187!                    tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) 
    187188!< jc tmp 
    188189                  END DO 
     
    194195         DO jn = n1,n2 
    195196            tabres(i1:i2,j1:j2,k1:k2,jn) =  tabres(i1:i2,j1:j2,k1:k2,jn) * e3t_0(i1:i2,j1:j2,k1:k2) & 
    196                                          & * tmask(i1:i2,j1:j2,k1:k2) 
     197              &                            * tmask(i1:i2,j1:j2,k1:k2) 
    197198         ENDDO 
    198199!< jc tmp 
     
    204205                     DO ji=i1,i2 
    205206                        IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 
    206                            ztb  = tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used 
     207                           ztb  =     tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used 
    207208                           ztnu = tabres(ji,jj,jk,jn) 
    208                            ztno = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs) 
    209                            tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) )  &  
    210                                      &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb) 
     209                           ztno =     tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs) 
     210                           tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) ) * tmask(ji,jj,jk)     & 
     211                             &                                                   /  e3t(ji,jj,jk,Kbb) 
    211212                        ENDIF 
    212213                     ENDDO 
Note: See TracChangeset for help on using the changeset viewer.