- Timestamp:
- 2015-11-30T12:48:01+01:00 (9 years ago)
- Location:
- branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD
- Files:
-
- 8 edited
-
trd_oce.F90 (modified) (1 diff)
-
trddyn.F90 (modified) (1 diff)
-
trdglo.F90 (modified) (3 diffs)
-
trdken.F90 (modified) (6 diffs)
-
trdmxl.F90 (modified) (7 diffs)
-
trdpen.F90 (modified) (2 diffs)
-
trdtra.F90 (modified) (1 diff)
-
trdvor.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trd_oce.F90
r5950 r5951 71 71 INTEGER, PUBLIC, PARAMETER :: jpdyn_bfri = 12 !: implicit bottom friction (ln_bfrimp=.TRUE.) 72 72 INTEGER, PUBLIC, PARAMETER :: jpdyn_ken = 13 !: use for calculation of KE 73 INTEGER, PUBLIC, PARAMETER :: jpdyn_spgflt = 14 !: filter contribution to surface pressure gradient (spg_flt)74 INTEGER, PUBLIC, PARAMETER :: jpdyn_spgexp = 15 !: explicit contribution to surface pressure gradient (spg_flt)75 73 ! 76 74 !!---------------------------------------------------------------------- -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trddyn.F90
r5950 r5951 113 113 CASE( jpdyn_spg ) ; CALL iom_put( "utrd_spg", putrd ) ! surface pressure gradient 114 114 CALL iom_put( "vtrd_spg", pvtrd ) 115 CASE( jpdyn_spgexp ); CALL iom_put( "utrd_spgexp", putrd ) ! surface pressure gradient (explicit)116 CALL iom_put( "vtrd_spgexp", pvtrd )117 CASE( jpdyn_spgflt ); CALL iom_put( "utrd_spgflt", putrd ) ! surface pressure gradient (filtered)118 CALL iom_put( "vtrd_spgflt", pvtrd )119 115 CASE( jpdyn_pvo ) ; CALL iom_put( "utrd_pvo", putrd ) ! planetary vorticity 120 116 CALL iom_put( "vtrd_pvo", pvtrd ) -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdglo.F90
r5950 r5951 4 4 !! Ocean diagnostics: global domain averaged tracer and momentum trends 5 5 !!===================================================================== 6 !! History : 1.0 ! 2004-08 (C. Talandier) New trends organization7 !! 3.5 ! 2012-02 (G. Madec) add 3D tracer zdf trend output using iom6 !! History : 1.0 ! 2004-08 (C. Talandier) New trends organization 7 !! 3.5 ! 2012-02 (G. Madec) add 3D tracer zdf trend output using iom 8 8 !!---------------------------------------------------------------------- 9 9 … … 19 19 USE trd_oce ! trends: ocean variables 20 20 USE phycst ! physical constants 21 USE ldftra _oce ! ocean active tracers: lateral physics22 USE ldfdyn _oce! ocean dynamics: lateral physics21 USE ldftra ! lateral diffusion: eddy diffusivity & EIV coeff. 22 USE ldfdyn ! ocean dynamics: lateral physics 23 23 USE zdf_oce ! ocean vertical physics 24 24 USE zdfbfr ! bottom friction … … 26 26 USE eosbn2 ! equation of state 27 27 USE phycst ! physical constants 28 ! 28 29 USE lib_mpp ! distibuted memory computing library 29 30 USE in_out_manager ! I/O manager -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90
r5950 r5951 13 13 USE oce ! ocean dynamics and tracers variables 14 14 USE dom_oce ! ocean space and time domain variables 15 USE sbc_oce ! surface boundary condition: ocean 15 16 USE zdf_oce ! ocean vertical physics variables 16 17 USE trd_oce ! trends: ocean variables 17 18 !!gm USE dynhpg ! hydrostatic pressure gradient 18 19 USE zdfbfr ! bottom friction 19 USE ldftra_oce ! ocean active tracers lateral physics 20 USE sbc_oce ! surface boundary condition: ocean 20 USE ldftra ! ocean active tracers lateral physics 21 21 USE phycst ! physical constants 22 22 USE trdvor ! ocean vorticity trends 23 23 USE trdglo ! trends:global domain averaged 24 USE trdmxl ! ocean active mixed layer tracers trends 24 USE trdmxl ! ocean active mixed layer tracers trends 25 ! 25 26 USE in_out_manager ! I/O manager 26 27 USE iom ! I/O manager library … … 96 97 nkstp = kt 97 98 DO jk = 1, jpkm1 98 bu (:,:,jk) = e1u(:,:) *e2u(:,:) * fse3u_n(:,:,jk)99 bv (:,:,jk) = e1v(:,:) *e2v(:,:) * fse3v_n(:,:,jk)99 bu (:,:,jk) = e1e2u(:,:) * fse3u_n(:,:,jk) 100 bv (:,:,jk) = e1e2v(:,:) * fse3v_n(:,:,jk) 100 101 r1_bt(:,:,jk) = 1._wp / ( e1e2t(:,:) * fse3t_n(:,:,jk) ) * tmask(:,:,jk) 101 102 END DO … … 117 118 ! 118 119 SELECT CASE( ktrd ) 119 CASE( jpdyn_hpg ) ; CALL iom_put( "ketrd_hpg", zke ) ! hydrostatic pressure gradient 120 CASE( jpdyn_spg ) ; CALL iom_put( "ketrd_spg", zke ) ! surface pressure gradient 121 CASE( jpdyn_spgexp ); CALL iom_put( "ketrd_spgexp", zke ) ! surface pressure gradient (explicit) 122 CASE( jpdyn_spgflt ); CALL iom_put( "ketrd_spgflt", zke ) ! surface pressure gradient (filter) 123 CASE( jpdyn_pvo ) ; CALL iom_put( "ketrd_pvo", zke ) ! planetary vorticity 124 CASE( jpdyn_rvo ) ; CALL iom_put( "ketrd_rvo", zke ) ! relative vorticity (or metric term) 125 CASE( jpdyn_keg ) ; CALL iom_put( "ketrd_keg", zke ) ! Kinetic Energy gradient (or had) 126 CASE( jpdyn_zad ) ; CALL iom_put( "ketrd_zad", zke ) ! vertical advection 127 CASE( jpdyn_ldf ) ; CALL iom_put( "ketrd_ldf", zke ) ! lateral diffusion 128 CASE( jpdyn_zdf ) ; CALL iom_put( "ketrd_zdf", zke ) ! vertical diffusion 120 CASE( jpdyn_hpg ) ; CALL iom_put( "ketrd_hpg", zke ) ! hydrostatic pressure gradient 121 CASE( jpdyn_spg ) ; CALL iom_put( "ketrd_spg", zke ) ! surface pressure gradient 122 CASE( jpdyn_pvo ) ; CALL iom_put( "ketrd_pvo", zke ) ! planetary vorticity 123 CASE( jpdyn_rvo ) ; CALL iom_put( "ketrd_rvo", zke ) ! relative vorticity (or metric term) 124 CASE( jpdyn_keg ) ; CALL iom_put( "ketrd_keg", zke ) ! Kinetic Energy gradient (or had) 125 CASE( jpdyn_zad ) ; CALL iom_put( "ketrd_zad", zke ) ! vertical advection 126 CASE( jpdyn_ldf ) ; CALL iom_put( "ketrd_ldf", zke ) ! lateral diffusion 127 CASE( jpdyn_zdf ) ; CALL iom_put( "ketrd_zdf", zke ) ! vertical diffusion 129 128 ! ! wind stress trends 130 131 132 133 134 135 136 137 138 139 140 141 142 129 CALL wrk_alloc( jpi, jpj, z2dx, z2dy, zke2d ) 130 z2dx(:,:) = un(:,:,1) * ( utau_b(:,:) + utau(:,:) ) * e1u(:,:) * e2u(:,:) * umask(:,:,1) 131 z2dy(:,:) = vn(:,:,1) * ( vtau_b(:,:) + vtau(:,:) ) * e1v(:,:) * e2v(:,:) * vmask(:,:,1) 132 zke2d(1,:) = 0._wp ; zke2d(:,1) = 0._wp 133 DO jj = 2, jpj 134 DO ji = 2, jpi 135 zke2d(ji,jj) = 0.5_wp * ( z2dx(ji,jj) + z2dx(ji-1,jj) & 136 & + z2dy(ji,jj) + z2dy(ji,jj-1) ) * r1_bt(ji,jj,1) 137 END DO 138 END DO 139 CALL iom_put( "ketrd_tau", zke2d ) 140 CALL wrk_dealloc( jpi, jpj , z2dx, z2dy, zke2d ) 141 CASE( jpdyn_bfr ) ; CALL iom_put( "ketrd_bfr", zke ) ! bottom friction (explicit case) 143 142 !!gm TO BE DONE properly 144 143 !!gm only valid if ln_bfrimp=F otherwise the bottom stress as to be recomputed at the end of the computation.... … … 162 161 ! ENDIF 163 162 !!gm end 164 163 CASE( jpdyn_atf ) ; CALL iom_put( "ketrd_atf", zke ) ! asselin filter trends 165 164 !! a faire !!!! idee changer dynnxt pour avoir un appel a jpdyn_bfr avant le swap !!! 166 165 !! reflechir a une possible sauvegarde du "vrai" un,vn pour le calcul de atf.... … … 184 183 ! CALL iom_put( "ketrd_bfri", zke2d ) 185 184 ! ENDIF 186 187 188 189 190 191 192 193 185 CASE( jpdyn_ken ) ; ! kinetic energy 186 ! called in dynnxt.F90 before asselin time filter 187 ! with putrd=ua and pvtrd=va 188 zke(:,:,:) = 0.5_wp * zke(:,:,:) 189 CALL iom_put( "KE", zke ) 190 ! 191 CALL ken_p2k( kt , zke ) 192 CALL iom_put( "ketrd_convP2K", zke ) ! conversion -rau*g*w 194 193 ! 195 194 END SELECT … … 263 262 ENDIF 264 263 ! ! allocate box volume arrays 265 IF 264 IF( trd_ken_alloc() /= 0 ) CALL ctl_stop('trd_ken_alloc: failed to allocate arrays') 266 265 ! 267 266 !!gm IF( .NOT. (ln_hpg_zco.OR.ln_hpg_zps) ) & 268 267 !!gm & CALL ctl_stop('trd_ken_init : only full and partial cells are coded for conversion rate') 269 268 ! 270 IF ( .NOT.lk_vvl ) THEN! constant volume: bu, bv, 1/bt computed one for all269 IF( .NOT.lk_vvl ) THEN ! constant volume: bu, bv, 1/bt computed one for all 271 270 DO jk = 1, jpkm1 272 bu (:,:,jk) = e1u(:,:) *e2u(:,:) * fse3u_n(:,:,jk)273 bv (:,:,jk) = e1v(:,:) *e2v(:,:) * fse3v_n(:,:,jk)271 bu (:,:,jk) = e1e2u(:,:) * fse3u_n(:,:,jk) 272 bv (:,:,jk) = e1e2v(:,:) * fse3v_n(:,:,jk) 274 273 r1_bt(:,:,jk) = 1._wp / ( e1e2t(:,:) * fse3t_n(:,:,jk) ) 275 274 END DO -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl.F90
r5950 r5951 22 22 USE trd_oce ! trends: ocean variables 23 23 USE trdmxl_oce ! ocean variables trends 24 USE ldftra _oce ! ocean active tracers lateral physics24 USE ldftra ! lateral diffusion: eddy diffusivity & EIV coeff. 25 25 USE zdf_oce ! ocean vertical physics 26 USE in_out_manager ! I/O manager27 26 USE phycst ! Define parameters for the routines 28 27 USE dianam ! build the name of file (routine) … … 30 29 USE zdfmxl ! mixed layer depth 31 30 USE zdfddm ! ocean vertical physics: double diffusion 32 USE ioipsl ! NetCDF library33 31 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 34 32 USE diadimg ! dimg direct access file format output 35 33 USE trdmxl_rst ! restart for diagnosing the ML trends 34 ! 35 USE in_out_manager ! I/O manager 36 USE ioipsl ! NetCDF library 36 37 USE prtctl ! Print control 37 38 USE restart ! for lrst_oce … … 73 74 !! * Substitutions 74 75 # include "domzgr_substitute.h90" 75 # include "ldftra_substitute.h90"76 76 # include "zdfddm_substitute.h90" 77 77 !!---------------------------------------------------------------------- … … 165 165 166 166 167 168 167 SELECT CASE( ktrd ) 168 CASE( jptra_npc ) ! non-penetrative convection: regrouped with zdf 169 169 !!gm : to be completed ! 170 170 ! IF( .... 171 171 !!gm end 172 173 172 CASE( jptra_zdfp ) ! iso-neutral diffusion: "pure" vertical diffusion 173 ! ! regroup iso-neutral diffusion in one term 174 174 tmltrd(:,:,jpmxl_ldf) = tmltrd(:,:,jpmxl_ldf) + ( tmltrd(:,:,jpmxl_zdf) - tmltrd(:,:,jpmxl_zdfp) ) 175 175 smltrd(:,:,jpmxl_ldf) = smltrd(:,:,jpmxl_ldf) + ( smltrd(:,:,jpmxl_zdf) - smltrd(:,:,jpmxl_zdfp) ) … … 801 801 END IF 802 802 803 IF( nn_cla == 1 ) CALL ctl_warn( ' You set n_cla = 1. Note that the Mixed-Layer diagnostics ', &804 & ' are not exact along the corresponding straits. ')805 806 803 ! ! allocate trdmxl arrays 807 804 IF( trd_mxl_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'trd_mxl_init : unable to allocate trdmxl arrays' ) … … 810 807 811 808 812 813 nkstp = nit000 - 1 ! current time step indicator initialization 809 nkstp = nit000 - 1 ! current time step indicator initialization 814 810 815 811 … … 851 847 IF( nn_ctls == 1 ) THEN 852 848 CALL ctl_opn( inum, 'ctlsurf_idx', 'OLD', 'UNFORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 853 READ ( inum ) nbol849 READ ( inum, * ) nbol 854 850 CLOSE( inum ) 855 851 END IF -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90
r5950 r5951 17 17 USE trd_oce ! trends: ocean variables 18 18 USE eosbn2 ! equation of state and related derivatives 19 USE ldftra _oce ! ocean active tracers lateral physics19 USE ldftra ! lateral diffusion: eddy diffusivity & EIV coeff. 20 20 USE zdfddm ! vertical physics: double diffusion 21 21 USE phycst ! physical constants 22 ! 22 23 USE in_out_manager ! I/O manager 23 24 USE iom ! I/O manager library … … 99 100 CALL wrk_alloc( jpi, jpj, z2d ) 100 101 z2d(:,:) = wn(:,:,1) * ( & 101 102 103 &) / fse3t(:,:,1)102 & - ( rab_n(:,:,1,jp_tem) + rab_pe(:,:,1,jp_tem) ) * tsn(:,:,1,jp_tem) & 103 & + ( rab_n(:,:,1,jp_sal) + rab_pe(:,:,1,jp_sal) ) * tsn(:,:,1,jp_sal) & 104 & ) / fse3t(:,:,1) 104 105 CALL iom_put( "petrd_sad" , z2d ) 105 106 CALL wrk_dealloc( jpi, jpj, z2d ) -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90
r5950 r5951 23 23 USE trdpen ! trends: Potential ENergy 24 24 USE trdmxl ! ocean active mixed layer tracers trends 25 USE ldftra_oce ! ocean active tracers lateral physics 25 USE ldftra ! ocean active tracers lateral physics 26 USE ldfslp 26 27 USE zdfddm ! vertical physics: double diffusion 27 28 USE phycst ! physical constants 29 ! 28 30 USE in_out_manager ! I/O manager 29 31 USE iom ! I/O manager library -
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90
r5950 r5951 20 20 USE sbc_oce ! surface boundary condition: ocean 21 21 USE phycst ! Define parameters for the routines 22 USE ldfdyn _oce! ocean active tracers: lateral physics22 USE ldfdyn ! ocean active tracers: lateral physics 23 23 USE dianam ! build the name of file (routine) 24 24 USE zdfmxl ! mixed layer depth 25 USE lbclnk ! ocean lateral boundary conditions (or mpp link)25 ! 26 26 USE in_out_manager ! I/O manager 27 27 USE ioipsl ! NetCDF library 28 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 28 29 USE lib_mpp ! MPP library 29 30 USE wrk_nemo ! Memory allocation … … 57 58 !! * Substitutions 58 59 # include "domzgr_substitute.h90" 59 # include "ldfdyn_substitute.h90"60 60 # include "vectopt_loop_substitute.h90" 61 61 !!----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.
