Ignore:
Timestamp:
2020-06-16T12:07:28+02:00 (3 months ago)
Author:
rblod
Message:

ticket #2129 : major corrections in domcfg

File:
1 edited

Legend:

Unmodified
Added
Removed
  • utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_dom_update.F90

    r13056 r13109  
    2222      IF( Agrif_Root() ) return 
    2323 
    24       CALL agrif_update_variable(bottom_level_id,locupdate=(/npt_copy,0/),procname = update_bottom_level) 
     24      CALL agrif_update_variable(bottom_level_id,procname = update_bottom_level) 
    2525      ! 
    2626      Agrif_UseSpecialValueInUpdate = .TRUE. 
     
    4141      ! 
    4242      !!----------------------------------------------------------------------  
     43      REAL(WP),DIMENSION(jpi,jpj) :: zk 
    4344      ! 
    4445      IF( before) THEN 
     
    4950         WHERE ( mbkt(i1:i2,j1:j2) .EQ. 0 ) 
    5051            ssmask(i1:i2,j1:j2) = 0. 
     52            mbkt(i1:i2,j1:j2)   = 1 
    5153         ELSEWHERE 
    5254            ssmask(i1:i2,j1:j2) = 1. 
    53          END WHERE            
     55         END WHERE  
     56         zk(:,:) = REAL(mbkt(:,:),wp); CALL lbc_lnk('update_bottom',zk,'T',1.); mbkt(:,:) = MAX(NINT(zk(:,:)),1) 
     57         CALL lbc_lnk('update_bottom',ssmask,'T',1.)           
    5458      ENDIF 
    5559      ! 
     
    7175            DO jj=j1,j2 
    7276               DO ji=i1,i2 
    73                    IF( mbkt(ji,jj) .LE. jk ) THEN 
     77                   IF( mbkt(ji,jj) .GE. jk ) THEN 
    7478                      tabres(ji,jj,jk) = e3t_0(ji,jj,jk) 
    7579                   ELSE 
     
    8387            DO jj=j1,j2 
    8488               DO ji=i1,i2 
    85                    IF( mbkt(ji,jj) .LE.jk ) THEN 
     89                   IF( mbkt(ji,jj) .GE. jk ) THEN 
    8690                      e3t_0(ji,jj,jk) = MAX(tabres(ji,jj,jk),MIN(e3zps_min,e3t_1d(jk)*e3zps_rat)) 
    8791                   ELSE 
     
    9195            END DO 
    9296         END DO 
     97 
     98         CALL lbc_lnk('update_e3t',e3t_0,'T',1.) 
    9399         ! 
    94100      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.