Changeset 9006


Ignore:
Timestamp:
2017-12-13T12:14:44+01:00 (3 years ago)
Author:
timgraham
Message:

Modification to updates

Location:
branches/2017/dev_METO_MERCATOR_2017_agrif/NEMOGCM/NEMO/NST_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_METO_MERCATOR_2017_agrif/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90

    r8998 r9006  
    309309               DO jj=j1,j2 
    310310                  DO ji=i1,i2 
    311                      tabres(ji,jj,jk,jn) = (tsn(ji,jj,jk,jn) * e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) & 
    312                                            * tmask(ji,jj,jk) * zrho_xy + (tmask(ji,jj,jk)-1)*999._wp 
     311                     tabres(ji,jj,jk,jn) = (tsn(ji,jj,jk,jn) * e3t_n(ji,jj,jk) ) & 
     312                                           * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
    313313                  END DO 
    314314               END DO 
     
    318318            DO jj=j1,j2 
    319319               DO ji=i1,i2 
    320                   tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e1e2t(ji,jj) * e3t_n(ji,jj,jk) * zrho_xy & 
     320                  tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t_n(ji,jj,jk) & 
    321321                                           + (tmask(ji,jj,jk)-1)*999._wp 
    322322               END DO 
     
    333333                  N_in = N_in + 1 
    334334                  tabin(jk,:) = tabres(ji,jj,jk,n1:n2-1)/tabres(ji,jj,jk,n2) 
    335                   h_in(N_in) = tabres(ji,jj,jk,n2)/e1e2t(ji,jj) 
     335                  h_in(N_in) = tabres(ji,jj,jk,n2) 
    336336               ENDDO 
    337337               N_out = 0 
     
    339339                  IF (tmask(ji,jj,jk) < -900) EXIT ! TODO: Will not work with ISF 
    340340                  N_out = N_out + 1 
    341                   h_out(N_out) = e3t_n(ji,jj,jk) !Parent grid scale factors. Could multiply by e1e2t here instead of division above 
     341                  h_out(N_out) = e3t_n(ji,jj,jk)  
    342342               ENDDO 
    343343               IF (N_in > 0) THEN !Remove this? 
  • branches/2017/dev_METO_MERCATOR_2017_agrif/NEMOGCM/NEMO/NST_SRC/agrif_top_update.F90

    r9005 r9006  
    8383                  DO ji=i1,i2 
    8484                     tabres(ji,jj,jk,jn) = (trn(ji,jj,jk,jn) * e3t_n(ji,jj,jk) ) & 
    85                                            * tmask(ji,jj,jk) * zrho_xy + (tmask(ji,jj,jk)-1)*999._wp 
     85                                           * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 
    8686                  END DO 
    8787               END DO 
Note: See TracChangeset for help on using the changeset viewer.