- Timestamp:
- 2020-08-20T12:10:55+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC
- Files:
-
- 2 added
- 3 edited
- 5 copied
-
diawri.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/diawri.F90) (3 diffs)
-
divhor.F90 (added)
-
domain.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/domain.F90) (1 diff)
-
dommsk.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/dommsk.F90) (5 diffs)
-
dynldf_lap_blp.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/dynldf_lap_blp.F90) (7 diffs)
-
oce.F90 (added)
-
sshwzv.F90 (copied) (copied from NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/OCE/DYN/sshwzv.F90) (2 diffs)
-
usrdef_nam.F90 (modified) (5 diffs)
-
usrdef_sbc.F90 (modified) (1 diff)
-
usrdef_zgr.F90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/diawri.F90
r13151 r13416 75 75 INTEGER :: nid_U, nz_U, nh_U, ndim_U, ndim_hU ! grid_U file 76 76 INTEGER :: nid_V, nz_V, nh_V, ndim_V, ndim_hV ! grid_V file 77 INTEGER :: nid_F, nz_F, nh_F, ndim_F, ndim_hF ! grid_F file 77 78 INTEGER :: nid_W, nz_W, nh_W ! grid_W file 78 79 INTEGER :: nid_A, nz_A, nh_A, ndim_A, ndim_hA ! grid_ABL file 79 80 INTEGER :: ndex(1) ! ??? 80 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hT, ndex_hU, ndex_hV 81 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hT, ndex_hU, ndex_hV, ndex_hF 81 82 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hA, ndex_A ! ABL 82 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_T, ndex_U, ndex_V 83 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_T, ndex_U, ndex_V, ndex_F 83 84 INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_bT 84 85 … … 138 139 CALL iom_put("e3u_0", e3u_0(:,:,:) ) 139 140 CALL iom_put("e3v_0", e3v_0(:,:,:) ) 141 CALL iom_put("e3f_0", e3f_0(:,:,:) ) 140 142 ! 141 143 !!st13 … … 173 175 CALL iom_put( "ssh" , ssh(:,:,Kmm) ) ! sea surface height 174 176 ENDIF 175 176 !!an 177 IF( iom_use("ht") ) & ! water column at t-point 177 ! 178 # if defined key_bvp 179 IF( iom_use("rpo") ) CALL iom_put( "rpo", rpo(:,:,1) ) 180 IF( iom_use("bmpt") ) CALL iom_put( "bmpt", bmpt(:,:,1) ) 181 IF( iom_use("bmpu") ) CALL iom_put( "bmpu", bmpu(:,:,1) ) 182 IF( iom_use("bmpv") ) CALL iom_put( "bmpv", bmpv(:,:,1) ) 183 #endif 184 ! 185 IF( iom_use("ht") ) & ! water column at t-point 178 186 CALL iom_put( "ht" , ht_0(:,:) + ssh(:,:,Kmm) ) 179 187 ! -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/domain.F90
r13151 r13416 154 154 CALL dom_zgr( ik_top, ik_bot ) ! Vertical mesh and bathymetry (return top and bottom ocean t-level indices) 155 155 156 CALL dom_msk( ik_top, ik_bot ) ! Masks 156 CALL dom_msk( ik_top, ik_bot ) ! Masks 157 ! Penalisation fields 158 ! 159 # if defined key_bvp 160 ! Multiply the scale factors by the porosity field 161 ! ------------------------------------------------- 162 e3t_0(:,:,:) = rpo(:,:,:) * e3t_0(:,:,:) 163 ! phiUe3U = mean_i(phiTe3T) 164 DO_3D_00_00(1,jpkm1) 165 e3u_0(ji,jj,jk) = 0.5_wp * ( e3t_0(ji,jj,jk) + e3t_0(ji+1,jj,jk) ) 166 e3v_0(ji,jj,jk) = 0.5_wp * ( e3t_0(ji,jj,jk) + e3t_0(ji,jj+1,jk) ) 167 e3f_0(ji,jj,jk) = 0.25_wp * ( e3t_0(ji,jj+1,jk) + e3t_0(ji+1,jj+1,jk) & 168 & + e3t_0(ji,jj ,jk) + e3t_0(ji+1,jj ,jk) ) 169 END_3D 170 ! copy the communication layer (kinf od periodic over the basin) 171 CALL lbc_lnk_multi( 'domain', e3t_0, 'T', 1., e3u_0, 'U', 1., & 172 & e3v_0, 'V', 1., e3f_0, 'F', 1., kfillmode=jpfillcopy ) 173 ! 174 ! why not (ca se fait une fois, ca ne coute pas si cher) 175 ! ou, WHERE( tmask(:,:,:) /= 0._wp ) e3t_0(:,:,:) = pry(:,:,:) * e3t_0(:,:,:) 176 ! phiU = mean_i(phiT) 177 !DO_3D_00_00 (1,jpkm1) 178 ! e3u_0(ji,jj,jk) = 0.5_wp * ( pry(ji,jj,jk) + pry(ji+1,jj,jk) ) * e3u_0(ji,jj,jk) 179 ! e3v_0(ji,jj,jk) = 0.5_wp * ( pry(ji,jj,jk) + pry(ji,jj+1,jk) ) * e3v_0(ji,jj,jk) 180 ! e3f_0(ji,jj,jk) = 0.25_wp * ( pry(ji,jj+1,jk) + pry(ji+1,jj+1,jk) & 181 ! & pry(ji,jj ,jk) + pry(ji+1,jj ,jk) ) * e3f_0(ji,jj,jk) 182 !END_3D 183 !CALL lbc_lnk_multi( 'domain', e3u_0, 'U', 1., e3v_0, 'V', 1., e3f_0, 'F', 1., kfillmode=jfillcopy ) 184 #endif 157 185 ! 158 186 ht_0(:,:) = 0._wp ! Reference ocean thickness -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/dommsk.F90
r13151 r13416 32 32 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 33 33 USE lib_mpp ! Massively Parallel Processing library 34 ! 35 USE usrdef_nam , ONLY : rn_abp, r1_abp, rn_fsp 34 36 35 37 IMPLICIT NONE … … 92 94 INTEGER :: iktop, ikbot ! - - 93 95 INTEGER :: ios, inum 94 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zwf ! 2D workspace 96 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zwf ! 2D workspace 97 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3 ! 3D workspace 95 98 !! 96 99 NAMELIST/namlbc/ rn_shlat, ln_vorlat … … 139 142 END_2D 140 143 ! 144 # if defined key_bvp 145 ! Add a land layer around the domain of porosity alpha 146 ! ----------------------------------------------------- 147 !!an to be added : lissage sur nn_cnp dx 148 WRITE(numout,*) 'dommsk rpo entering' 149 !zw3 (:,:,:) = tmask(:,:,:) 150 zw3 = tmask 151 rpo (:,:,:) = 1._wp 152 r1_rpo(:,:,:) = 1._wp 153 bmpt (:,:,:) = 0._wp 154 ! alternative pour hf dans les angles 155 WHERE(tmask == 0._wp) 156 rpo (:,:,:) = rn_abp ! used on T points 157 r1_rpo(:,:,:) = r1_abp ! 158 END WHERE 159 WRITE(numout,*) 'dommsk rpo loop' 160 ! futur : boucle 3D et test 2D sur les tmask(ji,jj,jk) pour choper la bathy 161 DO_2D_00_00 162 IF( tmask(ji,jj,1) < tmask(ji+1,jj+1,1) .OR. & 163 & tmask(ji,jj,1) < tmask(ji+1,jj-1,1) .OR. & 164 & tmask(ji,jj,1) < tmask(ji-1,jj-1,1) .OR. & 165 & tmask(ji,jj,1) < tmask(ji-1,jj+1,1) ) THEN 166 ! plus simple : tmask(ji) + tmask(ji+1)= 0.5 idem pour jj en .OR. 167 ! mais j'ai pas mon coin 168 ! pour smooth, boucle un peu similaire, si pénalisé alors son copain l'est aussi mais un peu moins 169 ! peut etre sur pry 170 ! to be added : pry = 0 dans la terre et r1_pry vaut 1 pour ne pas avoir de 1/0 dans certain cas 171 ! non, car e3t0 vaudra 0 alors 172 ! ! flat bottom 173 zw3 (ji,jj, 1:jpkm1) = 1._wp 174 !rpo (ji,jj, 1:jpkm1) = rn_abp ! used on T points 175 !r1_rpo(ji,jj, 1:jpkm1) = r1_abp ! 176 bmpt (ji,jj, 1:jpkm1) = rn_fsp ! 177 ENDIF 178 END_2D 179 tmask = zw3 180 CALL lbc_lnk_multi( 'dommsk', rpo, 'T', 1._wp, r1_rpo , 'T', 1._wp, & 181 & bmpt, 'T', 1._wp, kfillmode=jpfillcopy ) 182 WRITE(numout,*) 'dommsk rpo done' 183 #endif 184 ! 141 185 ! the following call is mandatory 142 186 ! it masks boundaries (bathy=0) where needed depending on the configuration (closed, periodic...) … … 157 201 END_3D 158 202 ENDIF 159 203 ! 160 204 ! Ocean/land mask at u-, v-, and f-points (computed from tmask) 161 205 ! ---------------------------------------- … … 174 218 END DO 175 219 CALL lbc_lnk_multi( 'dommsk', umask, 'U', 1., vmask, 'V', 1., fmask, 'F', 1. ) ! Lateral boundary conditions 176 220 221 # if defined key_bvp 222 ! Permeability applied on the inner boundary layer 223 ! -------------------------------------------- 224 WRITE(numout,*) 'dommsk bmp loop' 225 bmpu(:,:,:) = 0._wp 226 bmpv(:,:,:) = 0._wp 227 DO_2D_00_00 228 IF( ( umask(ji-1,jj+1,1) + umask(ji+1,jj+1,1) & ! U points 229 & + umask(ji-1,jj-1,1) + umask(ji+1,jj-1,1) ) == 3._wp) THEN 230 ! 231 bmpu(ji,jj, 1:jpkm1) = rn_fsp 232 ! 233 ENDIF 234 IF( ( vmask(ji-1,jj+1,1) + vmask(ji+1,jj+1,1) & ! V points 235 & + vmask(ji-1,jj-1,1) + vmask(ji+1,jj-1,1) ) == 3._wp) THEN 236 ! 237 bmpv(ji,jj, 1:jpkm1) = rn_fsp 238 ! 239 ENDIF 240 END_2D 241 WRITE(numout,*) 'dommsk bmp done' 242 CALL lbc_lnk_multi( 'dommsk', bmpu, 'U', 1._wp, bmpv , 'V', 1._wp, kfillmode=jpfillcopy ) 243 #endif 244 177 245 ! Ocean/land mask at wu-, wv- and w points (computed from tmask) 178 246 !----------------------------------------- -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/dynldf_lap_blp.F90
r13151 r13416 20 20 USE in_out_manager ! I/O manager 21 21 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 22 22 USE lib_mpp ! MPP library 23 ! 24 USE usrdef_nam , ONLY : nn_dynldf_lap_typ ! use laplacian parameter 25 ! 23 26 IMPLICIT NONE 24 27 PRIVATE … … 31 34 INTEGER, PUBLIC, PARAMETER :: np_dynldf_lap_symN = 3 ! symmetric laplacian (cartesian) 32 35 33 INTEGER, PUBLIC, PARAMETER :: ln_dynldf_lap_typ = 1 ! choose type of laplacian (ideally from namelist)36 !INTEGER, PUBLIC, PARAMETER :: nn_dynldf_lap_typ = 1 ! choose type of laplacian (ideally from namelist) 34 37 !!anSYM 35 38 !! * Substitutions … … 70 73 REAL(wp), DIMENSION(jpi,jpj) :: zcur, zdiv 71 74 REAL(wp), DIMENSION(jpi,jpj) :: zten, zshe ! tension (diagonal) and shearing (anti-diagonal) terms 75 ! 76 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ze3t,ze3u,ze3v, ze3f ! not penalised scale factors 72 77 !!---------------------------------------------------------------------- 73 78 ! … … 79 84 WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator, pass=', kpass 80 85 WRITE(numout,*) '~~~~~~~ ' 81 WRITE(numout,*) ' ln_dynldf_lap_typ = ', ln_dynldf_lap_typ82 SELECT CASE( ln_dynldf_lap_typ ) ! print the choice of operator86 WRITE(numout,*) ' nn_dynldf_lap_typ = ', nn_dynldf_lap_typ 87 SELECT CASE( nn_dynldf_lap_typ ) ! print the choice of operator 83 88 CASE( np_dynldf_lap_rot ) ; WRITE(numout,*) ' ==>>> div-rot laplacian' 84 89 CASE( np_dynldf_lap_sym ) ; WRITE(numout,*) ' ==>>> symmetric laplacian (covariant form)' 85 CASE( np_dynldf_lap_symN) ; WRITE(numout,*) ' ==>>> symmetric laplacian ( simpleform)'90 CASE( np_dynldf_lap_symN) ; WRITE(numout,*) ' ==>>> symmetric laplacian (cartesian form)' 86 91 END SELECT 87 92 ENDIF … … 91 96 ENDIF 92 97 ! 93 SELECT CASE( ln_dynldf_lap_typ ) 98 !!an Kbb because grad et LeapFrog explicit 99 !! ze3t(:,:,:) = e3t(:,:,:,Kbb) 100 !! ze3u(:,:,:) = e3u(:,:,:,Kmm) 101 !! ze3v(:,:,:) = e3v(:,:,:,Kmm) 102 !! ze3f(:,:,:) = e3f(:,:,:) 103 !!# if defined key_bvp 104 !! ! gradient and divergence not penalised 105 !! ze3t(:,:,:) = ze3t(:,:,:) * r1_rpo(:,:,:) 106 !! ze3u(:,:,:) = ze3u(:,:,:) * r1_rpo(:,:,:) 107 !! ze3v(:,:,:) = ze3v(:,:,:) * r1_rpo(:,:,:) 108 !! !ze3f(:,:,:) = ze3f(:,:,:) * r1_rpo(:,:,:) ! unused 109 !!#endif 110 SELECT CASE( nn_dynldf_lap_typ ) 94 111 ! 95 112 CASE ( np_dynldf_lap_rot ) !== Vorticity-Divergence form ==! … … 99 116 DO_2D_01_01 100 117 ! ! ahm * e3 * curl (computed from 1 to jpim1/jpjm1) 101 !!gm open question here : e3f at before or now ? probably now...102 118 !!gm note that ahmf has already been multiplied by fmask 103 119 zcur(ji-1,jj-1) = & … … 193 209 ! 194 210 CASE DEFAULT ! error 195 CALL ctl_stop('STOP','dyn_ldf_lap: wrong value for ln_dynldf_lap_typ' )211 CALL ctl_stop('STOP','dyn_ldf_lap: wrong value for nn_dynldf_lap_typ' ) 196 212 END SELECT 197 213 ! -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/sshwzv.F90
r13151 r13416 77 77 REAL(wp), DIMENSION(jpi,jpj,jpt), INTENT(inout) :: pssh ! sea-surface height 78 78 ! 79 INTEGER :: j k ! dummy loop index79 INTEGER :: ji, jj, jk ! dummy loop indices 80 80 REAL(wp) :: zcoef ! local scalar 81 81 REAL(wp), DIMENSION(jpi,jpj) :: zhdiv ! 2D workspace 82 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ze3t ! penalised scale factors 82 83 !!---------------------------------------------------------------------- 83 84 ! … … 98 99 CALL wad_lmt(pssh(:,:,Kbb), zcoef * (emp_b(:,:) + emp(:,:)), rDt, Kmm, uu, vv ) 99 100 ENDIF 100 101 ! 102 DO_3D_00_00( 1, jpkm1 ) 103 ze3t(ji,jj,jk) = e3t(ji,jj,jk,Kmm) 104 END_3D 105 # if defined key_bvp 106 ! divergence not penalised 107 ze3t(:,:,:) = ze3t(:,:,:) * r1_rpo(:,:,:) 108 #endif 101 109 CALL div_hor( kt, Kbb, Kmm ) ! Horizontal divergence 102 110 ! 103 111 zhdiv(:,:) = 0._wp 104 112 DO jk = 1, jpkm1 ! Horizontal divergence of barotropic transports 105 zhdiv(:,:) = zhdiv(:,:) + e3t(:,:,jk,Kmm) * hdiv(:,:,jk)113 zhdiv(:,:) = zhdiv(:,:) + ze3t(:,:,jk) * hdiv(:,:,jk) 106 114 END DO 107 115 ! ! Sea surface elevation time stepping -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/usrdef_nam.F90
r13005 r13416 42 42 ! ! = 3 symmetric laplacian (cartesian) 43 43 LOGICAL , PUBLIC :: ln_dynldf_lap_PM ! if True - apply the P.Marchand boundary condition on the laplacian 44 ! !!* penalisation *!! 45 REAL(wp), PUBLIC :: rn_abp ! alpha boundary parameter [-] 46 INTEGER , PUBLIC :: nn_cnp ! number of cell on which is smoothed the porosity (phi) [-] 47 REAL(wp), PUBLIC :: rn_fsp ! friction parameter 1/epsilon of the permeability [1/s] 48 ! 49 REAL(wp), PUBLIC :: r1_abp ! inverse alpha boundary parameter [-] 44 50 ! 45 51 !!---------------------------------------------------------------------- … … 74 80 & rn_f0 ,rn_beta, & ! coriolis parameter 75 81 & rn_modified_grav, rn_rfr , rn_tau, & ! reduced gravity, friction, wind 76 & nn_dynldf_lap_typ, ln_dynldf_lap_PM ! temporary parameter 82 & nn_dynldf_lap_typ, ln_dynldf_lap_PM, & ! temporary parameter 83 & rn_abp, nn_cnp, rn_fsp ! penalisation parameters 77 84 !!---------------------------------------------------------------------- 78 85 ! … … 97 104 kpi = ceiling(zbl / ze1 ) ! Global Domain size + ghost cells 98 105 kpj = ceiling(zbl / ze1 ) ! Global Domain size + ghost cells 99 100 !!an101 !#if defined key_agrif102 ! IF( .NOT. Agrif_Root() ) THEN103 ! kpi = nbcellsx + 2 + 2*nbghostcells104 ! kpj = nbcellsy + 2 + 2*nbghostcells105 ! ENDIF106 !#endif107 !!an108 106 ! 109 107 IF( rn_modified_grav /= 0._wp) grav = rn_modified_grav ! update the gravity 110 !111 ! !== Temporary namelist parameter ==!112 !113 ! ll_dynldf_lap_PM = ln_dynldf_lap_PM114 108 ! 115 109 kpk = jpkglo … … 117 111 kperio = 0 ! AM98 configuration : closed domain 118 112 ! 113 # if defined key_bvp 114 r1_abp = 1._wp / rn_abp 115 #endif 119 116 ! ! control print 120 117 IF(lwp) THEN … … 131 128 WRITE(numout,*) ' rotation angle chosen rn_theta = ', rn_theta, 'deg' 132 129 WRITE(numout,*) ' modified gravity rn_modified_grav = ', rn_modified_grav, 'm2/s' 133 #if defined key_agrif134 IF( Agrif_Root() ) THEN135 #endif136 WRITE(numout,*) ' jpiglo = Nx*nn_AM98 + 2*n_ghost jpiglo = ', kpi137 WRITE(numout,*) ' jpjglo = Nx*nn_AM98 + 2*n_ghost jpjglo = ', kpj138 #if defined key_agrif139 ENDIF140 #endif141 130 WRITE(numout,*) ' number of model levels jpkglo = ', kpk 142 131 WRITE(numout,*) ' ' -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/usrdef_sbc.F90
r13005 r13416 88 88 ! length of the domain : 2000km x 2000km 89 89 utau(ji,jj) = - ztauu * COS( rpi * gphiu(ji,jj) / 2000000._wp) 90 vtau(ji,jj) = - ztauv * COS( rpi * gphi u(ji,jj) / 2000000._wp)90 vtau(ji,jj) = - ztauv * COS( rpi * gphiv(ji,jj) / 2000000._wp) 91 91 END DO 92 92 END DO -
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/cfgs/AM98/MY_SRC/usrdef_zgr.F90
r13005 r13416 141 141 pdepw_1d(1) = 0._wp 142 142 pdept_1d(1) = 250._wp 143 143 ! 144 144 pdepw_1d(2) = 500._wp 145 145 pdept_1d(2) = 750._wp 146 147 146 ! 148 147 ! ! e3t and e3w from depth … … 194 193 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) 195 194 ! 196 ! rotated domain (45deg)197 !zylim0 = 100000._wp * SQRT( 2._wp ) / 3 ! dy*sqrt(2)/2 * 2/3198 !zylim1 = 2000000._wp - zylim0 ! 2000km - dy*sqrt(2)/2 * 2/3199 !zxlim0 = 100000._wp * SQRT( 2._wp ) / 3 ! dx*sqrt(2)/2 * 2/3200 !zxlim1 = 2000000._wp - zxlim0 ! 2000km - dx*sqrt(2)/2 * 2/3201 202 ! rotated domain203 !zylim0 = 10000._wp / REAL( nn_AM98, wp) ! dy/10204 !zylim1 = 2000000._wp + zylim0 ! 2000km + dy/10205 !zxlim0 = 10000._wp / REAL( nn_AM98, wp) ! dx/10206 !zxlim1 = 2000000._wp + zxlim0 ! 2000km + dx/10207 208 195 ! 209 196 zylim0 = 10000._wp ! +10km … … 211 198 zxlim0 = 10000._wp ! +10km 212 199 zxlim1 = 2010000._wp ! 2010km 213 200 ! 214 201 DO jj = 1, jpj 215 202 DO ji = 1, jpi … … 229 216 END DO 230 217 ! mask the lonely corners 231 DO jj = 1, jpj232 DO ji = 1, jpi218 DO jj = 2, jpim1 219 DO ji = 2, jpim1 233 220 zcoeff = k_top(ji+1,jj) + k_top(ji,jj+1) & 234 221 + k_top(ji-1,jj) + k_top(ji,jj-1) … … 239 226 END DO 240 227 END DO 241 ! k_bot(:,:) = NINT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere 242 ! 243 ! k_top(:,:) = MIN( 1 , k_bot(:,:) ) ! = 1 over the ocean point, =0 elsewhere 228 ! 244 229 ! 245 230 END SUBROUTINE zgr_msk_top_bot
Note: See TracChangeset
for help on using the changeset viewer.
