Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.F90
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.F90
r5082 r5087 3 3 ! 4 4 SUBROUTINE yamada_c(ngrid,timestep,plev,play & 5 &,pu,pv,pt,d_u,d_v,d_t,cd,q2,km,kn,kq,d_t_diss,ustar &6 &,iflag_pbl)5 ,pu,pv,pt,d_u,d_v,d_t,cd,q2,km,kn,kq,d_t_diss,ustar & 6 ,iflag_pbl) 7 7 USE dimphy, ONLY: klon, klev 8 8 USE print_control_mod, ONLY: prt_level … … 129 129 fsm(ri)=1.96*(0.1912-ri)*(0.2341-ri)/((1.-ri)*(0.2231-ri)) 130 130 fl(zzz,zl0,zq2,zn2)= & 131 &max(min(l0(ig)*kap*zlev(ig,k)/(kap*zlev(ig,k)+l0(ig)) &132 &,0.5*sqrt(q2(ig,k))/sqrt(max(n2(ig,k),1.e-10))) ,1.)131 max(min(l0(ig)*kap*zlev(ig,k)/(kap*zlev(ig,k)+l0(ig)) & 132 ,0.5*sqrt(q2(ig,k))/sqrt(max(n2(ig,k),1.e-10))) ,1.) 133 133 134 134 … … 242 242 DO ig=1,ngrid 243 243 zlev(ig,nlev)=zlay(ig,nlay) & 244 &+( zlay(ig,nlay) - zlev(ig,nlev-1) )244 +( zlay(ig,nlay) - zlev(ig,nlev-1) ) 245 245 ENDDO 246 246 !!!!!! <----
Note: See TracChangeset
for help on using the changeset viewer.