Changeset 10854
- Timestamp:
- 2019-04-09T12:31:34+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/AMM15_v3_6_STABLE_package_collate_coupling/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r10478 r10854 122 122 !! *** FUNCTION zdf_gls_alloc *** 123 123 !!---------------------------------------------------------------------- 124 ALLOCATE( mxln(jpi,jpj,jpk) , zwall(jpi,jpj,jpk) , & 125 & ustars2(jpi,jpj) , ustarb2(jpi,jpj) , & 124 ALLOCATE( mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) , & 126 125 & rsbc_tke1(jpi,jpj), rsbc_tke3(jpi,jpj) , & 127 & rsbc_psi1(jpi,jpj), rtrans(jpi,jpj), STAT= zdf_gls_alloc ) 126 & rsbc_psi1(jpi,jpj), rtrans(jpi,jpj), & 127 & ustars2(jpi,jpj), ustarb2(jpi,jpj) , STAT= zdf_gls_alloc ) 128 128 ! 129 129 IF( lk_mpp ) CALL mpp_sum ( zdf_gls_alloc ) … … 410 410 ! 411 411 ! One level below 412 en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1_default * ((zhsro(:,:)+fsdepw(:,:,2)) & 413 & / zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 412 en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1_default * ((zhsro(:,:)+fsdepw(:,:,2))/zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 414 413 en(:,:,2) = MAX(en(:,:,2), rn_emin ) 415 414 z_elem_a(:,:,2) = 0._wp … … 431 430 z_elem_a(:,:,2) = 0._wp 432 431 zkar(:,:) = (rl_sf + (vkarmn-rl_sf)*(1.-exp(-rtrans_default*fsdept(:,:,1)/zhsro(:,:)) )) 433 zflxs(:,:) = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) & 434 & * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 432 zflxs(:,:) = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 435 433 436 434 en(:,:,2) = en(:,:,2) + zflxs(:,:)/fse3w(:,:,2)
Note: See TracChangeset
for help on using the changeset viewer.