Changeset 6661 for branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90
- Timestamp:
- 2016-06-03T10:51:04+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90
r6630 r6661 27 27 28 28 PUBLIC zdf_mxl ! called by step.F90 29 PUBLIC zdf_mxl_tref ! called by asminc.F90 29 30 PUBLIC zdf_mxl_alloc ! Used in zdf_tke_init 30 31 … … 86 87 REAL(wp) :: zN2_c ! local scalar 87 88 INTEGER, POINTER, DIMENSION(:,:) :: imld ! 2D workspace 88 REAL(wp) :: t_ref ! Reference temperature89 REAL(wp) :: temp_c = 0.2 ! temperature criterion for mixed layer depth90 89 !!---------------------------------------------------------------------- 91 90 ! … … 141 140 CALL iom_put( "mldkz5" , hmld ) ! turbocline depth 142 141 ENDIF 142 143 IF(ln_ctl) CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmlp, clinfo2=' hmlp : ', ovlap=1 ) 144 ! 145 CALL wrk_dealloc( jpi,jpj, imld ) 146 ! 147 IF( nn_timing == 1 ) CALL timing_stop('zdf_mxl') 148 ! 149 END SUBROUTINE zdf_mxl 143 150 151 152 SUBROUTINE zdf_mxl_tref() 153 !!---------------------------------------------------------------------- 154 !! *** ROUTINE zdf_mxl_tref *** 155 !! 156 !! ** Purpose : Compute the mixed layer depth with temperature criteria. 157 !! 158 !! ** Method : The temperature-defined mixed layer depth is required 159 !! when assimilating SST in a 2D analysis. 160 !! 161 !! ** Action : hmld_tref 162 !!---------------------------------------------------------------------- 163 ! 164 INTEGER :: ji, jj, jk ! dummy loop indices 165 REAL(wp) :: t_ref ! Reference temperature 166 REAL(wp) :: temp_c = 0.2 ! temperature criterion for mixed layer depth 167 !!---------------------------------------------------------------------- 168 ! 169 ! Initialise array 170 IF( zdf_mxl_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'zdf_mxl_tref : unable to allocate arrays' ) 171 144 172 !For the AMM model assimiation uses a temperature based mixed layer depth 145 173 !This is defined here … … 162 190 ENDDO 163 191 ENDDO 164 165 IF(ln_ctl) CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmlp, clinfo2=' hmlp : ', ovlap=1 ) 166 ! 167 CALL wrk_dealloc( jpi,jpj, imld ) 168 ! 169 IF( nn_timing == 1 ) CALL timing_stop('zdf_mxl') 170 ! 171 END SUBROUTINE zdf_mxl 192 193 END SUBROUTINE zdf_mxl_tref 172 194 173 195 !!======================================================================
Note: See TracChangeset
for help on using the changeset viewer.