- Timestamp:
- 2017-06-06T15:55:44+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r8093 r8143 61 61 USE diurnal_bulk ! diurnal warm layer 62 62 USE cool_skin ! Cool skin 63 USE wrk_nemo ! working array64 63 65 64 IMPLICIT NONE … … 183 182 DO jj = 2, jpjm1 184 183 DO ji = fs_2, fs_jpim1 ! vector opt. 185 !!gm old186 !!gm BUG missing x 0.5187 zztmpx = ( bfrua(ji ,jj) * un(ji ,jj,mbku(ji ,jj)) &188 & + bfrua(ji-1,jj) * un(ji-1,jj,mbku(ji-1,jj)) )189 zztmpy = ( bfrva(ji, jj) * vn(ji,jj ,mbkv(ji,jj )) &190 & + bfrva(ji,jj-1) * vn(ji,jj-1,mbkv(ji,jj-1)) )191 z2d(ji,jj) = rau0 * SQRT( zztmpx * zztmpx + zztmpy * zztmpy ) * tmask(ji,jj,1)192 !!gm193 184 zztmp2 = ( ( rCdU_bot(ji+1,jj)+rCdU_bot(ji ,jj) ) * un(ji ,jj,mbku(ji ,jj)) )**2 & 194 185 & + ( ( rCdU_bot(ji ,jj)+rCdU_bot(ji-1,jj) ) * un(ji-1,jj,mbku(ji-1,jj)) )**2 & … … 196 187 & + ( ( rCdU_bot(ji,jj )+rCdU_bot(ji,jj-1) ) * vn(ji,jj-1,mbkv(ji,jj-1)) )**2 197 188 z2d(ji,jj) = zztmp * SQRT( zztmp2 ) * tmask(ji,jj,1) 198 !!gm new end199 189 ! 200 END DO201 END DO190 END DO 191 END DO 202 192 CALL lbc_lnk( z2d, 'T', 1. ) 203 193 CALL iom_put( "taubot", z2d ) … … 449 439 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 450 440 ! 451 REAL(wp), POINTER, DIMENSION(:,:) :: zw2d ! 2D workspace452 REAL(wp), POINTER, DIMENSION(:,:,:) :: zw3d ! 3D workspace441 REAL(wp), DIMENSION(jpi,jpj) :: zw2d ! 2D workspace 442 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d ! 3D workspace 453 443 !!---------------------------------------------------------------------- 454 444 ! 455 445 IF( nn_timing == 1 ) CALL timing_start('dia_wri') 456 446 ! 457 CALL wrk_alloc( jpi,jpj , zw2d ) 458 IF( .NOT.ln_linssh ) CALL wrk_alloc( jpi,jpj,jpk , zw3d ) 459 ! 460 ! Output the initial state and forcings 461 IF( ninist == 1 ) THEN 447 IF( ninist == 1 ) THEN !== Output the initial state and forcings ==! 462 448 CALL dia_wri_state( 'output.init', kt ) 463 449 ninist = 0 … … 467 453 ! ----------------- 468 454 469 ! local variable for debugging 470 ll_print = .FALSE. 455 ll_print = .FALSE. ! local variable for debugging 471 456 ll_print = ll_print .AND. lwp 472 457 … … 891 876 ENDIF 892 877 ! 893 CALL wrk_dealloc( jpi , jpj , zw2d )894 IF( .NOT.ln_linssh ) CALL wrk_dealloc( jpi , jpj , jpk , zw3d )895 !896 878 IF( nn_timing == 1 ) CALL timing_stop('dia_wri') 897 879 !
Note: See TracChangeset
for help on using the changeset viewer.