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_sponge.F90

    r10989 r11463  
    108108                  N_in = N_in + 1 
    109109                  tabin(jk,:) = tabres(ji,jj,jk,n1:n2-1) 
    110                   h_in(N_in) = tabres(ji,jj,jk,n2) 
     110                  h_in(N_in)  = tabres(ji,jj,jk,n2) 
    111111               END DO 
    112112               N_out = 0 
    113113               DO jk=1,jpk ! jpk of child grid 
    114114                  IF (tmask(ji,jj,jk) == 0) EXIT  
    115                   N_out = N_out + 1 
     115                  N_out     = N_out + 1 
    116116                  h_out(jk) = e3t(ji,jj,jk,Kmm) !Child grid scale factors. Could multiply by e1e2t here instead of division above 
    117117               ENDDO 
    118118               IF (N_in > 0) THEN 
    119                   h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
     119                  h_diff = SUM( h_out(1:N_out) ) - SUM( h_in(1:N_in) ) 
    120120                  tabres(ji,jj,k2,:) = tabres(ji,jj,k2-1,:) !what is this line for????? 
    121121                  DO jn=1,jptra 
     
    133133                  trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) - tabres_child(ji,jj,jk,1:jptra) 
    134134# else 
    135                   trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) - tabres(ji,jj,jk,1:jptra) 
     135                  trbdiff(ji,jj,jk,1:jptra) = tr(ji,jj,jk,1:jptra,Kbb) -       tabres(ji,jj,jk,1:jptra) 
    136136# endif 
    137137               ENDDO 
     
    153153                  DO ji = i1+1,i2-1 
    154154                     IF( .NOT. tabspongedone_trn(ji,jj) ) THEN  
    155                         tr(ji,jj,jk,jn,Krhs) = tr(ji,jj,jk,jn,Krhs) + (  ztu(ji,jj) - ztu(ji-1,jj  )     & 
    156                            &                                   + ztv(ji,jj) - ztv(ji  ,jj-1)  )  & 
    157                            &                                * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
     155                        tr(ji,jj,jk,jn,Krhs) = tr(ji,jj,jk,jn,Krhs) + (  ztu(ji,jj) - ztu(ji-1,jj  )       & 
     156                           &                                           + ztv(ji,jj) - ztv(ji  ,jj-1)  )    & 
     157                           &                                          * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
    158158                     ENDIF 
    159159                  END DO 
Note: See TracChangeset for help on using the changeset viewer.