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 10030 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OFF/dtadyn.F90 – NEMO

Ignore:
Timestamp:
2018-08-03T10:18:16+02:00 (6 years ago)
Author:
gm
Message:

#1911 (ENHANCE-04): RK3 branch - step II.3 remove e3uw_$ e3vw_$, except e3.w_0 and use only after e3 in dyn/trazdf

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OFF/dtadyn.F90

    r10009 r10030  
    332332         ! 
    333333         !                          !* BEFORE fields :  
    334          CALL ssh2e3_before               ! set:      hu , hv , r1_hu, r1_hv  
    335          !                                !      e3t, e3w, e3u, e3uw, e3v, e3vw        (from 1 to jpkm1) 
     334         CALL ssh2e3_before               ! set:       hu,  hv, r1_hu, r1_hv  
     335         !                                !      e3t, e3u, e3v,            (from 1 to jpkm1) 
     336         !                                !      e3w,                      (from 1 to jpk  ) 
    336337         ! 
    337338         !                                ! set jpk level one to the e3._0 values 
    338          e3t_b(:,:,jpk) = e3t_0(:,:,jpk)  ;   e3u_b(:,:,jpk) =  e3w_0(:,:,jpk)  ;   e3v_b(:,:,jpk) =  e3v_0(:,:,jpk) 
    339          e3w_b(:,:,jpk) = e3w_0(:,:,jpk)  ;  e3uw_b(:,:,jpk) = e3uw_0(:,:,jpk)  ;  e3vw_b(:,:,jpk) = e3vw_0(:,:,jpk) 
     339         e3t_b(:,:,jpk) = e3t_0(:,:,jpk)  ;  e3u_b(:,:,jpk) = e3u_0(:,:,jpk) 
     340                                             e3v_b(:,:,jpk) = e3v_0(:,:,jpk) 
    340341         ! 
    341342         !                          !* NOW fields :  
    342          CALL ssh2e3_now                  ! set: ht , hu , hv , r1_hu, r1_hv 
    343          !                                !      e3t, e3w, e3u, e3uw, e3v, e3vw, e3f   (from 1 to jpkm1) 
    344          !                                !      gdept_n, gdepw_n, gde3w_n 
    345 !!gm issue?   gdept_n, gdepw_n, gde3w_n never defined at jpk 
     343         CALL ssh2e3_now                  ! set:  ht,  hu,  hv, r1_hu, r1_hv 
     344         !                                !      e3t, e3u, e3v , e3f              (from 1 to jpkm1) 
     345         !                                !      e3w, gdept_n, gdepw_n, gde3w_n   (from 1 to jpk  ) 
    346346         ! 
    347347         !                                ! set one for all last level to the e3._0 value 
    348          e3t_n(:,:,jpk) = e3t_0(:,:,jpk)  ;   e3u_n(:,:,jpk) =  e3w_0(:,:,jpk)  ;   e3v_n(:,:,jpk) =  e3v_0(:,:,jpk) 
    349          e3w_n(:,:,jpk) = e3w_0(:,:,jpk)  ;  e3uw_n(:,:,jpk) = e3uw_0(:,:,jpk)  ;  e3vw_n(:,:,jpk) = e3vw_0(:,:,jpk) 
    350          e3f_n(:,:,jpk) = e3f_0(:,:,jpk) 
     348         e3t_n(:,:,jpk) = e3t_0(:,:,jpk)  ;  e3u_n(:,:,jpk) = e3u_0(:,:,jpk) 
     349         e3f_n(:,:,jpk) = e3f_0(:,:,jpk)  ;  e3v_n(:,:,jpk) = e3v_0(:,:,jpk) 
    351350         ! 
    352351         !                          !* AFTER fields : (last level for OPA, 3D required for AGRIF initialisation) 
    353          e3t_a(:,:,:) = e3t_n(:,:,:)   ;   e3u_a(:,:,:) = e3u_n(:,:,:)   ;   e3v_a(:,:,:) = e3v_n(:,:,:) 
     352         e3t_a(:,:,:) = e3t_n(:,:,:)   ;   e3u_a(:,:,:) = e3u_n(:,:,:) 
     353         e3w_a(:,:,:) = e3w_n(:,:,:)   ;   e3v_a(:,:,:) = e3v_n(:,:,:) 
    354354         ! 
    355355      ENDIF 
     
    370370                  END DO 
    371371                  nk_rnf(ji,jj) = jk 
    372                ELSEIF( h_rnf(ji,jj) == -1._wp   ) THEN   ;  nk_rnf(ji,jj) = 1 
    373                ELSEIF( h_rnf(ji,jj) == -999._wp ) THEN   ;  nk_rnf(ji,jj) = mbkt(ji,jj) 
     372               ELSEIF( h_rnf(ji,jj) == -1._wp   ) THEN   ;   nk_rnf(ji,jj) = 1 
     373               ELSEIF( h_rnf(ji,jj) == -999._wp ) THEN   ;   nk_rnf(ji,jj) = mbkt(ji,jj) 
    374374               ELSE 
    375375                  CALL ctl_stop( 'sbc_rnf_init: runoff depth not positive, and not -999 or -1, rnf value in file fort.999'  ) 
Note: See TracChangeset for help on using the changeset viewer.