Changeset 7753 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90
r7698 r7753 71 71 ! Constrain the thickness of the smallest category above himin 72 72 !---------------------------------------------------------------------- 73 !$OMP PARALLEL74 !$OMP DO schedule(static) private(jj,ji,rswitch)75 73 DO jj = 1, jpj 76 74 DO ji = 1, jpi … … 87 85 ! ice concentration should not exceed amax 88 86 !----------------------------------------------------- 89 !$OMP DO schedule(static) private(jj, ji) 90 DO jj = 1, jpj 91 DO ji = 1, jpi 92 at_i(ji,jj) = 0._wp 93 END DO 94 END DO 87 at_i(:,:) = 0._wp 95 88 DO jl = 1, jpl 96 !$OMP DO schedule(static) private(jj, ji) 97 DO jj = 1, jpj 98 DO ji = 1, jpi 99 at_i(ji,jj) = a_i(ji,jj,jl) + at_i(ji,jj) 100 END DO 101 END DO 89 at_i(:,:) = a_i(:,:,jl) + at_i(:,:) 102 90 END DO 103 91 104 92 DO jl = 1, jpl 105 !$OMP DO schedule(static) private(jj, ji)106 93 DO jj = 1, jpj 107 94 DO ji = 1, jpi … … 113 100 END DO 114 101 END DO 115 !$OMP END PARALLEL116 102 117 103 !--------------------- … … 120 106 IF ( nn_icesal == 2 ) THEN 121 107 DO jl = 1, jpl 122 !$OMP PARALLEL DO schedule(static) private(jj,ji,zsal,rswitch)123 108 DO jj = 1, jpj 124 109 DO ji = 1, jpi … … 149 134 ! Ice drift 150 135 !------------ 151 !$OMP PARALLEL DO schedule(static) private(jj, ji)152 136 DO jj = 2, jpjm1 153 137 DO ji = 2, jpim1 … … 164 148 CALL lbc_lnk( v_ice(:,:), 'V', -1. ) 165 149 !mask velocities 166 !$OMP PARALLEL 167 !$OMP DO schedule(static) private(jj, ji) 168 DO jj = 1, jpj 169 DO ji = 1, jpi 170 u_ice(ji,jj) = u_ice(ji,jj) * umask(ji,jj,1) 171 v_ice(ji,jj) = v_ice(ji,jj) * vmask(ji,jj,1) 172 END DO 173 END DO 150 u_ice(:,:) = u_ice(:,:) * umask(:,:,1) 151 v_ice(:,:) = v_ice(:,:) * vmask(:,:,1) 174 152 175 153 ! ------------------------------------------------- … … 177 155 ! ------------------------------------------------- 178 156 DO jl = 1, jpl 179 !$OMP DO schedule(static) private(jj, ji) 180 DO jj = 1, jpj 181 DO ji = 1, jpi 182 oa_i(ji,jj,jl) = oa_i(ji,jj,jl) + a_i(ji,jj,jl) * rdt_ice / rday ! ice natural aging 183 afx_thd(ji,jj) = afx_thd(ji,jj) + ( a_i(ji,jj,jl) - a_i_b(ji,jj,jl) ) * r1_rdtice 184 END DO 185 END DO 157 oa_i(:,:,jl) = oa_i(:,:,jl) + a_i(:,:,jl) * rdt_ice / rday ! ice natural aging 158 afx_thd(:,:) = afx_thd(:,:) + ( a_i(:,:,jl) - a_i_b(:,:,jl) ) * r1_rdtice 186 159 END DO 187 160 afx_tot = afx_thd + afx_dyn 188 161 189 !$OMP DO schedule(static) private(jj, ji)190 162 DO jj = 1, jpj 191 163 DO ji = 1, jpi … … 201 173 END DO 202 174 END DO 203 !$OMP END PARALLEL204 175 205 176 ! conservation test
Note: See TracChangeset
for help on using the changeset viewer.