- Timestamp:
- 2017-12-19T09:47:17+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90
r9019 r9125 24 24 USE iom ! I/O manager library 25 25 USE lib_mpp ! MPP library 26 USE wrk_nemo ! Memory allocation27 26 28 27 IMPLICIT NONE … … 70 69 ! 71 70 INTEGER :: jk ! dummy loop indices 72 REAL(wp), POINTER, DIMENSION(:,:) :: z2d ! 2D workspace73 REAL(wp), POINTER, DIMENSION(:,:,:):: zpe ! 3D workspace71 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: z2d ! 2D workspace 72 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zpe ! 3D workspace 74 73 !!---------------------------------------------------------------------- 75 74 ! 76 CALL wrk_alloc( jpi, jpj, jpk, zpe )77 75 zpe(:,:,:) = 0._wp 78 76 ! … … 96 94 CASE ( jptra_zad ) ; CALL iom_put( "petrd_zad", zpe ) ! vertical advection 97 95 IF( ln_linssh ) THEN ! cst volume : adv flux through z=0 surface 98 CALL wrk_alloc( jpi, jpj, z2d)96 ALLOCATE( z2d(jpi,jpj) ) 99 97 z2d(:,:) = wn(:,:,1) * ( & 100 98 & - ( rab_n(:,:,1,jp_tem) + rab_pe(:,:,1,jp_tem) ) * tsn(:,:,1,jp_tem) & … … 102 100 & ) / e3t_n(:,:,1) 103 101 CALL iom_put( "petrd_sad" , z2d ) 104 CALL wrk_dealloc( jpi, jpj,z2d )102 DEALLOCATE( z2d ) 105 103 ENDIF 106 104 CASE ( jptra_ldf ) ; CALL iom_put( "petrd_ldf" , zpe ) ! lateral diffusion … … 115 113 CASE ( jptra_atf ) ; CALL iom_put( "petrd_atf" , zpe ) ! asselin time filter (last trend) 116 114 !IF( ln_linssh ) THEN ! cst volume : ssh term (otherwise include in e3t variation) 117 ! CALL wrk_alloc( jpi, jpj, z2d)115 ! ALLOCATE( z2d(jpi,jpj) ) 118 116 ! z2d(:,:) = ( ssha(:,:) - sshb(:,:) ) & 119 117 ! & * ( dPE_dt(:,:,1) * tsn(:,:,1,jp_tem) & 120 118 ! & + dPE_ds(:,:,1) * tsn(:,:,1,jp_sal) ) / ( e3t_n(:,:,1) * pdt ) 121 119 ! CALL iom_put( "petrd_sad" , z2d ) 122 ! CALL wrk_dealloc( jpi, jpj,z2d )120 ! DEALLOCATE( z2d ) 123 121 !ENDIF 124 122 ! 125 123 END SELECT 126 124 ! 127 CALL wrk_dealloc( jpi, jpj, jpk, zpe )128 125 ! 129 126 END SUBROUTINE trd_pen
Note: See TracChangeset
for help on using the changeset viewer.