- Timestamp:
- 2015-11-30T20:55:41+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl.F90
r5621 r5956 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
Note: See TracChangeset
for help on using the changeset viewer.