Changeset 11285


Ignore:
Timestamp:
2019-07-18T11:48:27+02:00 (13 months ago)
Author:
dford
Message:

Merge in latest revision of dev_r5518_obs_oper_update.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_kd490/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r9463 r11285  
    115115      !!                ***  FUNCTION zdf_gls_alloc  *** 
    116116      !!---------------------------------------------------------------------- 
    117       ALLOCATE( mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) ,     & 
    118          &      ustars2(jpi,jpj) , ustarb2(jpi,jpj)   , STAT= zdf_gls_alloc ) 
     117      ALLOCATE( en(jpi,jpj,jpk),  mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) ,     & 
     118         &      avt_k (jpi,jpj,jpk) , avm_k (jpi,jpj,jpk),                    & 
     119         &      avmu_k(jpi,jpj,jpk) , avmv_k(jpi,jpj,jpk),                    & 
     120         &      ustars2(jpi,jpj), ustarb2(jpi,jpj)                      , STAT= zdf_gls_alloc ) 
    119121         ! 
    120122      IF( lk_mpp             )   CALL mpp_sum ( zdf_gls_alloc ) 
     
    322324      !  
    323325      ! One level below 
    324       en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2)) & 
    325           &            / zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
     326      en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2))/zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
    326327      en(:,:,2) = MAX(en(:,:,2), rn_emin ) 
    327328      z_elem_a(:,:,2) = 0._wp  
     
    344345      z_elem_a(:,:,2) = 0._wp 
    345346      zkar(:,:)       = (rl_sf + (vkarmn-rl_sf)*(1.-exp(-rtrans*fsdept(:,:,1)/zhsro(:,:)) )) 
    346       zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) & 
    347            &                      * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
     347      zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
    348348 
    349349      en(:,:,2) = en(:,:,2) + zflxs(:,:)/fse3w(:,:,2) 
Note: See TracChangeset for help on using the changeset viewer.