- Timestamp:
- 2011-11-04T10:08:18+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2787_MERCATOR2_tidalharm/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90
r3042 r3043 234 234 nhan = (nitend_han-nit000_han+1)/nstep_han ! Number of dumps used for analysis 235 235 236 NBINCO= 2*nb_ana236 ninco = 2*nb_ana 237 237 238 238 ksp = 0 … … 246 246 kun = kun + 1 247 247 ksp = ksp + 1 248 ISPARSE(ksp) = keq249 JSPARSE(ksp) = kun250 SPARSEVALUE(ksp)= &248 nisparse(ksp) = keq 249 njsparse(ksp) = kun 250 valuesparse(ksp)= & 251 251 +( MOD(jc,2) * ft(jh) * COS(ana_freq(jh)*ztime + vt(jh) + ut(jh)) & 252 252 +(1.-MOD(jc,2))* ft(jh) * SIN(ana_freq(jh)*ztime + vt(jh) + ut(jh))) … … 255 255 END DO 256 256 257 NBSPARSE=ksp257 nsparse=ksp 258 258 259 259 ! Elevation: … … 265 265 DO jc = 1,2 266 266 kun = kun + 1 267 TAB4(kun)=ana_temp(ji,jj,kun,1)267 tmp4(kun)=ana_temp(ji,jj,kun,1) 268 268 ENDDO 269 269 ENDDO … … 273 273 ! Fill output array 274 274 DO jh = 1, nb_ana 275 ana_amp(ji,jj,jh,1)= TAB7((jh-1)*2+1)276 ana_amp(ji,jj,jh,2)= TAB7((jh-1)*2+2)275 ana_amp(ji,jj,jh,1)=tmp7((jh-1)*2+1) 276 ana_amp(ji,jj,jh,2)=tmp7((jh-1)*2+2) 277 277 END DO 278 278 END DO … … 303 303 DO jc = 1,2 304 304 kun = kun + 1 305 TAB4(kun)=ana_temp(ji,jj,kun,2)305 tmp4(kun)=ana_temp(ji,jj,kun,2) 306 306 ENDDO 307 307 ENDDO … … 311 311 ! Fill output array 312 312 DO jh = 1, nb_ana 313 ana_amp(ji,jj,jh,1)= TAB7((jh-1)*2+1)314 ana_amp(ji,jj,jh,2)= TAB7((jh-1)*2+2)313 ana_amp(ji,jj,jh,1)=tmp7((jh-1)*2+1) 314 ana_amp(ji,jj,jh,2)=tmp7((jh-1)*2+2) 315 315 END DO 316 316 … … 337 337 DO jc = 1,2 338 338 kun = kun + 1 339 TAB4(kun)=ana_temp(ji,jj,kun,3)339 tmp4(kun)=ana_temp(ji,jj,kun,3) 340 340 ENDDO 341 341 ENDDO … … 345 345 ! Fill output array 346 346 DO jh = 1, nb_ana 347 ana_amp(ji,jj,jh,1)= TAB7((jh-1)*2+1)348 ana_amp(ji,jj,jh,2)= TAB7((jh-1)*2+2)347 ana_amp(ji,jj,jh,1)=tmp7((jh-1)*2+1) 348 ana_amp(ji,jj,jh,2)=tmp7((jh-1)*2+2) 349 349 END DO 350 350
Note: See TracChangeset
for help on using the changeset viewer.