Changeset 14223 for NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90
- Timestamp:
- 2020-12-19T11:22:45+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90
r13472 r14223 99 99 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace 100 100 INTEGER :: ji, jj 101 INTEGER :: igloi, igloj ! to be removed in the future, see comment bellow102 101 !!--------------------------------------------------------------------- 103 102 #if defined key_si3 … … 106 105 ! define unique value on each point. z2d ranging from 0.05 to -0.05 107 106 ! 108 ! WARNING: to keep compatibility with the trunk that was including periodocity into the input data,109 ! we must define z2d as bellow.110 ! Once we decide to forget trunk compatibility, we must simply define z2d as:111 !!$ DO_2D( 0, 0, 0, 0 )112 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) )113 !!$ END_2D114 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) )115 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) )116 107 DO_2D( 0, 0, 0, 0 ) 117 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0 _oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) )108 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 118 109 END_2D 119 110 utau_ice(:,:) = 0.1_wp + z2d(:,:)
Note: See TracChangeset
for help on using the changeset viewer.