- Timestamp:
- 2017-07-04T17:46:48+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/wrk_OMP_test_for_Silvia/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfddm.F90
r8055 r8279 95 95 !!gm and many acces in memory 96 96 97 DO jj = k_Jstr, k_Jend!== R=zrau = (alpha / beta) (dk[t] / dk[s]) ==!98 DO ji = k_Jstr, k_Iend97 DO jj = tnldj, tnlej !== R=zrau = (alpha / beta) (dk[t] / dk[s]) ==! 98 DO ji = tnldi, tnlei 99 99 zrw = ( gdepw_n(ji,jj,jk ) - gdept_n(ji,jj,jk) ) & 100 100 !!gm please, use e3w_n below … … 113 113 END DO 114 114 115 DO jj = k_Jstr, k_Jend!== indicators ==!116 DO ji = k_Jstr, k_Iend115 DO jj = tnldj, tnlej !== indicators ==! 116 DO ji = tnldi, tnlei 117 117 ! stability indicator: msks=1 if rn2>0; 0 elsewhere 118 118 IF( rn2(ji,jj,jk) + 1.e-12 <= 0. ) THEN ; zmsks(ji,jj) = 0._wp … … 145 145 ! ------------------ 146 146 ! Constant eddy coefficient: reset to the background value 147 DO jj = k_Jstr, k_Jend148 DO ji = k_Jstr, k_Iend147 DO jj = tnldj, tnlej 148 DO ji = tnldi, tnlei 149 149 zinr = 1._wp / zrau(ji,jj) 150 150 ! salt fingering
Note: See TracChangeset
for help on using the changeset viewer.