Changeset 6678


Ignore:
Timestamp:
2016-06-09T15:37:05+02:00 (4 years ago)
Author:
hliu
Message:

remove one bug in Wetting/Drying? part

Location:
branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r6603 r6678  
    10081008         !                                   ! Sum sea level 
    10091009         ssha(:,:) = ssha(:,:) + za1 * ssha_e(:,:) 
    1010          IF(ln_wd) ssha(:,:) = sshai(:,:) 
     1010         IF(ln_wd) ssha(:,:) = ssha(:,:) * (1 - wdmask(:,:)) + sshai(:,:) * wdmask(:,:) 
    10111011         !                                                 ! ==================== ! 
    10121012      END DO                                               !        end loop      ! 
  • branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/wet_dry.F90

    r6152 r6678  
    132132        
    133133        zflag  = 0 
    134         zdepwd = 50._wp   !maximum depth on which that W/D could possibly happen 
     134        zdepwd = 150._wp   !maximum depth on which that W/D could possibly happen 
    135135 
    136136        
     
    201201           
    202202                 IF(zdep1 > zdep2) THEN 
     203                   IF(jk1 .eq. 1) wdmask(ji,jj) = 1 
    203204                   zflag = 1 
    204                    wdmask(ji, jj) = 1 
    205205                   zcoef = ( ( zdep2 - rn_wdmin2 ) * ztmp - zzflxn * z2dt ) / ( zflxp(ji,jj) * z2dt ) 
    206206                   zcoef = max(zcoef, 0._wp) 
     
    209209                   IF(zflxu1(ji-1,jj) < 0._wp) zwdlmtu(ji-1,jj) = zcoef 
    210210                   IF(zflxv1(ji,  jj) > 0._wp) zwdlmtv(ji  ,jj) = zcoef 
    211                    IF(zflxv1(ji,jj-1) < 0._wp) zwdlmtv(ji-1,jj) = zcoef 
     211                   IF(zflxv1(ji,jj-1) < 0._wp) zwdlmtv(ji,jj-1) = zcoef 
    212212                 END IF 
    213213              END DO ! ji loop 
     
    332332              DO ji = 2, jpim1  
    333333         
    334                  wdmask(ji,jj) = 0 
     334                 !wdmask(ji,jj) = 0 
    335335                 IF(tmask(ji, jj, 1) < 0.5_wp) CYCLE  
    336336                 IF(bathy(ji,jj) > zdepwd) CYCLE 
     
    356356                   IF(zflxu1(ji-1,jj) < 0._wp) zwdlmtu(ji-1,jj) = zcoef 
    357357                   IF(zflxv1(ji,  jj) > 0._wp) zwdlmtv(ji  ,jj) = zcoef 
    358                    IF(zflxv1(ji,jj-1) < 0._wp) zwdlmtv(ji-1,jj) = zcoef 
     358                   IF(zflxv1(ji,jj-1) < 0._wp) zwdlmtv(ji,jj-1) = zcoef 
    359359                 END IF 
    360360              END DO ! ji loop 
Note: See TracChangeset for help on using the changeset viewer.