Changeset 7753 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90
r7698 r7753 83 83 IF( l_trddyn ) THEN ! temporary save of ta and sa trends 84 84 CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv ) 85 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji) 86 DO jk = 1, jpk 87 DO jj = 1, jpj 88 DO ji = 1, jpi 89 ztrdu(ji,jj,jk) = ua(ji,jj,jk) 90 ztrdv(ji,jj,jk) = va(ji,jj,jk) 91 END DO 92 END DO 93 END DO 85 ztrdu(:,:,:) = ua(:,:,:) 86 ztrdv(:,:,:) = va(:,:,:) 94 87 ENDIF 95 88 ! … … 98 91 .OR. nn_ice_embd == 2 ) THEN ! embedded sea-ice 99 92 ! 100 !$OMP PARALLEL DO schedule(static) private(jj, ji)101 93 DO jj = 2, jpjm1 102 94 DO ji = fs_2, fs_jpim1 ! vector opt. … … 108 100 IF( ln_apr_dyn .AND. .NOT.ln_dynspg_ts ) THEN !== Atmospheric pressure gradient (added later in time-split case) ==! 109 101 zg_2 = grav * 0.5 110 !$OMP PARALLEL DO schedule(static) private(jj, ji)111 102 DO jj = 2, jpjm1 ! gradient of Patm using inverse barometer ssh 112 103 DO ji = fs_2, fs_jpim1 ! vector opt. … … 124 115 CALL upd_tide( kt ) ! update tide potential 125 116 ! 126 !$OMP PARALLEL DO schedule(static) private(jj, ji)127 117 DO jj = 2, jpjm1 ! add tide potential forcing 128 118 DO ji = fs_2, fs_jpim1 ! vector opt. … … 138 128 zintp = REAL( MOD( kt-1, nn_fsbc ) ) / REAL( nn_fsbc ) 139 129 zgrau0r = - grav * r1_rau0 140 !$OMP PARALLEL 141 !$OMP DO schedule(static) private(jj, ji) 142 DO jj = 1, jpj 143 DO ji = 1, jpi 144 zpice(ji,jj) = ( zintp * snwice_mass(ji,jj) + ( 1.- zintp ) * snwice_mass_b(ji,jj) ) * zgrau0r 145 END DO 146 END DO 147 !$OMP DO schedule(static) private(jj, ji) 130 zpice(:,:) = ( zintp * snwice_mass(:,:) + ( 1.- zintp ) * snwice_mass_b(:,:) ) * zgrau0r 148 131 DO jj = 2, jpjm1 149 132 DO ji = fs_2, fs_jpim1 ! vector opt. … … 152 135 END DO 153 136 END DO 154 !$OMP END PARALLEL155 137 ! 156 138 CALL wrk_dealloc( jpi,jpj, zpice ) 157 139 ENDIF 158 140 ! 159 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji)160 141 DO jk = 1, jpkm1 !== Add all terms to the general trend 161 142 DO jj = 2, jpjm1 … … 177 158 ! 178 159 IF( l_trddyn ) THEN ! save the surface pressure gradient trends for further diagnostics 179 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji) 180 DO jk = 1, jpk 181 DO jj = 1, jpj 182 DO ji = 1, jpi 183 ztrdu(ji,jj,jk) = ua(ji,jj,jk) - ztrdu(ji,jj,jk) 184 ztrdv(ji,jj,jk) = va(ji,jj,jk) - ztrdv(ji,jj,jk) 185 END DO 186 END DO 187 END DO 160 ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 161 ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 188 162 CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 189 163 CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv )
Note: See TracChangeset
for help on using the changeset viewer.