Ignore:
Timestamp:
2015-01-07T19:03:53+01:00 (6 years ago)
Author:
hliu
Message:

upload the modifications for W/D based on r:4826

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r4826_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r4795 r5014  
    145145      fse3t_a(:,:,jpk) = e3t_0(:,:,jpk) 
    146146 
     147      IF(ln_wd) THEN 
     148        DO jj = 1, jpj 
     149          DO ji = 1, jpi 
     150            IF(e3t_0(ji,jj,1) < 0._wp) fse3t_a(ji,jj,1:2) = 0.5_wp * rn_wdmin1  
     151          ENDDO 
     152        ENDDO 
     153      END IF 
     154 
    147155      ! Reconstruction of all vertical scale factors at now and before time steps 
    148156      ! ============================================================================= 
     
    793801      CHARACTER(len=*), INTENT(in) ::   cdrw   ! "READ"/"WRITE" flag 
    794802      !! * Local declarations 
    795       INTEGER ::   jk 
     803      INTEGER ::   ji, jj, jk 
    796804      INTEGER ::   id1, id2, id3, id4, id5     ! local integers 
    797805      !!---------------------------------------------------------------------- 
     
    867875            fse3t_n(:,:,:) = e3t_0(:,:,:) 
    868876            sshn(:,:) = 0.0_wp 
     877 
     878            IF(ln_wd) THEN 
     879              DO jj = 1, jpj 
     880                DO ji = 1, jpi 
     881                  IF(e3t_0(ji,jj,1) < 0._wp) THEN 
     882                    fse3t_b(ji,jj,1:2) = 0.5_wp * rn_wdmin1  
     883                    fse3t_n(ji,jj,1:2) = 0.5_wp * rn_wdmin1  
     884                    sshb(ji,jj) = rn_wdmin1 - bathy(ji,jj) 
     885                    sshn(ji,jj) = rn_wdmin1 - bathy(ji,jj) 
     886                  ENDIF 
     887                ENDDO 
     888              ENDDO 
     889            END IF 
     890 
    869891            IF( ln_vvl_ztilde .OR. ln_vvl_layer) THEN 
    870892               tilde_e3t_b(:,:,:) = 0.0_wp 
Note: See TracChangeset for help on using the changeset viewer.