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 7210 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90 – NEMO

Ignore:
Timestamp:
2016-11-08T17:05:22+01:00 (7 years ago)
Author:
cbricaud
Message:

commit modification in CRS branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90

    r7207 r7210  
    223223     CALL crs_dom_e3( e1v, e2v, zfse3v, p_sfc_2d_crs=e1v_crs  , cd_type='V', p_mask=vmask, p_e3_crs=e3v_0_crs, p_e3_max_crs=e3v_max_0_crs) 
    224224 
    225      WHERE(e3t_max_0_crs == 0._wp) e3t_max_0_crs=r_inf 
    226      WHERE(e3u_max_0_crs == 0._wp) e3u_max_0_crs=r_inf 
    227      WHERE(e3v_max_0_crs == 0._wp) e3v_max_0_crs=r_inf 
    228      WHERE(e3w_max_0_crs == 0._wp) e3w_max_0_crs=r_inf 
     225     DO jk = 1, jpk 
     226        DO ji = 1, jpi_crs 
     227           DO jj = 1, jpj_crs 
     228              IF( e3t_max_0_crs(ji,jj,jk) == 0._wp ) e3t_max_0_crs(ji,jj,jk) = e3t_1d(jk) 
     229              IF( e3w_max_0_crs(ji,jj,jk) == 0._wp ) e3w_max_0_crs(ji,jj,jk) = e3w_1d(jk) 
     230              IF( e3u_max_0_crs(ji,jj,jk) == 0._wp ) e3u_max_0_crs(ji,jj,jk) = e3t_1d(jk) 
     231              IF( e3v_max_0_crs(ji,jj,jk) == 0._wp ) e3v_max_0_crs(ji,jj,jk) = e3t_1d(jk) 
     232           ENDDO 
     233        ENDDO 
     234     ENDDO 
    229235 
    230236#if defined key_vvl 
     
    281287     CALL crs_dom_ope( gdept_0, 'MAX', 'T', tmask, gdept_0_crs, p_e3=zfse3t, psgn=1.0 )  
    282288     CALL crs_dom_ope( gdepw_0, 'MAX', 'W', tmask, gdepw_0_crs, p_e3=zfse3w, psgn=1.0 ) 
     289 
     290     DO jk = 1, jpk 
     291        DO ji = 1, jpi_crs 
     292           DO jj = 1, jpj_crs 
     293              IF( gdept_0_crs(ji,jj,jk) .LE. 0._wp ) gdept_0_crs(ji,jj,jk) = gdept_1d(jk) 
     294              IF( gdepw_0_crs(ji,jj,jk) .LE. 0._wp ) gdepw_0_crs(ji,jj,jk) = gdepw_1d(jk) 
     295           ENDDO 
     296        ENDDO 
     297     ENDDO 
     298 
    283299#if defined key_vvl 
    284300     gdept_n_crs(:,:,:) = gdept_0_crs(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.