Changeset 13352


Ignore:
Timestamp:
2020-07-29T10:33:53+02:00 (13 days ago)
Author:
jchanut
Message:

#2222, #2129: correct depths used in linear vertical interpolation. NB: these are also used for initial state interpolation/extrapolation with identical vertical grids

Location:
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90

    r13351 r13352  
    883883                  z_in(1) = 0.5_wp * h_in(1) - zhtot + ht0_parent(ji,jj) 
    884884                  DO jk=2,N_in 
    885                      z_in(jk) = z_in(jk-1) + 0.5_wp * h_in(jk) 
     885                     z_in(jk) = z_in(jk-1) + 0.5_wp * (h_in(jk-1)+h_in(jk)) 
    886886                  END DO 
    887887 
     
    895895                  z_out(1) = 0.5_wp * h_out(1) - SUM(h_out(1:N_out)) + ht_0(ji,jj) 
    896896                  DO jk=2,N_out 
    897                      z_out(jk) = z_out(jk-1) + 0.5_wp * h_out(jk) 
     897                     z_out(jk) = z_out(jk-1) + 0.5_wp * (h_out(jk)+h_out(jk-1)) 
    898898                  END DO 
    899899 
     
    10341034                 z_in(1) = 0.5_wp * h_in(1) - zhtot + hu0_parent(ji,jj)  
    10351035                 DO jk=2,N_in 
    1036                     z_in(jk) = z_in(jk-1) + 0.5_wp * h_in(jk) 
     1036                    z_in(jk) = z_in(jk-1) + 0.5_wp * (h_in(jk)+h_in(jk-1)) 
    10371037                 END DO 
    10381038                      
     
    10461046                 z_out(1) = 0.5_wp * h_out(1) - SUM(h_out(1:N_out)) + hu_0(ji,jj) 
    10471047                 DO jk=2,N_out 
    1048                     z_out(jk) = z_out(jk-1) + 0.5_wp * h_out(jk)  
     1048                    z_out(jk) = z_out(jk-1) + 0.5_wp * (h_out(jk-1) + h_out(jk))  
    10491049                 END DO   
    10501050 
     
    11571157                  z_in(1) = 0.5_wp * h_in(1) - zhtot + hv0_parent(ji,jj) 
    11581158                  DO jk=2,N_in 
    1159                      z_in(jk) = z_in(jk-1) + 0.5_wp * h_in(jk) 
     1159                     z_in(jk) = z_in(jk-1) + 0.5_wp * (h_in(jk-1)+h_in(jk)) 
    11601160                  END DO 
    11611161 
     
    11691169                  z_out(1) = 0.5_wp * h_out(1) - SUM(h_out(1:N_out)) + hv_0(ji,jj) 
    11701170                  DO jk=2,N_out 
    1171                      z_out(jk) = z_out(jk-1) + 0.5_wp * h_out(jk) 
     1171                     z_out(jk) = z_out(jk-1) + 0.5_wp * (h_out(jk-1)+h_out(jk)) 
    11721172                  END DO 
    11731173  
  • NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_top_interp.F90

    r13351 r13352  
    140140                  z_in(1) = 0.5_wp * h_in(1) - zhtot + ht0_parent(ji,jj) 
    141141                  DO jk=2,N_in 
    142                      z_in(jk) = z_in(jk-1) + 0.5_wp * h_in(jk) 
     142                     z_in(jk) = z_in(jk-1) + 0.5_wp * (h_in(jk-1)+h_in(jk)) 
    143143                  END DO 
    144144 
     
    152152                  z_out(1) = 0.5_wp * h_out(1) - SUM(h_out(1:N_out)) + ht_0(ji,jj) 
    153153                  DO jk=2,N_out 
    154                      z_out(jk) = z_out(jk-1) + 0.5_wp * h_out(jk) 
     154                     z_out(jk) = z_out(jk-1) + 0.5_wp * (h_out(jk-1)+h_out(jk)) 
    155155                  END DO 
    156156 
Note: See TracChangeset for help on using the changeset viewer.