New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14019 – NEMO

Changeset 14019


Ignore:
Timestamp:
2020-12-02T18:50:19+01:00 (4 years ago)
Author:
jchanut
Message:

#2222, use right thicknesses for initial state interpolation from parent

File:
1 edited

Legend:

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

    r14013 r14019  
    10791079                     IF ( l_vremap ) THEN 
    10801080                        h_in(jk) = e3u0_parent(ji,jj,jk)  
    1081                         zhtot = zhtot + h_in(jk) 
     1081                     ELSE 
     1082                        IF (jk==N_in) THEN 
     1083                           h_in(jk) = hu0_parent(ji,jj) + ptab(ji,jj,k2,2) - zhtot 
     1084                        ELSE 
     1085                           h_in(jk) = ptab(ji,jj,jk,2)/(e2u(ji,jj)*zrhoy)  
     1086                        ENDIF 
    10821087                     ENDIF 
     1088                     zhtot = zhtot + h_in(jk) 
    10831089                     IF( h_in(jk) .GT. 0. ) THEN 
    10841090                     tabin(jk) = ptab(ji,jj,jk,1)/(e2u(ji,jj)*zrhoy*h_in(jk)) 
     
    12041210                     IF (l_vremap) THEN 
    12051211                        h_in(jk) = e3v0_parent(ji,jj,jk) 
    1206                         zhtot = zhtot + h_in(jk) 
     1212                     ELSE 
     1213                        IF (jk==N_in) THEN 
     1214                           h_in(jk) = hv0_parent(ji,jj) + ptab(ji,jj,k2,2) - zhtot 
     1215                        ELSE 
     1216                           h_in(jk) = ptab(ji,jj,jk,2)/(e1v(ji,jj)*zrhox)  
     1217                        ENDIF 
    12071218                     ENDIF 
     1219                     zhtot = zhtot + h_in(jk) 
    12081220                     IF( h_in(jk) .GT. 0. ) THEN 
    12091221                       tabin(jk) = ptab(ji,jj,jk,1)/(e1v(ji,jj)*zrhox*h_in(jk)) 
Note: See TracChangeset for help on using the changeset viewer.