Changeset 13295 for NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90
r13286 r13295 69 69 ! we must define z2d as bellow. 70 70 ! Once we decide to forget trunk compatibility, we must simply define z2d as: 71 !!$ DO_2D _00_0071 !!$ DO_2D( 0, 0, 0, 0 ) 72 72 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 73 73 !!$ END_2D 74 74 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 75 75 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 76 DO_2D _00_0076 DO_2D( 0, 0, 0, 0 ) 77 77 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 78 78 END_2D … … 81 81 pssh(:,:) = z2d(:,:) ! +/- 0.05 m 82 82 ! 83 DO_3D _00_00(1, jpkm1 )83 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 84 84 zfact = REAL(jk-1,wp) / REAL(jpk-1,wp) ! 0 to 1 to add a basic stratification 85 85 ! temperature choosen to lead to ~50% ice at the beginning if rn_thres_sst = 0.5
Note: See TracChangeset
for help on using the changeset viewer.